CN107656309B - A kind of method and device of determining fault throw - Google Patents
A kind of method and device of determining fault throw Download PDFInfo
- Publication number
- CN107656309B CN107656309B CN201710696640.4A CN201710696640A CN107656309B CN 107656309 B CN107656309 B CN 107656309B CN 201710696640 A CN201710696640 A CN 201710696640A CN 107656309 B CN107656309 B CN 107656309B
- Authority
- CN
- China
- Prior art keywords
- data
- coherence
- interval
- interest
- tomography
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000003325 tomography Methods 0.000 claims abstract description 60
- 238000012545 processing Methods 0.000 claims abstract description 45
- 238000005520 cutting process Methods 0.000 claims abstract description 44
- 230000001427 coherent effect Effects 0.000 claims abstract description 19
- 230000003595 spectral effect Effects 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 238000013507 mapping Methods 0.000 claims description 32
- 238000005070 sampling Methods 0.000 claims description 14
- 230000000149 penetrating effect Effects 0.000 claims 2
- 239000003086 colorant Substances 0.000 claims 1
- 238000013506 data mapping Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/65—Source localisation, e.g. faults, hypocenters or reservoirs
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application discloses a kind of method and device of determining fault throw.The method is provided with the seismic data of interval of interest, the described method includes: the frequency range based on the seismic data, spectral decomposition processing is carried out to the seismic data, multiple earthquakes frequency dividing data of the interval of interest are obtained, and obtains specified earthquake from the multiple earthquake frequency dividing data and divides data；Wherein, the fracture information that the specified earthquake frequency dividing data are used to reflect in the interval of interest；Coherent processing is carried out to the specified earthquake frequency dividing data, obtains coherence data；Wherein, the coherence data includes the coherence factor of multiple data points and the data point；Based on the coherence data, the cutting relationship between the tomography in the interval of interest and a karst beading is determined；According to the cutting relationship, the turn-off of the tomography is determined.The accuracy of determined fault throw can be improved in technical solution provided by the embodiments of the present application.
Description
Technical field
This application involves oil field prospecting seismic data interpretation technical field, in particular to a kind of method of determining fault throw
And device.
Background technique
The powerful pressure and tension generated in process of crustal movement is more than the intensity of rock stratum itself, causes rock stratum to occur disconnected
It splits, the sillar of time-parallel strata two sides is relatively moved along the plane of disruption and the structural configuration of formation, is properly termed as tomography.Fault development is early
Phase Diagenesis fluids form seam type hole reservoir along tomography flowing corrosion, become hydrocarbon pore volume, later period oil gas hangs down along tomography
To migration finally at hiding.Tomography and its near zone are usually the main distributed area of oil gas, are having for oil-gas exploration and well site deployment
Sharp target area.It can be seen that tomography with the formation of fractured-vuggy reservoir, the migration of oil gas at hide, distribution of oil gas etc. is close now
It is inseparable.Therefore, the research of fractured-vuggy reservoir interrupting layer is very important.And determining fault throw is study of fault rule
One of important evidences of fault features such as mould, fault sealing property, tomography associate structure.
Usually determining the main process of the method for fault throw in the prior art is: carrying out stratigraphic structure solution to seismic data
It releases, obtains corresponding structural map, tomography two sides stratigraphic structure data are obtained from structural map, according to tomography two sides stratigraphic structure number
According to determining fault throw.This method is mainly used for determining turn-off of the tomography two sides rock mass on vertical.However, being directed to fracture hole type
Carbonate Oil gas reservoir, internal tomography is based on strike-slip fault.Since the vertical changing of the relative positions of strike-slip fault is unobvious, cunning is walked
The turn-off of tomography walking sliding turn-off and being mainly horizontally oriented, the method for existing determining fault throw may be difficult based on construction
Figure accurately calculates the sliding turn-off of walking of strike-slip fault.
Summary of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device of determining fault throw, to improve determined tomography
The accuracy of turn-off.
In order to solve the above technical problems, the method and device that the embodiment of the present application provides a kind of determining fault throw is in this way
It realizes:
A kind of method of determining fault throw, is provided with the seismic data of interval of interest, which comprises
Based on the frequency range of the seismic data, spectral decomposition processing is carried out to the seismic data, obtains the mesh
Multiple earthquakes of interval divide data, and obtain specified earthquake from the multiple earthquake frequency dividing data and divide data；Wherein,
The fracture information that the specified earthquake frequency dividing data are used to reflect in the interval of interest；
Coherent processing is carried out to the specified earthquake frequency dividing data, obtains coherence data；Wherein, the coherence data includes
The coherence factor of multiple data points and the data point；
Based on the coherence data, the cutting relationship between the tomography in the interval of interest and a karst beading is determined；
According to the cutting relationship, the turn-off of the tomography is determined.
It is described to be based on the coherence data in preferred embodiment, determine tomography in the interval of interest and a karst beading it
Between cutting relationship, comprising:
Color mapping processing is carried out to the coherence data, obtains the color mapping treated coherence data；
Based on the color mapping treated coherence data, determine tomography in the interval of interest and a karst beading it
Between cutting relationship.
It is described that color mapping processing is carried out to the coherence data in preferred embodiment, after obtaining the color mapping processing
Coherence data, comprising:
For the coherence data, designated color component value is set；The designated color component value is according to the designated color
Type determines；
The new coherence data of the designated color component value will be provided with as the color mapping, and that treated is relevant
Data.
It is described that for the coherence data, designated color component value is set in preferred embodiment, comprising:
By coherence factor be 0 and coherence factor is type of the first numerical intervals between 0.4 according to the designated color
The first equal proportion is carried out to split point；By coherence factor be 0.5 and coherence factor is second value section between 1 according to described specified
The type of color carries out the second equal proportion and splits point；
The first numerical intervals and second equal proportion after being split point according to first equal proportion split point after the second number
It is worth section, is that numerical value is 0.5 in coherence factor and the coherence data between 0 to 0.4 by numerical value in the coherence data
Coherence factor between to 1 is mapped one by one with the designated color component value of the designated color type；And by the phase
Numerical value is the coherence factor between 0.4 to 0.5 in dry data, is reflected with the color component value of the black in the designated color type
It penetrates；Wherein, the color point for the white in coherence factor and the designated color type that numerical value is 0 and 1 in the coherence data
Magnitude mapping.
It is described based on the color mapping treated coherence data in preferred embodiment, it determines in the interval of interest
Cutting relationship between tomography and a karst beading, comprising:
Horizon slice processing is carried out to the color mapping treated coherence data, obtains the relevant number of specified horizon slice
According to；The specified horizon slice coherence data indicates the color mapping treated the relevant number of coherence data middle layer position section
According to；
Based on the specified horizon slice coherence data, the position of the interval of interest interrupting layer and a karst beading is identified；
According to the position of the interval of interest interrupting layer and a karst beading, the tomography and karst in the interval of interest are determined
Cutting relationship between a beading.
In preferred embodiment, it is described be based on the specified horizon slice coherence data, identify the interval of interest interrupting layer and
The position of a karst beading, comprising:
When the color component value in the specified horizon slice coherence data at the first sampling point position is less than or equal in advance
If when the first color threshold, using first sampling point position as the position of the tomography；Wherein, the position institute of the tomography
The graphic feature of composition is linear feature；
When the color component value in the specified horizon slice coherence data at the second sampling point position is greater than or equal in advance
If when the second color threshold, using second sampling point position as the position of the karst beading；Wherein, the karst beading
The graphic feature that is constituted of position be circular feature.
In preferred embodiment, the value range of default first color threshold is 0~50.
In preferred embodiment, the value range of default second color threshold is 200~256.
It is described according to the cutting relationship in preferred embodiment, determine the turn-off of the tomography, comprising:
According to the cutting relationship, the opposite shifting in the karst beading between the piecemeal figure of tomography two sides is determined
Dynamic distance；
Using the relative movement distance as the turn-off of the tomography.
A kind of device of determining fault throw, described device provide the seismic data of interval of interest；Described device includes: frequency
Spectral factorization processing module, Coherent processing module, cutting relationship determination module and turn-off determining module；Wherein,
The spectral decomposition processing module, for the frequency range based on the seismic data, to the seismic data into
Line frequency spectral factorization processing obtains multiple earthquakes frequency dividing data of the interval of interest, and from the multiple earthquake frequency dividing data
It obtains specified earthquake and divides data；Wherein, the fracture letter that the specified earthquake frequency dividing data are used to reflect in the interval of interest
Breath；
The Coherent processing module obtains coherence data for carrying out Coherent processing to the specified earthquake frequency dividing data；
Wherein, the coherence data includes the coherence factor of multiple data points and the data point；
The cutting relationship determination module, for being based on the coherence data, determine tomography in the interval of interest with
Cutting relationship between a karst beading；
The turn-off determining module, for determining the turn-off of the tomography according to the cutting relationship.
The embodiment of the present application provides a kind of method and device of determining fault throw, the frequency based on the seismic data
Range carries out spectral decomposition processing to the seismic data, obtains multiple earthquakes frequency dividing data of the interval of interest, and from institute
It states and obtains specified earthquake frequency dividing data in multiple earthquake frequency dividing data；Wherein, the specified earthquake frequency dividing data are for reflecting institute
State the fracture information in interval of interest；In this way, carrying out Coherent processing by dividing data to the specified earthquake, can be improved disconnected
The accuracy of identification of layer, so that the interval of interest can be accurately determined based on the coherence data obtained after the Coherent processing
In tomography and a karst beading between cutting relationship；Finally, the tomography can be accurately determined according to the cutting relationship
Turn-off.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of flow chart of the embodiment of the method for determining fault throw of the application；
Fig. 2 is the schematic diagram of coherence factor and color component value mapping in the embodiment of the present application；
Fig. 3 is the schematic diagram of the cutting relationship of the embodiment of the present application interrupting layer and a karst beading；
Fig. 4 is the flow chart for the Installation practice that the application determines fault throw.
Specific embodiment
The embodiment of the present application provides a kind of method and device of determining fault throw.
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
The embodiment of the present application provides a kind of method of determining fault throw.The method is provided with the earthquake of interval of interest
Data.
In the present embodiment, the interval of interest can be the underground interval in carbonate karst fracture cave reservoir.Institute
Stating interval of interest can be the still undetermined underground interval of fault throw.
In the present embodiment, the seismic data can be the knot of the earthquake frequency dividing data investigation of a variety of different frequencies
Fruit.The seismic data can be the set of the seismic data of multiple sampled points.The seismic data can be to be tieed up by main profile
The 3D data volume that degree, cross-track dimension and sampling time dimension are constituted.Three dimension pairwise orthogonals.The earthquake number
According to may include: amplitude at any sampling point position, and spectrum information corresponding with the seismic data.The frequency spectrum
Information may include the frequency range of the seismic data.The spectrum information can also include the dominant frequency of the seismic data.
Fig. 1 is a kind of flow chart of the embodiment of the method for determining fault throw of the application.As shown in Figure 1, the determination is disconnected
The method of fault throw, includes the following steps.
Step S101: the frequency range based on the seismic data carries out spectral decomposition processing to the seismic data, obtains
Multiple earthquakes to the interval of interest divide data, and obtain specified earthquake divider ratio from the multiple earthquake frequency dividing data
According to；Wherein, the fracture information that the specified earthquake frequency dividing data are used to reflect in the interval of interest.
In the present embodiment, the seismic data can be decomposed into pre- by the frequency range based on the seismic data
If the earthquake of the multi-frequency of spacing frequency divides data.For example, the frequency range of the seismic data can be 10~40 hertz
Hereby.The dominant frequency of the seismic data is 25 hertz.The preset interval frequency can be 5 hertz.It can be by the seismic data
It is decomposed into the earthquake that frequency is 10,15,20,25,30,35,40 hertz and divides data.Described each earthquake frequency dividing data can be with
Including single frequency content.The earthquake frequency dividing data are specifically as follows the 3D data volume of a certain frequency.
In the present embodiment, tuning thickness, which refers to, divides the minimum that data can be differentiated according to the earthquake of a certain frequency
Along thickness degree, usual value is a quarter that the earthquake of the frequency divides the wavelength of data.When a geologic body in interval of interest
Formation thickness close to the earthquake frequency dividing data of the frequency tuning thickness when, incident seismic wave at the geology body position and anti-
It penetrates seismic wave to interfere, amplitude enhancing, so that the geologic body can divide data in the earthquake comprising the frequency
Seismic data in be easily identified, for example, the seismic amplitude of the geologic body can be greater than other seismic datas amplitude.Due to
The tuning thickness of the earthquake frequency dividing data of different frequency is different, and therefore, dividing data according to the earthquake of different frequency can differentiate
The geologic body of different-thickness.Since the turn-off of strike-slip fault is smaller, the earthquake frequency dividing data of high frequency are chosen, can preferably be reflected
Fracture information in the interval of interest.Specified earthquake divider ratio is obtained in this way, can be divided in data from the multiple earthquake
According to.Wherein, the specified earthquake frequency dividing data can be 40 hertz.
Step S102: Coherent processing is carried out to the specified earthquake frequency dividing data, obtains coherence data；Wherein, the phase
Dry data include the coherence factor of multiple data points and the data point.
In the present embodiment, Coherent processing is carried out to the specified earthquake frequency dividing data, obtains coherence data, specifically may be used
To include, the designated layer position along the interval of interest can be calculated in the specified earthquake frequency dividing data in the designated layer position
It is default when window in specify coherence factor between adjacent seismic channel data, i.e. wave-form similarity.Wherein, the designated layer position can
Think any layer position in the interval of interest, it is described to specify adjacent seismic channel data that divide in data for the specified earthquake
Any group of adjacent seismic channel data at the designated layer position.The step of repeating the calculating coherence factor, available institute
State coherence data.
In the present embodiment, the coherence data can be the coherence factor by multiple data points and the data point
The three-dimensional coherent data volumes constituted.Wherein, the data point can be corresponded with the sampled point in the seismic data.
Step S103: it is based on the coherence data, determines cutting between tomography and the karst beading in the interval of interest
Cut relationship.
In the present embodiment, be based on the coherence data, determine tomography in the interval of interest and a karst beading it
Between cutting relationship, can specifically include, can to the coherence data carry out color mapping processing, the available color
Mapping treated coherence data.Based on the color mapping treated coherence data, can determine in the interval of interest
Tomography and a karst beading between cutting relationship.Wherein, the karst beading can indicate the karst in the interval of interest
Fracture hole.
In the present embodiment, color mapping processing is carried out to the coherence data, after obtaining the color mapping processing
Coherence data, can specifically include, can for the coherence data be arranged designated color component value.Wherein, the specified face
Color includes: white, grey and black.The designated color component value is determined according to the type of the designated color.For example, described
The type of designated color is 256 kinds, then the value of the designated color component value can be 0~255.It can will be provided with institute
The new coherence data of designated color component value is stated as the color mapping treated coherence data.
In the present embodiment, designated color component value is set for the coherence data, can specifically include, it can be by phase
Responsibility number is 0 and coherence factor is type progress first grade ratios of the first numerical intervals between 0.4 according to the designated color
Example is split point.It can be 0.5 by coherence factor and coherence factor is kind of the second value section between 1 according to the designated color
Class carries out the second equal proportion and splits point.Then, the first numerical intervals and described second etc. after being split point according to first equal proportion
Ratio split point after second value section, can by numerical value in the coherence data between 0 to 0.4 coherence factor and institute
Stating numerical value in coherence data is the coherence factor between 0.5 to 1, is carried out with the designated color component value of the designated color type
It maps one by one.It can be the coherence factor between 0.4 to 0.5 by numerical value in the coherence data, with the designated color type
In black color component value mapping.Wherein, the coherence factor and the specified face that numerical value is 0 and 1 in the coherence data
The color component value mapping of white in color type.For example, Fig. 2 is that coherence factor is reflected with color component value in the embodiment of the present application
The schematic diagram penetrated.
In the present embodiment, it based on the color mapping treated coherence data, determines in the interval of interest
Cutting relationship between tomography and a karst beading, can specifically include, can be to the color mapping treated coherence data
Horizon slice processing is carried out, specified horizon slice coherence data is obtained.Wherein, the specified horizon slice coherence data indicates institute
State color mapping treated the coherence data of coherence data middle layer position section.Based on the specified horizon slice coherence data,
It can identify the position of the interval of interest interrupting layer and a karst beading.According to the interval of interest interrupting layer and a karst beading
Position can determine the cutting relationship between the tomography in the interval of interest and a karst beading.
In the present embodiment, be based on the specified horizon slice coherence data, identify the interval of interest interrupting layer and
The position of a karst beading, can specifically include, when the face in the specified horizon slice coherence data at the first sampling point position
It, can be using first sampling point position as the position of the tomography when colouring component value is less than or equal to default first color threshold
It sets.Wherein, the graphic feature that the position of the tomography is constituted is linear feature.The value of default first color threshold
It may range from 0~50.When the color component value in the specified horizon slice coherence data at the second sampling point position is greater than
It, can be using second sampling point position as the position of the karst beading or when being equal to default second color threshold.Wherein,
The graphic feature that the position of the karst beading is constituted is circular feature.The value range of default second color threshold
It can be 200~256.
Step S104: according to the cutting relationship, the turn-off of the tomography is determined.
In the present embodiment, it according to the cutting relationship, determines the turn-off of the tomography, can specifically include: according to
The cutting relationship can determine the relative movement distance in the karst beading between the piecemeal figure of tomography two sides.
It can be using the relative movement distance as the turn-off of the tomography.
For example, Fig. 3 is the schematic diagram of the cutting relationship of the embodiment of the present application interrupting layer and a karst beading.In Fig. 3 black compared with
Thin lines indicate the position of strike-slip fault, and the annulus with one fixed width indicates the position of a karst beading in Fig. 3.Such as Fig. 3 institute
Show, after a karst beading are cut by strike-slip fault, due to the translation changing of the relative positions of strike-slip fault, is located at strike-slip fault two sides in a karst beading
Visuals be displaced, can be using mobile relative distance as the turn-off of strike-slip fault.
The embodiment of the method for the determining fault throw, based on the frequency range of the seismic data, to the earthquake number
According to spectral decomposition processing is carried out, multiple earthquakes frequency dividing data of the interval of interest are obtained, and from the multiple earthquake divider ratio
Specified earthquake frequency dividing data are obtained according to middle；Wherein, the specified earthquake frequency dividing data are used to reflect disconnected in the interval of interest
Split information；In this way, carrying out Coherent processing by dividing data to the specified earthquake, the accuracy of identification of tomography can be improved, from
And based on the coherence data obtained after the Coherent processing, the tomography and karst string in the interval of interest can be accurately determined
Cutting relationship between pearl；Finally, can accurately determine according to the cutting relationship and be located at tomography two in the karst beading
The distance of the visuals relative movement of side, and using the distance of the relative movement as the turn-off of the tomography.
Fig. 4 is the flow chart for the Installation practice that the application determines fault throw.Described device can provide interval of interest
Seismic data.As shown in figure 4, the apparatus may include: spectral decomposition processing module 100, is cut at Coherent processing module 200
Cut relationship determination module 300 and turn-off determining module 400.
The spectral decomposition processing module 100, can be used for the frequency range based on the seismic data, to the earthquake
Data carry out spectral decomposition processing, obtain multiple earthquakes frequency dividing data of the interval of interest, and divide from the multiple earthquake
Specified earthquake is obtained in data divides data；Wherein, the specified earthquake frequency dividing data are for reflecting in the interval of interest
Fracture information.
The Coherent processing module 200 can be used for carrying out Coherent processing to the specified earthquake frequency dividing data, obtain phase
Dry data；Wherein, the coherence data includes the coherence factor of multiple data points and the data point.
The cutting relationship determination module 300 can be used for determining in the interval of interest based on the coherence data
Cutting relationship between tomography and a karst beading.
The turn-off determining module 400 can be used for determining the turn-off of the tomography according to the cutting relationship.
The determining fault throw Installation practice is corresponding with the embodiment of the method for the determining fault throw, Ke Yishi
Now determine the technical solution of the embodiment of the method for fault throw, and the technical effect of adquisitiones embodiment.
In the 1990s, the improvement of a technology can be distinguished clearly be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow to be programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, programmable logic device
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, logic function determines device programming by user.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, designs and makes without asking chip maker
Dedicated IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " is patrolled
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but there are many kind, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed is most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only needing method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, so that it may it is readily available the hardware circuit for realizing the logical method process.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions for including in it can also be considered as in hardware component.Or
Even, can will be considered as realizing the device of various functions either the software module of implementation method can be Hardware Subdivision again
Structure in part.
Device that above-described embodiment illustrates, module can specifically realize by computer chip or entity, or by having certain
The product of function is planted to realize.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be embodied in the form of software products, in a typical configuration
In, calculating equipment includes one or more processors (CPU), input/output interface, network interface and memory.The computer is soft
Part product may include that some instructions are used so that a computer equipment (can be personal computer, server or network
Equipment etc.) execute method described in certain parts of each embodiment of the application or embodiment.The computer software product can
To be stored in memory, memory may include the non-volatile memory in computer-readable medium, random access memory
(RAM) and/or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer
The example of readable medium.Computer-readable medium includes that permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer readable instructions, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only memory
(ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory techniques, CD-ROM are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or
Other magnetic storage devices or any other non-transmission medium, can be used for storage can be accessed by a computing device information.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for device reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, clothes
Business device computer, handheld device or portable device, laptop device, multicomputer system, microprocessor-based system, set
Top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer including any of the above system or equipment
Distributed computing environment etc..
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and
Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's
Spirit.
Claims (8)
1. a kind of method of determining fault throw, which is characterized in that be provided with the seismic data of interval of interest, the method packet
It includes:
Based on the frequency range of the seismic data, spectral decomposition processing is carried out to the seismic data, obtains the target zone
Multiple earthquakes of section divide data, and obtain specified earthquake from the multiple earthquake frequency dividing data and divide data；Wherein, described
The fracture information that specified earthquake frequency dividing data are used to reflect in the interval of interest；
Coherent processing is carried out to the specified earthquake frequency dividing data, obtains coherence data；Wherein, the coherence data includes multiple
The coherence factor of data point and the data point；
Based on the coherence data, the cutting relationship between the tomography in the interval of interest and a karst beading is determined, comprising: right
The coherence data carries out color mapping processing, obtains the color mapping treated coherence data；It is reflected based on the color
Coherence data of penetrating that treated, determines the cutting relationship between the tomography in the interval of interest and a karst beading；
According to the cutting relationship, the turn-off of the tomography is determined, comprising: according to the cutting relationship, determine the karst string
Relative movement distance in pearl between the piecemeal figure of tomography two sides；Using the relative movement distance as the tomography
Turn-off.
2. a kind of method of determining fault throw according to claim 1, which is characterized in that described to the coherence data
Color mapping processing is carried out, the color mapping is obtained treated coherence data, comprising:
For the coherence data, designated color component value is set；The designated color component value is according to the type of the designated color
It determines；
The new coherence data of the designated color component value will be provided with as the color mapping treated coherence data.
3. a kind of method of determining fault throw according to claim 2, which is characterized in that described is the coherence data
Designated color component value is set, comprising:
By coherence factor be 0 and coherence factor is type progress of the first numerical intervals between 0.4 according to the designated color
First equal proportion is split point；By coherence factor be 0.5 and coherence factor is second value section between 1 according to the designated color
Type carry out the second equal proportion and split point；
The first numerical intervals and second equal proportion after being split point according to first equal proportion split point after second value area
Between, by numerical value in the coherence data be in coherence factor and the coherence data between 0 to 0.4 numerical value be 0.5 to 1 it
Between coherence factor, mapped one by one with the designated color component value of the designated color type；And by the relevant number
According to the coherence factor that middle numerical value is between 0.4 to 0.5, mapped with the color component value of the black in the designated color type；
Wherein, the color component value for the white in coherence factor and the designated color type that numerical value is 0 and 1 in the coherence data
Mapping.
4. a kind of method of determining fault throw according to claim 1, which is characterized in that described to be reflected based on the color
Coherence data of penetrating that treated, determines the cutting relationship between the tomography in the interval of interest and a karst beading, comprising:
Horizon slice processing is carried out to the color mapping treated coherence data, obtains specified horizon slice coherence data；
The specified horizon slice coherence data indicates the color mapping treated the coherence data of coherence data middle layer position section；
Based on the specified horizon slice coherence data, the position of the interval of interest interrupting layer and a karst beading is identified；
According to the position of the interval of interest interrupting layer and a karst beading, the tomography and a karst beading in the interval of interest are determined
Between cutting relationship.
5. a kind of method of determining fault throw according to claim 4, which is characterized in that described to be based on the specified edge
Layer slice coherence data, identifies the position of the interval of interest interrupting layer and a karst beading, comprising:
When the color component value in the specified horizon slice coherence data at the first sampling point position is less than or equal to default the
When one color threshold, using first sampling point position as the position of the tomography；Wherein, the position of the tomography is constituted
Graphic feature be linear feature；
When the color component value in the specified horizon slice coherence data at the second sampling point position is greater than or equal to default the
When second colors threshold value, using second sampling point position as the position of the karst beading；Wherein, the position of the karst beading
Constituted graphic feature is set as circular feature.
6. a kind of method of determining fault throw according to claim 5, which is characterized in that the default first color threshold
The value range of value is 0~50.
7. a kind of method of determining fault throw according to claim 5, which is characterized in that the default second color threshold
The value range of value is 200~256.
8. a kind of device of determining fault throw, which is characterized in that the seismic data of described device offer interval of interest；The dress
Set includes: spectral decomposition processing module, Coherent processing module, cutting relationship determination module and turn-off determining module；Wherein,
The spectral decomposition processing module carries out frequency to the seismic data for the frequency range based on the seismic data
Spectral factorization processing obtains multiple earthquakes frequency dividing data of the interval of interest, and obtains from the multiple earthquake frequency dividing data
Specified earthquake divides data；Wherein, the fracture information that the specified earthquake frequency dividing data are used to reflect in the interval of interest；
The Coherent processing module obtains coherence data for carrying out Coherent processing to the specified earthquake frequency dividing data；Its
In, the coherence data includes the coherence factor of multiple data points and the data point；
The cutting relationship determination module determines the tomography and karst in the interval of interest for being based on the coherence data
Cutting relationship between a beading；Wherein, the cutting relationship determination module is used to carry out at color mapping the coherence data
Reason, obtains the color mapping treated coherence data, and based on the color mapping treated coherence data, determines institute
State the cutting relationship between the tomography in interval of interest and a karst beading；
The turn-off determining module, for determining the turn-off of the tomography according to the cutting relationship；Wherein, the turn-off is true
Cover half block is used to be determined according to the cutting relationship opposite between the piecemeal figure of tomography two sides in the karst beading
Moving distance, and using the relative movement distance as the turn-off of the tomography.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710696640.4A CN107656309B (en) | 2017-08-15 | 2017-08-15 | A kind of method and device of determining fault throw |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710696640.4A CN107656309B (en) | 2017-08-15 | 2017-08-15 | A kind of method and device of determining fault throw |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107656309A CN107656309A (en) | 2018-02-02 |
| CN107656309B true CN107656309B (en) | 2019-09-06 |
Family
ID=61128572
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710696640.4A Active CN107656309B (en) | 2017-08-15 | 2017-08-15 | A kind of method and device of determining fault throw |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107656309B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107765318B (en) * | 2017-10-10 | 2019-06-11 | 中国石油天然气集团公司 | A kind of method and device of determining geologic horizon |
| CN109188518B (en) * | 2018-09-04 | 2019-06-28 | 中国矿业大学(北京) | The recognition methods of coal measure strata sandstone and system based on earthquake frequency splitting technology |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459801A (en) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | Coherence enhancement processing method used for recognizing fault |
| CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
| CN105425299A (en) * | 2015-11-17 | 2016-03-23 | 中国石油天然气股份有限公司 | Method and device for determining formation fracture distribution |
| CN105510964A (en) * | 2015-11-27 | 2016-04-20 | 中国石油大学(华东) | Seismic recognition method of low-order strike-slip faults in complex structural areas |
| CN106019370A (en) * | 2016-05-06 | 2016-10-12 | 中国石油天然气股份有限公司 | Information recommendation method and device |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1325935C (en) * | 2004-08-18 | 2007-07-11 | 中国石油天然气集团公司 | Superiority frequency band coherent handling method used for fine tomography explanation |
| CN105093300B (en) * | 2015-07-27 | 2017-10-17 | 中国石油天然气股份有限公司 | A kind of boundary recognition of geological body method and device |
| CN105334534A (en) * | 2015-10-21 | 2016-02-17 | 中国石油大学(华东) | Low order fault interpretation method based on construction mode guidance |
-
2017
- 2017-08-15 CN CN201710696640.4A patent/CN107656309B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459801A (en) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | Coherence enhancement processing method used for recognizing fault |
| CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
| CN105425299A (en) * | 2015-11-17 | 2016-03-23 | 中国石油天然气股份有限公司 | Method and device for determining formation fracture distribution |
| CN105510964A (en) * | 2015-11-27 | 2016-04-20 | 中国石油大学(华东) | Seismic recognition method of low-order strike-slip faults in complex structural areas |
| CN106019370A (en) * | 2016-05-06 | 2016-10-12 | 中国石油天然气股份有限公司 | Information recommendation method and device |
Non-Patent Citations (3)
| Title |
|---|
| "RGB分频技术在断块油藏断层识别中的应用";陈国飞 等;《复杂油气藏》;20150630;第8卷(第2期);第29-32、68页 |
| "分频相位分析技术在断层解释中的应用";王庆华 等;《工程地球物理学报》;20131130;第10卷(第6期);第896-899页 |
| "分频相干技术在复杂断裂解释中的应用";韩磊 等;《复杂油气藏》;20161231;第9卷(第4期);第16页左栏第1段-第21页右栏第1段、附图1-6 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107656309A (en) | 2018-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106324676B (en) | A kind of method and device of determining fault sealing property | |
| Neves et al. | Detection of potential fractures and small faults using seismic attributes | |
| CN106338763B (en) | A kind of determining seismic profile shows the method and device of data | |
| CN107942382B (en) | A kind of method and device removing sedimentation setting | |
| CN109100796A (en) | A kind of gas hydrates seismic data processing technique and device | |
| CA3137892C (en) | Identifying potential hydrocarbon traps in a subterranean region using recursive anisotropic erosion of seismic data | |
| CN106255902B (en) | Use the system and method for seismic data analysis geologic feature | |
| CN106353818B (en) | A kind of method and system of seismic data TEC time error correction | |
| CN107144878B (en) | A kind of fault recognizing method and device | |
| CN108252710A (en) | A kind of method and device of determining oil reservoir water saturation | |
| CN106873030B (en) | The determination method and apparatus of Fault Seal | |
| CN106896406A (en) | The method and apparatus that High-quality Reservoir is predicted based on impedance Domain Properties | |
| CN107656309B (en) | A kind of method and device of determining fault throw | |
| CN104375178B (en) | A kind of carbonate rock fractured cave method for predicting reservoir and device | |
| CN107179546B (en) | A kind of method and device of determining velocity field | |
| CN106970423B (en) | A kind of method and device of determining shale gas stratigraphic structure model | |
| Ngoc et al. | The application of seismic attributes for reservoir characterization in Pre-Tertiary fractured basement, Vietnam-Malaysia offshore | |
| CN108386185A (en) | A kind of method and device of determining sand body connectivity | |
| US10969507B2 (en) | Seismic channel geobody extraction for basin fairway simulation | |
| CN109143375A (en) | A kind of method and device in determining drilling well orientation | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| CN106845086B (en) | Formation pressure calculation method and device | |
| CN108710149A (en) | A kind of prediction technique and device of thin reservoir | |
| CN106772596B (en) | A kind of method and device of determining pre-stack time migration velocity field | |
| CN106908842B (en) | A kind of method and device of determining individual well sedimentary micro curve |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |