CN110050204A - A method of earthquake-capturing is improved using active ultralight seismic acquisition system - Google Patents
A method of earthquake-capturing is improved using active ultralight seismic acquisition system Download PDFInfo
- Publication number
- CN110050204A CN110050204A CN201780073625.2A CN201780073625A CN110050204A CN 110050204 A CN110050204 A CN 110050204A CN 201780073625 A CN201780073625 A CN 201780073625A CN 110050204 A CN110050204 A CN 110050204A
- Authority
- CN
- China
- Prior art keywords
- source
- receiver
- interest
- area
- trace
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000013508 migration Methods 0.000 claims abstract description 13
- 230000005012 migration Effects 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims description 9
- 239000004575 stone Substances 0.000 claims 1
- 239000002872 contrast media Substances 0.000 description 9
- 238000001514 detection method Methods 0.000 description 9
- 239000003153 chemical reaction reagent Substances 0.000 description 8
- 239000007789 gas Substances 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 239000012530 fluid Substances 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 229920000642 polymer Polymers 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000011084 recovery Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000009877 rendering Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000002105 nanoparticle Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011261 inert gas Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000024241 parasitism Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
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/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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, 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/282—Application of seismic models, synthetic seismograms
-
- 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
- G01V1/308—Time lapse or 4D effects, e.g. production related effects to the formation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3817—Positioning of seismic devices
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of methods for being determined to obtain source-receiver pair best orientation of seismic data, it include: first step, it identifies area-of-interest (32), the area-of-interest be early earthquake acquisition main body, so as to obtain the area-of-interest subsoil image；Second step obtains the seismic data obtained during the early earthquake acquisition of the area-of-interest within the interested time；Third step, using some or all of seismic data inverse migration, with each source-receiver of determination to the position of (31,34), each source-receiver contributes to the image of the subsoil of the area-of-interest within the interested time；Four steps obtains the source-receiver to the untreated trace of position (31,34)；5th step selects at least one best untreated trace from a plurality of untreated trace；And the 6th step, determine the source-receiver for corresponding to described at least one best untreated trace to position (31,34).
Description
This patent disclosure relates generally to land and sea area to acquire seismic data.
It is known from the prior art that the image that the geophysical structure of subsoil can be generated by acquiring seismic data.To the greatest extent
The position of natural resources of such as petroleum or natural gas etc cannot accurately be indicated by managing the image, but for art technology
The image may determine that these resource presence or absence really for personnel.Therefore, continuously generate this kind of image be explore subsoil with
It was found that an essential part for these natural resources processes.In order to generate these images, it is known in the state of the art have it is several
Kind method.For offshore exploration, it is possible to implement ship draws earthquake towed cable.This can be arranged relative to subsea horizontal with constant depth
A little towing cables.The ship also depicts one group of seismic origin that can generate seismic wave.These seimic wave propagations are to seabed and pass through sea
Bottom is to pass through it, until seismic wave encounters catoptric arrangement to reflect them.The earthquake wave direction oversea propagation of these reflections, Zhi Daoji
The seismic wave of these reflections is detected at the seismic sensor into towing cable.Based on the seismic data, therefore subsoil can be generated
Image.If being referred to as 2D image using single towing cable.If using the several towing cables drawn simultaneously by ship,
It is referred to as 3D rendering.When different time carries out 3D earthquake-capturing twice in the same area, 4D image is obtained.In time t0
First time acquisition is carried out, second is carried out in time t0+t and acquires, t can for example be equal to some months or several years.
Alternatively, the cable for being placed on seabed or sensor record seismic data, these technologies referred to as sea also can be used
Bottom cable (OBC), subsea node (OBN) or submarine earthquake meter (OBS).For example, in the case where OBC acquisition, multiple sensors
It links together in the cable.Then some in these cables are mounted in the region to be measured in seabed.It can also be by one
Or multiple sensors are incorporated in the underwater vehicle of seabed.It would be possible to be placed on seabed for the delivery vehicle of automatic Pilot,
To record seismic data.Then ship fetches the delivery vehicle so that the seismic data is transmitted to ship.
On land in the case where earthquake-capturing, the system for being exclusively used in 4D acquisition includes several sensors, for example, hydrophone or
Accelerometer, these hydrophones or accelerometer are placed on expectation and detect on the region of subsoil.These sensors and ground face contact.
It also placed several focus on the ground in region to be measured.Recording device is connected to sensor and for example in truck.Often
A focus includes the vibrator of variable number, and the quantity of vibrator is usually 1 to 5, and can also include local controller.
Central controller be can have to coordinate the launch time of focus.GPS system can be implemented, to be associated with the transmitting in source at any time
And data acquired in sensor.In the configuration, control focus is to generate seismic wave, and multiple sensors record is by oil
The wave that hiding and gas reservoir or other structures reflect.Seismic survey can be carried out in different time different moments, such as monthly or every year,
Subsoil to be imaged again, so that it is determined that the lasting variation that oil reservoir and gas reservoir occur at any time.
All these technologies can be used in monitoring production oil reservoir and gas reservoir.For these configurations, the purpose of 4D processing is logical
The difference for the seismic data that assessment different moments obtain is crossed to determine variation pattern and the position of rock physical property, is generally being grasped
Make to implement before mineral deposit, this is to refer to exploration (baseline exploration), and implement after having operated same mineral deposit, this is monitoring exploration.
Currently, the purpose of 4D exploration solution is to update the 3D acquisition of the subsoil in considered region completely.In other words, baseline is surveyed
It visits and monitoring exploration obtains a large amount of seismic data, to obtain complete 3D rendering.Therefore, come from the angle of repeatability
It sees, it is relatively less accurate that current techniques take a significant amount of time (about several thoughtful some months), because they use impermanency system
System, and generate great amount of cost.All of these factors taken together is all the obstacle of 4D exploration development.Therefore, the cost for not only needing 4D to explore
It is lower, also to make it faster more accurate.
Therefore, the detection that therefore is present invention simultaneously improves active ultralight seismic system especially for raising sensitivity horizontal.
Fig. 1 describes an a kind of example of active ultralight seismic acquisition system.Term is ultralight refer to relative to use it is hundreds of extremely
Millions of a source-receivers to traditional earthquake-capturing for, use source-receiver pair of limited quantity.Here it can be seen that
One single source 1 and receiver 2 are right.It can also realize source antenna and receiver antenna.Source antenna and receiver antenna are wrapped respectively
Aggregation containing source and receiver.For example, they pass through welding assembly.Source antenna expressively seems for single source, receiver antenna
Then expressively as being single receiver.Antenna source allows to the transmitting of focus wave, and receiver antenna allows to focus this
The reception of a little waves.This technology is known as wave beam control.In addition, this allows to create noise filter.Also illustrate water-bearing layer
3, tomography 4, gas reservoir 6, area-of-interest or point 7 and production site 5.Source-receiver pair quantity limits the sky so that traditional
Between noise suppressing method less be applicable in.These traditional spatial noise suppressing methods can cover low-producing influence sometimes, i.e.,
By the physical property put on the time span considered variation and change the seismic character less than 1%.In time span
It is especially short, such as about 3 months, one day or when a hour, it is even more so, because producing caused ground in these times
It is low to shake effect.It detects these variations and allows to predict upcoming influence, such as can usually be seen in a very long time
The variation of the saturation degree arrived and its relative influence.These influences disclose interested under the influence of variation in a short time
The early stage in region excites.
Therefore, for active ultralight seismic system, it is necessary to improve the detectivity in the short time.
In order to solve this problem, the purpose of the present invention is a kind of for being determined to obtain source-reception of seismic data
The method of the best orientation of device pair, comprising the following steps:
First step identifies area-of-interest, which is the main body of early earthquake acquisition, to obtain
State the image of the subsoil of area-of-interest；
Second step is obtained within the interested time and is obtained during the early earthquake acquisition of the area-of-interest
Seismic data；
Third step, it is described every with the position of each source-receiver pair of determination using the inverse migration of the seismic data
A source-receiver contributes to the image of the subsoil of the area-of-interest within the interested time；
- the four steps obtains the source-receiver to the untreated trace of position；
- the five step selects at least one best untreated trace from the untreated trace；
- the six step determines the best source-receiver for corresponding to the best untreated trace to position；
Preferably, the inverse migration of third step is carried out by octave coverage.
Preferably, the early earthquake acquisition is 4D acquisition.
Preferably, the 5th step further includes the optimal untreated trace of selection, allows to detect best 4D effect.
Preferably, the 5th step further includes using petroleum elastic model to select optimal untreated trace.
Preferably, the 5th step further includes considering surface factor.
Preferably, the 5th step further includes the reference zone for being considered as earthquake variation calibration.
Preferably, third step further include obtain the seismic wave of each source-receiver pair best slope and transmitting and
Receiving direction.
Preferably, this method further includes by source-receiver to the 7th step for being placed on the optimum position and described
Optimum position obtains the 8th step of new seismic data.
Preferably, the 8th step is based on the best slope of seismic wave and transmitting and receiving direction.
It is described below by reference to what attached drawing provided, the present invention will be best understood, in which:
- Fig. 1 shows an example of active ultralight seismic acquisition system；
Fig. 2 shows one embodiment of the present of invention；
- Fig. 3 shows another embodiment of the invention；
- Fig. 4 shows the example for selecting the realization of method of at least one best trace object of the invention；
- Fig. 5 shows the seismic wave across reference area.
- Fig. 6 shows the example for the method processing untreated trace of earthquake realized in through the invention.
It is well known in the art that chemical reagent, such as polymer are used during extracting the fluid being located in subsoil.Make
The yield of the rate of recovery can be made to improve 10% or more with chemical reagent.In the case where the secondary recovery of fluid, these changes are used
It is even more reasonable to learn reagent.These polymer are injected in the upstream in secondary recovery stage, i.e. another injection, to generate thrust dress
Set, this thrust device have with effect as piston type, then promote the thrust device.In addition, in most cases,
It is obtained and has handled 2D or 3D earthquake.This especially obtain following information, such as the seismic chart of rate pattern, migration
Picture, multiple wave pattern, earthquake be untreated and the trace of processing and the source and the receiver that are used in particular for existing 2D or 3D rendering
Accurate location and land data static models.By the way that following information can be obtained to being appreciated and understood by for subsoil, such as
The resilient nature of target rock and expected earthquake generate effect, i.e., under the influence of production resilient nature variation.To landform
Understanding also especially obtain following additional information, such as surface interruptions object (road, pipeline, well, platform, factory, live
Residence etc.) or not reproducible seismic noise generator.
The present invention is integrated in fluid injection especially for by earthquake contrast agent, and this earthquake contrast agent can be liquid
Or gas, to improve the earthquake-capturing for implementing ultralight seismic acquisition system.Most of ultralight seismic acquisition systems be intended to detect with
The relevant variation of the generation of fluid or gas.The chemical reagent of such as polymer is used to improve the rate of recovery of oil to be extracted.This
One embodiment of invention includes injection as chemical reagent and/or the earthquake contrast agent for the supplement for injecting fluid, thisly
Shake contrast agent is substantially enlarged or reduced seismic response.This allow to increase the seismic response of injection and therefore can to the greatest extent can
The arrival of the front for the fluid being injected into area to be monitored can be detected fastly.Most of the time, chemical reagent is located at front
Front, to form so-called chemical piston.This chemical reagent, liquid or gas are injected chemical internal piston, in chemistry
The above or below of piston.Preferably, the earthquake contrast agent must be neutral for chemical piston, i.e., it cannot be dropped
The property of low chemical reagent used.Earthquake contrast agent may include high density or low-density nano particle and/or resonant frequency and ground
Shake acquires compatible nano particle.Earthquake contrast agent can also contain inert gas, such as nitrogen.Earthquake contrast agent also may include polymerization
Object.Fig. 2 shows oil slick 21, wherein being injected water 25 and polymer 22 in oil slick 21 by syringe 24.Earthquake contrast agent 23
It is also coupled in the medium of polymer piston 22.
The present invention also aims to a kind of method, this method can improve the positioning of source-receiver pair, and can be with
Seismic processing, the i.e. record at the output of sensor are carried out to single untreated trace, or as the time is averaged, that is,
Increase during untreated trace or stack, to reduce non-repeatable noise or increasing under the frame of permanent or semipermanent system
Strong signal intensity.In one embodiment of the invention, identification needs the region detected to collect available information to implement this
Method.The first step is in interested region (also referred to as point) or nearby to use existing seismic data.The 2D, 3D, 4D
It shakes data or well and is used as information source, to know in surveyed subsoil, seismic wave especially in region to be detected
Path.Preferably, it checks the launch point centered on region to be detected, is deviated by removal, is i.e. removal launch point and receiver
The distance between it is preselected to execute first, the offset has too many noise for effectively detecting, i.e. signal-to-noise ratio is too low.So
After execute inverse migration, the i.e. converse digraph of seismic ray, available all or part earthquake data.Use this of area-of-interest
Each source of the available imaging for facilitating point-of-interest of return radiation figure and theoretical position X, Y and Z of receiver, can be with
Retrieve the time that area-of-interest appears in untreated trace.This time is known as the interested time.Pass through the step
It can also determine the arrival slope for having irradiated the seismic wave of area-of-interest.Fig. 3 shows the realization of the method.In inverse migration
The step of in, can see source 31, receiver 34, area-of-interest or point 32 and return radiation Figure 33 herein.Such as
In the case where complicated geology and the finer analysis of needs, inverse migration can also carry out in octave coverage.
Then the selection of at least one best trace is realized.This method is applied to contributive to the imaging of area-of-interest
Region, and include the steps that being selected from three parameters.First parameter is the effective contribution of source-receiver pair, i.e.,
The energy or amplitude of untreated record, the source-receiver are felt to the building in existing 2D, 3D seismic data or well is allowed to
The image in interest region, and there are also the global signal-to-noise ratio of trace.In fact, not only checking the row of trace in the region of interest
For, but also except area-of-interest check trace behavior, with detection low signal-to-noise ratio middle except area-of-interest.Trace
It can have high s/n ratio in the region of interest, but there is low signal-to-noise ratio in the outside of area-of-interest.Preferably, using with
Relevant second parameter of multidate information, i.e., existing 4D seismic data and/or expected production effect and its shake trace over the ground
Influence modeling.Preferably, using third parameter relevant to surface information, such as the obstacle in source or receiver positioning is prevented
Object, the strong noise generator to be avoided, a possibility that accessing the current source that must support or the network coverage.Before by considering
In the case that two parameters identify several optimum positions, final choice will be surface standard.
By during inverse migration step by the reality of theoretical position X, Y and Z of source and receiver and existing seismic data
Grid intersection is acquired to obtain the first parameter, restores untreated trace, i.e., effective survey of the subsoil obtained according to these positions
Magnitude.The temporal information on untreated trace is appeared in using these untreated traces and area-of-interest, it can be to latent
Qualitative and quantitative choosing is carried out in the subset of best trace.When trace can detect variation, which is optimal.Anti-
Temporal information is obtained in bias step, which is also referred to as the interested time.During the interested time, strong noise
Trace or the trace of altitude decay are dropped.On the contrary, remaining the trace with strong reflection or refraction, this trace is referred to as clear
Clean trace.In view of multiple that is visible on launch point and/or being obtained during handling existing seismic data influence so that
It can choose and not occur the trace of high magnification numbe during each interested time.Fig. 4 shows five traces 41,42,43,44
With 45.Axis 46 on axis of ordinates indicates offset, i.e. source-receiver distance.Trace 41 has too many noise and 42 He of trace
45 during the interested time excessive attenuation.Therefore, non-selected these three traces are (that is, trace 41,42 and 45).Another party
Face, trace 43 and 44 have good signal-to-noise ratio and the visible energy during the detection time of area-of-interest, therefore, it is considered that
The two traces 43 and 44 may be optimal, to be retained.Second parameter is related to multidate information, when 4D seismic data
When available, the multidate information can get.It is emerging using whether complete 4D image can influence considered sense with Detect change
Interesting region.If variation has influenced considered area-of-interest, by being obtained not to from basis and post-collection
Processing trace is analyzed, and can determine the trace for having best 4D effect for area-of-interest.These traces are to identification
It is this to change maximum trace of making contributions, therefore these traces are for detecting new 4D effect in the region of interest
It is optimal.As the supplement of existing 4D seismic trace, or when existing 4D trace is unavailable, mutually tied using with Dynamic deposition model
The rock elasticity model of the area-of-interest of conjunction can model synthesis earthquake-capturing, which estimates by producing
Earthquake sheet caused by raw.
Third parameter is related to the information in relation to surface factor, and the information of these surface factors is obtained from the first parameter or front two
The subset that the combination of a factor is identified can be such that the position of selected source-receiver pair intersects with surface information.This letter
Breath intersects the reality that on the one hand can verify these source-receiver pair positioning not in forbidden zone by inspection such as selected location
Border feasibility.This information, which is intersected, makes it also possible to consider that the high seismic noise that can be appeared between existing seismic data occurs
Device, such as in the road, gas compression unit or marine corridor built.
It can analyze the feasibility of active ultralight seismic survey by these steps.
Described method so far can be supplemented, to improve the seismic processing of active ultralight system, and
In variation of the Jiang Yuan-receiver to area-of-interest is detected when being permanently or semi-permanently placed on optimum position.This allows to
Detection is concentrated on the region of interest, without the need to build the complete image of subsoil.This is a kind of method of trace to trace,
It is middle respectively to each source-receiver to being observed and handled.In order to overcome the variation just occurred on the region of interest, it is assumed that
6 months about most in very short time span, can identify in subsoil square on the region of interest will not send out
At least one region for changing.The region is known as reference zone.The reference zone is related to ultralight detection system for obtaining
Ambient noise measured value.Reference zone is also used to correct the variation reached on this reference zone.It is located at ginseng to correct
The 4D effect of overlying regions is examined, wherein the variation at the reference zone by near surface becomes apparent, it is contemplated that the ground in the region
Ringing should must not change, and reference zone is used as calibration.Therefore can correct the region earthquake change and apply the correction with
Disclose the variation of area-of-interest.These variations are the variations of amplitude, journey time or Wavelet temporal.It can see in Fig. 5
Above-mentioned source 51, receiver 52, seismic wave 53, area-of-interest 54 and reference zone 55.Fig. 6, which is shown, is applied to this method not
One embodiment of trace is handled, wherein the method is subject of the present invention.Abscissa indicates the time as unit of day, indulges
Reference axis then indicates the difference of the traveling time of seismic wave whithin a period of time.The starting point of stroke is the seismic origin, followed by mineral deposit,
It is finally receiver.Upper curve corresponds to every day entry, and lower curve corresponds to first day derivative of record.It is observed at the 25th day
Significant changes are arrived, it means that geomechanics variation has occurred.This geomechanics variation can be such as saturation degree, pressure
Or temperature change.
In addition, the concept of cell or case no longer has in all senses when considering the single seismic trace of ultralight acquisition.
The spatial resolution of search coverage is considered as Fresnel region, and Fresnel region itself depends on frequency.Then can by frequency multiplication into
Row filtering, to improve the sensitivity of detection.By reducing the frequency of received data, increase in farther area-of-interest
Detectivity.To achieve it, ultralight acquisition system must be emitted and records several frequency ranges.For example,
For 5500 meters of the speed per second of area-of-interest, it is located at the Fresnel apart from 50,100,150,200 and 275 meters of dot center
Area corresponds respectively to 27,14,9,7 and 5 hertz of frequency, is used for even speed model and zero offset collection model.Therefore, in order to
The variation present in 5 hertz but being not present at 14 hertz is detected, this allows to be inferred to observe in the region of interest
Between 100 to 275 meters of the center of the variation interference distance area-of-interest arrived.
It on the other hand, can be for each source-receiver to estimation by using the information obtained in inverse migration step
The transmitting of seismic wave and reception slope, to illuminate area-of-interest.It in the transmission, can be by using several sources or a source day
Line is by the transmitting of wave or beam control focusing in the slope and optimum orientation for being used to illuminate area-of-interest.Therefore, it avoids
Many noise factors, such as diffraction, parasitism and multiple echo and parasitic reflection, and improve the threshold value of detectivity.?
It,, can be with by average one group of Proximity Sensor or by using receiver antenna by using three-component receiver in reception
Obtain average trace, which only includes earthquake information, earthquake information according to the direction defined in inverse migration step with
Slope reaches.The filtering allows to significantly reduce noise factor, especially multiple, and can increase the threshold value of detectivity.This
It needs to improve detection accuracy by reducing the search coverage in Fresnel region, so that the slope by adjusting wave beam control comes more
Focusing-detection well.Therefore, using single source-receiver pair, several different focusing can be obtained in the Fresnel region, this
It is identical from several different points are distinguished in this region.The two filtering can be realized independently or together, with further
Improve the threshold value of detectivity.It, can be by modifying transmitting according to frequency multiplication and receiving if carrying out inverse migration with octave coverage
Slope is further improved the filtering.Therefore, the present invention is it can be concluded that a kind of detection conclusion, on the one hand which considers sense emerging
On the other hand the variation that interesting region occurs at any time relative to reference zone considers the time change of these variations.
It can also implement the present invention to record complete trace, then the complete trace can be used for identifying and be located at initially
Other points for being referred to as chance above or below point, and can be for the focusing in each point adjustment Fresnel region.
Claims (10)
1. a kind of method for being determined to obtain source-receiver pair best orientation of seismic data, which is characterized in that packet
Include following steps:
First step identifies an area-of-interest (32) that the area-of-interest is the main body of early earthquake acquisition, to obtain
Obtain the image of the subsoil of the area-of-interest；
Second step obtains the ground obtained during the early earthquake acquisition of the area-of-interest within the interested time
Shake data；
Third step, using the inverse migration of the seismic data, position with each source-receiver of determination to (31,34), institute
It states each source-receiver and tribute has been made to the image of the subsoil of the area-of-interest within the interested time
It offers；
- the four steps obtains the source-receiver to the untreated trace (41,42,43,44,45) of position (31,34)；
- the five step selects at least one best untreated trace (43,44) from a plurality of untreated trace；
- the six step, determine correspond to best source-receiver of the best untreated trace (43,44) to position (31,
34)。
2. the method according to claim 1, wherein the inverse migration of third step is carried out by octave coverage.
3. method according to any of the preceding claims, which is characterized in that the early earthquake acquisition is 4D acquisition.
4. according to the method described in claim 3, it is characterized in that, the 5th step further includes the optimal untreated mark of selection
Line (43,44) allows to detect best 4D effect.
5. method according to any of the preceding claims, which is characterized in that the 5th step further includes using stone
Oily elastic model is to select optimal untreated trace (43,44).
6. method according to any of the preceding claims, which is characterized in that the 5th step further includes consideration table
Face factor.
7. method according to any of the preceding claims, which is characterized in that the 5th step further includes considering to use
Make the reference zone (55) of earthquake variation calibration.
8. method according to any of the preceding claims, which is characterized in that the third step further includes obtaining often
A source-receiver is to the best slope of the seismic wave of (31,34) and transmitting and receiving direction.
9. method according to any of the preceding claims, which is characterized in that the method also includes by source-receiver
7th step of the optimum position is placed on to (31,34) and obtains the 8th step of new seismic data in the optimum position
Suddenly.
10. according to the method described in claim 9, it is characterized in that, the 8th step is based on the best oblique of the seismic wave
Rate and transmitting and receiving direction.
Applications Claiming Priority (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1661842 | 2016-12-02 | ||
| FR1661842A FR3059700A3 (en) | 2016-12-02 | 2016-12-02 | METHOD OF IMPROVING THE ULTRALIGHT ACTIVE SEISMIC DETECTION SYSTEMS |
| FR1757024 | 2017-07-25 | ||
| FR1757024A FR3059786B1 (en) | 2016-12-02 | 2017-07-25 | METHOD FOR IMPROVING SEISMIC ACQUISITIONS USING ULTIMATE ACTIVE SEISMIC DETECTION SYSTEMS |
| PCT/EP2017/080582 WO2018099880A1 (en) | 2016-12-02 | 2017-11-28 | Method for improving seismic acquisitions utilising active ultralight seismic detection systems |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110050204A true CN110050204A (en) | 2019-07-23 |
| CN110050204B CN110050204B (en) | 2021-12-03 |
Family
ID=59974644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201780073625.2A Active CN110050204B (en) | 2016-12-02 | 2017-11-28 | Method for improving seismic acquisition using active ultra-light seismic detection system |
Country Status (9)
| Country | Link |
|---|---|
| US (1) | US11435492B2 (en) |
| EP (1) | EP3548930A1 (en) |
| JP (1) | JP7404069B2 (en) |
| CN (1) | CN110050204B (en) |
| CA (1) | CA3082926A1 (en) |
| FR (2) | FR3059700A3 (en) |
| RU (1) | RU2751573C2 (en) |
| SA (1) | SA519401832B1 (en) |
| WO (1) | WO2018099880A1 (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4188610A (en) * | 1977-08-29 | 1980-02-12 | Hydroacoustics, Inc. | Method of and apparatus for the generation and transmission of signals for echolocation and other signalling purposes, such as in geophysical exploration |
| US20100080081A1 (en) * | 2008-09-26 | 2010-04-01 | Providence technologies, Inc. | Method and apparatus for seismic exploration |
| CN102282481A (en) * | 2009-01-19 | 2011-12-14 | 兰德马克图形公司 | data acquisition and prestack migration based on seismic visibility analysis |
| CN102944897A (en) * | 2012-11-20 | 2013-02-27 | 成都晶石石油科技有限公司 | Correction method for sea well shock speed scissors difference based on standard reference layer |
| US20160162613A1 (en) * | 2014-12-09 | 2016-06-09 | Schlumberger Technology Corporation | Method for disposing seismic signal receivers for a seismic acquisition system |
| US20160161620A1 (en) * | 2013-06-28 | 2016-06-09 | Cgg Services Sa | System and method for estimating repeatability using base data |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6885918B2 (en) * | 2000-06-15 | 2005-04-26 | Geo-X Systems, Ltd. | Seismic monitoring and control method |
| FR2916859B1 (en) * | 2007-05-31 | 2009-08-21 | Cgg Services Sa | METHOD OF PROCESSING SEISMIC DATA |
| JP5446101B2 (en) | 2008-03-03 | 2014-03-19 | 株式会社ニコン | Electronics |
| US8339898B2 (en) * | 2008-05-25 | 2012-12-25 | Westerngeco L.L.C. | Processing seismic data using combined regularization and 4D binning |
| EP2376945A4 (en) * | 2008-12-17 | 2017-02-22 | Exxonmobil Upstream Research Company | System and method for performing time-lapse monitor surveying using sparse monitor data |
| US8553496B2 (en) * | 2010-02-09 | 2013-10-08 | Ion Geophysical Corporation | Seismic source separation |
| US10088588B2 (en) * | 2013-04-03 | 2018-10-02 | Cgg Services Sas | Device and method for stable least-squares reverse time migration |
| WO2016038466A1 (en) * | 2014-09-08 | 2016-03-17 | Cgg Services Sa | Seismic migration using an indexed matrix |
| JP2016085096A (en) * | 2014-10-24 | 2016-05-19 | 独立行政法人石油天然ガス・金属鉱物資源機構 | Data processing method, data processing device, and data processing program |
| EP3171203B1 (en) * | 2015-11-18 | 2019-01-02 | CGG Services SAS | Adaptive ensemble-based method and device for highly-nonlinear problems |
-
2016
- 2016-12-02 FR FR1661842A patent/FR3059700A3/en active Pending
-
2017
- 2017-07-25 FR FR1757024A patent/FR3059786B1/en active Active
- 2017-11-28 WO PCT/EP2017/080582 patent/WO2018099880A1/en unknown
- 2017-11-28 EP EP17822555.3A patent/EP3548930A1/en active Pending
- 2017-11-28 RU RU2019118875A patent/RU2751573C2/en active
- 2017-11-28 CN CN201780073625.2A patent/CN110050204B/en active Active
- 2017-11-28 US US16/465,669 patent/US11435492B2/en active Active
- 2017-11-28 CA CA3082926A patent/CA3082926A1/en active Pending
- 2017-11-28 JP JP2019529198A patent/JP7404069B2/en active Active
-
2019
- 2019-05-21 SA SA519401832A patent/SA519401832B1/en unknown
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4188610A (en) * | 1977-08-29 | 1980-02-12 | Hydroacoustics, Inc. | Method of and apparatus for the generation and transmission of signals for echolocation and other signalling purposes, such as in geophysical exploration |
| US20100080081A1 (en) * | 2008-09-26 | 2010-04-01 | Providence technologies, Inc. | Method and apparatus for seismic exploration |
| CN102282481A (en) * | 2009-01-19 | 2011-12-14 | 兰德马克图形公司 | data acquisition and prestack migration based on seismic visibility analysis |
| CN102944897A (en) * | 2012-11-20 | 2013-02-27 | 成都晶石石油科技有限公司 | Correction method for sea well shock speed scissors difference based on standard reference layer |
| US20160161620A1 (en) * | 2013-06-28 | 2016-06-09 | Cgg Services Sa | System and method for estimating repeatability using base data |
| US20160162613A1 (en) * | 2014-12-09 | 2016-06-09 | Schlumberger Technology Corporation | Method for disposing seismic signal receivers for a seismic acquisition system |
Also Published As
| Publication number | Publication date |
|---|---|
| US11435492B2 (en) | 2022-09-06 |
| WO2018099880A1 (en) | 2018-06-07 |
| FR3059700A3 (en) | 2018-06-08 |
| JP2019536050A (en) | 2019-12-12 |
| JP7404069B2 (en) | 2023-12-25 |
| EP3548930A1 (en) | 2019-10-09 |
| RU2019118875A (en) | 2021-01-14 |
| RU2751573C2 (en) | 2021-07-15 |
| SA519401832B1 (en) | 2022-10-25 |
| FR3059786B1 (en) | 2019-07-05 |
| CA3082926A1 (en) | 2018-06-07 |
| FR3059786A1 (en) | 2018-06-08 |
| US20190302299A1 (en) | 2019-10-03 |
| CN110050204B (en) | 2021-12-03 |
| RU2019118875A3 (en) | 2021-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Foti et al. | Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project | |
| AU2015218416B2 (en) | Impact assessment of marine seismic surveys | |
| EP3059615A1 (en) | Amplitude-versus-angle analysis for quantative interpretation | |
| CN103582826A (en) | Determining an indication of wavefield velocity | |
| EP2743734A2 (en) | Target-oriented 4D binning in common reflection point | |
| AU2013237741A1 (en) | Methods and systems for monitoring a petroleum reservoir | |
| US9759828B2 (en) | Determining a streamer position | |
| US20140297190A1 (en) | Monitoring of source signature directivity in seismic systems | |
| Zeng et al. | Turning a telecom fiber‐optic cable into an ultradense seismic array for rapid postearthquake response in an urban area | |
| EP2743737A2 (en) | Methods and systems for quality control of seismic illumination maps | |
| US20240069234A1 (en) | Seismic survey data access | |
| US20220026595A1 (en) | Systems and methods of performing velocity surveys using spaced source activation lines | |
| CN110050204A (en) | A method of earthquake-capturing is improved using active ultralight seismic acquisition system | |
| US20220236435A1 (en) | Low-Frequency Seismic Survey Design | |
| US20120081996A1 (en) | Determining Sea Conditions in Marine Seismic Spreads | |
| US9971053B2 (en) | Using crossline measurement data for an action relating to survey of a target structure | |
| Xie et al. | Integrating distributed acoustic sensing and computer vision for real-time seismic location of landslides and rockfalls along linear infrastructure | |
| Song | Locating Petroleum Sources Using Dsp Techniques | |
| Zarifi et al. | Development of a baseline seismic monitoring network for the Horda Platform region in support of CO2 storage at the Northern Lights project | |
| Share | Multi-Scale Imaging of Major Fault Zones in Southern California | |
| EA043815B1 (en) | LOW FREQUENCY SEISMIC SURVEY DESIGN | |
| Rentsch | A migration-type approach for the fast location of seismicity: theory and applications | |
| Dreger et al. | Rupture Process of the 2004 Parkfield Earthquake Utilizing Near-Fault Seismic Records and Implications for ShakeMap | |
| Tselentis et al. | Case History High-resolution passive seismic tomography for 3D velocity, Poisson’s ratio m, and P-wave quality QP in the Delvina hydrocarbon field, southern Albania | |
| Templeton et al. | Mapping Fracture Network Creation with Microseismicity During EGS Demonstrations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |