CN104698499A - Oil gas surveying method and device based on compression shielding layer seismic waves - Google Patents

Oil gas surveying method and device based on compression shielding layer seismic waves Download PDF

Info

Publication number
CN104698499A
CN104698499A CN201510153066.9A CN201510153066A CN104698499A CN 104698499 A CN104698499 A CN 104698499A CN 201510153066 A CN201510153066 A CN 201510153066A CN 104698499 A CN104698499 A CN 104698499A
Authority
CN
China
Prior art keywords
screen layer
layer
seismic
seismic event
wave energy
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
Application number
CN201510153066.9A
Other languages
Chinese (zh)
Other versions
CN104698499B (en
Inventor
田仁飞
王军
刘磊
李晶
班丽
孙渡
王玮
李忠新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sinopec Shengli Geological Scientific Reserch Institute
Chengdu Univeristy of Technology
Original Assignee
Sinopec Shengli Geological Scientific Reserch Institute
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinopec Shengli Geological Scientific Reserch Institute, Chengdu Univeristy of Technology filed Critical Sinopec Shengli Geological Scientific Reserch Institute
Priority to CN201510153066.9A priority Critical patent/CN104698499B/en
Publication of CN104698499A publication Critical patent/CN104698499A/en
Application granted granted Critical
Publication of CN104698499B publication Critical patent/CN104698499B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides an oil gas surveying method and device based on compression shielding layer seismic waves. The method comprises the steps that according to well seism data in a surveying zone, a shielding layer and a target layer are determined, the seism data comprise logging data and seismic waves; according to seismic wave energy values of the target layer and the shielding layer, the compressibility coefficient of the seismic waves of the shielding layer is determined; according to the compressibility coefficient, the seismic waves of the shielding layer are subjected to compression processing; and seismic waves of the target layer and the seismic waves of the shielding layer obtained after compressing processing are used for oil gas surveying. According to the method and device, the seismic waves of the shielding layer are compressed, meanwhile, the seismic waves of the target layer are kept unchanged, accordingly, the errors of oil gas surveying according to the seismic waves of the target layer can be reduced, and the success rate of oil gas surveying is improved.

Description

A kind of oil gas investigation method based on puncture mask layer seismic event and device
Technical field
The present invention relates to oil gas geophysical reconnaissance development field, in particular to a kind of oil gas investigation method based on puncture mask layer seismic event and device.
Background technology
At present, when carrying out seismic survey to oil gas, the seismic event that Water demand destination layer reflects Artificial Seismic Wave, to determine that whether this destination layer is for hydrocarbon zone.But because the earthquake wave intensity of destination layer reflection is very little, when there is the very strong screen layer of the seismic event of reflection in the stratum more than destination layer, the seismic event of destination layer can the seismic event of conductively-closed layer flood.Therefore how to weaken the earthquake wave intensity of screen layer, become urgent problem.
Current, provide a kind of method weakening the earthquake wave intensity of screen layer in correlation technique, comprising: by matching pursuit algorithm, wavelet decomposition is carried out to the seismic event that individual stratum is reflected, obtain multiple wavelet.Preset number maximum for amplitude in the wavelet an obtained wavelet is reconstructed, the waveform reconstructing and obtain is deducted by Artificial Seismic Wave, obtain the waveform that an amplitude is very little, using the seismic event of this waveform as screen layer, so weaken the earthquake wave intensity of screen layer, reduce its impact on the seismic event of destination layer.
But said method also needs the seismic event reconstructing destination layer in the process of carrying out wavelet decomposition and reconstruct, reconstruct can change the feature of the seismic event of destination layer, make the seismic event distortion of destination layer, the error causing oil gas to survey is very large, reduces the success ratio of oil gas exploration.
Summary of the invention
In view of this, the object of the embodiment of the present invention is to provide a kind of oil gas investigation method based on puncture mask layer seismic event and device, realize the seismic event of puncture mask layer, weaken the earthquake wave intensity of screen layer, keep the seismic event of destination layer constant simultaneously, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
First aspect, embodiments provide a kind of oil gas investigation method based on puncture mask layer seismic event, described method comprises:
According to the well shake data determination screen layer in survey area and destination layer, described well shake data comprise log data and seismic event;
The compressibility coefficient of the seismic event of described screen layer is determined according to the seismic wave energy value of described destination layer and described screen layer;
According to described compressibility coefficient, compression process is carried out to the seismic event of described screen layer;
The seismic event of the described screen layer after the seismic event of described destination layer and compression process is utilized to carry out oil gas exploration.
In conjunction with first aspect, embodiments provide the first possible implementation of first aspect, wherein, the described seismic wave energy value according to described destination layer and described screen layer determines the compressibility coefficient of the seismic event of described screen layer, comprising:
According to the seismic event of described screen layer, determine the seismic wave energy value of described screen layer;
According to the seismic event of described destination layer, determine the seismic wave energy value of described destination layer;
Calculate the ratio between the seismic wave energy value of described destination layer and the seismic wave energy value of described screen layer, described ratio is defined as the compressibility coefficient of the seismic event of described screen layer.
In conjunction with the first possible implementation of first aspect, embodiments provide the implementation that the second of first aspect is possible, wherein, the described seismic event according to described screen layer, determine the seismic wave energy value of described screen layer, comprising:
The predeterminated position of each waveform comprised by the seismic event of described screen layer by overall automatic tracing algorithm is coupled together, and obtains the lineups of described screen layer, and described predeterminated position comprises trough, crest or zero phase;
According to described lineups, from the seismic event of described screen layer, obtain the numerical value of the first preset number sampled point;
According to the numerical value of a described first preset number sampled point, calculate the seismic wave energy value of described screen layer.
In conjunction with first aspect, embodiments provide the third possible implementation of first aspect, wherein, described according to described compressibility coefficient to the seismic event of described screen layer carry out compression process, comprising:
The numerical value of the second preset number sampled point is obtained from the seismic event of described screen layer;
The numerical value of a described second preset number sampled point is reduced respectively described compressibility coefficient doubly.
In conjunction with first aspect, embodiments provide the 4th kind of possible implementation of first aspect, wherein, described according to described compressibility coefficient to the seismic event of described screen layer carry out compression process after, also comprise:
Obtain the seismic wave energy value of the rear described screen layer of compression, if the seismic wave energy value of described screen layer is more than or equal to predetermined threshold value after compression, then continue to perform the operation that the described seismic wave energy value according to described destination layer and described screen layer determines the compressibility coefficient of the seismic event of described screen layer.
Second aspect, embodiments provide a kind of oil gas surveying device based on puncture mask layer seismic event, described device comprises:
First determination module, for shaking data determination screen layer and destination layer according to the well in survey area, described well shake data comprise log data and seismic event;
Second determination module, for determining the compressibility coefficient of the seismic event of described screen layer according to the seismic wave energy value of described destination layer and described screen layer;
Compression module, for carrying out compression process according to described compressibility coefficient to the seismic event of described screen layer;
Exploration module, the seismic event of the described screen layer after processing for utilizing the seismic event of described destination layer and compression carries out oil gas exploration.
In conjunction with second aspect, embodiments provide the first possible implementation of second aspect, wherein, described second determination module comprises:
First determining unit, for the seismic event according to described screen layer, determines the seismic wave energy value of described screen layer;
Second determining unit, for the seismic event according to described destination layer, determines the seismic wave energy value of described destination layer;
Computing unit, for calculating the ratio between the seismic wave energy value of described destination layer and the seismic wave energy value of described screen layer, is defined as the compressibility coefficient of the seismic event of described screen layer by described ratio.
In conjunction with the first possible implementation of second aspect, embodiments provide the implementation that the second of second aspect is possible, wherein, described first determining unit comprises:
Connexon unit, the predeterminated position for each waveform comprised by the seismic event of described screen layer by overall automatic tracing algorithm is coupled together, and obtains the lineups of described screen layer, and described predeterminated position comprises trough, crest or zero phase;
Obtain subelement, for according to described lineups, from the seismic event of described screen layer, obtain the numerical value of the first preset number sampled point;
Computation subunit, for the numerical value according to a described first preset number sampled point, calculates the seismic wave energy value of described screen layer.
In conjunction with second aspect, embodiments provide the third possible implementation of second aspect, wherein, described compression module comprises:
Acquiring unit, for obtaining the numerical value of the second preset number sampled point in the seismic event from described screen layer;
Reducing unit, for reducing described compressibility coefficient respectively doubly by the numerical value of a described second preset number sampled point.
In conjunction with second aspect, embodiments provide the 4th kind of possible implementation of second aspect, wherein, described device also comprises:
Judge module, for obtaining the seismic wave energy value of the rear described screen layer of compression, if the seismic wave energy value of described screen layer is more than or equal to predetermined threshold value after compression, then triggers described second determination module and continue to perform the operation that the described seismic wave energy value according to described destination layer and described screen layer determines the compressibility coefficient of the seismic event of described screen layer.
In the method provided in the embodiment of the present invention and device, determine the compressibility coefficient of the seismic event of this screen layer according to the seismic wave energy value of survey area internal shield and destination layer; After utilizing this compressibility coefficient to carry out compression process to the seismic event of screen layer, carry out oil gas exploration again, owing to have compressed the seismic event of screen layer, weaken the earthquake wave intensity of screen layer, and do not change the seismic event of destination layer, therefore carry out oil gas exploration according to the seismic event of the seismic event of this destination layer and the screen layer after compressing, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
For making above-mentioned purpose of the present invention, feature and advantage become apparent, preferred embodiment cited below particularly, and coordinate appended accompanying drawing, be described in detail below.
Accompanying drawing explanation
In order to be illustrated more clearly in the technical scheme of the embodiment of the present invention, be briefly described to the accompanying drawing used required in embodiment below, be to be understood that, the following drawings illustrate only some embodiment of the present invention, therefore the restriction to scope should be counted as, for those of ordinary skill in the art, under the prerequisite not paying creative work, other relevant accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 shows a kind of oil gas investigation method process flow diagram based on puncture mask layer seismic event that the embodiment of the present invention 1 provides;
Fig. 2 A shows a kind of oil gas investigation method process flow diagram based on puncture mask layer seismic event that the embodiment of the present invention 2 provides;
Fig. 2 B shows a kind of destination layer well shake Fine calibration schematic diagram that the embodiment of the present invention 2 provides;
Fig. 2 C shows the seismic cross-section crossing Z1 well that the embodiment of the present invention 2 provides;
Fig. 2 D shows the lineups schematic diagram of the seismic event of the screen layer that the embodiment of the present invention 2 provides;
Fig. 2 E shows the seismic cross-section crossing Z1 well after the compression that the embodiment of the present invention 2 provides;
Fig. 3 A shows a kind of oil gas surveying device structural representation based on puncture mask layer seismic event that the embodiment of the present invention 3 provides;
Fig. 3 B shows another kind that the embodiment of the present invention 3 the provides oil gas surveying device structural representation based on puncture mask layer seismic event.
Embodiment
Below in conjunction with accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.The assembly of the embodiment of the present invention describing and illustrate in usual accompanying drawing herein can be arranged with various different configuration and design.Therefore, below to the detailed description of the embodiments of the invention provided in the accompanying drawings and the claimed scope of the present invention of not intended to be limiting, but selected embodiment of the present invention is only represented.Based on embodiments of the invention, the every other embodiment that those skilled in the art obtain under the prerequisite not making creative work, all belongs to the scope of protection of the invention.
When considering in correlation technique the earthquake wave intensity weakening screen layer, change the feature of the seismic event of destination layer, make the seismic event distortion of destination layer, and then cause the error of carrying out oil-gas exploration according to the seismic event of this destination layer very large, reduce the success ratio of oil-gas exploration.Based on this, embodiments provide a kind of method that the seismic event of screen layer is compressed and device, keep the seismic event of destination layer constant while compressing the seismic event of screen layer, thus reduce the error of oil-gas exploration, improve the success ratio of oil-gas exploration.Be described below by embodiment.
Embodiment 1
See Fig. 1, embodiments provide a kind of oil gas investigation method based on puncture mask layer seismic event, the method can be performed by the oil gas surveying device based on puncture mask layer seismic event.The method specifically comprises the following steps:
Step 101: according to the well shake data determination screen layer in survey area and destination layer, this well shake data comprise log data and seismic event, certainly, according to actual conditions, this well shake data can also comprise: geologic data and well data;
Step 102: according to the compressibility coefficient of the seismic event of the seismic wave energy value determination screen layer of destination layer and screen layer;
Step 103: compression process is carried out to the seismic event of screen layer according to compressibility coefficient;
Step 104: utilize the seismic event of the screen layer after the seismic event of destination layer and compression process to carry out oil gas exploration.
In the said method that the embodiment of the present invention provides, determine the compressibility coefficient of the seismic event of this screen layer according to the seismic wave energy value of survey area internal shield and destination layer; After utilizing this compressibility coefficient to carry out compression process to the seismic event of screen layer, carry out oil gas exploration again, owing to have compressed the seismic event of screen layer, weaken the earthquake wave intensity of screen layer, and do not change the seismic event of destination layer, therefore carry out oil gas exploration according to the seismic event of the seismic event of this destination layer and the screen layer after compressing, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
Wherein, above-mentioned seismic wave energy value can obtain according to the seismic event of corresponding geological stratification, and based on this, the compressibility coefficient of the seismic event of the above-mentioned seismic wave energy value determination screen layer according to destination layer and screen layer, comprising:
According to the seismic event of screen layer, determine the seismic wave energy value of screen layer;
The predeterminated position of each waveform first comprised by the seismic event of this screen layer by overall automatic tracing algorithm is coupled together, obtain the lineups of this screen layer, this predeterminated position can be trough, crest or zero phase, then according to these lineups, the numerical value of the first preset number sampled point is obtained from the seismic event of this screen layer, according to the numerical value of the first preset number sampled point, by following formula (1), calculate the seismic wave energy value of this screen layer;
B = 1 2 Σ j = 1 n y j 2 . . . ( 1 )
Wherein, B is the seismic wave energy value of screen layer, and n is the first preset number, and j is the sequence number of sampled point, y jfor the numerical value of a jth sampled point;
Wherein, overall automatic tracing algorithm can determine the lineups on all stratum of whole survey area rapidly, and is automatically coupled together according to crest, trough, positive zero phase or negative zero phase place by lineups.Wherein, when determining lineups according to crest, be the peak-peak of amplitude in each waveform is coupled together obtain lineups.When determining lineups according to trough, be the minimum peak of amplitude in each waveform coupled together to obtain lineups, the degree of accuracy of the lineups so obtained is very high.And the small patches of segmentation is dynamically connected certainly, can ensure that the lineups of the horizontal uncontinuities such as tomography are explained more rationally, and shorten the time determining lineups, improve the efficiency of oil gas exploration.After determining the lineups on each stratum, geologists only needs, according to research needs, to give rational geological Significance by these stratum.
According to the seismic event of destination layer, determine the seismic wave energy value of destination layer;
Particularly, in the seismic event of this destination layer, determine a waveform at random, this waveform is got the numerical value of a 3rd preset number sampled point equably.According to the numerical value of this 3rd preset number sampled point, by following formula (2), calculate the seismic wave energy value of this destination layer;
A = 1 2 Σ i = 1 m x i 2 . . . ( 2 )
Wherein, A is the seismic wave energy value of destination layer, and m is the 3rd preset number, and i is the sequence number of sampled point, x ibe the numerical value of i-th sampled point;
Calculate the ratio between the seismic wave energy value of destination layer and the seismic wave energy value of screen layer, ratio is defined as the compressibility coefficient of the seismic event of screen layer;
According to the seismic wave energy value of this destination layer and the seismic wave energy value of this screen layer, by following formula (3), calculate the compressibility coefficient of the seismic event of this screen layer.
C = A / B = ( Σ i = 1 m x i 2 ) / ( Σ j = 1 n y j 2 ) . . . ( 3 )
Wherein, C is the compressibility coefficient of the seismic event of screen layer.
Wherein, this compressibility coefficient is greater than 0 and is less than 1, utilizes this compressibility coefficient can the earthquake wave intensity of effective puncture mask layer, and the determination mode advantages of simple of this compressibility coefficient.
Wherein, according to compressibility coefficient, compression process is carried out to the seismic event of screen layer, comprising:
The numerical value of the second preset number sampled point is obtained from the seismic event of screen layer;
The numerical value of the second preset number sampled point is reduced respectively compressibility coefficient doubly.
Particularly, each waveform that the seismic event for this screen layer comprises, according to the lineups of this screen layer, is arranged in from this waveform the numerical value that the part of more than these lineups and part once obtain the second preset number sampled point equably.The numerical value of the second preset number sampled point is multiplied by this compressibility coefficient respectively, this waveform is reduced this compressibility coefficient doubly.Other each waveforms with this waveform, can be reduced this compressibility coefficient doubly, thus realize by this compressibility coefficient of the earthquake wave condensing of this screen layer doubly by other each waveforms that the seismic event for this screen layer comprises in the manner described above respectively.
Wherein, after compression process is carried out to the seismic event of this screen layer, weaken the earthquake wave intensity of this screen layer, thus eliminate this screen layer impact is flooded on destination layer.Do not change the earthquake wave amplitude on other stratum comprising destination layer, phase place and frequecy characteristic simultaneously, more can give prominence to the feature of the seismic event of destination layer, for the calculating of destination layer reservoir parameter provides reliable basis.
Further, because the earthquake wave intensity of the screen layer after compression may be still very strong, may be still very large on the impact of destination layer, so carry out after compression processes, also comprising to the seismic event of screen layer according to compressibility coefficient:
Obtain the seismic wave energy value of screen layer after compression, if the seismic wave energy value of screen layer is more than or equal to predetermined threshold value after compression, then continue to perform the operation of the compressibility coefficient of the seismic event of the seismic wave energy value determination screen layer according to destination layer and screen layer.
Wherein, this predetermined threshold value is that technician is obtained by great many of experiments, if the seismic wave energy value of this screen layer is less than this predetermined threshold value after compression, then show to produce a desired effect to the compression process of screen layer, the screen layer after compression is very little on the impact of destination layer.If the seismic wave energy value of this screen layer is still more than or equal to this predetermined threshold value after compression, then show that the screen layer after compressing is still very large on the impact of this destination layer, therefore need to utilize the rear seismic wave energy value of screen layer of compression and the seismic wave energy value of destination layer, again a compressibility coefficient is obtained, the seismic event of new compressibility coefficient to the screen layer after compression is utilized to compress further, until the seismic wave energy value of screen layer is less than this predetermined threshold value.
When the seismic wave energy value of the screen layer after compressing is less than this predetermined threshold value, the impact of this screen layer on destination layer diminishes, the seismic event of the screen layer after the seismic event of destination layer and compression can be utilized to process, in conjunction with the information such as geologic data and well data, predict that this destination layer is the probability of hydrocarbon zone, and predict the reservoir parameter such as position and storage capacity of this destination layer.
Wherein, in embodiments of the present invention, the strong amplitude impact of screen layer more simply can be eliminated efficiently by compressibility coefficient, and do not change the feature of the seismic event of other formation, more be conducive to the accurate calculating of reservoir parameter, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
Embodiment 2
See Fig. 2 A, embodiments provide a kind of oil gas investigation method based on puncture mask layer seismic event, the method can be performed by the oil gas surveying device based on puncture mask layer seismic event, and the method specifically comprises the following steps:
Step 201: according to the well shake data determination screen layer in survey area and destination layer, this well shake data comprise log data and seismic event;
Wherein, destination layer may be hydrocarbon zone, and destination layer is the formation at target locations that seismic prospecting will be studied.Screen layer is generally and causes the discontinuous unconformity surface of upper and lower Stratigraphic Time due to the stratum lacking a certain period between the stratum successively deposited, as above the Cretaceous System Tugulu group group stratum covered directly overlays on Jurassic three group stratum, work river, unify in intercalary delection Jurassic systerm upper system stratum, and the unconformity plane so formed is screen layer.The seismic event of screen layer generally has the features such as strong amplitude and low frequency.
Wherein, well shake data comprise seismic event and the log data on each stratum in survey area, the seismic event on stratum comprises multiple waveform, and log data comprises the data such as AC (Acoustic loggingCurve, acoustic logging) curve, density, current potential.In addition, well shake data can also comprise geologic data in survey area and well data etc., can also analyze by geologic data and well data to screen layer and destination layer.
In the seismic survey of oil gas, on earth's surface with manual method earthquake-wave-exciting, this seismic event is when propagating to depths, the earth's core, and each stratum can be reflected this seismic event.Wave detector is set on earth's surface, is detected the seismic event of each stratum reflection by this wave detector.Wave detector by the seismic input wave on each stratum that detects to terminal.The feature of terminal to the seismic event on each stratum is analyzed, utilize well to shake the method for demarcating and determine destination layer, then detect the seismic event whether existing near this destination layer and there is the features such as strong amplitude and low frequency, if exist, then the seismic event detected is defined as the seismic event of screen layer.
Wherein, the technology that well shake is demarcated is very ripe, does not repeat them here.
Wherein, the concrete operations detecting the seismic event of screen layer near destination layer can be as follows:
Default oscillator intensity and predeterminated frequency is stored in terminal.For with each waveform of this destination layer in the predeterminable range, terminal obtains oscillator intensity and the frequency of this waveform, the oscillator intensity of this waveform and default oscillator intensity is compared, and the frequency of this waveform and predeterminated frequency is compared.If the oscillator intensity of this waveform is greater than default oscillator intensity, and the frequency of this waveform is less than predeterminated frequency, then this waveform is defined as the seismic event of screen layer.For with this destination layer other each waveforms in the predeterminable range, can with this waveform, determine that whether other each waveforms are the seismic event of screen layer according to aforesaid operations.
Wherein, after determining screen layer and destination layer according to the operation of this step, the operation of 202-204 as follows can obtain compressibility coefficient, this compressibility coefficient is used for carrying out compression to the seismic event of screen layer and processes.
Step 202: according to the seismic event of screen layer, determines the seismic wave energy value of this screen layer;
Wherein, this step can realize in the operation of S1-S3 as follows, is specially:
S1: the predeterminated position of each waveform comprised by the seismic event of this screen layer by overall automatic tracing algorithm is coupled together, and obtains the lineups of this screen layer, this predeterminated position can be trough, crest or zero phase;
Wherein, lineups are the lines of the extreme value that in multiple waveforms of comprising of the seismic event on stratum, phase place is identical, as the line of the trough of each waveform, or the line of crest, or the line of zero phase.
Wherein, overall automatic tracing algorithm is the lineups being explained the seismic event on each stratum by the disposable batch of overall automatic tracing seismic stratigraphic interpretation.The well on each stratum in this survey area shake data are converted to geologic model, this geologic model are divided into the small patches presetting size one by one.For each stratum in this survey area, as this screen layer, by the mode that cost function minimization controls, according to the method extract layer position of global optimum from this screen layer that geologic model comprises, carry out stratigraphy and analysis of sequence structure again, find out the small patches at the predeterminated position place that phase place is identical in each waveform in this screen layer, each small patches found out is coupled together, just obtains the lineups of this screen layer.For other each stratum, with this screen layer, the lineups on other each stratum can be determined respectively according to aforesaid operations.
Wherein, overall automatic tracing algorithm can determine the lineups on all stratum of whole survey area rapidly, and is automatically coupled together according to crest, trough, positive zero phase or negative zero phase place by lineups.Wherein, when determining lineups according to crest, be the peak-peak of amplitude in each waveform is coupled together obtain lineups.When determining lineups according to trough, be the minimum peak of amplitude in each waveform coupled together to obtain lineups, the degree of accuracy of the lineups so obtained is very high.And the small patches of segmentation is dynamically connected certainly, can ensure that the lineups of the horizontal uncontinuities such as tomography are explained more rationally, and shorten the time determining lineups, improve the efficiency of oil gas exploration.After determining the lineups on each stratum, geologists only needs, according to research needs, to give rational geological Significance by these stratum.
S2: according to these lineups, obtains the numerical value of the first preset number sampled point from the seismic event of this screen layer;
Particularly, from the seismic event of this screen layer, a waveform is determined at random.According to these lineups, determine the portion waveshape being positioned at more than these lineups in this waveform, and determine the portion waveshape being positioned at below these lineups.The uniform numerical value obtaining the first preset number sampled point in this two parts waveform.
Wherein, obtain the numerical value of sampled point according to these lineups, the numerical value obtaining sampled point in the seismic event of screen layer equably can be guaranteed, thus make the seismic wave energy value that obtains according to the numerical value of sampled point more accurate.
S3: according to the numerical value of the first preset number sampled point, calculates the seismic wave energy value of this screen layer.
Particularly, according to the numerical value of the first preset number sampled point, by following formula (1), calculate the seismic wave energy value of this screen layer.
B = 1 2 Σ j = 1 n y j 2 . . . ( 1 )
Wherein, B is the seismic wave energy value of screen layer, and n is the first preset number, and j is the sequence number of sampled point, y jfor the numerical value of a jth sampled point.
Step 203: according to the seismic event of this destination layer, determines the seismic wave energy value of this destination layer;
Particularly, in the seismic event of this destination layer, determine a waveform at random, this waveform is got the numerical value of a 3rd preset number sampled point equably.According to the numerical value of this 3rd preset number sampled point, by following formula (2), calculate the seismic wave energy value of this destination layer.
A = 1 2 Σ i = 1 m x i 2 . . . ( 2 )
Wherein, A is the seismic wave energy value of destination layer, and m is the 3rd preset number, and i is the sequence number of sampled point, x ibe the numerical value of i-th sampled point.
Further, also first can determine the lineups of this destination layer, then according to the lineups of destination layer, from the seismic event of this destination layer, obtain a 3rd preset number sampled point equably.
Step 204: calculate the ratio between the seismic wave energy value of this destination layer and the seismic wave energy value of this screen layer, is defined as the compressibility coefficient of the seismic event of this screen layer by this ratio;
Wherein, according to the seismic wave energy value of this destination layer and the seismic wave energy value of this screen layer, by following formula (3), the compressibility coefficient of the seismic event of this screen layer is calculated.
C = A / B = ( Σ i = 1 m x i 2 ) / ( Σ j = 1 n y j 2 ) . . . ( 3 )
Wherein, C is the compressibility coefficient of the seismic event of screen layer.
Wherein, this compressibility coefficient is greater than 0 and is less than 1, utilizes this compressibility coefficient can the earthquake wave intensity of effective puncture mask layer.
Wherein, after obtaining compressibility coefficient by the operation of above-mentioned steps 202-204, the operation of 205 carrys out the seismic event of puncture mask layer as follows, to weaken the earthquake wave intensity of screen layer, avoids the seismic event of screen layer on the impact of the seismic event of destination layer.
Step 205: compression process is carried out to the seismic event of this screen layer according to this compressibility coefficient;
Particularly, each waveform that the seismic event for this screen layer comprises, according to the lineups of this screen layer, is arranged in from this waveform the numerical value that the part of more than these lineups and part once obtain the second preset number sampled point equably.The numerical value of the second preset number sampled point is multiplied by this compressibility coefficient respectively, this waveform is reduced this compressibility coefficient doubly.Other each waveforms with this waveform, can be reduced this compressibility coefficient doubly, thus realize by this compressibility coefficient of the earthquake wave condensing of this screen layer doubly by other each waveforms that the seismic event for this screen layer comprises in the manner described above respectively.
Wherein, after compression process is carried out to the seismic event of this screen layer, weaken the earthquake wave intensity of this screen layer, thus eliminate this screen layer impact is flooded on destination layer.Do not change the earthquake wave amplitude on other stratum comprising destination layer, phase place and frequecy characteristic simultaneously, more can give prominence to the feature of the seismic event of destination layer, for the calculating of destination layer reservoir parameter provides reliable basis.
Wherein, after by the operation of above-mentioned steps 205 compression process being carried out to the seismic event of this screen layer, also need the operation of as follows 206 to judge whether to produce a desired effect to the compression of screen layer, to determine whether to need to compress process further to the seismic event of this screen layer.
Step 206: the seismic wave energy value obtaining this screen layer after compressing, judges that after compressing, whether the seismic wave energy value of this screen layer is more than or equal to predetermined threshold value, if so, then returns and performs step 204, if not, then performs step 207;
Wherein, the operation obtaining the seismic wave energy value of this screen layer after compressing is identical with the operation of above-mentioned steps 202, does not repeat them here.
Wherein, this predetermined threshold value is that technician is obtained by great many of experiments, if the seismic wave energy value of this screen layer is less than this predetermined threshold value after compression, then show to produce a desired effect to the compression process of screen layer, the screen layer after compression is very little on the impact of destination layer.If the seismic wave energy value of this screen layer is still more than or equal to this predetermined threshold value after compression, then show that the screen layer after compressing is still very large on the impact of this destination layer, therefore need to return step 204, utilize the rear seismic wave energy value of screen layer of compression and the seismic wave energy value of destination layer, again a compressibility coefficient is obtained, the seismic event of new compressibility coefficient to the screen layer after compression is utilized to compress further, until the seismic wave energy value of screen layer is less than this predetermined threshold value.
Further, the seismic event of the screen layer after compression can also be put back into the position that in the seismic event on each stratum, this screen layer is corresponding, whether the ratio between the earthquake wave frequency of the screen layer then after the seismic event of evaluating objects layer and compression, the ratio between amplitude or the ratio between phase place reach default value, if reach default value, then stop the seismic event of this screen layer of compression, if do not reach default value, then continue to compress further according to the seismic event of method provided by the invention to the screen layer after compression.
Further, in embodiments of the present invention, when the seismic wave energy value of the screen layer after compressing is more than or equal to predetermined threshold value, the first preset number in all right set-up procedure 202 and the 3rd preset number in step 203, recalculate the seismic wave energy value of original screen layer and the seismic wave energy value of destination layer, obtain a new compressibility coefficient according to the seismic wave energy value of screen layer and destination layer, and again the seismic event of original screen layer is compressed according to this new compressibility coefficient.
Step 207: utilize the seismic event of the screen layer after the seismic event of destination layer and compression process to carry out oil gas exploration.
Wherein, utilize the seismic event of screen layer after the seismic event of destination layer and compression process, in conjunction with the information such as geologic data and well data, predict that this destination layer is the probability of hydrocarbon zone, and predict the reservoir parameters such as the position of this destination layer and storage capacity.
Wherein, in order to understand method provided by the invention more easily, be described in detail below in conjunction with figure.Destination layer well shake Fine calibration schematic diagram as shown in Figure 2 B, this schematic diagram comprises AC curve, stratum, theogram and crosses Z1 well seismic section.Wherein, the destination layer between J1s1 and the J1s21 of stratum is sound wave and densimetric curve in the well-log information utilizing somewhere Z1 well, and in conjunction with 3-D data volume, the well shake carrying out fine synthesis record is demarcated.As can be seen from Fig. 2 B, a stratum K1tg is had on the J1s1 of stratum, K1tg layer amplitude is very strong, in conjunction with geologic data and well data, K1tg layer is demarcated as screen layer, K1tg layer has very strong Seismic reflection character, has obvious energy barrier effect, causes the destination layer reflected energy between J1s1 and J1s21 to die down.Cross the seismic cross-section of Z1 well as that shown in fig. 2 c, can clearly be seen that the reflection strength of K1tg layer is very large, its earthquake wave intensity is more much larger than the earthquake wave intensity of the destination layer between J1s1 and J1s21, causes the earthquake wave amplitude of the destination layer between J1s1 and J1s21 very faint.
Wherein, the lineups schematic diagram of the seismic event of screen layer as shown in Figure 2 D, the seismic event of screen layer K1tg is divided into upper and lower two parts by the lineups of this screen layer K1tg, when obtaining sampled point in the seismic event of screen layer K1tg, from these upper and lower two parts, sampled point can be set uniformly, and obtain the numerical value of sampled point.Numerical value according to sampled point obtains compressibility coefficient, and after compressing according to the seismic event of compressibility coefficient to screen layer K1tg, see the seismic cross-section crossing Z1 well after the compression shown in Fig. 2 E, after comparison diagram 2C and Fig. 2 E can clearly be seen that compression, the earthquake wave intensity of screen layer K1tg has significantly died down, and the amplitude characteristic compressing the destination layer between rear J1s1 and J1s21 becomes more obvious.
Wherein, in embodiments of the present invention, the strong amplitude impact of screen layer more simply can be eliminated efficiently by compressibility coefficient, and do not change the feature of the seismic event of other formation, more be conducive to the accurate calculating of reservoir parameter, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
In embodiments of the present invention, the compressibility coefficient of the seismic event of this screen layer is determined according to the seismic wave energy value of survey area internal shield and destination layer; After utilizing this compressibility coefficient to carry out compression process to the seismic event of screen layer, carry out oil gas exploration again, owing to have compressed the seismic event of screen layer, weaken the earthquake wave intensity of screen layer, and do not change the seismic event of destination layer, therefore carry out oil gas exploration according to the seismic event of the seismic event of this destination layer and the screen layer after compressing, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
Embodiment 3
See Fig. 3 A, embodiments provide a kind of oil gas surveying device based on puncture mask layer seismic event, this device is for performing the above-mentioned oil gas investigation method based on puncture mask layer seismic event.This device specifically comprises:
First determination module 301, for shaking data determination screen layer and destination layer according to the well in survey area, this well shake data comprise log data and seismic event;
Second determination module 302, for the compressibility coefficient of the seismic event of the seismic wave energy value determination screen layer according to destination layer and screen layer;
Compression module 303, for carrying out compression process according to compressibility coefficient to the seismic event of screen layer;
Exploration module 304, the seismic event of the screen layer after processing for utilizing the seismic event of destination layer and compression carries out oil gas exploration.
Wherein, see Fig. 3 B, the second determination module 302 comprises:
First determining unit 3021, for the seismic event according to screen layer, determines the seismic wave energy value of screen layer;
Second determining unit 3022, for the seismic event according to destination layer, determines the seismic wave energy value of destination layer;
Computing unit 3023, for calculating the ratio between the seismic wave energy value of destination layer and the seismic wave energy value of screen layer, is defined as the compressibility coefficient of the seismic event of screen layer by ratio.
Wherein, see Fig. 3 B, the first determining unit 3021 comprises:
Connexon unit 30211, the predeterminated position for each waveform comprised by the seismic event of screen layer by overall automatic tracing algorithm is coupled together, and obtains the lineups of screen layer, and predeterminated position comprises trough, crest or zero phase;
Obtain subelement 30212, for according to lineups, from the seismic event of screen layer, obtain the numerical value of the first preset number sampled point;
Computation subunit 30213, for the numerical value according to the first preset number sampled point, calculates the seismic wave energy value of screen layer.
Wherein, see Fig. 3 B, compression module 303 comprises:
Acquiring unit 3031, for obtaining the numerical value of the second preset number sampled point from the seismic event of screen layer;
Reducing unit 3032, for reducing compressibility coefficient respectively doubly by the numerical value of the second preset number sampled point.
Further, see Fig. 3 B, this device also comprises:
Judge module 305, for obtaining the seismic wave energy value of the rear screen layer of compression, if the seismic wave energy value of screen layer is more than or equal to predetermined threshold value after compression, then trigger the operation that the second determination module 302 continues to perform the compressibility coefficient of the seismic event of the seismic wave energy value determination screen layer according to destination layer and screen layer.
In embodiments of the present invention, the compressibility coefficient of the seismic event of this screen layer is determined according to the seismic wave energy value of survey area internal shield and destination layer; After utilizing this compressibility coefficient to carry out compression process to the seismic event of screen layer, carry out oil gas exploration again, owing to have compressed the seismic event of screen layer, weaken the earthquake wave intensity of screen layer, and do not change the seismic event of destination layer, therefore carry out oil gas exploration according to the seismic event of the seismic event of this destination layer and the screen layer after compressing, reduce the error of oil gas exploration, improve the success ratio of oil gas exploration.
The oil gas surveying device based on puncture mask layer seismic event that the embodiment of the present invention provides can be the specific hardware on equipment or the software be installed on equipment or firmware etc.Those skilled in the art can be well understood to, and for convenience and simplicity of description, the specific works process of the system of aforementioned description, device and unit, all can with reference to the corresponding process in said method embodiment.
In several embodiments that the application provides, should be understood that disclosed apparatus and method can realize by another way.Device embodiment described above is only schematic, such as, the division of described unit, be only a kind of logic function to divide, actual can have other dividing mode when realizing, again such as, multiple unit or assembly can in conjunction with or another system can be integrated into, or some features can be ignored, or do not perform.Another point, shown or discussed coupling each other or direct-coupling or communication connection can be by some communication interfaces, and the indirect coupling of device or unit or communication connection can be electrical, machinery or other form.
The described unit illustrated as separating component or can may not be and physically separates, and the parts as unit display can be or may not be physical location, namely can be positioned at a place, or also can be distributed in multiple network element.Some or all of unit wherein can be selected according to the actual needs to realize the object of the present embodiment scheme.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, also can be that the independent physics of unit exists, also can two or more unit in a unit integrated.
If described function using the form of SFU software functional unit realize and as independently production marketing or use time, can be stored in a computer read/write memory medium.Based on such understanding, the part of the part that technical scheme of the present invention contributes to prior art in essence in other words or this technical scheme can embody with the form of software product, this computer software product is stored in a storage medium, comprising some instructions in order to make a computer equipment (can be personal computer, server, or the network equipment etc.) perform all or part of step of method described in each embodiment of the present invention.And aforesaid storage medium comprises: USB flash disk, portable hard drive, ROM (read-only memory) (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. various can be program code stored medium.
The above; be only the specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; change can be expected easily or replace, all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should described be as the criterion with the protection domain of claim.

Claims (10)

1. based on an oil gas investigation method for puncture mask layer seismic event, it is characterized in that, described method comprises:
According to the well shake data determination screen layer in survey area and destination layer, described well shake data comprise log data and seismic event;
The compressibility coefficient of the seismic event of described screen layer is determined according to the seismic wave energy value of described destination layer and described screen layer;
According to described compressibility coefficient, compression process is carried out to the seismic event of described screen layer;
The seismic event of the described screen layer after the seismic event of described destination layer and compression process is utilized to carry out oil gas exploration.
2. method according to claim 1, is characterized in that, the described seismic wave energy value according to described destination layer and described screen layer determines the compressibility coefficient of the seismic event of described screen layer, comprising:
According to the seismic event of described screen layer, determine the seismic wave energy value of described screen layer;
According to the seismic event of described destination layer, determine the seismic wave energy value of described destination layer;
Calculate the ratio between the seismic wave energy value of described destination layer and the seismic wave energy value of described screen layer, described ratio is defined as the compressibility coefficient of the seismic event of described screen layer.
3. method according to claim 2, is characterized in that, the described seismic event according to described screen layer, determines the seismic wave energy value of described screen layer, comprising:
The predeterminated position of each waveform comprised by the seismic event of described screen layer by overall automatic tracing algorithm is coupled together, and obtains the lineups of described screen layer, and described predeterminated position comprises trough, crest or zero phase;
According to described lineups, from the seismic event of described screen layer, obtain the numerical value of the first preset number sampled point;
According to the numerical value of a described first preset number sampled point, calculate the seismic wave energy value of described screen layer.
4. method according to claim 1, is characterized in that, described according to described compressibility coefficient to the seismic event of described screen layer carry out compression process, comprising:
The numerical value of the second preset number sampled point is obtained from the seismic event of described screen layer;
The numerical value of a described second preset number sampled point is reduced respectively described compressibility coefficient doubly.
5. the method according to any one of claim 1-4 claim, is characterized in that, described according to described compressibility coefficient to the seismic event of described screen layer carry out compression process after, also comprise:
Obtain the seismic wave energy value of the rear described screen layer of compression, if the seismic wave energy value of described screen layer is more than or equal to predetermined threshold value after compression, then continue to perform the operation that the described seismic wave energy value according to described destination layer and described screen layer determines the compressibility coefficient of the seismic event of described screen layer.
6. based on an oil gas surveying device for puncture mask layer seismic event, it is characterized in that, described device comprises:
First determination module, for shaking data determination screen layer and destination layer according to the well in survey area, described well shake data comprise log data and seismic event;
Second determination module, for determining the compressibility coefficient of the seismic event of described screen layer according to the seismic wave energy value of described destination layer and described screen layer;
Compression module, for carrying out compression process according to described compressibility coefficient to the seismic event of described screen layer;
Exploration module, the seismic event of the described screen layer after processing for utilizing the seismic event of described destination layer and compression carries out oil gas exploration.
7. device according to claim 6, is characterized in that, described second determination module comprises:
First determining unit, for the seismic event according to described screen layer, determines the seismic wave energy value of described screen layer;
Second determining unit, for the seismic event according to described destination layer, determines the seismic wave energy value of described destination layer;
Computing unit, for calculating the ratio between the seismic wave energy value of described destination layer and the seismic wave energy value of described screen layer, is defined as the compressibility coefficient of the seismic event of described screen layer by described ratio.
8. device according to claim 7, is characterized in that, described first determining unit comprises:
Connexon unit, the predeterminated position for each waveform comprised by the seismic event of described screen layer by overall automatic tracing algorithm is coupled together, and obtains the lineups of described screen layer, and described predeterminated position comprises trough, crest or zero phase;
Obtain subelement, for according to described lineups, from the seismic event of described screen layer, obtain the numerical value of the first preset number sampled point;
Computation subunit, for the numerical value according to a described first preset number sampled point, calculates the seismic wave energy value of described screen layer.
9. device according to claim 6, is characterized in that, described compression module comprises:
Acquiring unit, for obtaining the numerical value of the second preset number sampled point in the seismic event from described screen layer;
Reducing unit, for reducing described compressibility coefficient respectively doubly by the numerical value of a described second preset number sampled point.
10. the device according to any one of claim 6-9 claim, is characterized in that, described device also comprises:
Judge module, for obtaining the seismic wave energy value of the rear described screen layer of compression, if the seismic wave energy value of described screen layer is more than or equal to predetermined threshold value after compression, then triggers described second determination module and continue to perform the operation that the described seismic wave energy value according to described destination layer and described screen layer determines the compressibility coefficient of the seismic event of described screen layer.
CN201510153066.9A 2015-04-01 2015-04-01 A kind of oil gas investigation method and device based on puncture mask layer seismic wave Expired - Fee Related CN104698499B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510153066.9A CN104698499B (en) 2015-04-01 2015-04-01 A kind of oil gas investigation method and device based on puncture mask layer seismic wave

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510153066.9A CN104698499B (en) 2015-04-01 2015-04-01 A kind of oil gas investigation method and device based on puncture mask layer seismic wave

Publications (2)

Publication Number Publication Date
CN104698499A true CN104698499A (en) 2015-06-10
CN104698499B CN104698499B (en) 2017-10-31

Family

ID=53345817

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510153066.9A Expired - Fee Related CN104698499B (en) 2015-04-01 2015-04-01 A kind of oil gas investigation method and device based on puncture mask layer seismic wave

Country Status (1)

Country Link
CN (1) CN104698499B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107643539A (en) * 2016-07-21 2018-01-30 中国石油化工股份有限公司 A kind of method that strong screen layer is peeled off based on the analysis of coal seam seismic response features
CN109917458A (en) * 2019-04-16 2019-06-21 成都理工大学 Oil gas analyzing method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040122595A1 (en) * 2002-12-23 2004-06-24 Henri-Pierre Valero Methods and apparatus for processing acoustic waveforms received in a borehole
GB2426822A (en) * 2003-12-24 2006-12-06 Baker Hughes Inc Downhole library of master wavelets for vsp-while-drilling applications
CN102879825A (en) * 2012-09-18 2013-01-16 吉林大学 High power pulse noise detection and suppression method for vibroseis seismic data
CN103592684A (en) * 2013-10-21 2014-02-19 中国石油天然气集团公司 Massive seismic data compression method and device for preserving spatial attribute information
CN104375188A (en) * 2014-11-19 2015-02-25 中国石油天然气集团公司 Seismic wave transmission attenuation compensation method and device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040122595A1 (en) * 2002-12-23 2004-06-24 Henri-Pierre Valero Methods and apparatus for processing acoustic waveforms received in a borehole
GB2426822A (en) * 2003-12-24 2006-12-06 Baker Hughes Inc Downhole library of master wavelets for vsp-while-drilling applications
CN102879825A (en) * 2012-09-18 2013-01-16 吉林大学 High power pulse noise detection and suppression method for vibroseis seismic data
CN103592684A (en) * 2013-10-21 2014-02-19 中国石油天然气集团公司 Massive seismic data compression method and device for preserving spatial attribute information
CN104375188A (en) * 2014-11-19 2015-02-25 中国石油天然气集团公司 Seismic wave transmission attenuation compensation method and device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘炳杨 等: "《中国地球物理2012》", 31 October 2012 *
张军华 等: "强屏蔽层下弱反射储层特征分析及识别方法", 《特种油气藏》 *
李可恩 等: "高速屏蔽层下地震数据的采集", 《天然气地球科学》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107643539A (en) * 2016-07-21 2018-01-30 中国石油化工股份有限公司 A kind of method that strong screen layer is peeled off based on the analysis of coal seam seismic response features
CN109917458A (en) * 2019-04-16 2019-06-21 成都理工大学 Oil gas analyzing method and device

Also Published As

Publication number Publication date
CN104698499B (en) 2017-10-31

Similar Documents

Publication Publication Date Title
CN113759425B (en) Method and system for evaluating filling characteristics of deep paleo-karst reservoir stratum by well-seismic combination
US7872943B2 (en) Method of accounting for wavelet stretch in seismic data
CN113759424B (en) Karst reservoir filling analysis method and system based on spectral decomposition and machine learning
WO2019062655A1 (en) Method and device for determining thin interlayer
AU741349B2 (en) Method for characterizing subsurface petrophysical properties using linear shape attributes
WO2012139082A1 (en) Event selection in the image domain
CN106249299A (en) The lower thin layer weak reflected seismic energy restoration methods of strong reflection shielding and device
CN105652323A (en) Reservoir stratum prediction method
US20230058742A1 (en) Method and storage medium for quantitative reconstruction of paleowater depth based on milankovitch cycles
CN107576985B (en) A kind of method and apparatus of seismic inversion
CN104375178A (en) Carbonate rock crack hole reservoir prediction method and device
CN108254785B (en) Reservoir determination method and device
CN201747364U (en) Interval transit time curve reconstruction equipment
CN104698499A (en) Oil gas surveying method and device based on compression shielding layer seismic waves
RU2011148308A (en) METHOD FOR COMPREHENSIVE PROCESSING OF GEOPHYSICAL DATA AND TECHNOLOGICAL SYSTEM "LITOSCAN" FOR ITS IMPLEMENTATION
Infantino et al. Study of the spatial correlation of earthquake ground motion by means of physics-based numerical scenarios
CN106610507A (en) Method and device for removing ghost reflection interference
CN113219531A (en) Method and device for identifying gas-water distribution of tight sandstone
CN112711068A (en) Method and device for predicting effective oil and gas reservoir in sandstone
CN107576984B (en) The determination method and apparatus of window when during seismic amplitude attributes extraction
CN111025393B (en) Reservoir prediction method, device, equipment and medium for stratum containing thin coal seam
CN104111476A (en) Method and device for building formation velocity field
CN112130209A (en) Karst reservoir prediction method and device
CN107085236B (en) The determination method and apparatus of maximum offset
CN105989240A (en) Sand layer prediction method and apparatus

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
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: 20171031

Termination date: 20210401