CN1603860A - Method for seeking lithological anomalous body using orthogonal seismic trace data - Google Patents
Method for seeking lithological anomalous body using orthogonal seismic trace data Download PDFInfo
- Publication number
- CN1603860A CN1603860A CN 200410088559 CN200410088559A CN1603860A CN 1603860 A CN1603860 A CN 1603860A CN 200410088559 CN200410088559 CN 200410088559 CN 200410088559 A CN200410088559 A CN 200410088559A CN 1603860 A CN1603860 A CN 1603860A
- Authority
- CN
- China
- Prior art keywords
- seismic
- orthogonal
- seismic trace
- data
- amplitude
- 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 37
- 230000002547 anomalous effect Effects 0.000 title claims description 44
- 239000000284 extract Substances 0.000 claims description 11
- 239000004215 Carbon black (E152) Substances 0.000 claims description 9
- 229930195733 hydrocarbon Natural products 0.000 claims description 9
- 238000013508 migration Methods 0.000 claims description 6
- 230000005012 migration Effects 0.000 claims description 6
- 230000000052 comparative effect Effects 0.000 claims description 5
- 230000005856 abnormality Effects 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 125000001183 hydrocarbyl group Chemical group 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 239000003208 petroleum Substances 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 1
- 230000002159 abnormal effect Effects 0.000 abstract description 4
- 239000011435 rock Substances 0.000 abstract description 4
- 239000000463 material Substances 0.000 abstract 3
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000005755 formation reaction Methods 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 239000004576 sand Substances 0.000 description 5
- 239000010426 asphalt Substances 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000010363 phase shift Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000002405 diagnostic procedure Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Abstract
This invention relates to a gas oil field exploration method and belongs to oil exploration two-dimension earthquake material explanation technique field. The method comprises the following steps: first to figure out the cross earthquake channel; second to compare the explanation of two dimension earthquake; third to extract average vibration and form plane graph and display vibration abnormal area. It can identify the rock abnormal body according to regular excursion earthquake materials and crossing earthquake channel material vibration information difference and the longitudinal resolution rate of the two amounts to one-eighth lamda.
Description
Technical field
The present invention relates to the seismic exploration technique field of oil-gas exploration and development, particularly relate to the method that the conventional 2-d seismic data of a kind of usefulness is sought lithologic anomalous body.Conventional seismic data is converted to the orthogonal seismic trace data to seek lithologic anomalous body.This method belongs to petroleum prospecting 2-d seismic data interpretation technique field (IC: G01V 1/28).
Background technology
At present, utilize seismic exploration technique to seek lithologic anomalous body and mainly contain three kinds of methods: 1. seismic reflection pattern discrimination method: the SEQUENCE STRATIGRAPHIC model prediction sand body of setting up according to certain district can developable place and type, seismic phase analysis in combination then is according to the existence of seismic reflection structure identification particular lithologic body on the section; 2. earthquake-stratigraphic model diagnostic method: in structure, fracture complex area, when being difficult on the seismic section determine the horizontal change of tectonic structure and lithology, can set up geologic model earlier, import different parameters then and obtain theoretical seismic section, and with real seismic record compare with the checking geologic model correctness, and select reasonable earthquake interpretation scheme, to seek possible lithologic anomalous body; 3. seismic inversion stratigraphic section diagnostic method: this method computing machine treatment technology, to reflect that the earthquake reflection section of (interface) feature is directly changed into the wave impedance section of reflection subsurface formations feature, in conjunction with the on-the-spot actual well drilled result contrast that makes an explanation, thereby reach the purpose of seeking lithologic anomalous body.The thickness that these three kinds of methods all require lithologic anomalous body just has operability during more than or equal to 1/4 λ.And, on conventional seismic section or wave impedance section, can not produce distinguishable seismic response when the thickness of lithologic anomalous body during (the minimum thickness limit that vertically can differentiate of conventional earthquake is 1/4 λ, and λ is the seismic event wavelength) less than 1/4 λ.In the two-dimension earthquake exploration, echo amplitude is a very important parameters, and it is the important evidence of carrying out structure elucidation and oil-gas recognition (comprising bright spot, dim spot, flat spot etc.) on conventional seismic section.Seek lithologic anomalous body also based on the variation of amplitude parameter with orthogonal seismic trace.The superiority of orthogonal seismic trace be it can give prominence to a certain extent the reflection the unresolvable part of conventional 2-d seismic data (promptly being in the amplitude information of zero crossing location fuzzy).Therefore, orthogonal seismic trace and conventional two-dimensional migration seismic data combined compares explanation, and the comprehensive longitudinal frame of the two can reach 1/8 λ in theory.In real work, sometimes thickness is less than 1/4 λ but greater than having comprised the good important oil gas information such as sand body or hydrocarbon zone of preserving in the stratum of 1/8 λ, owing on conventional seismic section, can not discern and differentiate, thereby can miss possible lithologic deposit.
Summary of the invention
The present invention seeks to: a kind of method of seeking lithologic anomalous body with the orthogonal seismic trace data is provided, the thickness that makes anomalous body is less than 1/4 λ but more than or equal to 1/8 λ the time, still can carry out integrated interpretation in conjunction with conventional 2-d seismic data with orthogonal seismic trace, to discern the formed lithologic anomalous body of possible sand body or hydrocarbon zone.Overcome 2-d seismic data in lithologic anomalous body thickness unresolvable problem during less than 1/4 λ but more than or equal to 1/8 λ.This method can disclose abundanter formation information, can catch thickness less than 1/4 λ but more than or equal to hydrocarbon zone and the formed lithologic anomalous body of country rock of 1/8 λ, to find oil gas field.
Orthogonal seismic trace can be given prominence to the principle of differentiating fuzzy message in the conventional skew seismic data:
Routine is offset seismic data carries out 90 ° of phase-shift filterings when obtaining orthogonal seismic trace, seismic wavelet becomes asymmetric wavelet by zero-phase wavelet, and the difference of wavelet waveform is that orthogonal seismic trace can be given prominence to the key of differentiating fuzzy message in the conventional skew seismic data.Usually, a zero-phase wavelet is sinusoidal wave synthetic with many zero phases that equal peak amplitude is arranged.Keeping to change the shape of wavelet with the modification phase spectrum under the constant situation of spectral amplitude and frequency content.Wherein, linear phase is moved a constant time shift that is equal to wavelet, and the waveform of wavelet is constant.Constant phase is moved the shape that has then changed wavelet.Wherein 90 ° of phase shifts are converted to asymmetric wavelet (accompanying drawing 1 (a) (b)) with zero-phase wavelet.A real seismic trace through 90 ° of phase-shift filterings after, become its orthogonal seismic trace.Crest and trough in routine skew seismic data, vanishing intersects (zero crossing) on its orthogonal seismic trace, and the zero crossing in routine skew seismic data becomes crest and trough on its orthogonal seismic trace.Therefore orthogonal seismic trace can be given prominence to dim, the fuzzy amplitude information (accompanying drawing 2) that is in the zero crossing position in the conventional skew of the reflection seismic data.This is prerequisite and the basis of seeking lithologic anomalous body with orthogonal seismic trace.Orthogonal seismic trace and conventional skew seismic data are compared explanation, and the comprehensive longitudinal frame of the two can reach 1/8 λ in theory.
The technical scheme that the present invention is taked for its technical matters of solution is: seek the method for lithologic anomalous body based on the orthogonal seismic trace of conventional two-dimensional migration seismic trace, it is characterized in that this method mainly is made up of three parts:
1. calculate orthogonal seismic trace: open the menu that calculates orthogonal seismic trace, input two-dimension earthquake data obtain the orthogonal seismic trace data volume.
2. the two-dimension earthquake amplitude contrasts alternately: at the seismic profile interpretation window, section is shown by the variable density mode, the principle of selection or editor's colour code is to come amplitudes different on the representative profile with abundant colors, open two seismic profile interpretation windows, the two display parameter unanimity, seismic data volume is selected two-dimensional migration geological data and orthogonal seismic trace data respectively, on the orthogonal seismic trace section, select a phase place that has amplitude anomaly, contrast to two-dimensional quadrature seismic trace section as a seismic horizon, and is launched on this basis in the center line of this phase place and suitable position (when not existing unusually) thereof.
3. extract average amplitude one-tenth figure-demonstration amplitude anomaly district: open the menu that extracts average amplitude, the window of time interval unanimity is as the parameter of extracting average amplitude about the selection seismic horizon, the average amplitude that extracts is exported, in mapping softwares such as two fox softwares or Earthvision, the average amplitude value of this layer is carried out gridding, in planimetric map, demonstrate the amplitude anomaly district, this amplitude anomaly district is the distribution range of lithologic anomalous body, and whether this anomalous body is that the geological property of oil gas need be in conjunction with the well data, areal geology data or be aided with other means and comprehensively prove.
Further describe better scheme of the present invention below, illustrate again with three parts in landmark workstation seismic interpretation system, implementing this method:
1., calculate orthogonal seismic trace:orthogonal seismic trace is offset seismic data based on routine. is by landmark main menu (OpenWorks 2003 Command Menu) → Applications → PostStack/PAL → selection two-dimension earthquake work area → press Launch ... key → enter PostStack Family 2003{ two-dimension earthquake work area name } main menu → select input two-dimension earthquake data (segy form) or the two-dimension earthquake data (2v2 file) that generated at Input Data → create export file name (being the Vertical form of 9 characters) → Processes → Attributes → Instantaneous → Quadrature Trace → operation Run at Output Data; The orthogonal seismic trace of namely finishing after the end conventional migration technique seismic data (segy formatted file) calculates, or the shake data on the spot that existed (2v2 formatted file) are converted to orthogonal seismic trace data (2v2 formatted file).
2., the mutual comparative interpretation of two-dimension earthquake: enter SeisWorks2003 by landmark master menu (OpenWorks 2003Command Menu) → Applications → Seisworks → 2D, select the two-dimension earthquake work area.At the seismic profile interpretation window, seismic section shows by variable density, vertically between the optional 15.49cm/second-30.98cm/second of ratio, between the optional 32.28trace/cm-10.76trace/cm of grid scale, so that the different variations of amplitude on clearly showing the orthogonal seismic trace section and being offset the seismic trace section.It is principle (or select in the colour code dipazim) that the selection of colour code, editor show amplitudes different on the section with abundant colors as far as possible.Open two seismic profile interpretation windows, screen about placing respectively, the display parameter unanimity of the two.Seismic data volume selects to be offset geological data and orthogonal seismic trace data respectively.The amplitude of objective interval on conventional offset seismic profile and its orthogonal seismic trace section is observed, contrasted, on the orthogonal seismic trace section, select the phase place (crest or trough) that might have amplitude anomaly, and the medium line of this phase place as a seismic horizon h, suppose that its double-pass reflection time is the T0 millisecond, be approximately the T0+10 millisecond in the double-pass reflection time of the enterprising row amplitude equivalent beds position hl of conventional offset seismic profile (medium line of corresponding crest or trough) so.And these two layers are carried out comparative interpretation on section.If the two does not transversely have tangible amplitude difference, it is little to illustrate that then part that near the more conventional seismic trace of orthogonal seismic trace this place can give prominence to resolution and conventional seismic trace can be differentiated formation lithology horizontal change partly, if the amplitude of the two has evident difference, such as a certain position on offset seismic profile is the continuous reflection of a trough, and be the interrupted reflection of a trough in the relevant position of its orthogonal seismic trace section or disappear to a certain direction, we just think have amplitude anomaly on orthogonal seismic traces, can the more conventional seismic section of side light orthogonal seismic trace section can give prominence to the formation lithology of differentiating part and transversely exist variation.Only can't discern this lithologic body with offset seismic profile, and the orthogonal seismic trace section can give prominence to the reflection it.The center line of the anomalous amplitude on the orthogonal seismic trace section and the position (when unusually not existing) suitable with it are made an explanation as a seismic horizon, and on the two-dimension earthquake survey line, follow the trail of contrast, can demonstrate this unusual general plane scope at the SeisWorks/Map window.
3., extract the average amplitude parallel planes one-tenth figure that goes forward side by side: enter StratAmp 2003 master menus that extract average amplitude by landmark master menu (OpenWorks2003 Command Menu) → Data → Management → Seismic Project Manager → Horizons, select two-dimension earthquake work area and orthogonal seismic trace data, select the seismic horizon of input and the average amplitude layer position of name output again, calculate the input layer position average amplitude value of window (time window more little, Ke Hua amplitude anomaly distributive province is accurate more in the plane) sometime in each 1 millisecond of-15 milliseconds of scope up and down.Average amplitude is exported from workstation as a layer position, average amplitude value this layer in two fox softwares or Earthvision mapping software is carried out gridding, the plane distributing scope of amplitude can show abnormality, determine the geological property that this is unusual in conjunction with well data and areal geology data again, or be aided with additive method and verify.
Beneficial effect of the present invention is: in seismic prospecting, echo amplitude is a very important parameters.Layer thickness up and down more than or equal to the vertical distinguishable prerequisite of 1/4 λ earthquake under, on conventional seismic section, can form the wave impedance interface, show as crest or trough.And, on conventional seismic section, can not produce distinguishable seismic response when the thickness of lithology pinching body that exists in the stratum or hydrocarbon zone during greater than 1/8 λ but less than 1/4 λ.Though its vertical distinguishable thickness limit of simple orthogonal seismic trace also is 1/4 λ, because orthogonal seismic trace can be given prominence to the amplitude information of the conventional seismic trace zero crossing location fuzzy of reflection.Therefore, the two combines orthogonal seismic trace and conventional seismic trace, and its comprehensive longitudinal frame can reach 1/8 λ in theory.The amplitude of equivalent layer position in the horizontal change of orthogonal seismic trace layer amplitude and the conventional seismic data is compared explanation, (so-called " amplitude anomaly " is just for conventional offset seismic profile to find out amplitude anomaly on the orthogonal seismic trace section, the width of its phase place is not represented the thickness of possible lithologic anomalous body, but the power of its reflection can reflect the horizontal change of formation lithology difference up and down indirectly).Amplitude anomaly district on its plane is the distribution range of lithologic anomalous body.The present invention is a foundation with this principle, realizes determining the distribution range of lithologic anomalous body, the hydrocarbon-bearing pool that searching may exist.So this method has following effect:
1, discerns lithologic anomalous body according to the difference of amplitude information on routine skew seismic data and its orthogonal seismic trace data.Be in zero crossing location fuzzy, dim amplitude information because orthogonal seismic trace can be given prominence to the conventional seismic data of resolution, thereby can disclose abundanter formation information.The two comprehensive longitudinal frame can reach 1/8 λ in theory.When lithologic anomalous body thickness during less than 1/4 λ but greater than 1/8 λ, on conventional seismic section, can not produce distinguishable seismic response, reflect but information that can it is fuzzy, dim on the orthogonal seismic trace section is outstanding.The amplitude of relevant layers position on conventional offset seismic profile and its orthogonal seismic trace section is compared explanation, seek possible amplitude anomaly on the orthogonal seismic trace section, the amplitude anomaly district on its plane is exactly the distribution range of possible lithologic anomalous body.In real work, sometimes less than 1/4 λ but greater than may having comprised important oil gas information such as preserving sand body or hydrocarbon zone in the stratum of 1/8 λ, owing on conventional seismic section, can not discern and differentiate, thereby can miss possible lithologic deposit.Then overcome the conventional 2-d seismic data of simple usefulness in conjunction with the orthogonal seismic trace data and can not differentiate thickness less than 1/4 λ but greater than the deficiency on 1/8 λ stratum.
2, under the situation that conventional skew seismic data is only arranged, at favourable objective interval, utilize way of the present invention can carry out the search of extensive lithologic anomalous body, speed is fast, and is with low cost.
3, this method is easy to implement, and is easy to grasp, and simple to operate, practical value is big.To hydrocarbon zone thickness less than 1/4 λ but, have great importance greater than the exploration of the lithologic deposit of 1/8 λ.
Description of drawings
Figure below in the accompanying drawing 1 (a) is the waveform of zero phase seismic wavelet, and its last figure is the phase spectrum of this seismic wavelet.Figure below in the accompanying drawing 1 (b) is the waveform of asymmetric seismic wavelet, and its last figure is the phase spectrum of this seismic wavelet.Two width of cloth figure compare, the waveform of the seismic wavelet of the zero phase in the key diagram (a) and its phase spectrum, and this wavelet becomes asymmetric wavelet through quadrature filtering, and change has all taken place in its waveform and phase spectrum, becomes wavelet waveform and phase spectrum among the figure (b).
Accompanying drawing 2 is orthogonal seismic traces and shakes trace comparison figure on the spot.A is that orthogonal seismic trace can be given prominence to the part of differentiating real seismic trace among the figure; B is maximum positive amplitude; C is the minimal negative amplitude.
Accompanying drawing 3 is a section MD04-152 in Man Nan area, Tarim Basin, and this section is crossed 32 wells in the tower, and 34 wells are projection (apart from this survey line 5.7km) in the tower.It among the figure Silurian inside story reflection contrast situation, as can be seen in conventional offset seismic profile, in the tower the blue phase place (trough) of 32,34 wells be connect together (with this layer called after earthquake Tg42 reflection horizon, be in the center line of this blueness trough, in order to show the variation of blue phase place, do not mark concrete position.It is the equivalent beds position in earthquake Tg43 reflection horizon on the orthogonal seismic trace section, and, from seismic reflection two-way time Tg43 ≈ Tg42-10 millisecond), there is not any difference, and on its orthogonal seismic trace section, 32 wells are in the edge of blue phase place in the tower, and this well sandstone section under Silurian is found the show of oil and gas of 1.1m, through fixed at the enterprising rower of orthogonal seismic trace section with VSP speed, its position is just on blue phase place.34 wells corresponding display segment of 32 well Silurian in tower does not have show of oil and gas in the tower, the position that is equivalent to 34 well projections in the tower on the orthogonal seismic trace section of MD04-152 does not have the abnormal show of blue phase place yet, from the real brill data of two mouthfuls of wells, 32 well show of oil and gas sections are delta river mouth dam deposition in the tower, rerum natura is better, 34 wells are delta-plain facies in the corresponding well section in the tower, and rerum natura is relatively poor.Can think, be to cause the main cause of two mouthfuls of well abnormal show differences on the orthogonal seismic trace section in the physical difference of Silurian inside story interval and the difference of oil-gas possibility.
Accompanying drawing 4 is full south regional Silurian inside story amplitude anomaly flat distribution maps.This figure is Silurian inside story earthquake Tg43 reflection horizon 5 milliseconds of average amplitude plane variation diagrams up and down, and the scope of darker curve trap is the amplitude anomaly district of Silurian inside story among the figure, also is the scope of the favourable lithologic anomalous body of Silurian inside story.The area of this exceptions area is 1400 square kilometres.Can find out obviously that from figure 32 wells are in the edge in this amplitude anomaly district in the tower, 34 wells then are in outside this exceptions area in the tower.
Embodiment
Embodiment: use based on orthogonal seismic trace 32 wellblocks in tower of routine skew seismic data and found Silurian inside story amplitude anomaly district, the i.e. distribution range of lithologic anomalous body.
32,34 wellblock Silurian are the monocline by north orientation south lifting in the tower of Tarim Basin, the arch structure agensis.Nineteen ninety-five, to degrade trap with Silurian asphalt sandstones be target to 32 wells in the tower of this district probing, and this well has been found 1.1 meters show of oil and gas in Silurian 3791.71-3792.81 well section.Wherein, 0.77 meter of oil immersion, 0.33 meter at oil stain.Drilled 34 wells in the tower in the high-order bit that this well southeast 5km Silurian asphalt sandstones is degraded trap again in 1 year, the result does not find any show of oil and gas.From the results of drilling of two mouthfuls of wells, Silurian is degraded trap and is lacked effective cap rock.At the updip direction of MD04-152 survey line Silurian, Silurian is finally degraded totally by superstratum-Carboniferous system, and Carboniferous system bottom does not have effective cap rocks such as mud stone, and the confirmation Silurian asphalt sandstones is degraded trap and can not be formed effective hydrocarbon trap.The show of oil and gas of 32 wells should be a part that is subjected to the formed lithologic deposit of lithologic character controlling in the tower.Though on area deposition angle and geology aspect, this district is had higher rating, on conventional seismic section, be difficult to find and restricted the exploration progress in this district for the target of probing.Disposed 31 two-dimension earthquake surveys line in 2004 in this district, show to implement three Silurian inside story stratigraphic onlap trap made from old survey line of the past, through conventional skew seismic data in 2004 is made an explanation into figure, confirm that the trap of three Silurian inside story onlaps does not all exist.Afterwards, the author has identified the amplitude anomaly district (area is 1400 square kilometres) of Silurian inside story with orthogonal seismic trace data 32 well Silurian downdip directions in tower, be the distribution range of Silurian inside story lithologic anomalous body, it is the favourable exploration target area of Silurian inside story lithologic deposit.Actual mechanical process is:one, calculate orthogonal seismic trace:orthogonal seismic trace is offset seismic data based on routine. is by landmark main menu (OpenWorks 2003 Command Menu) → Applications → PostStack/PAL → selection mn2d two-dimension earthquake work area → press Launch ... key → enter PostStack Family 2003{mn2d} main menu → select the two-dimension earthquake data file mn2d.2v2 that input generated → create export file name quadratur → Processes → Attributes → Instantaneous → Quadrature Trace → operation Run at Output Data at Input Data finishes namely the shake data on the spot that existed (2v2 formatted file) to be converted to orthogonal seismic trace data (quadratur.2v2 formatted file) afterwards.
Two, the mutual comparative interpretation of two-dimension earthquake: advance SeisWorks2003 by landmark master menu (OpenWorks 2003Command Menu) → Applications → Seisworks → 2D, select mn2d two-dimension earthquake work area.At the seismic profile interpretation window, seismic section shows that by variable density vertically ratio is selected 23.23cm/second, and grid scale is selected 32.28trace/cm, and the colour code that the local area section shows is selected dipazim.Open two seismic profile interpretation windows, screen about placing respectively, the display parameter unanimity of the two.Seismic data volume selects to be offset geological data mn2d.2v2 and orthogonal seismic trace data quadratur.2v2 respectively.The amplitude of objective interval on conventional offset seismic profile and its orthogonal seismic trace section is observed, contrasted, select a phase place trough (its seismic reflection layer is Tg43) that has amplitude anomaly in the reflection of orthogonal seismic trace section Silurian inside story, this layer double-pass reflection time of 32 Jing Chu in tower is 2900 milliseconds, and the double-pass reflection time in the layer position (seismic reflection layer is Tg42) that the enterprising row amplitude of conventional offset seismic profile contrasts is 2910 milliseconds.And these two layers are carried out comparative interpretation on section.The latter on the conventional offset seismic profile of MD04-152 in the tower in 32 wells and the tower position of 34 wells (projection) be a continuous trough reflection, and degraded by the Carboniferous system at updip direction.And be an interrupted trough reflection in the respective layer position of MD04-152 orthogonal seismic trace (Tg43), then there is tangible amplitude anomaly with respect to conventional migrated section, updip direction in Silurian, 32 wells are in this unusual marginal position in the tower, and then there is not amplitude anomaly in projection place of 34 wells in being equivalent to tower.We think the amplitude anomaly that exists on orthogonal seismic trace, can the more conventional seismic section of side light orthogonal seismic trace section can give prominence to the formation lithology of differentiating part and transversely exist variation.Only can't discern this lithologic body with offset seismic profile, and the orthogonal seismic trace section can give prominence to the reflection it.The center line of the anomalous amplitude on the orthogonal seismic trace section is made an explanation as a seismic horizon (earthquake Tg43 reflection horizon), and on the two-dimension earthquake survey line, follow the trail of contrast, can demonstrate the distribution situation of this amplitude anomaly on two-dimentional survey line at the SeisWorks/Map window.
Three, extract the average amplitude parallel planes one-tenth figure that goes forward side by side: enter StratAmp 2003 master menus that extract average amplitude by landmark master menu (OpenWorks2003 Command Menu) → Data → Management → Seismic Project Manager → Horizons, select two-dimension earthquake work area and orthogonal seismic trace data quadratur.2v2, select the earthquake Tg43 reflection horizon of input and the average amplitude layer position Tg43_5 of name output (being Tg43 reflection horizon 5 milliseconds average amplitude up and down) again, calculate input layer position each average amplitude value of 5 milliseconds up and down.This average amplitude is exported from workstation as a layer position, average amplitude value the Tg43_5 layer in the Earthvision mapping software is carried out gridding, the plane distributing scope of amplitude can show abnormality, determine the geological property that this is unusual in conjunction with well data and areal geology data again, or be aided with the lithologic deposit that additive method verifies that 32 wellblocks may be missed in the tower of Tarim Basin.
On MD04-152 orthogonal seismic trace section, 32 well VSP speed are demarcated in tower, have comprised 1.1 meters show of oil and gas section in the amplitude anomaly of this well place Silurian inside story, and the reflection that further specifies this unusual and oil gas or well preserve sand body matches.Through on the orthogonal seismic trace section, this amplitude anomaly being followed the trail of, find that this exceptions area contour area in the plane reaches 1400km
2, be likely a large-scale lithologic deposit of Silurian inside story to have good exploration prospect.
Claims (9)
1. one kind is applied to method petroleum exploration field, seek lithologic anomalous body based on the orthogonal seismic trace of conventional two-dimensional migration seismic trace, seek zone thickness less than 1/4 λ (λ is the seismic event wavelength) but more than or equal to the lithologic anomalous body of 1/8 λ, find lithologic anomalous body to find possible hydrocarbon-bearing pool, this method goes up in workstation seismic interpretation system (software) and realizes, it is characterized in that: this method mainly is made up of three parts
1. calculate orthogonal seismic trace: open the menu that calculates orthogonal seismic trace, input two-dimension earthquake data obtain the orthogonal seismic trace data volume;
2. the two-dimension earthquake amplitude contrasts alternately: at the seismic profile interpretation window, section is shown by the variable density mode, the principle of selection or editor's colour code is to come amplitudes different on the representative profile with abundant colors, open two seismic profile interpretation windows, the two display parameter unanimity, seismic data volume is selected two-dimensional migration geological data and orthogonal seismic trace data respectively, on the orthogonal seismic trace section, select a phase place that has amplitude anomaly, contrast to two-dimensional quadrature seismic trace section as a seismic horizon, and is launched on this basis in the center line of this phase place and suitable position (when not existing unusually) thereof;
3. extract average amplitude one-tenth figure-demonstration amplitude anomaly district: open the menu that extracts average amplitude, the window of time interval unanimity is as the parameter of extracting average amplitude about the selection seismic horizon, the average amplitude that extracts is exported, in mapping softwares such as two fox softwares or Earthvision, the average amplitude value of this layer is carried out gridding, in planimetric map, demonstrate the amplitude anomaly district, this amplitude anomaly district is the distribution range of lithologic anomalous body, determines the oily attribute of this exceptions area then in conjunction with well data and areal geology data.
2, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, it is characterized in that: described calculating orthogonal seismic trace is conventional skew seismic data, by landmark master menu (OpenWorks 2003 Command Menu) → Applications → PostStack/PAL → selection two-dimension earthquake work area → press Launch ... key → enter PostStack Family 2003{ two-dimension earthquake work area name } master menu → select input two-dimension earthquake data (segy form) → create export file name (being the Vertical form of 9 characters) → Processes → Attributes → Instantaneous → Quadrature Trace → operation Run at Output Data at Input Data, the orthogonal seismic trace of promptly finishing after the end routine skew seismic data (segy formatted file) calculates, and converts orthogonal seismic trace data (2v2 formatted file) to.
3, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, it is characterized in that: the mutual comparative interpretation of described two-dimension earthquake is to enter SeisWorks2003 by landmark master menu (OpenWorks 2003Command Menu) → Applications → Seisworks → 2D, select the two-dimension earthquake work area, at the seismic profile interpretation window, seismic section shows by variable density, vertically between the optional 15.49cm/second-30.98cm/second of ratio, between the optional 32.28trace/cm-10.76trace/cm of grid scale, open two seismic profile interpretation windows, screen about placing respectively, the display parameter unanimity of the two, seismic data volume selects to be offset geological data and orthogonal seismic trace data respectively, the amplitude of objective interval on conventional offset seismic profile and its orthogonal seismic trace section is observed, contrast, selecting one on the orthogonal seismic trace section exists the phase place (crest or trough) of amplitude anomaly that the center line of the anomalous amplitude on the orthogonal seismic trace section and the position (when unusually not existing) suitable with it are made an explanation as a seismic horizon, and on the two-dimension earthquake survey line, follow the trail of contrast, can demonstrate this unusual general plane scope at the SeisWorks/Map window.
4, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, it is characterized in that: the described extraction average amplitude parallel planes of going forward side by side becomes figure to enter StratAmp 2003 master menus that extract average amplitude by landmark master menu (OpenWorks2003 Command Menu) → Data → Management → Seismic Project Manager → Horizons, select two-dimension earthquake work area and orthogonal seismic trace data, select the seismic horizon of input and the average amplitude layer position of name output again, calculate the average amplitude value of 2 milliseconds to 30 milliseconds of input layer positions (promptly this layer is up and down respectively between 1 millisecond-15 milliseconds) window in certain in a period of time, average amplitude is exported from workstation as a layer position, average amplitude value this layer in two fox softwares or Earthvision mapping software is carried out gridding, the plane distributing scope of the amplitude that shows abnormality.
5, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, and it is characterized in that: the two-dimension earthquake data are the skew geological data of segy form.
6, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, it is characterized in that: the skew seismic data volume that the two-dimension earthquake data are explained for the conventional two-dimension earthquake of confession that loads the back generation.
7, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, and it is characterized in that: the amplitude anomaly on the orthogonal seismic trace section is a crest.
8, orthogonal seismic trace as claimed in claim 1 is sought the method for lithologic anomalous body, and it is characterized in that: the amplitude anomaly on the orthogonal seismic trace section is a trough.
9, seek the method for lithologic anomalous body as the described orthogonal seismic trace of claim 1-8, it is characterized in that: when extracting the average amplitude value, can select the time window of input layer position, its scope is between 2 milliseconds-30 milliseconds.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200410088559 CN1267749C (en) | 2004-11-09 | 2004-11-09 | Method for seeking lithological anomalous body using orthogonal seismic trace data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200410088559 CN1267749C (en) | 2004-11-09 | 2004-11-09 | Method for seeking lithological anomalous body using orthogonal seismic trace data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1603860A true CN1603860A (en) | 2005-04-06 |
| CN1267749C CN1267749C (en) | 2006-08-02 |
Family
ID=34667157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200410088559 Expired - Fee Related CN1267749C (en) | 2004-11-09 | 2004-11-09 | Method for seeking lithological anomalous body using orthogonal seismic trace data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1267749C (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201408B (en) * | 2006-12-12 | 2010-06-09 | 中国石油天然气集团公司 | Method for constructing complex 2D seismic data flat spot control image |
| CN106772586A (en) * | 2017-02-20 | 2017-05-31 | 长江大学 | A kind of disguised fracture detection method based on seismic signal singularity |
| CN111665562A (en) * | 2020-07-20 | 2020-09-15 | 西南石油大学 | High-precision seismic sequence dividing method |
-
2004
- 2004-11-09 CN CN 200410088559 patent/CN1267749C/en not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201408B (en) * | 2006-12-12 | 2010-06-09 | 中国石油天然气集团公司 | Method for constructing complex 2D seismic data flat spot control image |
| CN106772586A (en) * | 2017-02-20 | 2017-05-31 | 长江大学 | A kind of disguised fracture detection method based on seismic signal singularity |
| CN106772586B (en) * | 2017-02-20 | 2019-04-02 | 长江大学 | A kind of concealment fracture detection method based on seismic signal singularity |
| CN111665562A (en) * | 2020-07-20 | 2020-09-15 | 西南石油大学 | High-precision seismic sequence dividing method |
| CN111665562B (en) * | 2020-07-20 | 2022-03-01 | 西南石油大学 | High-precision seismic sequence dividing method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1267749C (en) | 2006-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107817535B (en) | The determination method and apparatus of short lap | |
| Bouvier et al. | Three-dimensional seismic interpretation and fault sealing investigations, Nun River Field, Nigeria | |
| Lashin et al. | Characterization of the Qishn sandstone reservoir, Masila Basin–Yemen, using an integrated petrophysical and seismic structural approach | |
| CA2114424C (en) | Method and apparatus for petroleum and gas exploration | |
| CN104678438B (en) | A kind of CO2Time lapse seismic data CO in geological storage2The method of forecast of distribution | |
| Wang et al. | Structural characterization and hydrocarbon prediction for the SB5M strike-slip fault zone in the Shuntuo Low Uplift, Tarim Basin | |
| Halbouty | Geologic significance of Landsat data for 15 giant oil and gas fields | |
| CN104049283A (en) | Method for identification and prediction of fore-land basin extremely-thick conglomerate body | |
| Li | Geophysical exploration technology: Applications in lithological and stratigraphic reservoirs | |
| CN107861152A (en) | A kind of Mesozoic Strata fine correlation method | |
| White et al. | Sedimentology, statistics, and flow behavior for a tide-influenced deltaic sandstone, Frontier Formation, Wyoming, United States | |
| Karstens et al. | Deep-seated focused fluid migration as indicator for hydrocarbon leads in the East Shetland Platform, North Sea Province | |
| Zhao et al. | Relationship between fractures, stress, strike-slip fault and reservoir productivity, China Shunbei oil field, Tarim Basin | |
| Imam et al. | Mapping the geological structures in the Ras El Ush field (Gulf of Suez, Egypt), based on seismic interpretation and 3D modeling techniques | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| Al-Khazraji et al. | Reducing the risk and residual errors by uncertainty analysis of depth conversion in Nasiriya Oil Field, Iraq | |
| CN1267749C (en) | Method for seeking lithological anomalous body using orthogonal seismic trace data | |
| Ogbe et al. | Reservoir hydrocarbon prospectivity and productivity evaluations of sands S-600 and S-700 of Fega field, onshore Niger Delta Basin, Nigeria | |
| CN105093266A (en) | Seismogeological comprehensive description method for stratigraphic overlap line | |
| Jiang et al. | Recognizing the internal structure of normal faults in clastic rocks and its impact on hydrocarbon migration: A case study from Nanpu Depression in the Bohai Bay Basin, China | |
| Lyons et al. | Seismic Exploration for Stratigraphic Traps: Geophysical Exploration Methods | |
| Maleki et al. | Application of seismic attributes in structural study and fracture analysis of DQ oil field, Iran | |
| RU2145101C1 (en) | Method for estimation of service properties of gas-oil pool | |
| CN113093275B (en) | Method and device for improving drilling success rate of curved-flow river and curved-flow river delta oilfield | |
| Khatri | Controls on fracture network characteristics of the middle member of the Bakken Formation, Elm Coulee Field, Williston Basin, USA |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20060802 |