US20150177400A1 - Detecting and correcting changes in signal polarity for seismic data processing - Google Patents
Detecting and correcting changes in signal polarity for seismic data processing Download PDFInfo
- Publication number
- US20150177400A1 US20150177400A1 US14/409,515 US201314409515A US2015177400A1 US 20150177400 A1 US20150177400 A1 US 20150177400A1 US 201314409515 A US201314409515 A US 201314409515A US 2015177400 A1 US2015177400 A1 US 2015177400A1
- Authority
- US
- United States
- Prior art keywords
- seismic
- stack
- data
- traces
- polarity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000012545 processing Methods 0.000 title claims description 33
- 238000001514 detection method Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 109
- 230000008569 process Effects 0.000 claims description 21
- 230000007774 longterm Effects 0.000 claims description 5
- 238000010795 Steam Flooding Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000013500 data storage Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims 1
- 230000007246 mechanism Effects 0.000 abstract description 3
- 230000001131 transforming effect Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 41
- 238000003860 storage Methods 0.000 description 28
- 230000015572 biosynthetic process Effects 0.000 description 15
- 238000005755 formation reaction Methods 0.000 description 15
- 238000003672 processing method Methods 0.000 description 8
- 239000012530 fluid Substances 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 238000013507 mapping Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/366—Seismic filtering by correlation of seismic signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/42—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/123—Passive source, e.g. microseismics
- G01V2210/1234—Hydrocarbon reservoir, e.g. spontaneous or induced fracturing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/129—Source location
- G01V2210/1299—Subsurface, e.g. in borehole or below weathering layer or mud line
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1429—Subsurface, e.g. in borehole or below weathering layer or mud line
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/322—Trace stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
Definitions
- Reservoir characterization may be accomplished in a variety of ways for modeling behavior of fluids within a reservoir under different sets of circumstances and for finding optimal production techniques to maximize production.
- Seismic and microseismic surveys may be used for many applications, including for characterizing: structure, lithology, fractures, and fluid distribution in a reservoir.
- One example of an application for seismic or microseismic surveying is in hydraulic fracturing operations, wherein microseismic monitoring can be used to track the propagation of a hydraulic fracture through a formation. Seismic events can be detected, located and used to detect propagation of fluids/fractures) in a formation.
- Software may provide modeling, survey design, microseismic event detection and location (which may be optionally automated), uncertainty analysis, data integration, and visualization for interpretation.
- Seismic events are acoustic events generated by a force, such as a force generated by airguns, vibroseis systems, dynamite, etc.
- Microseismic events are elastic events generated by rock movement. Microseismic events may be generated during hydraulic fracturing of a formation as well as during other oilfield operations such as fluid production; water, gas or steam flooding; or formation compaction.
- a seismic source may be used to induce seismic waves in the earth.
- the sources may be natural subterranean events or induced subterranean events, such as propagation of hydraulics induced fractures through an earth formation.
- the seismic waves may propagate through the earth and may be transmitted, reflected, and diffracted by formations or discontinuities in the earth.
- the seismic waves may be detected by a plurality of sensors that may be positioned at the surface or within the earth, i.e., in shallow pilot wells or down boreholes in the earth. Each sensor may be used to monitor and record the seismic wavefield.
- the data received and recorded by a sensor is collectively referred to as a trace.
- the collection of traces from one or more sensors may be processed to gain information about the earth's subsurface or stored for later processing.
- Seismic and microseismic events can be characterized by a moment tensor that describes a unique radiation pattern having a polarity component and an amplitude component of the compressional and shear seismic energy radiated from the source of the seismic event.
- a large-aperture seismic array may be useful to observe changes in polarity and amplitude of the energy radiated. Inversion of seismic moment tensors can provide a way to characterize microseismic events to gain an understanding of stresses and strains in a subterranean field, including an understanding of the orientation and propagation of fractures in the field.
- Non-linear event location methods may involve selection and time picking of discrete microseismic arrivals for each of a plurality of seismic detectors and mapping to visually locate the source of the measured microseismic energy.
- the discrete time picks for each seismic detector must correspond to the same arrival of either a “P” or “S” wave and measure an arrival originating from a unique microseismic event.
- Imaging approaches relating to detection and location of seismic events may involve summing the signals recorded by a seismic array. However, when there are changes in the polarity of the moment tensor, signals having opposing polarity cancel each other out during stacking computations, rather than summing constructively. As such, in a summing technique information may be lost due to polarity changes.
- Feedback on seismic events occurring in the formation can be used to plan various stages of wellsite operations.
- Feedback on seismic events determined in real-time enables operators of wellsite procedures to intervene, direct or redirect the operations during the process to optimize the procedure and/or avoid undesired results.
- Some feedback methods involve processing microseismic event locations by mapping microseismic arrival times and polarization information into three-dimensional space through the use of modeled travel times and/or ray paths. Travel time look-up tables may be generated by modeling seismic data/an earth formation for a given velocity model.
- One mapping method for seismic/microseismic interpretation is commonly known as the “Non-Linear Event Location” method. Non-linear event location has also been used to locate macroseismic events such as earthquakes. Additional information on the topic can be found in U.S. Pat. No. 7,391,675 to Drew, incorporated herein in its entirety for all purposes, and the references listed therein.
- a method for seismic data processing by nonlinear stacking methods.
- the method can include receiving a seismic data signal comprising a plurality of traces; deriving a non-linear stack as a function of the seismic data signal; and detecting a polarity reversal in at least one of the plurality of traces based on the non-linear stack.
- a method can include disposing a sensor network in a seismic survey geometry; obtaining seismic survey data with the sensor network; correcting for move-out for a given horizon in the seismic survey data; segmenting the seismic survey data into a plurality of local groups; calculating a non-linear stack from the seismic survey data for each one of the local groups; locating a polarity reversal in the seismic survey data within the non-linear stack for each local group based on the non-linear stack for the particular group; and reversing the phase for the located polarity reversal in the seismic survey data within the non-linear stack for each local group.
- a system for acquiring and processing the seismic data.
- the system can include a seismic acquisition network of seismic sensors configured to acquire a seismic survey comprising: a plurality of data traces and a data processing apparatus.
- the data processing apparatus receives the plurality of data traces acquired in a formation from the seismic acquisition apparatus; computes a non-linear stack as a function of the seismic data signal; detects a polarity reversal in at least one of the plurality of traces based on the non-linear stack; and outputs a phase reversal to correct the trace in which the polarity reversal is detected.
- a system for acquiring seismic data and processing the seismic data.
- the system can include a data processing apparatus and a seismic acquisition network of seismic sensors configured to acquire a seismic survey.
- the data processing apparatus having a processor and a data storage device can perform a data processing method disclosed above for polarity detection and correction.
- FIG. 1 illustrates a wellsite system, in accordance with one or more implementations of various techniques described herein.
- FIG. 2 shows a first seismic data processing method work flow, in accordance with one or more implementations of various techniques described herein.
- FIG. 3 shows a plot of a gathered seismic signal of eight seismic data traces, i.e., a synthetic signal plus real noise, on 9 channels, plotted versus time (milliseconds).
- FIG. 4 shows a plot of the seismic signal of FIG. 3 with specific quantities of interest before a first iteration of the processing method of the present disclosure.
- FIG. 5 shows a plot of the seismic signal of FIG. 3 with specific quantities of interest after three iterations.
- FIG. 6 shows the data traces of FIG. 3 windowed by the AGCed STA/LTA of the phase-weighted stack, where “AGCed STA/LTA” stands for Automatic Gain Controlled ratio of the short term average (STA) to the long term average (LTA).
- AGCed STA/LTA stands for Automatic Gain Controlled ratio of the short term average (STA) to the long term average (LTA).
- FIG. 7 shows plots of correlation coefficients between channels 1 - 8 and the AGCed STA/LTA of the phase-weighted stack of the seismic data of FIG. 3 .
- FIG. 8 shows a plot of a final phase-weighted stack after three iterations of processing are complete.
- FIG. 9 shows a plot of the seismic data as a function of the AGCed STA/LTA of the final phase-weighted stack after three iterations are complete.
- FIG. 10 shows a plot of the correlation coefficients by channel after four iterations are complete.
- FIG. 11 shows a seismic data processing method work flow, in accordance with one or more implementations of various techniques described herein.
- FIGS. 12A and 12B show plots of a gathered seismic signal of many seismic data traces, i.e., a synthetic signal plus real noise.
- FIG. 12A shows the synthetic signal with added noise
- FIG. 12B shows the synthetic signal isolated from the noise.
- FIG. 13 shows a magnified view of all the traces from the area highlighted in a rectangle on FIG. 12 .
- FIG. 13A shows the highlighted portion of synthetic signal with added noise
- FIG. 13B shows the highlighted portion of the synthetic signal isolated from the noise.
- FIG. 14 shows a further magnified view of all the traces from the area highlighted in a rectangle on FIG. 13 .
- FIG. 14A shows the further magnified portion of the synthetic signal with added noise while FIG. 14B shows the further magnified portion of the synthetic signal isolated from the noise.
- FIG. 15A shows a plot of ideal correlation coefficients for representative groups of traces for the signal of FIG. 12 while FIG. 15B shows a plot of the actual correlation coefficients computed for the same representative groups.
- FIG. 16 illustrates a computer network into which implementations of various technologies described herein may be implemented.
- the embodiments may be described as a process which is depicted as a flowchart, a flow diagram, a data flow diagram, a structure diagram, or a block diagram. Although a flowchart may describe the operations as a sequential process, many of the operations can be performed in parallel or concurrently. In addition, the order of the operations may be re-arranged.
- a process is terminated when its operations are completed, but could have additional steps not included in the figure.
- a process may correspond to a method, a function, a procedure, a subroutine, a subprogram, etc. When a process corresponds to a function, its termination corresponds to a return of the function to the calling function or the main function.
- the term “storage medium” may represent one or more devices for storing data, including read only memory (ROM), random access memory (RAM), magnetic RAM, core memory, magnetic disk storage mediums, optical storage mediums, flash memory devices and/or other machine readable mediums for storing information.
- ROM read only memory
- RAM random access memory
- magnetic RAM magnetic RAM
- core memory magnetic disk storage mediums
- optical storage mediums flash memory devices and/or other machine readable mediums for storing information.
- computer-readable medium includes, but is not limited to portable or fixed storage devices, optical storage devices, wireless channels and various other mediums capable of storing, containing or carrying instruction(s) and/or data.
- embodiments may be implemented by hardware, software, firmware, middleware, microcode, hardware description languages, or any combination thereof.
- the program code or code segments to perform the necessary tasks may be stored in a machine readable medium such as storage medium.
- a processor(s) may perform the necessary tasks.
- a code segment may represent a procedure, a function, a subprogram, a program, a routine, a subroutine, a module, a software package, a class, or any combination of instructions, data structures, or program statements.
- a code segment may be coupled to another code segment or a hardware circuit by passing and/or receiving information, data, arguments, parameters, or memory contents. Information, arguments, parameters, data, etc. may be passed, forwarded, or transmitted via any suitable means including memory sharing, message passing, token passing, network transmission, etc.
- the disclosure relates to methods for data processing seismic survey data; including, but not limited to, microseismic data.
- Various methods utilize the amplitude component of acquired seismic survey data, including, but not limited to, stacking, migration, velocity analysis and attribute analysis.
- Various methods utilize the amplitude component of microseismic survey data; see, e.g., utilizing the energy onset SPE 95513, entitled “Automated Microseismic Event Detection and Location by Continuous Spatial Mapping,” Drew et al, 2005, and utilizing signal envelope, “The Source-Scanning Algorithm: Mapping the Distribution of Seismic Sources in Time and Space,” Kao, H., et al, 2004, Geophys. J.
- a semblance migration, or other stacking method stacks the entire signal, inclusive of both amplitude and phase.
- the present disclosure is directed to methods utilizing, by comparison, various non-linear stacking methods and systems that can be used for detecting and correcting signal polarity changes. This may be of use, for instance, in detecting polarity changes that occur over a focal mechanism of the seismic event, detecting polarity changes that occur within a seismic offset gather due to AVO (amplitude versus offset) variations, detecting polarity changes that occur in the recording system, or detecting polarity changes attributed to other reasons, and correcting for such polarity reversals.
- FIG. 1 illustrates a wellsite system 100 in accordance with one or more implementations of various techniques described herein.
- the operation may also be referred to herein as the hydraulic fracturing operation, though various additional production operations are contemplated herein as well.
- the operation may be conducted in concert with an active seismic survey in order to improve the effectiveness of the operation.
- the wellsite system 100 may optionally include an operations mechanism 102 for applying an operation such as pumping a fluid (such as a fracturing fluid or proppant) into or applying pressure or seismic forces to a well bore 104 disposed in a hydrocarbon reservoir 108 .
- the wellsite system 100 may also include a surface seismic array 112 , a buried seismic array 116 , a borehole seismic array 118 , or any combination thereof.
- a seismic source 114 may be included in the system 100 .
- the hydrocarbon reservoir 108 may be within a subsurface formation 110 , such as one of sandstone, carbonate, or chalk. Fractures 106 caused by, for example, the operation may be present in the formation 110 . The fractures 106 may improve the flow of hydrocarbons to the wellbore 104 .
- Receivers may be used to collect seismic/microseismic, monitoring data; and can be deployed in a number of ways: inserted near the depth of investigation (i.e., at the depth of a fracture) in a borehole (either the well being monitored or a nearby well); placed on the surface; buried relatively shallowly compared to formation depths (for example, at a range of about 1 to 300 feet); and any combination thereof.
- the surface seismic array 112 , the buried seismic array 116 , and the borehole seismic array 118 may include standard seismic receiver arrays used in seismic surveying and may include geophones, receivers, or other seismic sensing equipment configured to acquire seismic or microseismic data.
- the surface seismic array 112 may be positioned on the surface as a network of sensors.
- the buried seismic array 116 may be buried at any depth. In an embodiment, the buried seismic array 116 may be buried at a depth that is shallow relative to the depth of the formation 110 such as, but not limited to, from a range of about 1 to 300 feet.
- the borehole seismic array 118 may be placed permanently in the borehole, for example, or conveyed to any depth via a wireline cable, coiled tubing or the like.
- the seismic source 114 may be a standard seismic source used in seismic surveying, such as a vibroseis or dynamite.
- the seismic source 114 may be located on the surface or in a borehole.
- the seismic source 114 , the surface seismic array 112 , the buried seismic array 116 and/or the borehole seismic array 118 may be used to perform a seismic survey during the operation.
- the wellsite system 100 may also include a data processing apparatus 120 in a truck (or other convenient location, for instance, such as a vehicle, shelter or building) with computer equipment, discussed further below, to receive data acquired during a seismic survey by the surface seismic array 112 , the buried seismic array 116 , the borehole seismic array 118 , or any combination thereof according to a survey geometry.
- a data connection may be present between the surface seismic array 112 , the buried seismic array 116 , the borehole seismic array 118 and the data processing apparatus 120 ; and the data connection may be a wired connection or a wireless connection.
- a satellite 122 may be included, providing a wireless data connection between the surface seismic array 112 , the buried seismic array 116 , and the borehole seismic array 118 (or any combination thereof) to computer equipment, discussed further below, located at the data processing location.
- the manner of telemetry between the various system components is intended to be inclusive of all wired and wireless communications known and developed in the future; it is not the focus of the present disclosure.
- the data from the surface seismic array 112 , the buried seismic array 116 , and/or the borehole seismic array 118 may be pre-processed (for example, by the data processing apparatus 120 ) before the methods described in this patent are applied to the data.
- the pre-processing might include any one or more of the following steps in any order: filtering, transforming, interpolating, combinations of 3-component signals (for instance XYZ sensors), and modifying the time-index; the steps can be applied to either a single or several traces together, either deterministic or data-dependent.
- the specific individual processes may include, but are not limited to, frequency filtering, time resampling, deconvolution, signal rotation or re-orientation, noise-attenuation, FK filtering, tau-p filtering, time-variant move-out corrections, static-time corrections, modelled-time corrections, interpolation of traces, interpolation of sample values, and the like.
- the data as processed according to the methods described herein may be pre-processed upon receipt from the surface seismic array 112 , the buried seismic array 116 , and the borehole seismic array 118 .
- a seismic survey may be used to improve the effectiveness of a fracturing operation. For example, by performing a seismic survey during the fracturing operation, it may be possible to identify where in the formation 110 the fractures 106 are induced. In other implementations, seismic surveys may be used to improve the effectiveness of a fluid production operation; the effectiveness of a water, gas or steam flooding operation; and the effectiveness of a formation compaction operation.
- FIG. 2 shows a seismic data processing method work flow 250 , in accordance with one or more implementations of various techniques described herein.
- the method 250 begins with receiving a seismic data signal comprising a plurality of data traces at 200 .
- the method proceeds with computing a non-linear stack based on the seismic data signal at 202 .
- any non-linear or linear stacking method could be applied in accordance with the present processing method, including: phase normalized stack, an N th root stack, double-phase-weighted stack, double-phase-weighted N th root stack, real phasor stack, squared real phasor stack, or some combination of these and other non-linear stacks; some of which are described in an application filed concurrently herewith as “Methods and Systems for Seismic Data Processing by Nonlinear Stacking Methods” (Atty. Docket Number IS12.2232.1). Some of the non-linear stacks are described briefly below: N th Root Stack
- a product may be calculated by an equation according to the form:
- y NR ( t,n ) sgn[u ( t,n )]
- sgn( ⁇ circumflex over (x) ⁇ n (t)) is the sign or signum function (preserving the sign of x n (t) when the N th power is applied), and the summation or Nth root stack (see, e.g., Kanasewich et al., 1973) may be calculated by an equation according to the following form:
- subscript i denotes the different receiver signals in the data
- sgn[x i (t)] represents the sign or signum function (preserving the sign of x i (t) when the N th power is applied)
- L is the number of traces used in the stack
- phase-weighted stack see, e.g., Schimmel and Paulssen, 1997, can be calculated according to the form:
- ⁇ x(t) ⁇ is defined as a Hilbert transform of the signal.
- a Hilbert transform is an analytical transform discussed in further detail in U.S. Pat. No. 6,748,330, which is hereby incorporated in its entirety by reference for all purposes.
- the stack may be referred to as “phase weighted” as each sample of the linear stack will be weighted by the coherency of its instantaneous phases.
- phase weighted Nth root stack may be derived according to an equation of the form:
- y P (t) is the phase stack derived in 202 and y NR (t,n) is the n th root stack derived in 204 ; the resulting phase weighted Nth root stack y PWNR (t) is a product of the two.
- a double-phase-weighted stack may be governed by a function of the form:
- a double-phase-weighted Nth root stack may be governed by an equation of the form:
- y DPWNR ( t,m,n ) [ y DP ( t )] m y NR ( t,n ) eq. 9
- the normalized data may be referred to as real-phasor (RP) x RP (t).
- RP real-phasor
- in this context, can be derived from an analytic signal x A (t), defined as in eq. 5, or:
- x RP ⁇ ( t ) x ⁇ ( t ) ⁇ x A ⁇ ( t ) ⁇ eq . ⁇ 10
- Eq. 10 is a definition of real-phasor
- Eq. 11 shows an identity or an alternative definition, which is defined by a phase function only. Depending on the data availability, the real-phasor may be computed by eq. 10 or 11.
- the transformed signal can be represented by a function of the form:
- y DPWSRP ( t,m ) [ y DP ( t )] m y SRR ( t ) eq. 14
- y DPWSRP (t,m), y DP (t), and y SRP (t) are defined above.
- the linear stack has the form of:
- Lin is the number of folds in the linear stack. It is included here for comparison.
- the method proceeds with computing a function that emphasises, or “detects,” a particular time interval (or time window) containing a signal, such as the “short term average long term average ratio” (referred to herein as STA/LTA) of the non-linear stack (which was computed in 202 ) at 204 .
- STA/LTA short term average long term average ratio
- y i can represent the non-linear stack from 202 .
- the number of points in a short-term window can be represented as ns
- the number of points in a long-term window can be represented as nl, with nl>ns.
- the average energies in the short and long term windows preceding the time index i can be represented in the form:
- the ratio STA/LTA can be directly calculated.
- the method of detection of the signal itself is not the focus of the present disclosure, however, as there are many detection methods, there are many forms of STA/LTA in different references and in common use (e.g., in the definition of the sum limits).
- the detection function we use the detection function to illustrate certain time windows that may be most likely to contain signal.
- the STA/LTA function or other detection function can, for instance, be used in a cross multiplication for each time sample with the stacked function, represented by an expression such as:
- the STA/LTA ratio is from stacked function y(t) and can be replaced by any detection function and any non-linear stacking function, as described above.
- the detection function increases in value for time windows containing coherent signal in the stacked function; and provides a further enhancement of the stacked signals by multiplication of the functions when compared to the background noise.
- the window containing the signal could also be defined by calculation, prediction, or inspection of the data.
- the method proceeds with defining a window about each input trace (from 200 ) using the STA/LTA at 206 , or other detection method.
- the STA/LTA may also be processed by automatic gain control and applied to define a window framing about each trace.
- the method proceeds with cross-correlating the windowed input traces (from 206 ) with the non-linear stack (from 204 ) at 208 .
- the method proceeds with detecting, or locating, the trace amongst the cross-correlated traces that contains the largest correlation error and reversing the polarity of the trace having the largest correlation error at 210 , such that upon stacking (or summing), portions with reversed polarity do not cancel, but constructively add, resulting in more accurate detection of the seismic signal amongst the noise.
- the method iterates at 212 to repeat blocks 202 through 210 for each of the plurality of traces until polarity reversals in the plurality of traces have been detected, or until the channels for the plurality of traces correlate to within a target error threshold.
- the linear or non-linear stacking ( 202 ) is used when the signal-to-noise ratio is low in many seismic surveys. But if the acquired data is of high quality and the SNR is above a desired threshold, then the stacking step ( 202 ) may be skipped and a polarity detection function (e.g. 204 ) may applied to the data.
- a polarity detection function e.g. 204
- FIG. 3 illustrates an example, showing a plot of a seismic signal of eight seismic data traces, i.e., a synthetic signal plus real noise, on 9 channels, plotted versus time (milliseconds).
- the SNR is ⁇ 8 dB.
- Four of the eight channels have had their signal components reversed in polarity for the purpose of illustrating the method of FIG. 2 .
- channels 2 , 4 , 5 , and 6 specifically show a polarity reversal.
- FIG. 4 shows a plot of the seismic signal of FIG. 3 with specific quantities of interest before a first iteration of the processing method of FIG. 2 is applied.
- Channel ⁇ 1 shows the total data signal (i.e., the synthetic signal combined with a real noise signal), plotted as a function of time.
- Channel 0 shows the synthetic signal, isolated for clarity, as a function of time.
- Channel 1 shows one form of non-linear stack, specifically, the phase stack, as a function of time.
- Channel 2 shows a linear stack of the eight channels.
- Channel 3 shows the phase-weighted stack (such as may be computed in block 202 of FIG. 2 ) and
- Channel 4 shows the STA/LTA of the phase-weighted stack (such as the one computed in block 204 of FIG. 2 ).
- FIG. 5 shows a plot of the seismic signal of FIG. 3 with the quantities of interest shown in FIG. 4 after three iterations.
- Channel ⁇ 1 shows the total data signal (i.e., the synthetic signal combined with a real noise signal), plotted as a function of time.
- Channel 0 shows the synthetic signal, isolated for clarity, as a function of time.
- Channel 1 shows one form of non-linear stack, that is, the phase stack, as a function of time after 3 iterations of the process of FIG. 2 .
- Channel 2 shows a linear stack of the signal of Channel ⁇ 1 as a function of time after 3 iterations of the process of FIG. 2 .
- Channel 3 shows the phase-weighted stack (such as the one computed in block 202 of FIG.
- Channel 4 shows the STA/LTA of the phase-weighted stack (such as the one computed in block 204 of FIG. 2 ) after 3 iterations of the process of FIG. 2 .
- the result of 3 iterations through the method is that the SNR improves and the stack results are more similar to the synthetic data signal.
- the polarity reversals in channels 2 , 4 , and 5 have been detected and switched relative to the polarities shown in FIG. 3 .
- FIG. 6 shows various data traces windowed by the automatic gain controlled (AGCed) STA/LTA (channel 4 from FIG. 5 ) after 3 iterations of processing according to a method of the present disclosure.
- Channels 1 - 8 show the windowed input traces
- channel 9 shows the windowed phase-weighted stack (that is, a product of channel 3 and channel 4 from FIG. 5 ).
- AGCed automatic gain controlled
- FIG. 7 shows plots of correlation coefficients between channels 1 - 8 and the AGCed STA/LTA of the phase-weighted stack of the seismic data of FIG. 3 , confirming that only channel 6 still has reversed polarity after 3 iterations of the process of FIG. 2 , because each of the other channels (which have been corrected) correlates positively with the synthetic signal.
- FIG. 8 shows a plot of a final phase-weighted stack after four iterations of processing are complete; all four polarity reversals as corrected by the process are shown.
- Channel ⁇ 1 shows the total data signal (i.e., the synthetic signal combined with a real noise signal) plotted over time.
- Channel 0 shows the synthetic signal isolated for clarity over time.
- Channel 1 shows one form of non-linear stack, that is, the phase stack, as a function of time, after all iterations of the process of FIG. 2 , such that all four polarity reversals have been detected.
- Channel 2 shows a linear stack of all eight channels as a function of time after 4 iterations of the process of FIG. 2 .
- Channel 3 shows the phase-weighted stack (such as that computed in block 202 of FIG. 2 ) after 4 iterations of the process of FIG. 2
- Channel 4 shows the STA/LTA of the phase-weighted stack (such as that computed in block 204 of FIG. 2 ) after 4 iterations of the process of FIG. 2 .
- the result of all four iterations through the method is that the SNR improves further, still, and the stack results are more closely correlated to the synthetic data signal.
- FIG. 9 shows a plot of the seismic data as a function of the automatic gain controlled STA/LTA of the final phase-weighted stack after four iterations are complete. Specifically, the polarity reversal present for channel 6 in FIG. 6 is now corrected (having been detected).
- FIG. 10 shows a plot of the correlation coefficients by channel after four iterations are complete. All have positive values and confirms that the non-linear stack positively correlates to the original synthetic signal for all channels.
- the method shown in FIG. 2 can automatically detect polarity changes in traces and reverse them. Once the polarity reversals are detected and corrected, weak microseimic events may be more effectively detected. As mentioned above, even though the example was generated using phase-weighted stack, any other non-linear stacks may also be used.
- the traces may be further processed using any desired processes.
- the data may be processed for determining a parameter for a subterranean property.
- the data may be used for detecting or determining a microseismic event in a production operation.
- the microseismicity gives rise to an event, such as a fracturing operation; or a water, gas or steam flooding operation; may also be evaluated.
- the data with much improved SNR may also be used for creating an image of a subterranean structure.
- FIG. 11 shows another seismic data processing method work flow in accordance with one or more implementations of various techniques described herein.
- This method may be used to detect polarity reversals along a seismic horizon under poor SNR conditions.
- This method may be used as an automatic fault plane solution (e.g. beach ball diagram) on surface arrays; or automatic polarity reversal correction to enable semblance-based processing.
- the method begins with disposing a sensor network in a seismic survey geometry and dividing data from a surface sensor network into local groups at 1100 , such that data may be stacked within the local groups to increase the signal to noise ratio.
- the particular group size may in some embodiments be a parameter set by the user, and may depend on the particular problem. In an embodiment, the particular group size may be selected in such a way, based on data signal polarities, that the data in each group has signal that can be stacked constructively, added, or combined in some manner.
- the method includes obtaining seismic survey data with the sensor network at 1101 .
- the method further includes segmenting the seismic survey data a second time at 1102 into segments.
- Segmenting refers to associating specific portions of the seismic data (or seismic gather) to one of the two polarities that exist.
- segmenting can refer to labeling each trace (or part of each trace around a seismic horizon) with a positive polarity or a negative polarity.
- the segment size may be as small as one.
- the method further includes selecting a horizon among the obtained seismic survey data at 1103 .
- Moveout can be defined as relating to the difference in the arrival times or travel times of a wave measured by receivers at two different offset locations.
- One form of moveout is caused by the separation between the source and the receiver.
- the method proceeds with correcting for move-out for the selected horizon in the seismic survey data at 1104 .
- Various techniques of correcting for move-out are known, and any of such techniques could be employed at 1104 .
- the method proceeds with calculating a non-linear stack of the data associated with the horizon, which has been corrected for move-out at 1106 for each local group.
- the method proceeds with locating a polarity reversal in the seismic survey data based on the non-linear stack (calculated in 1106 for each local group) at 1108 , which may be performed using the LTA/STA (as described above) or any other detection method.
- the method further includes correcting the polarity of at least one segment of the seismic survey data according to the location of the polarity reversal at 1110 .
- FIG. 12A shows a plot of a gathered seismic signal of many seismic data traces, i.e., a synthetic signal plus real noise.
- a synthetic signal plus real noise We arbitrarily show only one trace in 5 from a large range of traces and display them at a convenient trace spacing to show the variation in traces, amplitude, and plotting offset (in meters) versus time (milliseconds).
- the signal to noise ratio is ⁇ 18 dB.
- FIG. 12B shows the synthetic signal isolated from noise, with the polarity reversed along certain segments, plotting offset (in meters) versus time (in milliseconds).
- FIG. 12B A first polarity reversal occurs at just before 900 m and a second polarity reversal occurs at roughly 1200 m, and a third at 1650 m, shown in FIG. 12B .
- a rectangle highlights the section containing several polarity changes, aligning with where the local groups of the sensor network are and with the manner in which the data are segmented into local groups and stacked within the local groups, as described above with reference to FIG. 11 .
- FIGS. 13A and B show magnified views of all the traces from the areas highlighted in a rectangle on FIGS. 12A and 12B in order to highlight a single polarity reversal. Whereas, FIG. 13A shows the synthetic signal with noise added, and FIG. 13B shows the synthetic signal isolated from noise.
- FIGS. 14A and 14B show a further magnified view of the area highlighted in the rectangle on FIG. 13 , showing just the one segment of data containing the polarity reversal at about 880 m.
- FIG. 14B shows the synthetic signal with noise added
- FIG. 14A shows the synthetic signal without the noise added.
- FIG. 15A shows a plot of ideal correlation coefficients for representative groups of traces for the signal of FIG. 12A
- FIG. 15B shows a plot of the actual correlation coefficients computed for the same representative groups.
- the data is segmented into local groups (as described above), phase-weighted N th root stack applied and the cross-correlation computed according to the method of FIG. 11 .
- the polarity can be handled so that the non-linear stack is cumulatively added, and polarity reversals do not mathematically cancel out portions of the data.
- FIG. 16 illustrates a computing system 1600 , into which implementations of various technologies described herein may be implemented.
- the computing system 1600 may include one or more system computers (processors) 1630 , which may be implemented as any conventional personal computer or server.
- processors system computers
- FIG. 16 illustrates a computing system 1600 , into which implementations of various technologies described herein may be implemented.
- the computing system 1600 may include one or more system computers (processors) 1630 , which may be implemented as any conventional personal computer or server.
- system computers processors
- FIG. 16 illustrates a computing system 1600 , into which implementations of various technologies described herein may be implemented.
- the computing system 1600 may include one or more system computers (processors) 1630 , which may be implemented as any conventional personal computer or server.
- HTTP hypertext transfer protocol
- the system computer 1630 may be in communication with disk storage devices 1629 , 1631 , and 1633 , which may be external hard disk storage devices. It is contemplated that disk storage devices 1629 , 1631 , and 1633 are conventional hard disk drives, and as such, will be implemented by way of a local area network or by remote access. Of course, while disk storage devices 1629 , 1631 , and 1633 are illustrated as separate devices, a single disk storage device may be used to store any and all the program instructions, measurement data, and results as desired.
- seismic data from the receivers may be stored in disk storage device 1631 .
- the system computer 1630 may retrieve the appropriate data from the disk storage device 1631 to process seismic data according to program instructions that correspond to implementations of various technologies described herein.
- the program instructions may be written in a computer programming language, such as C++, Java and the like.
- the program instructions may be stored in a computer-readable medium, such as program disk storage device 1633 .
- Such computer-readable media may include computer storage media and communication media.
- Computer storage media may include volatile and non-volatile media, and removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules or other data.
- Computer storage media may further include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid state memory technology, CD-ROM, digital versatile disks (DVD), or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the system computer 1630 .
- RAM random access memory
- ROM read-only memory
- EPROM erasable programmable read-only memory
- EEPROM electrically erasable programmable read-only memory
- flash memory or other solid state memory technology
- CD-ROM compact discs
- DVD digital versatile disks
- magnetic cassettes magnetic tape
- magnetic disk storage magnetic disk storage devices
- Communication media may embody computer-readable instructions, data structures, program modules or other data in a modulated data signal, such as a carrier wave or other transport mechanism. It may include any information delivery media.
- modulated data signal may mean a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal.
- communication media may include wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above may also be included within the scope of computer readable media.
- system computer 1630 may present output primarily onto graphics display 1627 , or alternatively via printer 1628 .
- the system computer 1630 may store the results of the methods described above on disk storage 1629 , for later use and further analysis.
- the keyboard 1626 and the pointing device (e.g., a mouse, trackball, or the like) 1625 may be provided with the system computer 1630 to enable interactive operation.
- the system computer 1630 may be located at a data center remote from the survey region.
- the system computer 1630 may be in communication with the receivers (either directly or via a recording unit, not shown), to receive signals indicative of the reflected seismic energy. These signals, after conventional formatting and other initial processing, may be stored by the system computer 1630 as digital data in the disk storage 1631 for subsequent retrieval and processing in the manner described above. While FIG. 16 illustrates the disk storage 1631 as directly connected to the system computer 1630 , it is also contemplated that the disk storage device 1631 may be accessible through a local area network or by remote access.
- disk storage devices 1629 , 1631 are illustrated as separate devices for storing input seismic data and analysis results, the disk storage devices 1629 , 1631 may be implemented within a single disk drive (either together with or separately from program disk storage device 1633 ), or in any other conventional manner as will be fully understood by one of skill in the art having reference to this specification.
- Complex wavelets are a more general way to describe the signal, and encompasses both the complex time domain representation (wavelet is a “delta function” in time domain), frequency domain representation (wavelet is a “sine wave” in time domain) and other classes of “wavelets” that to a greater or lesser degree are compact in the time domain or frequency domain.
- the concept of the complex “wavelet” representation of the signal includes the concept of amplitude and phase.
- the methods and systems of the present disclosure may be fully automated and able to operate continuously in time for monitoring, detecting, and locating seismic and microseismic events from on or multiple seismic detectors.
- Methods of the present invention utilize the contributions from one or multiple seismic detectors.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Utilizing the phase component of a moment tensor for a seismic data signal, isolated from the amplitude component, by automatically detecting polarity changes that occur over a focal mechanism of the seismic event, and correcting for such polarity reversals. Transforming seismic (including microseismic) signals as recorded by one or more seismic detectors to enhance detection of arrivals. The transforms enable the generation of an image, or map, representative of the likelihood that there was a source of seismic energy occurring at a given point in time at a particular point in time.
Description
- This application is related to and claims priority from the following U.S. Provisional patent applications:
- Application Ser. No. 61/663,403 filed on 22 Jun. 2012, (Atty. Docket Number IS12.2232-1), titled “Seismic data processing using a phase weighted Nth root stack”;
- Application Ser. No. 61/663,416 filed on 22 Jun. 2012, (Atty. Docket Number IS12.2232-2), titled “Detecting and correcting changes in signal polarity for seismic data processing”; and
- Application Ser. No. 61/663,449 filed on 22 Jun. 2012, (Atty. Docket Number IS12.2232-3), titled “Seismic data processing by nonlinear stacking methods”; the disclosures of which are incorporated by reference herein in their entirety for all purposes.
- Reservoir characterization may be accomplished in a variety of ways for modeling behavior of fluids within a reservoir under different sets of circumstances and for finding optimal production techniques to maximize production. Seismic and microseismic surveys may be used for many applications, including for characterizing: structure, lithology, fractures, and fluid distribution in a reservoir. One example of an application for seismic or microseismic surveying is in hydraulic fracturing operations, wherein microseismic monitoring can be used to track the propagation of a hydraulic fracture through a formation. Seismic events can be detected, located and used to detect propagation of fluids/fractures) in a formation. Software may provide modeling, survey design, microseismic event detection and location (which may be optionally automated), uncertainty analysis, data integration, and visualization for interpretation.
- Seismic events are acoustic events generated by a force, such as a force generated by airguns, vibroseis systems, dynamite, etc. Microseismic events are elastic events generated by rock movement. Microseismic events may be generated during hydraulic fracturing of a formation as well as during other oilfield operations such as fluid production; water, gas or steam flooding; or formation compaction.
- In seismic surveys, a seismic source may be used to induce seismic waves in the earth. By contrast, in a microseismic survey, the sources may be natural subterranean events or induced subterranean events, such as propagation of hydraulics induced fractures through an earth formation. In both cases, the seismic waves may propagate through the earth and may be transmitted, reflected, and diffracted by formations or discontinuities in the earth. The seismic waves may be detected by a plurality of sensors that may be positioned at the surface or within the earth, i.e., in shallow pilot wells or down boreholes in the earth. Each sensor may be used to monitor and record the seismic wavefield. The data received and recorded by a sensor is collectively referred to as a trace. The collection of traces from one or more sensors may be processed to gain information about the earth's subsurface or stored for later processing.
- Seismic and microseismic events, both naturally occurring and induced, can be characterized by a moment tensor that describes a unique radiation pattern having a polarity component and an amplitude component of the compressional and shear seismic energy radiated from the source of the seismic event. A large-aperture seismic array may be useful to observe changes in polarity and amplitude of the energy radiated. Inversion of seismic moment tensors can provide a way to characterize microseismic events to gain an understanding of stresses and strains in a subterranean field, including an understanding of the orientation and propagation of fractures in the field.
- Non-linear event location methods may involve selection and time picking of discrete microseismic arrivals for each of a plurality of seismic detectors and mapping to visually locate the source of the measured microseismic energy. However, to successfully and accurately locate the microseismic event, the discrete time picks for each seismic detector must correspond to the same arrival of either a “P” or “S” wave and measure an arrival originating from a unique microseismic event.
- Imaging approaches relating to detection and location of seismic events may involve summing the signals recorded by a seismic array. However, when there are changes in the polarity of the moment tensor, signals having opposing polarity cancel each other out during stacking computations, rather than summing constructively. As such, in a summing technique information may be lost due to polarity changes.
- Feedback on seismic events occurring in the formation can be used to plan various stages of wellsite operations. Feedback on seismic events determined in real-time enables operators of wellsite procedures to intervene, direct or redirect the operations during the process to optimize the procedure and/or avoid undesired results.
- Some feedback methods involve processing microseismic event locations by mapping microseismic arrival times and polarization information into three-dimensional space through the use of modeled travel times and/or ray paths. Travel time look-up tables may be generated by modeling seismic data/an earth formation for a given velocity model. One mapping method for seismic/microseismic interpretation is commonly known as the “Non-Linear Event Location” method. Non-linear event location has also been used to locate macroseismic events such as earthquakes. Additional information on the topic can be found in U.S. Pat. No. 7,391,675 to Drew, incorporated herein in its entirety for all purposes, and the references listed therein.
- This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
- In an embodiment, a method is disclosed for seismic data processing by nonlinear stacking methods. The method can include receiving a seismic data signal comprising a plurality of traces; deriving a non-linear stack as a function of the seismic data signal; and detecting a polarity reversal in at least one of the plurality of traces based on the non-linear stack.
- In another embodiment, a method is disclosed where method can include disposing a sensor network in a seismic survey geometry; obtaining seismic survey data with the sensor network; correcting for move-out for a given horizon in the seismic survey data; segmenting the seismic survey data into a plurality of local groups; calculating a non-linear stack from the seismic survey data for each one of the local groups; locating a polarity reversal in the seismic survey data within the non-linear stack for each local group based on the non-linear stack for the particular group; and reversing the phase for the located polarity reversal in the seismic survey data within the non-linear stack for each local group.
- In one embodiment, a system is disclosed for acquiring and processing the seismic data. The system can include a seismic acquisition network of seismic sensors configured to acquire a seismic survey comprising: a plurality of data traces and a data processing apparatus. The data processing apparatus receives the plurality of data traces acquired in a formation from the seismic acquisition apparatus; computes a non-linear stack as a function of the seismic data signal; detects a polarity reversal in at least one of the plurality of traces based on the non-linear stack; and outputs a phase reversal to correct the trace in which the polarity reversal is detected.
- In yet another embodiment, a system is disclosed for acquiring seismic data and processing the seismic data. The system can include a data processing apparatus and a seismic acquisition network of seismic sensors configured to acquire a seismic survey. The data processing apparatus having a processor and a data storage device can perform a data processing method disclosed above for polarity detection and correction.
- Embodiments of a method and system for seismic data processing for detecting and correcting changes in signal polarity are described with reference to the following figures. Like numbers are used throughout the figures to reference like features and components.
-
FIG. 1 illustrates a wellsite system, in accordance with one or more implementations of various techniques described herein. -
FIG. 2 shows a first seismic data processing method work flow, in accordance with one or more implementations of various techniques described herein. -
FIG. 3 shows a plot of a gathered seismic signal of eight seismic data traces, i.e., a synthetic signal plus real noise, on 9 channels, plotted versus time (milliseconds). -
FIG. 4 shows a plot of the seismic signal ofFIG. 3 with specific quantities of interest before a first iteration of the processing method of the present disclosure. -
FIG. 5 shows a plot of the seismic signal ofFIG. 3 with specific quantities of interest after three iterations. -
FIG. 6 shows the data traces ofFIG. 3 windowed by the AGCed STA/LTA of the phase-weighted stack, where “AGCed STA/LTA” stands for Automatic Gain Controlled ratio of the short term average (STA) to the long term average (LTA). -
FIG. 7 shows plots of correlation coefficients between channels 1-8 and the AGCed STA/LTA of the phase-weighted stack of the seismic data ofFIG. 3 . -
FIG. 8 shows a plot of a final phase-weighted stack after three iterations of processing are complete. -
FIG. 9 shows a plot of the seismic data as a function of the AGCed STA/LTA of the final phase-weighted stack after three iterations are complete. -
FIG. 10 shows a plot of the correlation coefficients by channel after four iterations are complete. -
FIG. 11 shows a seismic data processing method work flow, in accordance with one or more implementations of various techniques described herein. -
FIGS. 12A and 12B show plots of a gathered seismic signal of many seismic data traces, i.e., a synthetic signal plus real noise. For illustrating the present method,FIG. 12A shows the synthetic signal with added noise whereasFIG. 12B shows the synthetic signal isolated from the noise. -
FIG. 13 shows a magnified view of all the traces from the area highlighted in a rectangle onFIG. 12 . For illustrating the method ofFIG. 11 ,FIG. 13A shows the highlighted portion of synthetic signal with added noise whereasFIG. 13B shows the highlighted portion of the synthetic signal isolated from the noise. -
FIG. 14 shows a further magnified view of all the traces from the area highlighted in a rectangle onFIG. 13 . For illustrating the method ofFIG. 11 ,FIG. 14A shows the further magnified portion of the synthetic signal with added noise whileFIG. 14B shows the further magnified portion of the synthetic signal isolated from the noise. -
FIG. 15A shows a plot of ideal correlation coefficients for representative groups of traces for the signal ofFIG. 12 whileFIG. 15B shows a plot of the actual correlation coefficients computed for the same representative groups. -
FIG. 16 illustrates a computer network into which implementations of various technologies described herein may be implemented. - In the following description, numerous details are set forth to provide an understanding of the present disclosure. However, it will be understood by those skilled in the art that the present disclosure may be practiced without these details and that numerous variations or modifications from the described embodiments are possible.
- Specific details are given in the following description to provide a thorough understanding of the embodiments. However, it will be understood by one of ordinary skill in the art that the embodiments maybe practiced without these specific details. For example, circuits may be shown in block diagrams in order not to obscure the embodiments in unnecessary detail. In other instances, well-known circuits, processes, algorithms, structures, and techniques may be shown without unnecessary detail in order to avoid obscuring the embodiments.
- Also, it is noted that the embodiments may be described as a process which is depicted as a flowchart, a flow diagram, a data flow diagram, a structure diagram, or a block diagram. Although a flowchart may describe the operations as a sequential process, many of the operations can be performed in parallel or concurrently. In addition, the order of the operations may be re-arranged. A process is terminated when its operations are completed, but could have additional steps not included in the figure. A process may correspond to a method, a function, a procedure, a subroutine, a subprogram, etc. When a process corresponds to a function, its termination corresponds to a return of the function to the calling function or the main function.
- Moreover, as disclosed herein, the term “storage medium” may represent one or more devices for storing data, including read only memory (ROM), random access memory (RAM), magnetic RAM, core memory, magnetic disk storage mediums, optical storage mediums, flash memory devices and/or other machine readable mediums for storing information. The term “computer-readable medium” includes, but is not limited to portable or fixed storage devices, optical storage devices, wireless channels and various other mediums capable of storing, containing or carrying instruction(s) and/or data.
- Furthermore, embodiments may be implemented by hardware, software, firmware, middleware, microcode, hardware description languages, or any combination thereof. When implemented in software, firmware, middleware or microcode, the program code or code segments to perform the necessary tasks may be stored in a machine readable medium such as storage medium. A processor(s) may perform the necessary tasks. A code segment may represent a procedure, a function, a subprogram, a program, a routine, a subroutine, a module, a software package, a class, or any combination of instructions, data structures, or program statements. A code segment may be coupled to another code segment or a hardware circuit by passing and/or receiving information, data, arguments, parameters, or memory contents. Information, arguments, parameters, data, etc. may be passed, forwarded, or transmitted via any suitable means including memory sharing, message passing, token passing, network transmission, etc.
- The disclosure relates to methods for data processing seismic survey data; including, but not limited to, microseismic data. Various methods utilize the amplitude component of acquired seismic survey data, including, but not limited to, stacking, migration, velocity analysis and attribute analysis. Various methods utilize the amplitude component of microseismic survey data; see, e.g., utilizing the energy onset SPE 95513, entitled “Automated Microseismic Event Detection and Location by Continuous Spatial Mapping,” Drew et al, 2005, and utilizing signal envelope, “The Source-Scanning Algorithm: Mapping the Distribution of Seismic Sources in Time and Space,” Kao, H., et al, 2004, Geophys. J. Int., 157, 589-594. A semblance migration, or other stacking method, stacks the entire signal, inclusive of both amplitude and phase. The present disclosure is directed to methods utilizing, by comparison, various non-linear stacking methods and systems that can be used for detecting and correcting signal polarity changes. This may be of use, for instance, in detecting polarity changes that occur over a focal mechanism of the seismic event, detecting polarity changes that occur within a seismic offset gather due to AVO (amplitude versus offset) variations, detecting polarity changes that occur in the recording system, or detecting polarity changes attributed to other reasons, and correcting for such polarity reversals.
-
FIG. 1 illustrates awellsite system 100 in accordance with one or more implementations of various techniques described herein. The operation may also be referred to herein as the hydraulic fracturing operation, though various additional production operations are contemplated herein as well. In thewellsite system 100, the operation may be conducted in concert with an active seismic survey in order to improve the effectiveness of the operation. Thewellsite system 100 may optionally include anoperations mechanism 102 for applying an operation such as pumping a fluid (such as a fracturing fluid or proppant) into or applying pressure or seismic forces to a well bore 104 disposed in ahydrocarbon reservoir 108. Thewellsite system 100 may also include a surfaceseismic array 112, a buriedseismic array 116, a boreholeseismic array 118, or any combination thereof. Aseismic source 114 may be included in thesystem 100. - The
hydrocarbon reservoir 108 may be within asubsurface formation 110, such as one of sandstone, carbonate, or chalk.Fractures 106 caused by, for example, the operation may be present in theformation 110. Thefractures 106 may improve the flow of hydrocarbons to thewellbore 104. - Receivers (i.e., geophones) may be used to collect seismic/microseismic, monitoring data; and can be deployed in a number of ways: inserted near the depth of investigation (i.e., at the depth of a fracture) in a borehole (either the well being monitored or a nearby well); placed on the surface; buried relatively shallowly compared to formation depths (for example, at a range of about 1 to 300 feet); and any combination thereof.
- The surface
seismic array 112, the buriedseismic array 116, and the boreholeseismic array 118 may include standard seismic receiver arrays used in seismic surveying and may include geophones, receivers, or other seismic sensing equipment configured to acquire seismic or microseismic data. The surfaceseismic array 112 may be positioned on the surface as a network of sensors. The buriedseismic array 116 may be buried at any depth. In an embodiment, the buriedseismic array 116 may be buried at a depth that is shallow relative to the depth of theformation 110 such as, but not limited to, from a range of about 1 to 300 feet. The boreholeseismic array 118 may be placed permanently in the borehole, for example, or conveyed to any depth via a wireline cable, coiled tubing or the like. - The
seismic source 114 may be a standard seismic source used in seismic surveying, such as a vibroseis or dynamite. Theseismic source 114 may be located on the surface or in a borehole. Theseismic source 114, the surfaceseismic array 112, the buriedseismic array 116 and/or the boreholeseismic array 118 may be used to perform a seismic survey during the operation. - The
wellsite system 100 may also include adata processing apparatus 120 in a truck (or other convenient location, for instance, such as a vehicle, shelter or building) with computer equipment, discussed further below, to receive data acquired during a seismic survey by the surfaceseismic array 112, the buriedseismic array 116, the boreholeseismic array 118, or any combination thereof according to a survey geometry. A data connection may be present between the surfaceseismic array 112, the buriedseismic array 116, the boreholeseismic array 118 and thedata processing apparatus 120; and the data connection may be a wired connection or a wireless connection. - Alternatively, for data processing at a data processing location remote to the wellsite, a
satellite 122 may be included, providing a wireless data connection between the surfaceseismic array 112, the buriedseismic array 116, and the borehole seismic array 118 (or any combination thereof) to computer equipment, discussed further below, located at the data processing location. The manner of telemetry between the various system components is intended to be inclusive of all wired and wireless communications known and developed in the future; it is not the focus of the present disclosure. - In an embodiment, the data from the surface
seismic array 112, the buriedseismic array 116, and/or the boreholeseismic array 118 may be pre-processed (for example, by the data processing apparatus 120) before the methods described in this patent are applied to the data. The pre-processing might include any one or more of the following steps in any order: filtering, transforming, interpolating, combinations of 3-component signals (for instance XYZ sensors), and modifying the time-index; the steps can be applied to either a single or several traces together, either deterministic or data-dependent. The specific individual processes may include, but are not limited to, frequency filtering, time resampling, deconvolution, signal rotation or re-orientation, noise-attenuation, FK filtering, tau-p filtering, time-variant move-out corrections, static-time corrections, modelled-time corrections, interpolation of traces, interpolation of sample values, and the like. Thus, the data as processed according to the methods described herein may be pre-processed upon receipt from the surfaceseismic array 112, the buriedseismic array 116, and the boreholeseismic array 118. - In one implementation, a seismic survey may be used to improve the effectiveness of a fracturing operation. For example, by performing a seismic survey during the fracturing operation, it may be possible to identify where in the
formation 110 thefractures 106 are induced. In other implementations, seismic surveys may be used to improve the effectiveness of a fluid production operation; the effectiveness of a water, gas or steam flooding operation; and the effectiveness of a formation compaction operation. -
FIG. 2 shows a seismic data processing method work flow 250, in accordance with one or more implementations of various techniques described herein. The method 250 begins with receiving a seismic data signal comprising a plurality of data traces at 200. The method proceeds with computing a non-linear stack based on the seismic data signal at 202. Though examples presented below are directed to using a phase-weighted stack, for practical purposes, any non-linear or linear stacking method could be applied in accordance with the present processing method, including: phase normalized stack, an Nth root stack, double-phase-weighted stack, double-phase-weighted Nth root stack, real phasor stack, squared real phasor stack, or some combination of these and other non-linear stacks; some of which are described in an application filed concurrently herewith as “Methods and Systems for Seismic Data Processing by Nonlinear Stacking Methods” (Atty. Docket Number IS12.2232.1). Some of the non-linear stacks are described briefly below: Nth Root Stack - For example, in a non-linear stacking method referred to as the Nth root stack, a product may be calculated by an equation according to the form:
-
y NR(t,n)=sgn[u(t,n)]|u(t,n)|n eq. 1 - where sgn({circumflex over (x)}n(t)) is the sign or signum function (preserving the sign of xn(t) when the Nth power is applied), and the summation or Nth root stack (see, e.g., Kanasewich et al., 1973) may be calculated by an equation according to the following form:
-
- wherein subscript i denotes the different receiver signals in the data, sgn[xi(t)] represents the sign or signum function (preserving the sign of xi(t) when the Nth power is applied), L is the number of traces used in the stack, and the power n is a number larger than or equal to 1 that amplifies relatively small amplitude signals. For example, calculating an Nth root stack with n=1 results in a linear stack, and n>1 leads to a reduction in the signal variance.
- In another embodiment, a phase-weighted stack, see, e.g., Schimmel and Paulssen, 1997, can be calculated according to the form:
-
- is the phase stack derived from the phase of the analytical signal
- wherein {x(t)} is defined as a Hilbert transform of the signal. A Hilbert transform is an analytical transform discussed in further detail in U.S. Pat. No. 6,748,330, which is hereby incorporated in its entirety by reference for all purposes. The stack may be referred to as “phase weighted” as each sample of the linear stack will be weighted by the coherency of its instantaneous phases.
- As described in the related application filed concurrently herewith cited above (Atty. Docket Number IS12.2232.1), the phase weighted Nth root stack may be derived according to an equation of the form:
-
y PWNR(t,m,n)=[y P(t)]m y NR(t,n) eq. 6 - where yP(t) is the phase stack derived in 202 and yNR(t,n) is the nth root stack derived in 204; the resulting phase weighted Nth root stack yPWNR(t) is a product of the two.
- Also described in the related application filed concurrently herewith cited above (Atty. Docket Number IS12.2232.1), and analogous to the phase weighted stack described above, a double-phase-weighted stack may be governed by a function of the form:
-
- is the double-phase stack derived from the phase of the analytical signal of eq. 5.
- Also described in the related application filed concurrently herewith cited above (Atty. Docket Number IS12.2232.1), and analogous to the phase weighted stack described above, a double-phase-weighted Nth root stack may be governed by an equation of the form:
-
y DPWNR(t,m,n)=[y DP(t)]m y NR(t,n) eq. 9 - where yDP(t) is the double phase stack and yNR(t,n) is the nth root stack; the result represents a product of the two.
- The normalized data may be referred to as real-phasor (RP) xRP(t). The real-phasor (RP) and the envelope function |xA(t)|, in this context, can be derived from an analytic signal xA(t), defined as in eq. 5, or:
-
- where xA(t) and φx(t) are defined in eq. 5. Eq. 10 is a definition of real-phasor, while Eq. 11 shows an identity or an alternative definition, which is defined by a phase function only. Depending on the data availability, the real-phasor may be computed by eq. 10 or 11.
- From real phasor, we may define several more non-linear stacks; mathematically, the stacked functions may be described by:
-
- where subscript i denotes the different receiver signals in the data and xRP(t) represents the normalized signal in the manner described above. These non-linear stacks in eq. 12 and 13 are referred to as real-phasor and squared real-phasor stacks, respectively.
- In yet another non-linear stack, referred to as a double phase weighted squared real-phasor stack, the transformed signal can be represented by a function of the form:
-
y DPWSRP(t,m)=[y DP(t)]m y SRR(t) eq. 14 - where yDPWSRP(t,m), yDP(t), and ySRP(t) are defined above.
- The linear stack has the form of:
-
- wherein Lin is the number of folds in the linear stack. It is included here for comparison.
- In an embodiment, the method proceeds with computing a function that emphasises, or “detects,” a particular time interval (or time window) containing a signal, such as the “short term average long term average ratio” (referred to herein as STA/LTA) of the non-linear stack (which was computed in 202) at 204. Various known methods for detection such as computation of the STA/LTA may be employed, and the details of which are not related to present novelty. In an embodiment, mathematically, yi can represent the non-linear stack from 202. The number of points in a short-term window can be represented as ns, and the number of points in a long-term window can be represented as nl, with nl>ns. Then, the average energies in the short and long term windows preceding the time index i can be represented in the form:
-
- from which the ratio STA/LTA can be directly calculated. The method of detection of the signal itself is not the focus of the present disclosure, however, as there are many detection methods, there are many forms of STA/LTA in different references and in common use (e.g., in the definition of the sum limits). In the present disclosure we use the detection function to illustrate certain time windows that may be most likely to contain signal. The STA/LTA function or other detection function can, for instance, be used in a cross multiplication for each time sample with the stacked function, represented by an expression such as:
-
- where the STA/LTA ratio is from stacked function y(t) and can be replaced by any detection function and any non-linear stacking function, as described above. The detection function increases in value for time windows containing coherent signal in the stacked function; and provides a further enhancement of the stacked signals by multiplication of the functions when compared to the background noise. The window containing the signal could also be defined by calculation, prediction, or inspection of the data.
- In an embodiment, the method proceeds with defining a window about each input trace (from 200) using the STA/LTA at 206, or other detection method. In an embodiment, the STA/LTA may also be processed by automatic gain control and applied to define a window framing about each trace.
- The method proceeds with cross-correlating the windowed input traces (from 206) with the non-linear stack (from 204) at 208.
- The method proceeds with detecting, or locating, the trace amongst the cross-correlated traces that contains the largest correlation error and reversing the polarity of the trace having the largest correlation error at 210, such that upon stacking (or summing), portions with reversed polarity do not cancel, but constructively add, resulting in more accurate detection of the seismic signal amongst the noise.
- The method iterates at 212 to repeat
blocks 202 through 210 for each of the plurality of traces until polarity reversals in the plurality of traces have been detected, or until the channels for the plurality of traces correlate to within a target error threshold. - Not all the steps in the above example 250 are necessary. For example, the linear or non-linear stacking (202) is used when the signal-to-noise ratio is low in many seismic surveys. But if the acquired data is of high quality and the SNR is above a desired threshold, then the stacking step (202) may be skipped and a polarity detection function (e.g. 204) may applied to the data.
-
FIG. 3 illustrates an example, showing a plot of a seismic signal of eight seismic data traces, i.e., a synthetic signal plus real noise, on 9 channels, plotted versus time (milliseconds). The SNR is −8 dB. Four of the eight channels have had their signal components reversed in polarity for the purpose of illustrating the method ofFIG. 2 . As shown inFIG. 3 ,channels -
FIG. 4 shows a plot of the seismic signal ofFIG. 3 with specific quantities of interest before a first iteration of the processing method ofFIG. 2 is applied. Channel −1 shows the total data signal (i.e., the synthetic signal combined with a real noise signal), plotted as a function of time.Channel 0 shows the synthetic signal, isolated for clarity, as a function of time.Channel 1 shows one form of non-linear stack, specifically, the phase stack, as a function of time.Channel 2 shows a linear stack of the eight channels.Channel 3 shows the phase-weighted stack (such as may be computed inblock 202 ofFIG. 2 ) andChannel 4 shows the STA/LTA of the phase-weighted stack (such as the one computed inblock 204 ofFIG. 2 ). -
FIG. 5 shows a plot of the seismic signal ofFIG. 3 with the quantities of interest shown inFIG. 4 after three iterations. Channel −1 shows the total data signal (i.e., the synthetic signal combined with a real noise signal), plotted as a function of time.Channel 0 shows the synthetic signal, isolated for clarity, as a function of time.Channel 1 shows one form of non-linear stack, that is, the phase stack, as a function of time after 3 iterations of the process ofFIG. 2 .Channel 2 shows a linear stack of the signal of Channel −1 as a function of time after 3 iterations of the process ofFIG. 2 .Channel 3 shows the phase-weighted stack (such as the one computed inblock 202 ofFIG. 2 ) after 3 iterations of the process ofFIG. 2 , andChannel 4 shows the STA/LTA of the phase-weighted stack (such as the one computed inblock 204 ofFIG. 2 ) after 3 iterations of the process ofFIG. 2 . The result of 3 iterations through the method is that the SNR improves and the stack results are more similar to the synthetic data signal. Specifically, the polarity reversals inchannels FIG. 3 . -
FIG. 6 shows various data traces windowed by the automatic gain controlled (AGCed) STA/LTA (channel 4 fromFIG. 5 ) after 3 iterations of processing according to a method of the present disclosure. Channels 1-8 show the windowed input traces, andchannel 9 shows the windowed phase-weighted stack (that is, a product ofchannel 3 andchannel 4 fromFIG. 5 ). As shown, of the four channels with reversed polarities originally, only channel 6 remains reversed, while the other 3 channels with polarity reversals have been reversed upon detection of the polarity reversal within each channel. -
FIG. 7 shows plots of correlation coefficients between channels 1-8 and the AGCed STA/LTA of the phase-weighted stack of the seismic data ofFIG. 3 , confirming thatonly channel 6 still has reversed polarity after 3 iterations of the process ofFIG. 2 , because each of the other channels (which have been corrected) correlates positively with the synthetic signal. -
FIG. 8 shows a plot of a final phase-weighted stack after four iterations of processing are complete; all four polarity reversals as corrected by the process are shown. Channel −1 shows the total data signal (i.e., the synthetic signal combined with a real noise signal) plotted over time.Channel 0 shows the synthetic signal isolated for clarity over time.Channel 1 shows one form of non-linear stack, that is, the phase stack, as a function of time, after all iterations of the process ofFIG. 2 , such that all four polarity reversals have been detected.Channel 2 shows a linear stack of all eight channels as a function of time after 4 iterations of the process ofFIG. 2 .Channel 3 shows the phase-weighted stack (such as that computed inblock 202 ofFIG. 2 ) after 4 iterations of the process ofFIG. 2 , andChannel 4 shows the STA/LTA of the phase-weighted stack (such as that computed inblock 204 ofFIG. 2 ) after 4 iterations of the process ofFIG. 2 . The result of all four iterations through the method is that the SNR improves further, still, and the stack results are more closely correlated to the synthetic data signal. -
FIG. 9 shows a plot of the seismic data as a function of the automatic gain controlled STA/LTA of the final phase-weighted stack after four iterations are complete. Specifically, the polarity reversal present forchannel 6 inFIG. 6 is now corrected (having been detected). -
FIG. 10 shows a plot of the correlation coefficients by channel after four iterations are complete. All have positive values and confirms that the non-linear stack positively correlates to the original synthetic signal for all channels. - The method shown in
FIG. 2 can automatically detect polarity changes in traces and reverse them. Once the polarity reversals are detected and corrected, weak microseimic events may be more effectively detected. As mentioned above, even though the example was generated using phase-weighted stack, any other non-linear stacks may also be used. - Once the polarity changes in traces are detected and corrected, the traces may be further processed using any desired processes. For example, the data may be processed for determining a parameter for a subterranean property. The data may be used for detecting or determining a microseismic event in a production operation. The microseismicity gives rise to an event, such as a fracturing operation; or a water, gas or steam flooding operation; may also be evaluated. The data with much improved SNR may also be used for creating an image of a subterranean structure.
-
FIG. 11 shows another seismic data processing method work flow in accordance with one or more implementations of various techniques described herein. This method may be used to detect polarity reversals along a seismic horizon under poor SNR conditions. This method may be used as an automatic fault plane solution (e.g. beach ball diagram) on surface arrays; or automatic polarity reversal correction to enable semblance-based processing. - The method begins with disposing a sensor network in a seismic survey geometry and dividing data from a surface sensor network into local groups at 1100, such that data may be stacked within the local groups to increase the signal to noise ratio. The particular group size may in some embodiments be a parameter set by the user, and may depend on the particular problem. In an embodiment, the particular group size may be selected in such a way, based on data signal polarities, that the data in each group has signal that can be stacked constructively, added, or combined in some manner. The method includes obtaining seismic survey data with the sensor network at 1101.
- The method further includes segmenting the seismic survey data a second time at 1102 into segments. Segmenting, here, refers to associating specific portions of the seismic data (or seismic gather) to one of the two polarities that exist. In other words, segmenting can refer to labeling each trace (or part of each trace around a seismic horizon) with a positive polarity or a negative polarity. There can be 2 or more segments depending on the nature of the data. In an embodiment, the segment size may be as small as one.
- The method further includes selecting a horizon among the obtained seismic survey data at 1103.
- Moveout can be defined as relating to the difference in the arrival times or travel times of a wave measured by receivers at two different offset locations. One form of moveout is caused by the separation between the source and the receiver. The method proceeds with correcting for move-out for the selected horizon in the seismic survey data at 1104. Various techniques of correcting for move-out are known, and any of such techniques could be employed at 1104.
- The method proceeds with calculating a non-linear stack of the data associated with the horizon, which has been corrected for move-out at 1106 for each local group. The method proceeds with locating a polarity reversal in the seismic survey data based on the non-linear stack (calculated in 1106 for each local group) at 1108, which may be performed using the LTA/STA (as described above) or any other detection method. The method further includes correcting the polarity of at least one segment of the seismic survey data according to the location of the polarity reversal at 1110.
-
FIG. 12A shows a plot of a gathered seismic signal of many seismic data traces, i.e., a synthetic signal plus real noise. We arbitrarily show only one trace in 5 from a large range of traces and display them at a convenient trace spacing to show the variation in traces, amplitude, and plotting offset (in meters) versus time (milliseconds). The signal to noise ratio is −18 dB.FIG. 12B shows the synthetic signal isolated from noise, with the polarity reversed along certain segments, plotting offset (in meters) versus time (in milliseconds). A first polarity reversal occurs at just before 900 m and a second polarity reversal occurs at roughly 1200 m, and a third at 1650 m, shown inFIG. 12B . A rectangle highlights the section containing several polarity changes, aligning with where the local groups of the sensor network are and with the manner in which the data are segmented into local groups and stacked within the local groups, as described above with reference toFIG. 11 .FIGS. 13A and B show magnified views of all the traces from the areas highlighted in a rectangle onFIGS. 12A and 12B in order to highlight a single polarity reversal. Whereas,FIG. 13A shows the synthetic signal with noise added, andFIG. 13B shows the synthetic signal isolated from noise. In turn,FIGS. 14A and 14B show a further magnified view of the area highlighted in the rectangle onFIG. 13 , showing just the one segment of data containing the polarity reversal at about 880 m.FIG. 14B shows the synthetic signal with noise added, andFIG. 14A shows the synthetic signal without the noise added. -
FIG. 15A shows a plot of ideal correlation coefficients for representative groups of traces for the signal ofFIG. 12A , whereasFIG. 15B shows a plot of the actual correlation coefficients computed for the same representative groups. The data is segmented into local groups (as described above), phase-weighted Nth root stack applied and the cross-correlation computed according to the method ofFIG. 11 . By not only detecting polarity changes, but correctly identifying wherein the data the polarity changes occur, the polarity can be handled so that the non-linear stack is cumulatively added, and polarity reversals do not mathematically cancel out portions of the data. -
FIG. 16 illustrates acomputing system 1600, into which implementations of various technologies described herein may be implemented. Thecomputing system 1600 may include one or more system computers (processors) 1630, which may be implemented as any conventional personal computer or server. However, those skilled in the art will appreciate that implementations of various technologies described herein may be practiced in other computer system configurations, including hypertext transfer protocol (HTTP) servers, hand-held devices, multiprocessor systems, microprocessor-based or programmable consumer electronics, network PCs, minicomputers, mainframe computers, and the like. - The
system computer 1630 may be in communication withdisk storage devices disk storage devices disk storage devices - In one implementation, seismic data from the receivers may be stored in
disk storage device 1631. Thesystem computer 1630 may retrieve the appropriate data from thedisk storage device 1631 to process seismic data according to program instructions that correspond to implementations of various technologies described herein. The program instructions may be written in a computer programming language, such as C++, Java and the like. The program instructions may be stored in a computer-readable medium, such as programdisk storage device 1633. Such computer-readable media may include computer storage media and communication media. Computer storage media may include volatile and non-volatile media, and removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules or other data. Computer storage media may further include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid state memory technology, CD-ROM, digital versatile disks (DVD), or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by thesystem computer 1630. - Communication media may embody computer-readable instructions, data structures, program modules or other data in a modulated data signal, such as a carrier wave or other transport mechanism. It may include any information delivery media. The term “modulated data signal” may mean a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal. By way of example, and not limitation, communication media may include wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above may also be included within the scope of computer readable media.
- In one implementation, the
system computer 1630 may present output primarily ontographics display 1627, or alternatively viaprinter 1628. Thesystem computer 1630 may store the results of the methods described above ondisk storage 1629, for later use and further analysis. Thekeyboard 1626 and the pointing device (e.g., a mouse, trackball, or the like) 1625 may be provided with thesystem computer 1630 to enable interactive operation. - The
system computer 1630 may be located at a data center remote from the survey region. Thesystem computer 1630 may be in communication with the receivers (either directly or via a recording unit, not shown), to receive signals indicative of the reflected seismic energy. These signals, after conventional formatting and other initial processing, may be stored by thesystem computer 1630 as digital data in thedisk storage 1631 for subsequent retrieval and processing in the manner described above. WhileFIG. 16 illustrates thedisk storage 1631 as directly connected to thesystem computer 1630, it is also contemplated that thedisk storage device 1631 may be accessible through a local area network or by remote access. Furthermore, whiledisk storage devices disk storage devices - Note that while the present disclosure describes the concepts as applied to a time domain representation of the signal, it is intended that the same principals apply if treating the signal in other forms including complex wavelets. Complex wavelets are a more general way to describe the signal, and encompasses both the complex time domain representation (wavelet is a “delta function” in time domain), frequency domain representation (wavelet is a “sine wave” in time domain) and other classes of “wavelets” that to a greater or lesser degree are compact in the time domain or frequency domain. In all cases the concept of the complex “wavelet” representation of the signal includes the concept of amplitude and phase.
- Additionally, the methods and systems of the present disclosure may be fully automated and able to operate continuously in time for monitoring, detecting, and locating seismic and microseismic events from on or multiple seismic detectors. Methods of the present invention utilize the contributions from one or multiple seismic detectors.
- While the disclosure has been disclosed with respect to a limited number of embodiments, those skilled in the art, having the benefit of this disclosure, will appreciate numerous modifications and variations therefrom. While the disclosure has been described in the context of applications in downhole tools, the apparatus of the disclosure can be used in many applications.
- Although a few example embodiments have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the example embodiments without materially departing from this disclosure. Accordingly, all such modifications are intended to be included within the scope of this disclosure as defined in the following claims. In the claims, means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not simply structural equivalents, but also equivalent structures. Thus, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wooden parts together, whereas a screw employs a helical surface, in the environment of fastening wooden parts, a nail and a screw may be equivalent structures. It is the express intention of the applicant not to invoke 35 U.S.C. §112,
paragraph 6 for any limitations of any of the claims herein, except for those in which the claim expressly uses the words ‘means for’ together with an associated function. - The preferred aspects and embodiments were chosen and described in order to best explain the principles of the invention and its practical application. The preceding description is intended to enable others skilled in the art to best utilize the invention in various aspects and embodiments and with various modifications as are suited to the particular use contemplated. The description may be implemented in any seismic or microseismic measurement system, and particularly for microseismic for hydraulic fracture monitoring. In addition, the methods may be programmed and saved as a set of instructions, that, when executed, perform the methods described herein. It is intended that the scope of the invention be defined by the following claims.
Claims (16)
1. A computer implemented method for seismic data processing by detecting one or more polarity reversals, comprising:
receiving a seismic data signal comprising a plurality of traces;
stacking the traces if the signal-to-noise ratio is below a level; and
applying a detection function to the stacked traces to detect polarity reversal in at least one of the plurality of traces based on the stacked traces.
2. The method according to claim 1 , further comprising:
iterating through each of the plurality of traces to detect polarity reversals.
3. The method according to claim 1 , further comprising:
reversing phase for the at least one of the plurality of traces in which polarity reversal is detected.
4. The method according to claim 1 , wherein detecting the polarity reversal further comprises defining a windowed trace by multiplying each trace of the seismic data signal with a temporal window function.
5. The method according to claim 4 , wherein detecting the polarity reversal further comprises correlating a stack and each windowed trace.
6. The method according to claim 4 , wherein detecting the polarity reversal further comprises defining the windowed trace based on a short term average long term average ratio (STA/LTA) (1204).
7. The method according to claim 4 , wherein detecting the polarity reversal further comprises:
correlating each windowed trace with one of (1) a subset of the other windowed traces in the seismic data set and (2) each of the remaining windowed traces in the seismic data set; and
identifying at least one trace with the largest correlation error containing a polarity reversal.
8. The method according to claim 4 , wherein detecting the polarity reversal further comprises:
correlating (1) the stack and (2) each windowed trace; and
identifying at least one trace with the largest correlation error containing a polarity reversal.
9. The method according to claim 1 , further comprising:
segmenting the seismic data signal into one or more segments according to positive polarity and negative polarity;
dividing the seismic survey data into a plurality of local groups;
for each of the local groups, calculating a non-linear stack from the seismic survey data;
for each of the local groups, locating a polarity reversal in the seismic survey data based on the non-linear stack for the local group; and
reversing the phase for the segment containing the polarity reversal in the seismic survey data.
10. The method according to claim 1 , further comprising
disposing a sensor network in a seismic survey geometry, each sensor in the network adapted to obtain one of the plurality of traces comprising the seismic data signal; and
for a given horizon in the seismic survey data, correcting for move-out.
11. The method according to claim 1 , wherein the stacking comprises nonlinear stacking.
12. The method according to claim 11 , wherein a non-linear stack produced by the non-linear stacking comprises one of a phase stack, a phase normalized stack, a double phase stack, a phase weighted Nth root stack, a double phase weighted Nth root stack, an Nth root stack, a phase weighted stack, a double phase weighted stack, a real phasor, a squared real phasor, a double phase weighted squared real-phasor stack or a combination thereof.
13. The method according to claim 11 , wherein the seismic data signal comprises a microseismic data signal.
14. The method according to claim 11 , further comprising at least one of:
determining a parameter for a subterranean property;
detecting a microseismic event;
evaluating an event giving rise to microseismic event;
evaluating a fracturing process;
evaluating a water, gas or steam flooding operation; and
obtaining an image of a subterranean structure.
15. A data processing apparatus processing the seismic data, comprising:
a processor; and
a data storage device communicating with the computer, wherein the data storage device contains seismic data and program instructions;
wherein program instructions causes the computer to perform a method as in claims 1 -9, 11-14 to detect polarity reversal.
16. A system for processing the seismic data, comprising:
a seismic acquisition network of seismic sensors configured to acquire a seismic survey comprising a plurality of data traces; and
a data processing apparatus of claim 15 .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/409,515 US20150177400A1 (en) | 2012-06-22 | 2013-06-21 | Detecting and correcting changes in signal polarity for seismic data processing |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261663449P | 2012-06-22 | 2012-06-22 | |
US201261663403P | 2012-06-22 | 2012-06-22 | |
US201261663416P | 2012-06-22 | 2012-06-22 | |
US14/409,515 US20150177400A1 (en) | 2012-06-22 | 2013-06-21 | Detecting and correcting changes in signal polarity for seismic data processing |
PCT/IB2013/055103 WO2013190512A2 (en) | 2012-06-22 | 2013-06-21 | Detecting and correcting changes in signal polarity for seismic data processing |
Publications (1)
Publication Number | Publication Date |
---|---|
US20150177400A1 true US20150177400A1 (en) | 2015-06-25 |
Family
ID=49769629
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/409,515 Abandoned US20150177400A1 (en) | 2012-06-22 | 2013-06-21 | Detecting and correcting changes in signal polarity for seismic data processing |
US14/397,196 Active 2035-10-21 US10641916B2 (en) | 2012-06-22 | 2013-06-21 | Processing seismic data by nonlinear stacking |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/397,196 Active 2035-10-21 US10641916B2 (en) | 2012-06-22 | 2013-06-21 | Processing seismic data by nonlinear stacking |
Country Status (3)
Country | Link |
---|---|
US (2) | US20150177400A1 (en) |
EP (2) | EP2864820A4 (en) |
WO (2) | WO2013190513A2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017048884A1 (en) * | 2015-09-16 | 2017-03-23 | Schlumberger Technology Corporation | Bayseian polarity determination |
CN107843924A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | Utilize the seismic source location of moving constraint at the beginning of P ripples and focal mechanism joint inversion method |
US10338245B2 (en) * | 2013-12-05 | 2019-07-02 | Cgg Services Sas | Methods and systems of detecting a microseismic event using an iterative non-linear inversion algorithm |
US11327188B2 (en) * | 2018-08-22 | 2022-05-10 | Saudi Arabian Oil Company | Robust arrival picking of seismic vibratory waves |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013190513A2 (en) | 2012-06-22 | 2013-12-27 | Schlumberger Technology B.V. | Processing seismic data by nonlinear stacking |
CN103995290B (en) * | 2014-06-03 | 2016-11-02 | 山东科技大学 | A kind of microseism P ripple seismic phase first arrival automatic pick method |
CN104122584B (en) * | 2014-08-08 | 2017-05-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for determining directionality according to seismic data |
CN104155688A (en) * | 2014-08-13 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | High precision weighted stack method |
WO2016120814A2 (en) * | 2015-01-28 | 2016-08-04 | University Of The Witwatersrand, Johannesburg | Method of determining the location and depth of magnetic sources in the earth |
CN106168679B (en) * | 2015-05-18 | 2018-03-09 | 中国石油化工股份有限公司 | Seismic acquisition records the processing method of polarity |
CN109477904B (en) | 2016-06-22 | 2020-04-21 | 休斯敦大学*** | Nonlinear signal comparison and high resolution measurement of seismic or acoustic dispersion |
US10295685B2 (en) | 2017-04-06 | 2019-05-21 | Saudi Arabian Oil Company | Generating common image gather using wave-field separation |
US11016212B2 (en) | 2017-04-11 | 2021-05-25 | Saudi Arabian Oil Company | Compressing seismic wavefields in three-dimensional reverse time migration |
US11275190B2 (en) | 2018-05-16 | 2022-03-15 | Saudi Arabian Oil Company | Generating diffraction images based on wave equations |
CN110967744A (en) * | 2018-09-28 | 2020-04-07 | 中国石油化工股份有限公司 | Oil-gas identification method based on biphase medium theory and storage medium |
US11681043B2 (en) * | 2019-09-03 | 2023-06-20 | Saudi Arabian Oil Company | Diffraction imaging using pseudo dip-angle gather |
US11313988B2 (en) | 2019-12-13 | 2022-04-26 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction imaging |
US11402529B2 (en) | 2020-01-09 | 2022-08-02 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging |
US11467303B2 (en) | 2020-03-09 | 2022-10-11 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using a post-stack seismic diffraction imaging condition |
US11320557B2 (en) | 2020-03-30 | 2022-05-03 | Saudi Arabian Oil Company | Post-stack time domain image with broadened spectrum |
US11656378B2 (en) | 2020-06-08 | 2023-05-23 | Saudi Arabian Oil Company | Seismic imaging by visco-acoustic reverse time migration |
US11513500B2 (en) * | 2020-10-09 | 2022-11-29 | Halliburton Energy Services, Inc. | Method for equipment control |
US11899147B2 (en) | 2021-08-12 | 2024-02-13 | Saudi Arabian Oil Company | Method and system for seismic denoising using omnifocal reformation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4516206A (en) * | 1982-10-21 | 1985-05-07 | Mcevilly Thomas V | Post-processing of seismic parameter data based on valid seismic event determination |
US20110002194A1 (en) * | 2008-05-22 | 2011-01-06 | Matthias Imhof | Method For Geophysical and Geological Interpretation of Seismic Volumes In The Domains of Depth, Time, and Age |
Family Cites Families (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4740928A (en) * | 1984-10-09 | 1988-04-26 | Amoco Corporation | Sonic logging system |
US5197039A (en) | 1988-03-29 | 1993-03-23 | Shell Oil Company | Methods for processing seismic data |
US5073875A (en) * | 1990-11-15 | 1991-12-17 | Amoco Corporation | Method of enhancing geophysical data |
US6061298A (en) * | 1998-06-08 | 2000-05-09 | Pgs Tensor, Inc. | Method of and system for processing multicomponent seismic data |
GB2348003B (en) | 1999-03-19 | 2001-02-07 | Geco Prakla | Seismic data processing method for data acquired using overlapping vibratory sweeps |
US6363327B1 (en) * | 2000-05-02 | 2002-03-26 | Chroma Graphics, Inc. | Method and apparatus for extracting selected feature information and classifying heterogeneous regions of N-dimensional spatial data |
US6678207B2 (en) | 2001-05-25 | 2004-01-13 | Exxonmobil Upstream Research Company | Trapped water bottom multiple and peg-leg multiple suppression for ocean bottom seismic data |
KR100412097B1 (en) * | 2001-07-26 | 2003-12-24 | 주식회사 산하이앤씨 | Method, system and apparatus for processing seismic data using the time-varying optimum offset concept |
US6748330B2 (en) | 2002-04-10 | 2004-06-08 | Schlumberger Technology Corporation | Method and apparatus for anisotropic vector plane wave decomposition for 3D vertical seismic profile data |
WO2005098731A2 (en) | 2004-03-29 | 2005-10-20 | German Peter T | Systems and methods to determine elastic properties of materials |
US7042214B2 (en) | 2004-04-13 | 2006-05-09 | The Board Of Trustees Of The Leland Stanford Junior University | Non-linear symmetric sweep spectral-spatial RF pulses for MR spectroscopy |
US7525873B1 (en) * | 2004-07-23 | 2009-04-28 | Bf Corporation North America Inc. | Seismic inversion of conditioned amplitude spectra |
US7391675B2 (en) * | 2004-09-17 | 2008-06-24 | Schlumberger Technology Corporation | Microseismic event detection and location by continuous map migration |
EP1815272B1 (en) | 2004-11-24 | 2021-07-28 | Paradigm Geophysical Ltd. | System and method for fault identification |
US8368399B2 (en) * | 2005-05-05 | 2013-02-05 | New York University | Mode-scanning excitation magnetic resonance imaging method and system |
WO2008053436A1 (en) * | 2006-10-31 | 2008-05-08 | Philips Intellectual Property & Standards Gmbh | Mri rf encoding using multiple transmit coils |
WO2008076191A2 (en) | 2006-12-15 | 2008-06-26 | Exxonmobil Upstream Research Company | Identification and suppression of multiples in ocean bottom seismic data |
US8116167B2 (en) | 2008-06-12 | 2012-02-14 | Chevron U.S.A. Inc. | Method and system for generating a beam of acoustic energy from a borehole, and applications thereof |
US9354343B2 (en) | 2009-03-09 | 2016-05-31 | Ion Geophysical Corporation | Declination compensation for seismic survey |
US7978563B2 (en) | 2009-08-18 | 2011-07-12 | Microseismic, Inc. | Method for passive seismic emission tomography including polarization correction for source mechanism |
FR2960304B1 (en) | 2010-05-19 | 2012-09-14 | Cggveritas Services Sa | PASSIVE MONITORING METHOD OF SEISMIC EVENTS |
FR2961316A1 (en) * | 2010-06-10 | 2011-12-16 | Cggveritas Services Sa | PROCESS FOR PROCESSING SEISMIC MARINE DATA |
WO2013190513A2 (en) | 2012-06-22 | 2013-12-27 | Schlumberger Technology B.V. | Processing seismic data by nonlinear stacking |
-
2013
- 2013-06-21 WO PCT/IB2013/055104 patent/WO2013190513A2/en active Application Filing
- 2013-06-21 EP EP13807168.3A patent/EP2864820A4/en not_active Withdrawn
- 2013-06-21 US US14/409,515 patent/US20150177400A1/en not_active Abandoned
- 2013-06-21 WO PCT/IB2013/055103 patent/WO2013190512A2/en active Application Filing
- 2013-06-21 US US14/397,196 patent/US10641916B2/en active Active
- 2013-06-21 EP EP13806301.1A patent/EP2864819A4/en not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4516206A (en) * | 1982-10-21 | 1985-05-07 | Mcevilly Thomas V | Post-processing of seismic parameter data based on valid seismic event determination |
US20110002194A1 (en) * | 2008-05-22 | 2011-01-06 | Matthias Imhof | Method For Geophysical and Geological Interpretation of Seismic Volumes In The Domains of Depth, Time, and Age |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10338245B2 (en) * | 2013-12-05 | 2019-07-02 | Cgg Services Sas | Methods and systems of detecting a microseismic event using an iterative non-linear inversion algorithm |
WO2017048884A1 (en) * | 2015-09-16 | 2017-03-23 | Schlumberger Technology Corporation | Bayseian polarity determination |
EP3350626A4 (en) * | 2015-09-16 | 2019-08-21 | Services Petroliers Schlumberger | Bayseian polarity determination |
US11035967B2 (en) | 2015-09-16 | 2021-06-15 | Schlumberger Technology Corporation | Bayesian polarity determination |
CN107843924A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | Utilize the seismic source location of moving constraint at the beginning of P ripples and focal mechanism joint inversion method |
US11327188B2 (en) * | 2018-08-22 | 2022-05-10 | Saudi Arabian Oil Company | Robust arrival picking of seismic vibratory waves |
Also Published As
Publication number | Publication date |
---|---|
WO2013190512A2 (en) | 2013-12-27 |
WO2013190512A3 (en) | 2014-04-10 |
EP2864820A2 (en) | 2015-04-29 |
EP2864819A2 (en) | 2015-04-29 |
EP2864819A4 (en) | 2015-06-24 |
US10641916B2 (en) | 2020-05-05 |
US20150112601A1 (en) | 2015-04-23 |
EP2864820A4 (en) | 2015-06-17 |
WO2013190513A3 (en) | 2014-04-10 |
WO2013190513A2 (en) | 2013-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20150177400A1 (en) | Detecting and correcting changes in signal polarity for seismic data processing | |
Daley et al. | Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle, Alabama | |
Behm et al. | Retrieval of local surface wave velocities from traffic noise–an example from the La Barge basin (Wyoming) | |
RU2457513C2 (en) | Methods and systems for processing microseismic data | |
Eisner et al. | Uncertainties in passive seismic monitoring | |
Chadwick et al. | Quantitative analysis of time-lapse seismic monitoring data at the Sleipner CO 2 storage operation | |
US10324211B2 (en) | Seismic spectral balancing | |
Binder et al. | Modeling the seismic response of individual hydraulic fracturing stages observed in a time-lapse distributed acoustic sensing vertical seismic profiling survey | |
US20180203144A1 (en) | Interferometric Microseismic Imaging Methods and Apparatus | |
US20130215716A1 (en) | Integrated Passive and Active Seismic Surveying Using Multiple Arrays | |
CA2818790C (en) | Seismic trace attribute | |
US9234977B2 (en) | Processing collected survey data | |
Lambert et al. | Low‐frequency microtremor anomalies at an oil and gas field in Voitsdorf, Austria | |
US20190094397A1 (en) | Surface detection and location of microseismic events and earthquakes without the use of a velocity model | |
US20160334528A1 (en) | Systems and methods for characterizing subterranean formations utilizing azimuthal data | |
RU2737846C2 (en) | System for installing ground-based seismic sensors with pairs of adjacent multicomponent seismic sensors at an average distance of at least twenty meters | |
Titov et al. | Modeling and interpretation of scattered waves in interstage distributed acoustic sensing vertical seismic profiling survey | |
US20210208298A1 (en) | Measurement of in situ rock formation properties using surface seismic sources and downhole receivers | |
CN106383360A (en) | Seismic data based submarine sediment temperature-pressure field recognition method | |
Behm | Feasibility of borehole ambient noise interferometry for permanent reservoir monitoring | |
Freudenreich et al. | Microseismic feasibility studies–assessing the probability of success of monitoring projects | |
Raub et al. | Seismic‐wave propagation in shallow layers at the GONAF‐Tuzla Site, Istanbul, Turkey | |
Lebedev et al. | Accuracy in Estimating the Parameters of a Layered Medium Using Coherent Vector Reception of a Rayleigh Surface Wave | |
Mizuno et al. | Real-time passive seismic monitoring using DAS—Today's solutions and remaining challenges | |
Wilczynski et al. | Drill‐bit position monitoring using seismic‐while‐drilling data; numerical and field examples from Sweden |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SCHLUMBERGER TECHNOLOGY CORPORATION, TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:OZBEK, ALI;PROBERT, ANTHONY;RAYMER, DANIEL GORDON;SIGNING DATES FROM 20150115 TO 20150119;REEL/FRAME:035569/0432 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |