WO2010080366A1 - Hydrocarbon detection with passive seismic data - Google Patents
Hydrocarbon detection with passive seismic data Download PDFInfo
- Publication number
- WO2010080366A1 WO2010080366A1 PCT/US2009/067913 US2009067913W WO2010080366A1 WO 2010080366 A1 WO2010080366 A1 WO 2010080366A1 US 2009067913 W US2009067913 W US 2009067913W WO 2010080366 A1 WO2010080366 A1 WO 2010080366A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- data
- earthquake
- passive
- source
- subsurface region
- Prior art date
Links
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. analysis, for interpretation, for correction
-
- 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. analysis, for interpretation, for correction
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- 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
Definitions
- This invention pertains generally to hydrocarbon exploration and, more particularly, to inferring the presence or absence of hydrocarbons in the subsurface.
- the invention is a method for using seismic data generated by passive seismic sources such as earthquake activity in conjunction with other geophysical data for petroleum exploration.
- the invention is a method for hydrocarbon detection in a subsurface region using passive-source seismic data in conjunction with at least one other type of geophysical data, comprising: (a) obtaining passive-source seismic survey data for the subsurface region, wherein the survey's receivers were seismometers located generally over said subsurface region and spaced for hydrocarbon prospecting, said seismometers being suitable for global seismology earthquake detection, and wherein said passive-source seismic data include at least one event identified in an earthquake table and estimated to have a dominant frequency at the subsurface region in a range of 0 to 8 Hz based on said at least one event's magnitude and distance from the subsurface region;
- step (d) using the physical properties model to predict hydrocarbon potential for the subsurface region.
- a computer i.e. a data processing machine, programmed in accordance with the invention description that follows.
- Fig. 1 is a schematic plot illustrating the frequency content of different types of data
- Fig. 2 is a flow chart showing the use of passive seismic and other geophysical data to derive subsurface rock physical property models for hydrocarbon exploration;
- Fig. 3 shows a flow chart of the joint-inversion implementation of the invention in which a traveltime tomography analysis of the passive seismic data is combined with information from other type(s) of geophysical data;
- Fig. 4 shows a flow chart of the joint-inversion implementation of the invention in which reflection seismic data generated from the passive seismic data by cross- correlation are inverted jointly other type(s) of geophysical data;
- Fig. 5 is a schematic illustration of the method of synthesizing reflection seismic data from the passive seismic data using the cross-correlation technique
- Fig. 6 is a schematic diagram showing a problem layout for joint inversion
- Fig. 7 is a schematic diagram illustrating how joint inversion computational grids may be updated using gradients
- Fig. 8 shows a synthetic receiver function waveform as result of deconvolving the vertical component seismogram from a horizontal component seismogram for the simple velocity model shown in the inset;
- Figs. 9A-D show schematic diagrams of typical seismic phases expected on receiver function waveforms
- Fig. 10 is a flow chart showing basic steps in one embodiment of the present inventive method.
- Fig. 11 is a graph showing computed ray coverage from regional and teleseismic events at a producing field in the western United States;
- Fig. 12 is a plot of earthquake magnitude versus earthquake frequency content expressed by the dominant frequency;
- Fig. 13 is a map view display of travel time anomalies relative to a ID earth model
- Fig. 14 shows receiver function results that reveal reflectors in a shallow crustal section where incoming P-waves converted to S-waves.
- FIG. 1 shows approximately where the various types of geophysical data fit into the picture.
- the solid line curve 11 represents active seismic data.
- the approximate useful frequency range for controlled source electromagnetic (“CSEM”) data is represented by curve 12, and curve 13 shows the frequency range for gravity and NMO (normal moveout) velocity data.
- CSEM controlled source electromagnetic
- the present invention uses passive source seismic data, primarily earthquake data, to fill this gap, and teaches how to exploit these data.
- Passive seismic data have potential to fill the gap because such data can be rich in low frequency information.
- Earthquakes including tremors, volcanos, and disturbances of all magnitudes, and tidal waves are examples of passive seismic sources, which can be defined as seismic signals generated by natural sources or by man-made sources not employed purposefully in a seismic data acquisition program.
- Man-made disturbances may be present as noise because it may not be possible to entirely exclude incidental events such as road construction reverberations from the desired earthquake data.
- the curve 14 in Fig. 1 indicates the approximate frequency range that earthquakes can provide. Curve 14 is obviously very advantageously positioned relative to gap 10. It is particularly obvious from Fig. 1 how passive-source seismic data complements active source seismic data in terms of frequency bandwidth.
- CSEM magneto-telluric or gravity
- the earthquakes of most interest to this invention are large events (generally larger than magnitude 3.0) that are likely to occur far away from the region of interest (hundreds to thousands of kilometers away).
- Figure 12 shows a calculation of earthquake frequency content versus earthquake magnitude. The larger the earthquake, the more low- frequency energy that tends to be emitted. Moreover, the earth acts as a filter, which means that the further away an earthquake occurs, the more strongly the high-frequencies will be attentuated. Consequently, large earthquakes that occur far away will arrive at a recording station with lower frequency content, i.e. a lower dominant frequency, than earthquakes of a similar or lesser magnitude that occur nearer to the recording station.
- Earthquakes greater than magnitude 3.0 or 4.0 will likely contain energy within the desired frequency range for this invention, i.e, within the frequency gap 10, although the stress drop is what ultimately determines the amount of energy emanated.
- the earthquakes required in this invention to provide low frequency information differ from micro-earthquakes which are used to monitor gas and oil production. Micro-earthquakes are due to slip on tiny little faults and therefore have a typical frequency content of 10's to 100's of Hz, which can possibly contribute to a better result but cannot help to fill in the low- frequency gap indicated in Fig. 1. [0026] Earthquakes provide lower frequency vibrations than human operated seismic sources because of the magnitude of the disturbance, which increases the low frequency content.
- a magnitude 1.0 earthquake is roughly equivalent to the energy release from an explosion of about 70 pounds of TNT (a mid-sized construction site blast), which is larger than most explosions set off for seismic survey purposes.
- a magnitude 2.0 earthquake is similar to an explosion of 1 metric ton of TNT
- a magnitude 4.0 earthquake is approximately equivalent to the energy release from an explosion of 1,000 tons of TNT (a small nuclear blast).
- Other examples include the 1980 Mount St. Helens eruption, which released the energy equivalent of a magnitude 7.8 earthquake (or just over 500 megatons of TNT), and the 1883 Krakatoa eruption, which released the energy equivalent of a magnitude 8.5 earthquake (or about 5.6 gigatons of TNT). More information about natural sources can be obtained from the U.S.
- Passive seismic data have been used for hydrocarbon exploration and reservoir monitoring and characterization, by techniques that include seismic trace cross-correlation or interferometric techniques, which convert transmission seismic data to reflection seismic data; see for example Claerbout, "Synthesis of a Layered Medium From its Acoustic Transmission Response," Geophys. 33, 264 (1968); and Wapenaar, "Retrieving the Elastodynamic Green's Function of an Arbitrary Inhomogeneous Medium by Cross Correlation,” Phys. Rev. Lett. 93, 254301 (2004).
- microearthquake or micro-seismic or micro-tremor
- the term microearthquake will refer to low magnitude, naturally occurring seismic events originating below the survey seismometers, i.e. in or below the survey target region.
- Micro-seismicity is generally taken to be Mw ⁇ 0. (Earthquake magnitudes are expressed on a logarithmic scale, and thus very small earthquakes have negative magnitudes.)
- a hydrocarbon exploration case study describes 40 stations covering an area of 3000 km 2 that recorded for 10 months (i.e., station spacing of at least 10 km).
- a total of 900 earthquakes were recorded of which half (450) were used for tomographic inversion because those 450 events were (1) located within the seismic network and (2) had at least 20 P- and S-wave arrivals and (3) had estimated location errors of less than 1 km.
- the inversion grid was 2 km by 4 km.
- a second case study involved recording 200 micro-earthquakes with magnitudes between -0.5 and 3.0, locating the hypocenters, and determining the velocity structure via tomographic inversion. Station spacing was 500 m in this example.
- passive seismic data can be used as (1) a reconnaissance tool in order to optimize the cost of 3D conventional seismic surveys and (2) as a complementary method for reinterpretation of 2D regional seismic data or reprocessing of 2D seismic data by using the velocity model derived from the passive seismic data.
- Tselentis et al. describe a micro-earthquake network of 70 stations at 500 m spacing that recorded 220 local micro-earthquake events. These events were used to determine a ID velocity model. That model was then used as the initial model for a nonlinear inversion of
- the present invention in some of its embodiments is a method for using passive source earthquake data in conjunction with at least one other geophysical datatype for hydrocarbon exploration.
- the passive data must be acquired with a dense survey design, i.e. receiver spacing suitable for hydrocarbon exploration.
- Receiver spacing for hydrocarbon exploration is typically less than 1 km, preferably less than 250 m.
- Techniques developed within the earthquake seismology community are modified for exploring targets that are smaller than such large scale targets as the Earth's crust.
- the present invention derives earth rock property model(s) by jointly using passive seismic data together with at least one other type of non-seismic geophysical data such as active source seismic reflection data, controlled source electromagnetic (CSEM), magnetotelluric (MT), magnetic and gravity or gravity gradiometry data.
- CSEM controlled source electromagnetic
- MT magnetotelluric
- the present invention provides a way to close the low frequency gap or lacuna (between 0 and ⁇ 5 to 8 Hz) which plagues active source seismic techniques currently available and to improve joint inversion of multiple geophysical data.
- Figure 11 shows that contrary to traditional belief, earthquake data hold promise for imaging subsurface structure on a petroleum exploration scale, because they do not travel nearly as vertically as has been assumed and have better coverage than had been assumed, particularly if such data are collected for an extended period, e.g., six months.
- the inventors deployed an array of seismometers (-250 m spacing) over a producing field in the western United States for about 6 months time and recorded seismic data continuously at a rate of 100 samples per second. Preliminary analysis of these data find that relative travel time measurements can indeed be made (Fig. 13) and that the measured relative delays (or advances) are likely related to shallow velocity structure.
- 131 indicates travel times that are late relative to a ID earth model, and 132 indicates early travel times. Receiver function analysis was performed on some of these data.
- a first step (step 21 in the flow chart of
- Fig. 2 includes a passive seismic survey specifically designed for hydrocarbon exploration.
- This design preferably includes a station spacing much smaller (e.g., ⁇ 1 km) than typically used in earthquake seismology work (generally 10's of km).
- the smaller spacing and continuous or triggered recording over a long period of time (step 22) is intended to allow resolution of shallow targets.
- the present inventive method can make full use of this dataset, including earthquakes of all distances (local, regional, teleseismic) all source (occuring naturally or man-made sources not employed purposefully) and wave (body and surface waves) types, and all magnitudes (including microseismic scale).
- Standard techniques used in earthquake seismological work can then be employed (step 23) such as but not limited to tomographic approaches, receiver functions, or the interferometric techniques, modified where necessary to target a smaller scale appropriate for hydrocarbon exploration.
- step 22 of the present inventive method it may be preferable in performing step 22 of the present inventive method to conduct continuous recording over a period of time, such as several months up to years, in order to gather a sufficient number of "sources".
- all events that are detectable by the survey seismographs are recorded, from microseismicity (very small earthquakes in the surrounding area of the receiver) to larger earthquakes thousands of kilometers away.
- There is an "indistinct" categorization in earthquake reporting where earthquakes are sorted by their source-receiver distance: local (within the vicinity of the recording station), regional (up to 3000 km distance between source and receiver), and teleseismic (30-90 degrees distance). In the latter case, the seismic rays traverse through the Earth's mantle.
- Table 1 shows information published by the USGS, which demonstrates the frequency of occurrence of earthquakes (http ://neic .us gs . go v/neis/eqlists/eqstats . html) .
- Earthquakes are defined by their origin times and hypocenters (i.e., x,y,z location of the event). Catalogs of earthquakes occurring anywhere in the world with a magnitude Mw > 5.5 are maintained by scientists funded by the National Science Foundation and made publicly available (published annually in Physics of the Earth and Planetary Interiors and posted on the internet). Published papers describing the Centroid Moment Tensor (CMT) methodology for locating and describing earthquakes include: Dziewonski and Woodhouse, "An Experiment in the Systematic Study of Global Seismicity: Centroid Moment Tensor Solutions for 201 Moderate and Large Earthquakes of 1981," J. Geophysi.
- Permanent arrays of global seismometers are used to locate earthquakes occurring all over the world.
- the procedure for locating an earthquake is to observe the times at which the wavefront passes each seismic station.
- the location, depth and origin time of the earthquake are inferred typically via an inversion of the observed travel times.
- the process usually begins with an initial guess of the hypocenter, and then the hypocenter location is iteratively updated until a final location is found that best predicts the observed arrival times at the stations in the seismic network.
- An example of data from an earthquake catalog showing magnitudes, times and locations of a series of earthquakes is given in Table 2.
- earthquake catalogs provide more information on the fault that ruptured, for example its length, angle, and direction of slip.
- seismometer data from step 22 is examined, and likely earthquake events are identified. Time segments from the seismometer records corresponding to these selected seismic events are collected, and in some embodiments of the invention, the remainder of the passive data, which would include microearthquake data, is ignored. An earthquake catalog is then searched to identify the specific earthquake corresponding to each selected data segment.
- Source information particularly geographical location, depth, and time, may be obtained from the earthquake catalog or similar source as explained above.
- standard source location techniques can be used to locate or relocate events that occur by studying the arrival times at permanent seismic stations that continuously record seismic data around the world. Source location is not required if the interferometric and/or other ambient noise techniques are used to process the passive seismic data.
- the passive data are processed simultaneously (step 26) with other geophysical data (such as active seismic reflection data, gravity, magnetic, CSEM, MT. etc.) to achieve a better answer than would be obtained by inverting any of the individual geophysical data types alone.
- Simultaneous processing means deriving information from two or more sets of geophysical data by setting up a system of equations and solving them either as a single, mathematical optimization problem or by sequentially solving each equation by alternating which parameters are fixed and which are solved for.
- the simultaneous processing of step 26 may be joint inversion.
- Inversion as defined herein means an automated procedure for obtaining an earth model that is consistent with the measured data and is typically implemented as an optimization problem in which the mismatch between measured and forward synthesized data is minimized.
- manual forward-modeling iterative approaches the user updates the model to reduce the mismatch between predicted and measured data.
- Such manual forward-modeling, performed simultaneously on passive seismic data plus at least one other data type is an alternative way to perform step 26.
- Inversion techniques for individual data types have been used by the oil industry for many years. They are commonly used to derive subsurface rock properties from active seismic reflection data (Jin, et.
- Joint inversion of multiple data is a relatively new concept in the oil industry, but common in global seismology where different types of seismic waves are often analyzed together (i.e., body waves and surface waves).
- the data misfits for multiple data types are simultaneously optimized using hypothesized relations between the relevant physical properties such as P-wave velocity and conductivity. Any missing information content in one geophysical data used by the inversion can potentially be filled with information from a distinct geophysical data type in order to reduce the uncertainty in the rock property models that are predicted.
- Standard methods for integrating multiple geophysical data typically rely on human-directed interpretation processes, in which a person develops a hypothesis of the earth structure based on one of the data types (generally the seismic reflection data) and then determines whether structures inferred from the other geophysical data (e.g., gravity, CSEM, MT, etc.) are consistent with that initial guess.
- a manual iteration between the various data types ensues until a compromise of some sort has been reached. This sort of approach is extremely time -intensive and becomes increasingly difficult to carry out as additional data types are considered.
- Strack et al. describe schematically a type of joint inversion of time domain controlled source electromagnetic data with (active-source) seismic data in U.S. Patent No. 7,328,107. No examples or test results are presented.
- Patent Application Publication US 2007/0294036 Strack and Allegar disclose a method for acquiring and interpreting electroseismic and seismoelectric data in which joint inversion of the two data types is proposed. Joint inversion with electromagnetic data is also proposed. No test results are presented.
- the passive data were comprised solely of small, local earthquakes near to the area of interest, the earthquake times and locations were unknown. Therefore, they had to be inferred during the inversion process which introduced velocity errors because of trade -off between the velocities and source locations. Because the locations and times of the active data were known, the addition of those traveltimes to the passive traveltime inversion helped reduce this trade-off and therefore improved the crustal velocity estimates. In contrast, the present inventive method uses larger earthquakes whose locations are known (contained in earthquake catalogs) so that the trade-off referred to above is not an issue. In the present invention, the joint inversion of active and passive seismic datasets is used for an entirely different purpose, i.e. to help reduce the frequency gap that plagues active seismic inversions, as reflected by the requirement for earthquakes that fill the 0 - 8 hz frequency range.
- the passive seismic data can be used jointly with other data types (step 26) in a number of ways.
- Fig. 3 shows a flow chart in which the travel time information from the passive seismic data is used in a joint-inversion implementation of the invention through travel time tomography.
- the passive seismic data are processed.
- desired arrivals are identified.
- the arrival times for the identified arrivals are picked.
- the initial subsurface models including the seismic velocity model(s), and models for other parameters (e.g. conductivity) whose values affect the remaining data types, are built.
- other geophysical data used in the joint inversion are predicted using the initial models.
- the locations and times of occurrence for the sources in the passive seismic data are estimated using the initial seismic velocity model(s).
- the arrival travel times of those sources in step 306 are predicted based on the current seismic velocity model(s).
- the total data misfit between the predicted and the measured geophysical data is calculated.
- the data misfit for the passive seismic is a measure of the travel time difference between the arrival times computed based on the current velocity model(s) and the picked travel times from the passive seismic data.
- the total data misfit is compared against a preset value. If the total data misfit is less than the preset value (or other stopping condition is reached), the inversion process is stopped, and the current subsurface model(s) is the final result of the inversion.
- the subsurface models are updated at step 310, and the inversion process returns to steps 305 and 306 for another cycle.
- This implementation of joint inversion using passive data is preferred when the passive seismic signals come from earthquakes, in which the desired seismic mode arrivals can be picked.
- Fig. 4 is a flow chart showing another implementation of the invention involving joint inversion using passive seismic data in step 26.
- cross-correlation see Claerbout (1968) op. cit. or Wapenaar (2004) op. cit.
- the passive seismic data are processed.
- the cross-correlation technique is applied to the passive seismic data to generate quasi-shot gathers. Either passive noise (Artman, op. cit.
- This method converts the passive data into a number of "quasi-shot gathers," where each receiver location becomes an effective "shot.”
- the initial subsurface models, including the seismic velocity model(s) are built.
- the geophysical data other than the passive seismic data to be used in the joint inversion are predicted using the initial models.
- passive seismic data corresponding to the quasi-shot gathers from step 402 are predicted using the current velocity model(s).
- the total data misfit between the predicted and the measured geophysical data is calculated.
- the total data misfit is compared against a preset value.
- the inversion process is stopped, and the current subsurface models are the final result of the inversion. If the total data misfit is greater than the preset value, the subsurface models are updated at step 408, and the inversion process repeats beginning with steps 404 and 405.
- a joint inversion technique that may be used in preferred embodiments of the present invention as part of step 26 will be described in more detail next. However, it is important to understand that any joint inversion technique may be used in the present invention to invert earthquake data simultaneously with one or more other types of geophysical data as in step 26.
- the advantages of the joint inversion technique described herein are not limited to instances where earthquake data are one type of data being jointly inverted, although that application is central to this patent application.
- earthquake data are one type of data being jointly inverted
- a fundamental aspect of the present joint inversion technique is the recognition that rocks are systems of grains and pores (dry or fluid filled) contained within a matrix. This system links the physical properties of those rocks (e.g., velocity, resistivity, and density) to their structural composite properties (e.g., porosity, lithology, and water/hydrocarbon saturation).
- This linkage provides a means to complement information missing from one set of measurements (i.e., active source seismic) with information from another set of measurements (i.e., in particular the CSEM/MT/Gravity). It also provides a means for constructing the medium averages for grid transformations that are required in efficient finite difference computations.
- Figures 6 and 7 are schematic diagrams that summarize the following discussion of an inventive joint inversion engine.
- Equation (1) an error functional such as described by Equation (1) can be defined for each data type separately, followed by a multi-objective formulation of the joint inversion problem.
- Multi-objective optimization tools provide the capability to solve problems where several objective functions are optimized simultaneously (Steuer, “Multiple Criteria Optimization, Theory, Computations and Applications,” John Wiley & Sons (1986); Das, et al., “Normal-Boundary Intersection: A New Method for Generating the Pareto Surface in Nonlinear Multicriteria Optimization Problems," SIAM Journal on Optimization 8, 631-657 (1998)).
- the total joint error functional can be defined as a vector functional:
- ⁇ mal [ ⁇ x (m) W 1 (Jn) ... Wn (m)Y (2) where n is the number of data types involved in the joint inversion formulation.
- Possible approaches to solve such a problem include minimizing a weighted linear or nonlinear combination of the components of the vector function, or minimizing any one component and constraining the remaining to be within certain limits.
- the latter approach makes the inverse problem a constrained optimization problem, where the constraints define a set or sets of combinations of model parameters that are allowed in the inversion search space. Other constraints on the model parameters can be imposed to ensure that the inversion results are physically meaningful, or could be mathematical in nature - to make the inversion numerically stable.
- the model parameters are the lithology (e.g., shale volume content in a clastic rock or dolomite or Calcium Carbonate content in a carbonate rock), the water saturation of the pore fluid (Sw) , the porosity ( ⁇ ), and a second lithology parameter if required (Vuth), respectively, are defined on a fine scale grid.
- Rock physics equations are used to transform the model properties into the physical properties of conductivity ( ⁇ ) , elastic moduli (C ⁇ ) and density (p) , which are then used to solve Maxwell's equations, the elastic (or acoustic) wave equation or the Kirchhoff equation as a high-frequency approximation, and the gravitational potential equation.
- It includes the gradient in the data fitting error and the smoothing part for one of the physical properties ( ⁇ ,C kl ,or p) , and transforms the physical property from a macro- to a dmicro dmicro micro-grid, and computes the derivative with respect to one of the rock properties drock prop (Vsh, Sw, ⁇ , or Vhth) using an appropriate rock physics model such as the Xu- White model (Xu and White, "A Physical Model for Shear-Wave Velocity Prediction," Geophysical Prospecting 44, 687-717 (1996).
- the quality of the starting model may be improved by jointly inverting
- the Earthquake/CSEM/MT/Gravity information along with the low-frequency seismic travel time information (NMO velocities).
- the Earthquake/CSEM/MT/Gravity inversion produces macro-model P-wave and S-wave velocity volumes that rely on more than just the travel time information, but instead include information that is extracted from the non-seismic data (CSEM/MT/gravity) via the rock physics relations.
- various optimization techniques such as differential semblance optimization may be used.
- the active source seismic forward calculation may be performed as follows.
- the travel-times are computed via ray-tracing through the macro-model interval velocity model.
- the offset-dependent reflectivity of the pre-stack seismic gathers are inverted to find the elastic parameter variations using a convolutional or other appropriate approximation.
- the DSO penalty is applied in order to find the velocity model that flattens all data in the gathers.
- the wavelet scaling parameter is updated (if there are multiples in the data that can be used to constrain this information). This set of computations form an inner loop within the larger framework of the joint inversion that can be repeated many times until convergence is attained.
- the active source seismic piece has been separated into two pieces, i.e. the separation of dependence of seismic data non-linearly on the macromodel velocity information which controls traveltime information and dependence of the seismic data more or less linearly on the variations in elastic properties which control reflection amplitudes. Because of this separation, there are two separate error gradients to be computed and propagated to the rock physics grid.
- the first set relates to the macro-model Vp (or both Vp and Vs) parameter. It is derived from the travel-time part of the computation and is similar to the error gradients computed in equation 3 for the CSEM/MT/Gravity data fits.
- the second set of error gradients is with respect to the AVO piece:
- simulation grid for the AVO piece is on the same fine scale as the inversion grid, no up-scaling/down-scaling between these two is required.
- the amplitude information from the seismic is separated from the travel time piece of the seismic.
- the amplitude information is not included in the joint inversion objective function until after the low-frequency parts of the model have converged.
- the implied velocity model is fixed (for raytracing) and just the amplitude information from the seismic is inverted with the other datatypes (e.g., gravity, CSEM, MT, earthquake) until the solution for that modified set of datatypes has converged.
- step 21 involves design of a passive seismic monitoring survey, which is designed to fit the scale of hydrocarbon exploration.
- a survey of the present invention targets the hydrocarbon exploration case, and the survey design will therefore not only take the shallower target into account, but also the sources and its associated azimuths and distances that can be expected.
- the receiver spacing will be that suitable for hydrocarbon prospecting, i.e. the same as might be used in active source seismic hydrocarbon surveys, the receivers must be earthquake seismometers or similar acoustic detectors with frequency response extending down to about 1 Hz instead of the typical hydrophone or geophone used in active source seismic work which have frequency response down to only 5-6 Hz.
- the types of seismometer(s) used for data acquisition are selected according to the frequency ranges of the passive seismic signals that need to be recorded.
- the frequency range is related to the source-receiver distance, i.e. generally events of larger distance between the event and the recording station have lower frequency content.
- Earthquakes occurring at teleseismic distances contain the lowest frequency content, on the order of 1- 3 Hz, regional events (distances less than 30° or -3000 km), and local earthquakes (within 100 km of the recording station) produce signals of highest frequency content (since the high frequencies are not filtered out by a long propagation path), which can overlap with the active seismic frequency band.
- the receiver surface coverage area and spacing is preferably designed according to the frequency range and subsurface area to be covered. For example, receivers of low-frequency response may be placed more sparsely than those of high-frequency response.
- the required receiver spacing may be determined as a function of the maximum desired frequency, maximum dip of the arriving wavefront, and subsurface velocity, in a manner that is known to those familiar with the technical field of the active-source seismic data acquisition. For example, the receiver spacing is typically made proportional to the inverse of frequency.
- the receiver spacing for passive seismic acquisition (except for very local micro-earthquakes) is usually larger than the normal active seismic survey due to the lower frequency content of the passive seismic data one tries to acquire.
- this passive monitoring system is used, on a continuous or triggered basis, to record arrivals from all sources, such as earthquakes and man-made activities not intended as seismic sources.
- sources such as earthquakes and man-made activities not intended as seismic sources.
- an ocean bottom cable system or independent sensor packages with geophones capable of recording signals down to 0.5 Hz could be deployed for an extended time period in the survey area.
- a dense array of broadband receivers would be recording continuously or triggered over a predetermined period of time.
- a long continuous deployment time allows for recording of more data from more sources, providing a larger signal-to-noise ratio in the processed passive seismic data, whereas a triggering mechanism could be used to cause recording to begin only when an event is detected; reducing the amount of memory required to store the data.
- the data are processed to extract the signals from the continuous or triggered recording.
- the processing is according to the manner in which the passive seismic data are going to be used in step 26.
- Some embodiments of the invention use signals of all passive sources, i.e. natural and man-made not employed purposefully in a seismic acquisition program (such as sources from mining activities, etc.), all magnitudes (to as small as microseismicity), all distance ranges (including but not limited to local, regional, and teleseismic), and all wave types (body and surface waves). All these sources provide pieces of information, which together allow the frequency range to be expanded to a broad band between 0 Hz and 10 Hz.
- the extracted seismic signals may be used jointly with other type(s) of geophysical data to obtain subsurface structure and rock properties.
- both passive seismic data and the additional geophysical data are jointly used to derive subsurface images and rock property models, either by using joint inversion techniques or by applying inversion and/or imaging methods on the individual datasets sequentially.
- this step is implemented within a joint or simultaneous inversion, in which a subsurface model is obtained by minimizing the difference between the predicted and measured data.
- Alternative implementations are also possible. For example, one can derive the subsurface rock properties and/or structural image from the passive seismic data first, and then use the other geophysical data to refine the results.
- the geophysical data other than the passive seismic data can for example be either one of or a combination of active seismic data, CSEM, MT, gravity, magnetic, and well-log data.
- the joint-inversion can be performed in either the time or frequency domain.
- Different aspects of the passive seismic data can be used with appropriate techniques in the joint inversion process.
- the arrival times in the passive seismic data can be used to constrain, via travel time tomography, the subsurface seismic wave velocity model in the joint inversion; or reflection seismic data similar to active seismic data can be generated from the passive seismic data using the cross correlation technique (see Artman, op.cit. (2006)).
- the passive source, receiver function techniques (e.g., Langston and Phinney, op. cit.) will next be extended to the smaller scales of interest to hydrocarbon exploration.
- receiver functions calculated for a multi-component, active-source dataset. That reference is "Receiver Function Method in Reflection Seismology,” Edme and Singh, Geophysical Prospecting 56, 327-340 (2008).
- the paper discusses the application of the receiver function technique on wide-angle multi-component active seismic data, in an attempt to highlight the P-to-S converted phases, which can then be used for amplitude vs. slowness or offset analysis.
- the present invention uses a survey design with a significantly smaller spacing between the receivers than is normally employed in global seismology.
- the increased density of receivers is designed to allow for target depths typical of hydrocarbon exploration to be resolved.
- a long deployment with continuous recording will usually be necessary and will supply data from events of all distance ranges, which will provide an increased signal-to-noise ratio in the processing and imaging of the passive source seismic data.
- the receiver function technique removes the effects of source and propagation; its result is a waveform that reflects the medium in the vicinity of the recording station.
- Receiver functions are calculated by deconvolving (or cross-correlating) the vertical component seismogram from the horizontal component seismogram.
- the expected seismic phases are the P-to-S converted phase, and its reverberations PpPs and PsPs + PpSs; however, the technique can also be modified to highlight S to P converted phases.
- Figures 9A-D show a simplified illustration of the seismic phases of the P-to-S converted phase (9A) and its reverberations (PpPs in 9B, PsPs in 9C, and PpSs in 9D) that can be expected from the receiver function analysis.
- the differential arrival times can be used to construct the depth of the reflector of interest, and can be integrated with other geophysical datasets, such as active source seismic datasets, gravity, CSEM, MT, etc.. If traveltime information from multicomponent active source datasets is available, the differential arrival times from the receiver function method can be used along with the active source dataset. In addition to using the traveltime information independently, use can be made of the amplitude information from the receiver functions, which can be accomplished by a full waveform inversion. In this approach, initial models of the earth and the source are assumed, which allows computation of initial data by solving the appropriate wave equation with a sufficiently broadband source spectrum.
- these initial models may be modified in an inversion algorithm which is adjusted for both the earth medium and source parameters. This procedure can be accomplished in either the time or frequency domain.
- the main advantage of an error function realization in the frequency domain is that specific missing frequency values or low signal-to-noise frequency values will not affect good data values as they would with an error function in the time domain.
- the present inventive method can use events from the local and regional distance range for the deconvolution.
- the advantage of including these events in a closer vicinity to the recording station is the opportunity to increase the frequency band to a broad range of ⁇ 1 Hz up to frequencies comparable to these of active seismic exploration, and thus reduce the frequency gap 10 shown in Figure 1 and discussed previously herein.
- the difficulty in adding these events is that the assumptions made for teleseismic receiver functions may not be valid any longer, since the deconvolution may not completely remove the effects of the source-time function.
- the present invention can deal with this problem by (1) calibrating the earth response with some man-made sources (such as mining activities or other active sources) which have simple or known source-time functions, and (2) pursuing a full waveform inversion approach, which will enable determination of the remaining source effects.
- man-made sources such as mining activities or other active sources
- the objective is a technique for recovering subsurface geophysical parameters such as P- and S-wave velocity, density, elastic anisotropy coefficients ( ⁇ and ⁇ ), P- and S- wave attenuation quality factors, etc. based upon broadband measurements of multiple- component surface (or near surface) vector (i.e. three component), acoustic (or visco- acoustic) or elastic (or visco-elastic) displacement, velocity, as well as acceleration and surface (or near surface) elastic (or visco-elastic) tensor (i.e.
- broadband multi-component earthquake data nine component stress (collectively called broadband multi-component earthquake data) and forming multi-channel ratios of the broadband multi-component earthquake data (in the spirit of receiver functions) to reduce the unknown effects of the broadband earthquake source, aspects of earth close to the detector, and some intrinsic aspects of the detection system itself as well.
- the resulting multi-channel ratios of the broadband multi-component earthquake data will obey equations that are obtained from the appropriate acoustic, elastic, visco-acoustic, visco-elastic, wave equations. The complexity of the specific equations selected will depend upon survey objectives, data coverage and many other factors.
- models may be one, two, two and one-half or three dimensional and may involve acoustic, elastic or visco-elastic wave propagation and be isotropic, VTI or contain completely general twenty-one (complex) moduli plus density varying in all three dimensions.
- the multi-channel ratios of the broadband multi-component earthquake data preferably satisfy:
- the measurements should be sufficiently broadband i.e. the measurements should contain useful signal content down to very small frequencies and be without significant gaps up to some maximum frequency of interest.
- the number, quality, vector or tensor component and spatial distribution of the complete set of useful observations should effectively enable determination of both the medium properties to the level desired and the surviving unknown parameters describing the earthquake source dependence in the multi-channel ratios.
- Clearly ratio data from a single distant earthquake source will determine at best only certain averages of the earth volume between the earthquake source and the observation points along the paths taken by P-waves in the acoustic approximation, or P-wave and S-wave rays in the elastic or visco-elastic approximation or the qP-wave, qSV-wave and qSH-wave (quasi Compressional wave, quasi Shear Vertical wave and quasi Shear Horizontal wave in the one dimensional VTI model, etc.
- Data analysis or inversion in this embodiment of the present invention then involves a least-squares inversion based upon an error or objective function (including possible smoothing and penality terms) that contains the weighted and squared difference between measured and computed ratios constructed from suitably time windowed real and computed data.
- an error or objective function including possible smoothing and penality terms
- frequency domain error functions are more robust in the presence of missing frequencies. Again the time windowing process is preferred to ensure that only one significant earthquake (or other passive) source is responsible for a significant portion of the events observed in the windowed data.
- the earthquake data is assumed to be broadband, convergence of the least-squares gradient-based family of algorithms to false local minima is not foreseen as a problem (unlike the active-source reflection seismic inversion problem which unavoidedly suffers from contamination with local false minima in the absence of a nearly perfect initial velocity model).
- an initial earth model and an initial source model are assumed, and initial data are computed using a sufficiently broadband source spectrum with arbitrary amplitude using finite difference or other methods in the time or frequency domain as required by the selected numerical approach.
- derivatives of the computed data are obtained in order to allow the inversion algorithm to adjust both earth medium parameters (velocities, densities etc) and earthquake source parameters to decrease the initial fit error.
- the method of alternation between earth medium and earthquake source parameters widely employed in gradient-based optimization methods with non-homogeneous parameters may be employed.
- One alternative is the sub-space iteration method which groups inversion unknowns into physically non-heterogeneous classes to avoid the pitfalls of stiff optimization in which variables of vastly different numerical scale appear in the normal equations.
- the iterative inversion procedure continues adjusting medium and source parameters until a satisfactory fit to the measured and windowed data is achieved. Subsequent model testing can be used to help determine well constrained versus poorly constrained parameters.
- Receiver function results obtained either via an inversion method such as is described above or by the standard approach of manual iterative forward modeling can then be used to create a low- frequency or background model for AVO inversion of active source seismic reflection data in order to obtain estimates of P and S-impedance or even rock properties such as the porosity, lithology and fluid content of subsurface rocks (e.g., Saltzer, Leading Edge 24, 732-736 (2005)).
- the scattered wave results can be used to improve the velocity model used for depth migration of active source seismic reflection data.
- the passively recorded seismic data can be jointly inverted with other geophysical data such as active source seismic reflection data, CSEM, MT, gravity and gravity gradiometry, in order to obtain more robust estimates of elastic and/or rock properties.
- CSEM active source seismic reflection data
- MT gravity and gravity gradiometry
- the passively recorded seismic data can be jointly inverted with other geophysical data such as active source seismic reflection data, CSEM, MT, gravity and gravity gradiometry, in order to obtain more robust estimates of elastic and/or rock properties.
- the event information gathered from the prior step is used to cut the recordings into time windows, with one event per window. For each event identified, all receiver data are collected and an "event gather” is created, and noise is suppressed by filtering and stacking (step 45). All event gathers are collected, and at step 46 receiver ratios ("receiver functions") are calculated. A tomographic inversion is calculated (step 47) using the receiver functions, which has to include a P-wave starting model, determined from another data source (such as an active seismic source experiment).
- an initial velocity model is constructed based on the tomographic inversion and possibly other data, and a full waveform inversion for the initial model and source parameters is performed until the maximum resolution answer is achieved, i.e. until the inversion bandwidth is extended to an upper frequency limit that is high enough to be deemed acceptable.
- the results are verified for geologic reasonableness, and step 48 is iterated as needed, after adjusting the velocity model and/or the source parameters.
Abstract
Description
Claims
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200980154238.7A CN102272631B (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
EA201170920A EA026650B1 (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
CA2747426A CA2747426C (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
AU2009335964A AU2009335964B2 (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
BRPI0923960-0A BRPI0923960B1 (en) | 2009-01-09 | 2009-12-14 | METHODS FOR HYDROCARBON DETECTION IN A SUBSUPERFACE REGION AND TO PRODUCE HYDROCARBON FROM A SUBSUPERFACE REGION |
EP09837859.9A EP2386065A4 (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
US13/140,749 US8923094B2 (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14369509P | 2009-01-09 | 2009-01-09 | |
US61/143,695 | 2009-01-09 | ||
US26769509P | 2009-12-08 | 2009-12-08 | |
US61/267,695 | 2009-12-08 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2010080366A1 true WO2010080366A1 (en) | 2010-07-15 |
WO2010080366A8 WO2010080366A8 (en) | 2010-12-16 |
Family
ID=42316710
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2009/067913 WO2010080366A1 (en) | 2009-01-09 | 2009-12-14 | Hydrocarbon detection with passive seismic data |
Country Status (8)
Country | Link |
---|---|
US (1) | US8923094B2 (en) |
EP (1) | EP2386065A4 (en) |
CN (1) | CN102272631B (en) |
AU (1) | AU2009335964B2 (en) |
BR (1) | BRPI0923960B1 (en) |
CA (1) | CA2747426C (en) |
EA (1) | EA026650B1 (en) |
WO (1) | WO2010080366A1 (en) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012135225A2 (en) * | 2011-03-30 | 2012-10-04 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
WO2012030486A3 (en) * | 2010-08-30 | 2012-12-06 | Microseismic, Inc. | Method for detection of subsurface seismic events in vertically transversely isotropic media |
WO2012173718A1 (en) | 2011-06-17 | 2012-12-20 | Exxonmobil Upstream Research Company | Domain freezing in joint inversion |
WO2013012470A1 (en) * | 2011-07-21 | 2013-01-24 | Exxonmobil Upstream Research Company | Adaptive weighting of geophysical data types in joint inversion |
WO2012160331A3 (en) * | 2011-05-23 | 2013-03-07 | Imperial Innovations Limited | A method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits |
WO2013016011A3 (en) * | 2011-07-28 | 2013-08-15 | Bp Corporation North America Inc. | Field correlation for real-time passive seismic surveillance |
US8584519B2 (en) | 2010-07-19 | 2013-11-19 | Halliburton Energy Services, Inc. | Communication through an enclosure of a line |
US8633700B1 (en) | 2013-03-05 | 2014-01-21 | Hunt Energy Enterprises, Llc | Sensors for passive electroseismic and seismoelectric surveying |
US8873334B2 (en) | 2013-03-05 | 2014-10-28 | Hunt Energy Enterprises, L.L.C. | Correlation techniques for passive electroseismic and seismoelectric surveying |
CN104536047A (en) * | 2015-01-04 | 2015-04-22 | 奥菲(北京)石油技术有限公司 | Technology for applying microseism to rock burst monitoring |
WO2015102498A1 (en) * | 2014-01-03 | 2015-07-09 | Octio As | System and method for processing microseismic data |
WO2015124961A3 (en) * | 2014-02-24 | 2015-11-26 | Cgg Services Sa | Combination of controlled and uncontrolled seismic data |
US9213116B2 (en) | 2011-07-19 | 2015-12-15 | Halliburton Energy Services, Inc. | System and method for moment tensor migration imaging |
EP2715603A4 (en) * | 2011-06-02 | 2016-07-13 | Exxonmobil Upstream Res Co | Joint inversion with unknown lithology |
US9823373B2 (en) | 2012-11-08 | 2017-11-21 | Halliburton Energy Services, Inc. | Acoustic telemetry with distributed acoustic sensing system |
WO2017220918A1 (en) | 2016-06-23 | 2017-12-28 | Storengy | Method for characterising the underlying ground of a region using passive seismic signals, and corresponding system |
WO2018093476A1 (en) * | 2016-11-18 | 2018-05-24 | Landmark Graphics Corporation | Boundary adjustment of vertical seismic sections of three-dimensional seismic surveys to remove data gaps |
CN109653725A (en) * | 2018-09-13 | 2019-04-19 | 山东鼎维石油科技有限公司 | A layer water flooding degree log interpretation method is stored up based on sedimentary micro and the mixed of rock phase |
WO2020073048A1 (en) * | 2018-10-05 | 2020-04-09 | Es Xplore, L.L.C. | Passive electroseismic surveying |
WO2020099838A1 (en) * | 2018-11-14 | 2020-05-22 | Bp Exploration Operating Company Limited | Passive seismic imaging |
CN112444891A (en) * | 2020-11-20 | 2021-03-05 | 核工业北京地质研究院 | Volcanic rock type uranium polymetallic ore deep ore body detection method |
CN113466939A (en) * | 2021-07-20 | 2021-10-01 | 北京市水电物探研究所 | Micromotion exploration method and micromotion exploration system |
CN115327620A (en) * | 2021-05-11 | 2022-11-11 | 中国石油化工股份有限公司 | Microseism joint time difference stacking positioning method |
Families Citing this family (80)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9285502B2 (en) * | 2009-12-08 | 2016-03-15 | Chevron U.S.A. Inc. | System and method for lacunarity analysis |
US20110218735A1 (en) * | 2010-03-05 | 2011-09-08 | Baker Hughes Incorporated | Real-Time Lithology and Mineralogy Interpretation |
GB2479172B (en) * | 2010-03-31 | 2016-02-10 | Statoil Petroleum As | Estimating interval velocities |
US8358561B2 (en) * | 2010-04-13 | 2013-01-22 | Spectraseis Ag | Bayesian DHI for seismic data |
US8498845B2 (en) * | 2010-04-21 | 2013-07-30 | Exxonmobil Upstream Research Company | Method for geophysical imaging |
US8694299B2 (en) * | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
US9195783B2 (en) * | 2010-08-16 | 2015-11-24 | Exxonmobil Upstream Research Company | Reducing the dimensionality of the joint inversion problem |
US9134442B2 (en) * | 2010-12-16 | 2015-09-15 | Bp Corporation North America Inc. | Seismic acquisition using narrowband seismic sources |
US8700372B2 (en) * | 2011-03-10 | 2014-04-15 | Schlumberger Technology Corporation | Method for 3-D gravity forward modeling and inversion in the wavenumber domain |
US20130085731A1 (en) * | 2011-10-04 | 2013-04-04 | Michele De Stefano | Methods and systems for multiple-domain inversion of collected data |
US9194967B2 (en) * | 2011-11-22 | 2015-11-24 | Global Ambient Seismic, Inc. | Tomographic imaging of fracture-fault permeability zones during drilling operations |
US9075158B2 (en) * | 2011-11-22 | 2015-07-07 | Global Microseismic Services, Inc. | Using a drill bit as a seismic source for SET velocity analysis |
US9091787B2 (en) * | 2011-11-28 | 2015-07-28 | Westerngeco L.L.C. | Separation of simultaneous source data |
US8664586B2 (en) * | 2011-12-08 | 2014-03-04 | Saudi Arabian Oil Company | Super-resolution formation fluid imaging |
US8937279B2 (en) | 2011-12-08 | 2015-01-20 | Saudi Arabian Oil Company | Super-resolution formation fluid imaging with contrast fluids |
KR101413752B1 (en) * | 2012-03-13 | 2014-07-01 | 서울대학교산학협력단 | Seismic imaging system using cosine transform in logarithmic axis |
CN102819038B (en) * | 2012-07-27 | 2013-07-10 | 中国石油天然气股份有限公司 | Method and system for identifying oil and water in carbonate rocks |
EP2932304A4 (en) | 2012-12-14 | 2015-12-30 | Landmark Graphics Corp | Methods and systems for seismic modeling using multiple seismic source types |
US20140236487A1 (en) * | 2013-02-21 | 2014-08-21 | Westerngeco L.L.C. | Methods and computing systems for processing seismic data |
US20140254321A1 (en) * | 2013-03-08 | 2014-09-11 | Chevron U.S.A. Inc. | Methods and systems for determining clathrate presence and saturation using simulated well logs |
US9121965B2 (en) | 2013-03-11 | 2015-09-01 | Saudi Arabian Oil Company | Low frequency passive seismic data acquisition and processing |
US9952340B2 (en) | 2013-03-15 | 2018-04-24 | General Electric Company | Context based geo-seismic object identification |
CN104755961B (en) * | 2013-03-15 | 2018-05-04 | 雪佛龙美国公司 | Wave beam inverting based on Monte Carlo back projection |
US10261215B2 (en) * | 2013-04-02 | 2019-04-16 | Westerngeco L.L.C. | Joint inversion of geophysical attributes |
EP2992360B1 (en) * | 2013-04-29 | 2019-09-18 | Westerngeco LLC | Deghosting with adaptive operators |
EA032008B1 (en) * | 2013-10-28 | 2019-03-29 | Бипи Корпорейшн Норд Америка Инк. | Two stage seismic velocity model generation |
CA2931805C (en) * | 2013-12-05 | 2022-07-12 | Cgg Services Sa | Methods and systems of detecting a microseismic event using an iterative non-linear inversion algorithm |
US10345463B2 (en) * | 2014-02-21 | 2019-07-09 | Ion Geophysical Corporation | Methods and systems for using known source events in seismic data processing |
US9852373B2 (en) * | 2014-06-02 | 2017-12-26 | Westerngeco L.L.C. | Properties link for simultaneous joint inversion |
US20150362623A1 (en) * | 2014-06-12 | 2015-12-17 | Westerngeco, Llc | Joint inversion of attributes |
US10386514B2 (en) | 2014-07-24 | 2019-08-20 | Conocophillips Company | Target-oriented process for estimating fracture attributes from seismic data |
US10401516B2 (en) * | 2014-10-24 | 2019-09-03 | Westerngeco L.L.C. | Travel-time objective function for full waveform inversion |
PL230219B1 (en) * | 2014-10-30 | 2018-10-31 | Inst Technik Innowacyjnych Emag | Method and the system for the analysis of the geological structure and relative stresses in the layers situated over the mining headings in the deep mines |
WO2016099747A1 (en) | 2014-12-18 | 2016-06-23 | Exxonmobil Upstream Research Company | Scalable scheduling of parallel iterative seismic jobs |
TWI541528B (en) * | 2014-12-24 | 2016-07-11 | 財團法人國家實驗研究院 | Method of multi-triggering |
CA2973278C (en) * | 2015-01-14 | 2023-09-19 | Statoil (Beijing) Technology Services Co., Ltd. | Method, system and non-transitory computer-readable medium for forming a seismic image of a geological structure |
US10126450B2 (en) * | 2015-02-18 | 2018-11-13 | Conocophillips Company | Black hole boundary conditions |
US20160266275A1 (en) * | 2015-03-10 | 2016-09-15 | Schlumberger Technology Corporation | Methods for estimating formation parameters |
CA2983269A1 (en) * | 2015-04-30 | 2016-11-03 | Schlumberger Canada Limited | Method to estimate water saturation in electromagnetic measurements |
US10338242B2 (en) | 2015-05-20 | 2019-07-02 | Conocophillips Company | Surface wave tomography using sparse data acquisition |
US10520643B2 (en) * | 2015-10-20 | 2019-12-31 | Pgs Geophysical As | Geophysical inversion using sparse modeling |
US10677948B2 (en) | 2016-03-04 | 2020-06-09 | General Electric Company | Context based bounded hydrocarbon formation identification |
US11454102B2 (en) * | 2016-05-11 | 2022-09-27 | Baker Hughes, LLC | Methods and systems for optimizing a drilling operation based on multiple formation measurements |
EP3458883A1 (en) * | 2016-05-20 | 2019-03-27 | Exxonmobil Research And Engineering Company | Shape-based geophysical parameter inversion |
CN106094053A (en) * | 2016-06-07 | 2016-11-09 | 中国石油天然气股份有限公司 | A kind of hydrocarbon source rock life residence strength calculation method and device |
US11209566B2 (en) | 2016-06-22 | 2021-12-28 | University Of Houston System | Nonlinear signal comparison and high-resolution measurement of seismic or acoustic wave dispersion |
US10067252B2 (en) | 2016-07-25 | 2018-09-04 | Chevron U.S.A. Inc. | Methods and systems for identifying a clathrate deposit |
TWI661214B (en) * | 2016-11-29 | 2019-06-01 | National Applied Research Laboratories | On-site earthquake early warning system and method thereof accommodating automatic site effect calibration |
CA3043310C (en) * | 2016-12-02 | 2022-07-26 | Ratnanabha SAIN | Method for estimating petrophysical properties for single or multiple scenarios from several spectrally variable seismic and full wavefield inversion products |
GB201700399D0 (en) * | 2017-01-10 | 2017-02-22 | Reeves Wireline Tech Ltd | Improved method of and apparatus for carrying out acoustic well logging |
CN106886047A (en) * | 2017-02-28 | 2017-06-23 | 中国地质大学(北京) | A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio |
US10317558B2 (en) | 2017-03-14 | 2019-06-11 | Saudi Arabian Oil Company | EMU impulse antenna |
US10416335B2 (en) | 2017-03-14 | 2019-09-17 | Saudi Arabian Oil Company | EMU impulse antenna with controlled directionality and improved impedance matching |
US10330815B2 (en) | 2017-03-14 | 2019-06-25 | Saudi Arabian Oil Company | EMU impulse antenna for low frequency radio waves using giant dielectric and ferrite materials |
CN107102355B (en) * | 2017-04-27 | 2017-12-05 | 吉林大学 | The parallel Marchenko imaging methods of low-frequency reconfiguration |
US10705240B2 (en) * | 2017-05-11 | 2020-07-07 | Saudi Arabian Oil Company | Capacitive electromagnetic formation surveillance using passive source |
US11180975B2 (en) * | 2017-05-31 | 2021-11-23 | Schlumberger Technology Corporation | Geologic structural model generation |
CN107622324A (en) * | 2017-09-01 | 2018-01-23 | 燕山大学 | A kind of robust environmental economy dispatching method for considering more microgrid energy interactions |
US10365393B2 (en) | 2017-11-07 | 2019-07-30 | Saudi Arabian Oil Company | Giant dielectric nanoparticles as high contrast agents for electromagnetic (EM) fluids imaging in an oil reservoir |
CN108169817B (en) * | 2017-12-15 | 2019-11-15 | 中国石油集团川庆钻探工程有限公司 | Hydrocarbon index inverting Fluid Identification Method is combined in well logging earthquake |
CN108107478B (en) * | 2017-12-25 | 2019-06-04 | 湖南科技大学 | Magnetotelluric synchronizing detection and real time inversion method and system |
US11112515B2 (en) * | 2018-03-30 | 2021-09-07 | Bp Corporation North America Inc. | Seismic velocity derived hydrocarbon indication |
CL2018001759A1 (en) * | 2018-06-27 | 2018-12-07 | Tecnoexplora Spa | System and method of three-dimensional identification of geological alterations for mining prospecting and determination of rock quality. |
WO2020009752A1 (en) * | 2018-07-02 | 2020-01-09 | Exxonmobil Upstream Research Company | Full wavefield inversion with an image-gather-flatness constraint |
SG11202012669VA (en) | 2018-07-31 | 2021-02-25 | Exxonmobil Upstream Res Co | Fluid saturation model for petrophysical inversion |
US11181653B2 (en) | 2018-09-14 | 2021-11-23 | Exxonmobil Upstream Research Company | Reservoir characterization utilizing ReSampled seismic data |
US11009617B2 (en) * | 2019-02-20 | 2021-05-18 | Saudi Arabian Oil Company | Method for fast calculation of seismic attributes using artificial intelligence |
AU2020325168A1 (en) | 2019-08-06 | 2022-03-03 | ExxonMobil Technology and Engineering Company | Petrophysical inversion with machine learning-based geologic priors |
CN110579804B (en) * | 2019-10-10 | 2020-12-25 | 中国石油化工股份有限公司 | Diffusion filtering method under structure tensor trace constraint based on absolute square gradient |
US11493658B2 (en) * | 2019-12-03 | 2022-11-08 | China Petroleum & Chemical Corporation | Computer-implemented method and system employing nonlinear direct prestack seismic inversion for poisson impedance |
CN111190225B (en) * | 2020-01-13 | 2021-05-18 | 吉林大学 | Rock physical mechanism and earthquake prediction method for salt-space shale oil reservoir |
KR102528547B1 (en) * | 2020-02-07 | 2023-05-04 | 포항공과대학교 산학협력단 | Apparatus and method for early detecting earthquakes |
US11662493B2 (en) * | 2020-04-08 | 2023-05-30 | ExxonMobil Technology and Engineering Comany | Methodology for enhancing properties of geophysical data with deep learning networks |
CN111650669B (en) * | 2020-05-27 | 2022-02-08 | 吉林大学 | Hydraulic fracturing seismoelectric signal combined processing method |
US11905818B2 (en) * | 2020-12-08 | 2024-02-20 | Halliburton Energy Services, Inc. | Deep learning methods for wellbore pipe inspection |
US11892590B2 (en) | 2021-02-11 | 2024-02-06 | Saudi Arabian Oil Company | Tight reservoir fluid characterization by using sonic slowness well logs |
CN113049202B (en) * | 2021-03-08 | 2022-07-12 | 中国地震局工程力学研究所 | Local weighted regression correction method and system for acceleration integral displacement |
CN113219532B (en) * | 2021-04-26 | 2022-08-16 | 中国石油大学(华东) | Full waveform inversion speed modeling method based on low-frequency data compensation |
CN113391351B (en) * | 2021-06-10 | 2024-02-13 | 中国地质科学院 | Method for extracting mine collection area structure based on passive source seismic wave field analysis |
US11965996B2 (en) * | 2022-02-28 | 2024-04-23 | Saudi Arabian Oil Company | Generating low frequency models for seismic waveform inversion in formation regions with limited control wells |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000023821A1 (en) * | 1998-10-16 | 2000-04-27 | Strm, Llc | Method for 4d permeability analysis of geologic fluid reservoirs |
US20080002523A1 (en) * | 2006-06-09 | 2008-01-03 | Spectraseis Ag | VH Reservoir Mapping |
US20080021655A1 (en) * | 2006-06-30 | 2008-01-24 | Saenger Erik H | Signal Integration Measure for Seismic Data |
WO2008042081A1 (en) * | 2006-09-28 | 2008-04-10 | Exxonmobil Upstream Research Company | Iterative inversion of data from simultaneous geophysical sources |
US7424367B2 (en) * | 2004-05-27 | 2008-09-09 | Exxonmobil Upstream Research Company | Method for predicting lithology and porosity from seismic reflection data |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5200928A (en) | 1991-11-07 | 1993-04-06 | Chevron Research And Technology Company | Method for using mode converted P- to S- wave data to delineate an anomalous geologic structure |
US6128580A (en) | 1998-04-17 | 2000-10-03 | Bp Amoco Corporation | Converted-wave processing in many-layered anisotropic media |
WO2000060379A1 (en) * | 1999-04-02 | 2000-10-12 | Conoco, Inc. | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
US7164619B2 (en) | 2001-01-25 | 2007-01-16 | Johan Robertsson | Method of processing marine seismic data and a method of seismic surveying |
US6473696B1 (en) | 2001-03-13 | 2002-10-29 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
US7219115B2 (en) | 2002-03-28 | 2007-05-15 | Council Of Scientific & Industrial Research | Method for stable and accurate solution for an ill-conditioned system of normal equations associated with a prediction filter whose output is defined as a time advanced version of the input |
US6739165B1 (en) | 2003-02-05 | 2004-05-25 | Kjt Enterprises, Inc. | Combined surface and wellbore electromagnetic measurement system and method for determining formation fluid properties |
GB0319201D0 (en) | 2003-08-15 | 2003-09-17 | Westerngeco Ltd | Multi-component seismic data processing |
EP1692545B1 (en) * | 2003-10-24 | 2012-08-15 | Shell Internationale Research Maatschappij B.V. | Time-lapse seismic survey of a reservoir region |
WO2005085909A1 (en) | 2004-02-24 | 2005-09-15 | Kjt Enterprises, Inc. | Combined surface and wellbore electromagnetic measurement system and method for determining formation fluid properties |
US7486589B2 (en) | 2006-02-09 | 2009-02-03 | Schlumberger Technology Corporation | Methods and apparatus for predicting the hydrocarbon production of a well location |
US7328107B2 (en) | 2006-04-28 | 2008-02-05 | Kjt Enterprises, Inc. | Integrated earth formation evaluation method using controlled source electromagnetic survey data and seismic data |
US7620534B2 (en) | 2006-04-28 | 2009-11-17 | Saudi Aramco | Sound enabling computerized system for real time reservoir model calibration using field surveillance data |
US7340348B2 (en) | 2006-06-15 | 2008-03-04 | Kjt Enterprises, Inc. | Method for acquiring and interpreting seismoelectric and electroseismic data |
US20080075920A1 (en) | 2006-09-21 | 2008-03-27 | Yung-Hung Chih | Aluminum alloy housing structure and method making the same |
US20080162050A1 (en) * | 2006-12-28 | 2008-07-03 | Peter Harris | Method for interpreting seismic data and controlled source electromagnetic data to estimate subsurface reservoir properties |
WO2008087505A2 (en) | 2007-01-20 | 2008-07-24 | Spectraseis Ag | Time reverse reservoir localization |
US7724608B2 (en) | 2007-07-20 | 2010-05-25 | Wayne Simon | Passive reflective imaging for visualizing subsurface structures in earth and water |
-
2009
- 2009-12-14 CA CA2747426A patent/CA2747426C/en active Active
- 2009-12-14 WO PCT/US2009/067913 patent/WO2010080366A1/en active Application Filing
- 2009-12-14 EA EA201170920A patent/EA026650B1/en not_active IP Right Cessation
- 2009-12-14 EP EP09837859.9A patent/EP2386065A4/en not_active Withdrawn
- 2009-12-14 US US13/140,749 patent/US8923094B2/en active Active
- 2009-12-14 CN CN200980154238.7A patent/CN102272631B/en active Active
- 2009-12-14 AU AU2009335964A patent/AU2009335964B2/en active Active
- 2009-12-14 BR BRPI0923960-0A patent/BRPI0923960B1/en active IP Right Grant
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000023821A1 (en) * | 1998-10-16 | 2000-04-27 | Strm, Llc | Method for 4d permeability analysis of geologic fluid reservoirs |
US7424367B2 (en) * | 2004-05-27 | 2008-09-09 | Exxonmobil Upstream Research Company | Method for predicting lithology and porosity from seismic reflection data |
US20080002523A1 (en) * | 2006-06-09 | 2008-01-03 | Spectraseis Ag | VH Reservoir Mapping |
US20080021655A1 (en) * | 2006-06-30 | 2008-01-24 | Saenger Erik H | Signal Integration Measure for Seismic Data |
WO2008042081A1 (en) * | 2006-09-28 | 2008-04-10 | Exxonmobil Upstream Research Company | Iterative inversion of data from simultaneous geophysical sources |
Non-Patent Citations (1)
Title |
---|
See also references of EP2386065A4 * |
Cited By (55)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9003874B2 (en) | 2010-07-19 | 2015-04-14 | Halliburton Energy Services, Inc. | Communication through an enclosure of a line |
US8584519B2 (en) | 2010-07-19 | 2013-11-19 | Halliburton Energy Services, Inc. | Communication through an enclosure of a line |
CN103109207A (en) * | 2010-08-30 | 2013-05-15 | 麦克罗地震探测公司 | Method for detection of subsurface seismic events in vertically transversely isotropic media |
WO2012030486A3 (en) * | 2010-08-30 | 2012-12-06 | Microseismic, Inc. | Method for detection of subsurface seismic events in vertically transversely isotropic media |
US9513395B2 (en) | 2010-08-30 | 2016-12-06 | Microseismic, Inc. | Method for detection of subsurface seismic events in vertically transversely isotropic media |
US10378316B2 (en) | 2011-03-30 | 2019-08-13 | Es Xplore, L.L.C. | Method and system for passive electroseismic surveying |
GB2504419A (en) * | 2011-03-30 | 2014-01-29 | Hunt Energy Entpr Llc | Method and system for passive electroseismic surveying |
WO2012135225A2 (en) * | 2011-03-30 | 2012-10-04 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
US8478533B2 (en) | 2011-03-30 | 2013-07-02 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
US8494776B2 (en) | 2011-03-30 | 2013-07-23 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
US8347658B2 (en) | 2011-03-30 | 2013-01-08 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
US8560241B2 (en) | 2011-03-30 | 2013-10-15 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
US8566037B2 (en) | 2011-03-30 | 2013-10-22 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
US8566036B2 (en) | 2011-03-30 | 2013-10-22 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
AU2015271919B1 (en) * | 2011-03-30 | 2016-02-04 | Es Xplore, L.L.C. | Method and system for passive electroseismic surveying |
US8589079B2 (en) | 2011-03-30 | 2013-11-19 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
WO2012135225A3 (en) * | 2011-03-30 | 2012-11-22 | Hunt Energy Enterprises, Llc | Method and system for passive electroseismic surveying |
CN103562749A (en) * | 2011-05-23 | 2014-02-05 | 帝国创新有限公司 | A method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits |
WO2012160331A3 (en) * | 2011-05-23 | 2013-03-07 | Imperial Innovations Limited | A method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits |
EP2715603A4 (en) * | 2011-06-02 | 2016-07-13 | Exxonmobil Upstream Res Co | Joint inversion with unknown lithology |
US9453929B2 (en) | 2011-06-02 | 2016-09-27 | Exxonmobil Upstream Research Company | Joint inversion with unknown lithology |
WO2012173718A1 (en) | 2011-06-17 | 2012-12-20 | Exxonmobil Upstream Research Company | Domain freezing in joint inversion |
EP2721478A4 (en) * | 2011-06-17 | 2015-12-02 | Exxonmobil Upstream Res Co | Domain freezing in joint inversion |
US9213116B2 (en) | 2011-07-19 | 2015-12-15 | Halliburton Energy Services, Inc. | System and method for moment tensor migration imaging |
WO2013012470A1 (en) * | 2011-07-21 | 2013-01-24 | Exxonmobil Upstream Research Company | Adaptive weighting of geophysical data types in joint inversion |
EP2734866A4 (en) * | 2011-07-21 | 2016-02-17 | Exxonmobil Upstream Res Co | Adaptive weighting of geophysical data types in joint inversion |
US9494711B2 (en) | 2011-07-21 | 2016-11-15 | Garrett M Leahy | Adaptive weighting of geophysical data types in joint inversion |
US9348045B2 (en) | 2011-07-28 | 2016-05-24 | Bp Corporation North America Inc. | Field correlation for real-time passive seismic surveillance |
WO2013016011A3 (en) * | 2011-07-28 | 2013-08-15 | Bp Corporation North America Inc. | Field correlation for real-time passive seismic surveillance |
EA026290B1 (en) * | 2011-07-28 | 2017-03-31 | Бп Корпорейшн Норт Америка Инк. | Method of seismic surveillance above a region of the subsurface of the earth |
US9823373B2 (en) | 2012-11-08 | 2017-11-21 | Halliburton Energy Services, Inc. | Acoustic telemetry with distributed acoustic sensing system |
US10203427B2 (en) | 2013-03-05 | 2019-02-12 | Es Xplore, L.L.C. | Correlation techniques for passive electroseismic and seismoelectric surveying |
US8873334B2 (en) | 2013-03-05 | 2014-10-28 | Hunt Energy Enterprises, L.L.C. | Correlation techniques for passive electroseismic and seismoelectric surveying |
US8633700B1 (en) | 2013-03-05 | 2014-01-21 | Hunt Energy Enterprises, Llc | Sensors for passive electroseismic and seismoelectric surveying |
US9759838B2 (en) | 2013-03-05 | 2017-09-12 | Es Xplore, L.L.C. | Correlation techniques for passive electroseismic and seismoelectric surveying |
US10670754B2 (en) | 2014-01-03 | 2020-06-02 | Octio As | System and method for processing microseismic data |
WO2015102498A1 (en) * | 2014-01-03 | 2015-07-09 | Octio As | System and method for processing microseismic data |
WO2015124961A3 (en) * | 2014-02-24 | 2015-11-26 | Cgg Services Sa | Combination of controlled and uncontrolled seismic data |
US11169292B2 (en) | 2014-02-24 | 2021-11-09 | Sercel Sas | Combination of controlled and uncontrolled seismic data |
CN104536047B (en) * | 2015-01-04 | 2018-08-21 | 奥菲(北京)石油技术有限公司 | A kind of application process of microseism in bump monitoring |
CN104536047A (en) * | 2015-01-04 | 2015-04-22 | 奥菲(北京)石油技术有限公司 | Technology for applying microseism to rock burst monitoring |
WO2017220918A1 (en) | 2016-06-23 | 2017-12-28 | Storengy | Method for characterising the underlying ground of a region using passive seismic signals, and corresponding system |
GB2570811B (en) * | 2016-11-18 | 2021-11-03 | Landmark Graphics Corp | Boundary adjustment of vertical seismic sections of three-dimensional seismic surveys to remove data gaps |
WO2018093476A1 (en) * | 2016-11-18 | 2018-05-24 | Landmark Graphics Corporation | Boundary adjustment of vertical seismic sections of three-dimensional seismic surveys to remove data gaps |
GB2570811A (en) * | 2016-11-18 | 2019-08-07 | Landmark Graphics Corp | Boundary adjustment of vertical seismic sections of three-dimensional seismic surveys to remove data gaps |
US10768326B2 (en) | 2016-11-18 | 2020-09-08 | Landmark Graphics Corporation | Boundary adjustment of vertical seismic sections of three-dimensional seismic surveys to remove data gaps |
CN109653725A (en) * | 2018-09-13 | 2019-04-19 | 山东鼎维石油科技有限公司 | A layer water flooding degree log interpretation method is stored up based on sedimentary micro and the mixed of rock phase |
WO2020073048A1 (en) * | 2018-10-05 | 2020-04-09 | Es Xplore, L.L.C. | Passive electroseismic surveying |
WO2020099838A1 (en) * | 2018-11-14 | 2020-05-22 | Bp Exploration Operating Company Limited | Passive seismic imaging |
US11079507B2 (en) | 2018-11-14 | 2021-08-03 | Bp Exploration Operating Company Limited | Passive seismic imaging |
CN112444891A (en) * | 2020-11-20 | 2021-03-05 | 核工业北京地质研究院 | Volcanic rock type uranium polymetallic ore deep ore body detection method |
CN112444891B (en) * | 2020-11-20 | 2023-11-14 | 核工业北京地质研究院 | Volcanic type uranium polymetallic ore deep ore body detection method |
CN115327620A (en) * | 2021-05-11 | 2022-11-11 | 中国石油化工股份有限公司 | Microseism joint time difference stacking positioning method |
CN113466939A (en) * | 2021-07-20 | 2021-10-01 | 北京市水电物探研究所 | Micromotion exploration method and micromotion exploration system |
CN113466939B (en) * | 2021-07-20 | 2024-03-01 | 北京市水电物探研究所 | Micro-motion exploration method and micro-motion exploration system |
Also Published As
Publication number | Publication date |
---|---|
CA2747426C (en) | 2017-05-23 |
EP2386065A4 (en) | 2017-12-06 |
CN102272631A (en) | 2011-12-07 |
AU2009335964A1 (en) | 2011-07-28 |
EA201170920A1 (en) | 2012-02-28 |
WO2010080366A8 (en) | 2010-12-16 |
US20110255371A1 (en) | 2011-10-20 |
EA026650B1 (en) | 2017-05-31 |
BRPI0923960B1 (en) | 2019-07-30 |
CN102272631B (en) | 2015-03-25 |
AU2009335964B2 (en) | 2015-05-14 |
EP2386065A1 (en) | 2011-11-16 |
BRPI0923960A2 (en) | 2016-01-19 |
CA2747426A1 (en) | 2010-07-15 |
US8923094B2 (en) | 2014-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8923094B2 (en) | Hydrocarbon detection with passive seismic data | |
Wapenaar et al. | Tutorial on seismic interferometry: Part 2—Underlying theory and new advances | |
Prieux et al. | Multiparameter full waveform inversion of multicomponent ocean-bottom-cable data from the Valhall field. Part 1: Imaging compressional wave speed, density and attenuation | |
Warner et al. | Anisotropic 3D full-waveform inversion | |
Prieux et al. | On the footprint of anisotropy on isotropic full waveform inversion: the Valhall case study | |
Veeken | Seismic stratigraphy and depositional facies models | |
US10267937B2 (en) | Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data | |
AU2010201504B2 (en) | Method for calculation of seismic attributes from seismic signals | |
EP3359982B1 (en) | Seismic sensor orientation | |
Harmankaya et al. | Locating near-surface scatterers using non-physical scattered waves resulting from seismic interferometry | |
WO2018075738A1 (en) | Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data | |
Harding et al. | Velocity structure near IODP Hole U1309D, Atlantis Massif, from waveform inversion of streamer data and borehole measurements | |
EP4042210A1 (en) | Determining properties of a subterranean formation using an acoustic wave equation with a reflectivity parameterization | |
Jaiswal et al. | Seismic imaging of the Naga thrust using multiscale waveform inversion | |
Bian et al. | Evaluation of multi-scale full waveform inversion with marine vertical cable data | |
Greenwood et al. | Characterization and imaging of a hydrothermally active near-vertical fault zone in crystalline rocks based on hydrophone VSP data | |
Jreij et al. | Improving point-sensor image resolution with distributed acoustic sensing at Brady’s enhanced geothermal system | |
Barros et al. | Basement depths in the Parecis Basin (Amazon) with receiver functions from small local earthquakes in the Porto dos Gaúchos seismic zone | |
Jang et al. | Application of Reverse Time Migration to Faults Imaging in Rakhine Basin, Myanmar | |
Tylor-Jones et al. | Processing Essentials | |
Jang et al. | Research Article Application of Reverse Time Migration to Faults Imaging in Rakhine Basin, Myanmar | |
Seeni et al. | Processing and preliminary interpretation of the ultra high-density full-azimuth 3D seismic survey, Dukhan field, Qatar | |
Sæther | Seismic Forward Modeling of Deltaic Sequences | |
Guo et al. | Seismic evidence for velocity heterogeneity along~ 40 Ma old oceanic crustal segment formed at the slow-spreading Mid-Atlantic Ridge in the equatorial Atlantic Ocean from full waveform inversion of ocean bottom seismic data | |
Zong et al. | Examining the significance of the converted-modes beneath an ultra-high-velocity caprock: Physical modeling and a Gulf of Mexico salt proximity VSP survey |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200980154238.7 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09837859 Country of ref document: EP Kind code of ref document: A1 |
|
DPE1 | Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101) | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2747426 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13140749 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009335964 Country of ref document: AU |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2009335964 Country of ref document: AU Date of ref document: 20091214 Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009837859 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 201170920 Country of ref document: EA |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: PI0923960 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: PI0923960 Country of ref document: BR Kind code of ref document: A2 Effective date: 20110708 |