CN107238863A - A kind of method and device for judging arch structure authenticity - Google Patents
A kind of method and device for judging arch structure authenticity Download PDFInfo
- Publication number
- CN107238863A CN107238863A CN201710407665.8A CN201710407665A CN107238863A CN 107238863 A CN107238863 A CN 107238863A CN 201710407665 A CN201710407665 A CN 201710407665A CN 107238863 A CN107238863 A CN 107238863A
- Authority
- CN
- China
- Prior art keywords
- distortion
- arch structure
- domain
- depth
- piecewise fitting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/643—Horizon tracking
Abstract
The embodiment of the present application discloses a kind of method and device for judging arch structure authenticity.Methods described includes：The time-domain seismic profile data in purpose work area are obtained, carrying out layer thickness division to purpose work area obtains multiple layers of position, determines the time-domain horizon texture cross-sectional data of multiple layers of position；The time-domain horizon texture cross-sectional data of multiple layers of position is converted into Depth Domain horizon texture cross-sectional data；It is determined that distinguishing corresponding distortion part with multiple layers of position and distinguishing corresponding part of not distorting with multiple layers of position；It is determined that the distortion area of distortion part, and distance of the part away from preset reference line that do not distort；Determine a plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；Based on default criterion and relevant parameter, the arch structure authenticity of construction interval corresponding with a plurality of piecewise fitting straight line difference is judged.The judgement to arch structure authenticity in the stratum of work area can be realized.
Description
Technical field
The application is related to petroleum geology exploration technical field, more particularly to a kind of method for judging arch structure authenticity and
Device.
Background technology
During petroleum geology exploration, researcher usually requires first to obtain in the time-domain or Depth Domain in research work area
Reflection seismic waves data, then using obtain reflection seismic waves data to research work area carry out subsurface geological structure explanation,
To set up the subsurface geological structure model in research work area, so as to instruct follow-up actual geology investigation and prospecting.
For the work area with complex geological structure, because spread speed of the seismic wave in the stratum in the work area may be sent out
Raw vertical and horizontal change, frequently can lead to obtained reflection seismic waves data, such as TWT seismic profile data, show
Go out some false arch structure phenomenons, for example, the construction of the false anticline with antiform form.These false arch structures
Phenomenon may mislead understanding of the researcher to the subsurface geological structure in the work area.Therefore, need one kind badly and judge arch structure
The method of authenticity, so that researcher is accurately explained subsurface geological structure.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device for judging arch structure authenticity, to realize to work
The judgement of arch structure authenticity in area stratum.
In order to solve the above technical problems, the embodiment of the present application provides a kind of method and device for judging arch structure authenticity
It is realized in：
A kind of method for judging arch structure authenticity, including：
The time-domain seismic profile data in purpose work area are obtained, carrying out layer thickness division to the purpose work area obtains multiple layers
Position, based on the time-domain seismic profile data, determines the time-domain horizon texture cross-sectional data of the multiple layer position；
The time-domain horizon texture cross-sectional data of the multiple layer position is converted into Depth Domain horizon texture cross-sectional data；
Based on the Depth Domain horizon texture cross-sectional data, it is determined that with the multiple layer position respectively it is corresponding distortion part with
And with the multiple layer position corresponding part of not distorting respectively；The distortion area of the distortion part is determined, and it is described not abnormal
Become distance of the part away from preset reference line；
According to the distortion area and the distance of the part away from preset reference line that do not distort, determine that a plurality of piecewise fitting is straight
Line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；A segmentation in a plurality of piecewise fitting straight line
Fitting a straight line is corresponding with a construction interval in the purpose work area；One construction interval includes at least two layers of position；
Based on default criterion and the relevant parameter, judge and the corresponding structure of a plurality of piecewise fitting straight line difference
Make the arch structure authenticity of interval.
In preferred scheme, the time-domain horizon texture cross-sectional data by the multiple layer position is converted to Depth Domain layer position
Structural profile data, including：
According to default interval velocity, the time-domain horizon texture cross-sectional data of the multiple layer position is converted into the depth
Domain horizon texture cross-sectional data.
In preferred scheme, the default interval velocity is the conforming layer velocity field that speed is preset value.
Do not distort the distance away from preset reference line of part in preferred scheme, described in the determination, including：
According to the depth of stratum of the part of not distorting, and the depth of stratum residing for the preset reference line, institute is determined
State distance of the distortion part away from preset reference line.
Do not distort the depth of stratum of part in preferred scheme, described in the basis, and residing for the preset reference line
Depth of stratum, it is determined that the distance of the part away from preset reference line that do not distort, including：
, will part of not distorting corresponding with this layer of position when part of not distorting corresponding with a layer position is with plane-parallel
Depth of stratum subtract the absolute value after the depth of stratum residing for the preset reference line, be used as layer bit line corresponding with this layer of position
Do not distort distance of the part away from preset reference line；Or,
, will the portion that do not distort corresponding with this layer of position when part of not distorting corresponding with a layer position is not parallel with horizontal plane
Point and intersecting vertical range of all intersection points away from the preset reference line in distortion part average value, as corresponding with this layer of position
Do not distort part the distance away from preset reference line.
In preferred scheme, the distortion area for determining the distortion part, including：
According to the depth of stratum of the distortion part and the depth of stratum of the part of not distorting, the distortion part is determined
Distortion area.
In preferred scheme, it is described according to it is described distortion area and it is described do not distort part the distance away from preset reference line, really
Fixed a plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference, including：
Linear fit processing is carried out to the distortion area and the distance of the part away from preset reference line that do not distort, it is determined that
The a plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；The relevant parameter
Including：Slope and linear fit coefficient.
In preferred scheme, the default criterion includes：
When with the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences be more than or equal to 0.3
When, the arch structure with the construction interval of the piecewise fitting line correspondences is real；Or,
When the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences is less than 0.3 and slope
When absolute value is less than 0.01, the arch structure with the construction interval of the piecewise fitting line correspondences is false.
In preferred scheme, the default criterion also includes：
When with the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences be more than or equal to 0.3
And slope be more than 0 when, with the piecewise fitting line correspondences construction interval arch structure type for slip arch structure；Or
Person,
When with the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences be more than or equal to 0.3
And slope be less than 0 when, with the piecewise fitting line correspondences construction interval arch structure type for growth arch structure.
A kind of device for judging arch structure authenticity, described device includes：Time-domain layer position cross-sectional data determining module,
Time and depth transfer module, distance and area determining module, fitting a straight line determining module and arch structure authenticity judge module；Its
In,
The time-domain layer position cross-sectional data determining module, the time-domain seismic profile data for obtaining purpose work area,
Layer thickness division is carried out to the purpose work area and obtains multiple layers of position, based on the time-domain seismic profile data, is determined described many
The time-domain horizon texture cross-sectional data of individual layer position；
The time and depth transfer module, for the time-domain horizon texture cross-sectional data of the multiple layer position to be converted into depth
Domain horizon texture cross-sectional data；
The distance and area determining module, for based on the Depth Domain horizon texture cross-sectional data, it is determined that with it is described
Multiple layers of position distinguish corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively；Determine the distortion
Partial distortion area, and the distance of the part away from preset reference line that do not distort；
The fitting a straight line determining module, for according to it is described distortion area and it is described do not distort part away from preset reference line
Distance, determine a plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；It is described
A piecewise fitting straight line is corresponding with a construction interval in the purpose work area in a plurality of piecewise fitting straight line；One construction
Interval includes at least two layers of position；
The arch structure authenticity determining module, for based on default criterion and the relevant parameter, judge with
The a plurality of piecewise fitting straight line distinguishes the arch structure authenticity of corresponding construction interval.
In preferred scheme, the time-domain layer position cross-sectional data determining module includes：Seismic profile data acquisition module, layer
Position division module and horizon texture cross-sectional data determining module；Wherein,
The seismic profile data acquisition module, the time-domain seismic profile data for obtaining purpose work area；
The layer thickness division module, multiple layers of position are obtained for carrying out layer thickness division to the purpose work area；
The horizon texture cross-sectional data determining module, for based on the time-domain seismic profile data, it is determined that described
The time-domain horizon texture cross-sectional data of multiple layers of position.
In preferred scheme, the distance includes with area determining module：Layer position distortion determining module, distortion area determine mould
Block and do not distort apart from determining module；Wherein,
Layer position distortion determining module, for based on the Depth Domain horizon texture cross-sectional data, it is determined that with it is described many
Individual layer position respectively corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively；
The distortion area determining module, the distortion area for determining the distortion part；Specifically, according to it is described not
The depth of stratum of distortion part, and the depth of stratum residing for the preset reference line, it is determined that the part of not distorting is away from default
The distance of reference line；
It is described not distort apart from determining module, for determining the distance of the part away from preset reference line that do not distort；Specifically
Ground, according to the depth of stratum of the distortion part and the depth of stratum of the part of not distorting, determines the abnormal of the distortion part
Variable area.
The embodiment of the present application provides a kind of method and device for judging arch structure authenticity, can be to the purpose work
Area carries out layer thickness division and obtains multiple floor positions, based on the time-domain seismic profile data, it may be determined that the multiple layer position
Time-domain horizon texture cross-sectional data；The time-domain horizon texture cross-sectional data of the multiple layer position can be converted to Depth Domain
Horizon texture cross-sectional data；Based on the Depth Domain horizon texture cross-sectional data, it may be determined that right respectively with the multiple layer position
The distortion part answered and with the multiple layer position corresponding part of not distorting respectively；Determine the distortion face of the distortion part
Product, and the distance of the part away from preset reference line that do not distort；According to the distortion area and the part of not distorting away from pre-
If the distance of reference line, determine that a plurality of piecewise fitting straight line and associate corresponding with a plurality of piecewise fitting straight line difference are joined
Number；A piecewise fitting straight line can be relative with a construction interval in the purpose work area in a plurality of piecewise fitting straight line
Should, a plurality of piecewise fitting straight line can be corresponding with multiple construction intervals；One construction interval can include at least two layers of position；
Based on default criterion and the relevant parameter, it can be determined that with the corresponding techonosphere of a plurality of piecewise fitting straight line difference
The arch structure authenticity of section.Based on time-domain seismic profile data, the application method can be realized to anticline in the stratum of work area
Construct the judgement of authenticity.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart for the embodiment of the method for judging arch structure authenticity of the application；
Fig. 2 is the signal of the time-domain horizon texture cross-sectional data of the multiple layers of position in purpose work area in the embodiment of the present application
Figure；
Fig. 3 is the signal of the Depth Domain horizon texture cross-sectional data of the multiple layers of position in purpose work area in the embodiment of the present application
Figure；
Fig. 4 is the distortion area with different layers position corresponding distortion part respectively in the embodiment of the present application, and with difference
Distinguish the schematic diagram of corresponding distance of the part away from preset reference line that do not distort in layer position；
Fig. 5 is the schematic diagram of piecewise fitting straight line in the embodiment of the present application；
Fig. 6 is the composition structure chart for the device embodiment that the application judges arch structure authenticity；
Fig. 7 is the composition knot that the application judges time-domain layer position cross-sectional data in the device embodiment of arch structure authenticity
Composition；
Fig. 8 is the composition knot that the application judges time-domain layer position cross-sectional data in the device embodiment of arch structure authenticity
Composition.
Embodiment
The embodiment of the present application provides a kind of method and device for judging arch structure authenticity.
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
Fig. 1 is a kind of flow chart for the embodiment of the method for judging arch structure authenticity of the application.The judgement anticline structure
The method for making authenticity, comprises the following steps.
Step S101：The time-domain seismic profile data in purpose work area are obtained, layer thickness division is carried out to the purpose work area
Multiple layers of position are obtained, based on the time-domain seismic profile data, the time-domain horizon texture section of the multiple layer position are determined
Data.
Along the direction with horizontal plane, the purpose work area can be divided into multiple layers of position.Can between each layer of position
There is certain formation thickness.
By way of seismic prospecting and data acquisition, the time-domain seismic profile number in the purpose work area can be obtained
According to.Seismic horizon tracking can be carried out to the time-domain seismic profile data, it may be determined that first layer in the multiple layer position
TWT of the position in any layer position section horizontal position.The first layer position can be any layer in the multiple layer position
Position.According to first layer position any layer position section horizontal position TWT, it may be determined that the multiple layer position
Time-domain horizon texture cross-sectional data.The time-domain horizon texture cross-sectional data can include：Layer position section horizontal level
And TWT.The layer position section horizontal level is corresponded with the TWT.For example, Fig. 2 is that the application is real
Apply the schematic diagram of the time-domain horizon texture cross-sectional data of the multiple layers of position in purpose work area in example.Abscissa is the layer in Fig. 2
Position section horizontal level, ordinate is the TWT.In Fig. 2 white line represent in time-domain with the multiple layer
The corresponding layer bit line of position difference.Heavy black line represents the position of fault in time-domain in Fig. 2.
Step S102：The time-domain horizon texture cross-sectional data of the multiple layer position is converted into Depth Domain horizon texture to cut open
Face data.
Specifically, according to default interval velocity, the time-domain horizon texture cross-sectional data of the multiple layer position can be turned
It is changed to the Depth Domain horizon texture cross-sectional data.The Depth Domain horizon texture cross-sectional data can include：Layer position section water
Prosposition is put and depth of stratum.The layer position section horizontal level is corresponded with the depth of stratum.The default interval velocity
It can be the conforming layer velocity field that speed is preset value.The preset value can be 5,000 meter per seconds.
In one embodiment, following formula can be used by the layer position section horizontal position of a layer position
TWT is converted to the depth of stratum of this layer of position section horizontal position：
In formula (1), D is the depth of stratum of this layer of position section horizontal position, and t is this layer of position section horizontal position
TWT, V is the interval velocity that speed is the preset value.
Step S103：Based on the Depth Domain horizon texture cross-sectional data, it is determined that corresponding respectively with the multiple layer position
Distortion part and with the multiple layer position corresponding part of not distorting respectively；The distortion area of the distortion part is determined, with
And the distance of the part away from preset reference line that do not distort.
The preset reference line can be the horizontal line of a depth in the Depth Domain horizon texture section.The level
Line and plane-parallel.
In one embodiment, based on the Depth Domain horizon texture cross-sectional data, it is determined that dividing with the multiple layer position
Not corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively.Usual layer position can include many
Individual sampled point.One sampled point can be corresponded with layer position section horizontal level, a TWT.Based on described
Depth Domain horizon texture cross-sectional data, it may be determined that the corresponding layer position section horizontal level of any sampled point and double in layer position
When journey is travelled.Can set up using layer position section horizontal level be abscissa and TWT as ordinate coordinate system.At this
In coordinate system, a layer position is a layer bit line, and this layer of bit line can be made up of multiple sampled points.If first adopts in this layer of position
The slope for the straight line that sampling point is respectively constituted with two sampled points adjacent from first sampled point in this layer of position is different, can be by
First sampled point as this layer of position lopsided part；Otherwise, can using the first sampling point position described in this layer of position as
The not lopsided part of this layer of position.First sampled point can be any sampled point in this layer of position.
In one embodiment, according to the depth of stratum of the part of not distorting, and residing for the preset reference line
Depth of stratum, it may be determined that it is described do not distort part the distance away from preset reference line.
Further, do not distort the depth of stratum of part described in the basis, and the ground residing for the preset reference line
Layer depth, it is determined that the distance of the part away from preset reference line that do not distort.It can specifically include：When corresponding not with a layer position
When distortion part is with plane-parallel, the depth of stratum for not distorting part corresponding with this layer of position can be subtracted to the default ginseng
The absolute value after the depth of stratum residing for line is examined, distance of the part away from preset reference line that do not distort corresponding with this layer of position is used as.
Or, when part of not distorting corresponding with a layer position is not parallel with horizontal plane, can it not distorted corresponding with this layer of position
The average value of intersecting distance of all intersection points away from the preset reference line in part and distortion part, as corresponding with this layer of position
Do not distort distance of the part away from preset reference line.
For example, Fig. 3 is the Depth Domain horizon texture cross-sectional data of the multiple layers of position in purpose work area in the embodiment of the present application
Schematic diagram.Abscissa is layer position section horizontal level in Fig. 3, and ordinate is the depth of stratum.Black lines table in Fig. 3
Show the layer bit line corresponding with the multiple layer position difference in Depth Domain.Heavy black line represents the tomography in Depth Domain in Fig. 3
Position.1~20, B1~B4 is respectively different layers position of the purpose work area in Depth Domain in Fig. 3.Dotted line represents institute in Fig. 3
State preset reference line.
In one embodiment, according to the depth of stratum of the distortion part and the ground layer depth of the part of not distorting
Degree, it may be determined that the distortion area of the distortion part.
For example, Fig. 4 is the distortion area of distortion part corresponding with different layers position difference in the embodiment of the present application, Yi Jiyu
Different layers position respectively it is corresponding do not distort part the distance away from preset reference line schematic diagram.Fig. 4 grey areas St1And St14
The distortion area of distortion part corresponding with Fig. 3 middle levels position 1 and 14 difference of layer position is represented respectively.Ht in Fig. 41And Ht14Difference table
Show and distinguish corresponding distance of the part away from the preset reference line that do not distort with Fig. 3 middle levels position 1 and layer position 14, Hr represents described pre-
If the depth of stratum residing for reference line.
Step S104：According to the distortion area and the distance of the part away from preset reference line that do not distort, determine a plurality of
Piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference.
Linear fit processing can be carried out to the distortion area and the distance of the part away from preset reference line that do not distort,
Can determine a plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference.It is described
A piecewise fitting straight line is corresponding with a construction interval in the purpose work area in a plurality of piecewise fitting straight line.A plurality of segmentation
Fitting a straight line can be corresponding with multiple construction intervals.Wherein, a construction interval can include at least two layers of position.Described point
Section fitting a straight line can be associated with the distance and the distortion area of the part of not distorting away from preset reference line.
The relevant parameter can include：Slope and linearly dependent coefficient.The linearly dependent coefficient can be used to indicate that
Linear correlation degree between the distance and the distortion area of the part of not distorting away from preset reference line.It can use following
Formula calculates the linearly dependent coefficient with a piecewise fitting line correspondences：
In formula (1), R represents the linearly dependent coefficient with a piecewise fitting line correspondences, xiRepresent the piecewise fitting
Do not distort distance of the part away from preset reference line, y for i-th layer of position of line correspondencesiRepresent the i-th of the piecewise fitting line correspondences
The distortion area of individual layer position, n represents the layer position number of the piecewise fitting line correspondences.
For example, Fig. 5 is the schematic diagram of piecewise fitting straight line in the embodiment of the present application.In Fig. 5 abscissa for it is described it is vertical away from
From unit is km, and ordinate is the distortion area, and unit is sq-km.As shown in figure 5, pair with Fig. 3 in multiple layers of position
Corresponding layer bit line does not distort vertical height of the part away from preset reference line respectively, and corresponding respectively with multiple layers of position in Fig. 3
Layer bit line distortion part distortion area carry out linear fit processing, four piecewise fitting straight lines can be obtained.Further,
Intervals can be constructed by multiple layer thickness divisions as four using in Fig. 3 according to four piecewise fitting straight lines.Four construction intervals
Respectively suprastructure section, middle construction interval, infrastructure section and bottom construct interval.Wherein, piecewise fitting straight line can be with
Corresponding with a construction interval, four piecewise fitting straight lines can be corresponding with four construction intervals respectively.Suprastructure section
It is made up of layer position 1~4, middle construction interval is made up of layer position 5~14, infrastructure section is made up of layer position 15~19, bottom techonosphere
Section is made up of 20 and B1~B4 of layer position.
Step S105：Based on default criterion and the relevant parameter, judge and a plurality of piecewise fitting straight line point
The arch structure authenticity of not corresponding construction interval.
Specifically, the default criterion can include：When with the relevant parameter of a piecewise fitting line correspondences
When the absolute value of linearly dependent coefficient is more than or equal to 0.3, the arch structure with the construction interval of the piecewise fitting line correspondences
Can be real；When linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences absolute value be less than 0.3 and
When the absolute value of slope is less than 0.01, the arch structure with the construction interval of the piecewise fitting line correspondences can be untrue
's.Relevant parameter based on default criterion, with the first piecewise fitting line correspondences, it can be determined that straight with the first piecewise fitting
The arch structure authenticity of the corresponding construction interval of line.The first piecewise fitting straight line can be straight for a plurality of piecewise fitting
Any piecewise fitting straight line in line.
In one embodiment, the default criterion can also include：When with a piecewise fitting line correspondences
Relevant parameter in linearly dependent coefficient absolute value be more than or equal to 0.3 and slope be more than 0 when, with the piecewise fitting straight line pair
The arch structure type for the construction interval answered can be slippage arch structure；When associating ginseng with piecewise fitting line correspondences
In number the absolute value of linearly dependent coefficient be more than or equal to 0.3 and slope be less than 0 when, the construction with the piecewise fitting line correspondences
The arch structure type of interval can be growth arch structure.
For example, the linearly dependent coefficient of the corresponding piecewise fitting straight line of suprastructure section in Fig. 5 is 0.8857, slope
For -1.3048.The linearly dependent coefficient absolute value that the corresponding piecewise fitting straight line of interval is constructed in Fig. 5 is 0.7969, slope for-
0.1126.The linearly dependent coefficient absolute value of the corresponding piecewise fitting straight line of infrastructure section in Fig. 5 is 0.9989, and slope is
0.993.The linearly dependent coefficient absolute value of the corresponding piecewise fitting straight line of bottom construction interval in Fig. 5 is 0.0748, and slope is
0.0045., can be respectively to suprastructure section, middle construction interval, infrastructure section and bottom techonosphere according to default criterion
The authenticity of the arch structure of section is judged, it may be determined that the arch structure of suprastructure section and middle construction interval is true
, and arch structure type is growth anticline construction, the arch structure of infrastructure section is also true that, and arch structure class
Type is slips arch structure, and the arch structure of bottom construction interval is false.
The embodiment of the method for judging arch structure authenticity, can carry out layer thickness division to the purpose work area and obtain
Multiple layers of position, based on the time-domain seismic profile data, it may be determined that the time-domain horizon texture section of the multiple layer position
Data；The time-domain horizon texture cross-sectional data of the multiple layer position can be converted to Depth Domain horizon texture cross-sectional data；
Based on the Depth Domain horizon texture cross-sectional data, it may be determined that with the multiple layer position respectively it is corresponding distortion part and with
The multiple layer position corresponding part of not distorting respectively；Determine the distortion area of the distortion part, and the portion that do not distort
Divide the distance away from preset reference line；According to the distortion area and the distance of the part away from preset reference line that do not distort, it is determined that
A plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；The a plurality of piecewise fitting
A piecewise fitting straight line can be corresponding with a construction interval in the purpose work area in straight line, a plurality of piecewise fitting straight line
Can be corresponding with multiple construction intervals；One construction interval can at least include two layer positions；Based on default criterion and
The relevant parameter, it can be determined that the arch structure for constructing interval corresponding with a plurality of piecewise fitting straight line difference is true
Property.Based on time-domain seismic profile data, the application method can realize the judgement to arch structure authenticity in the stratum of work area.
Fig. 6 is the composition structure chart for the device embodiment that the application judges arch structure authenticity.As shown in fig. 6, described
Judging the device of arch structure authenticity can include：Time-domain layer position cross-sectional data determining module 100, time and depth transfer module
200th, distance and area determining module 300, fitting a straight line determining module 400 and arch structure authenticity judge module 500.
The time-domain layer position cross-sectional data determining module 100, the time-domain earthquake that can be used for obtaining purpose work area is cutd open
Face data, multiple layers of position can be obtained by carrying out layer thickness division to the purpose work area, based on the time-domain seismic profile data,
The time-domain horizon texture cross-sectional data of the multiple layer position can be determined.
The time and depth transfer module 200, can be used for turning the time-domain horizon texture cross-sectional data of the multiple layer position
It is changed to Depth Domain horizon texture cross-sectional data.
The distance and area determining module 300, can be used for being based on the Depth Domain horizon texture cross-sectional data, it is determined that
With the multiple layer position respectively corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively；Can be true
The distortion area of the fixed distortion part, and the distance of the part away from preset reference line that do not distort.
The fitting a straight line determining module 400, can be used for according to the distortion area and the part of not distorting away from pre-
If the distance of reference line, determine that a plurality of piecewise fitting straight line and associate corresponding with a plurality of piecewise fitting straight line difference are joined
Number；A piecewise fitting straight line can be relative with a construction interval in the purpose work area in a plurality of piecewise fitting straight line
Should；One construction interval can include at least two layers of position.
The arch structure authenticity judge module 500, can be used for based on default criterion and the relevant parameter,
It may determine that the arch structure authenticity for constructing interval corresponding with a plurality of piecewise fitting straight line difference.
Fig. 7 is the composition knot that the application judges time-domain layer position cross-sectional data in the device embodiment of arch structure authenticity
Composition.As shown in fig. 7, time-domain layer position cross-sectional data determining module 100 can include in Fig. 6：Seismic profile data acquisition mould
Block 110, layer thickness division module 120 and horizon texture cross-sectional data determining module 130.
The seismic profile data acquisition module 110, can be used for the time-domain seismic profile data for obtaining purpose work area.
The layer thickness division module 120, can be used for obtaining multiple layers of position to purpose work area progress layer thickness division.
The horizon texture cross-sectional data determining module 130, can be used for being based on the time-domain seismic profile data, really
The time-domain horizon texture cross-sectional data of fixed the multiple layer position.
Fig. 8 is the composition knot that the application judges time-domain layer position cross-sectional data in the device embodiment of arch structure authenticity
Composition.As shown in figure 8, distance can include with area determining module 300 in Fig. 6：Layer position distortion determining module 310, distortion face
Accumulate determining module 320 and do not distort apart from determining module 330.
The layer position distortion determining module 310, can be used for being based on the Depth Domain horizon texture cross-sectional data, it is determined that with
The multiple layer position respectively corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively.
The distortion area determining module 320, is determined for the distortion area of the distortion part.Specifically, may be used
With the depth of stratum for part of not distorted according to, and the depth of stratum residing for the preset reference line, it is determined that described not abnormal
Become distance of the part away from preset reference line.
It is described not distort apart from determining module 330, be determined for it is described do not distort part away from preset reference line away from
From.Specifically, it can be determined described abnormal according to the depth of stratum and the depth of stratum of the part of not distorting of the distortion part
Become the distortion area of part.
The method of the device embodiment for judging arch structure authenticity and the judgement arch structure authenticity is implemented
Example is corresponding, it is possible to achieve the present processes embodiment, and the technique effect of adquisitiones embodiment.
In the 1990s, for a technology improvement can clearly distinguish be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (for the improvement of 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 is 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, PLD
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make
Special IC chip 2.Moreover, nowadays, substitution manually makes IC chip, and this programming is also used instead mostly
" logic compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development
Seemingly, and the source code before compiling also handy specific programming language is write, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, 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 are 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 to method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
Controller can be implemented in any suitable manner, for example, controller can take such as microprocessor or processing
Device and storage can by the computer of the computer readable program code (such as software or firmware) of (micro-) computing device
Read medium, gate, switch, application specific integrated circuit (Application Specific Integrated Circuit,
ASIC), the form of programmable logic controller (PLC) and embedded microcontroller, the example of controller includes but is not limited to following microcontroller
Device：ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320, are deposited
Memory controller is also implemented as a part for the control logic of memory.
It is also known in the art that in addition to realizing controller in pure computer readable program code mode, it is complete
Controller can be caused with gate, switch, application specific integrated circuit, programmable by the way that method and step is carried out into programming in logic entirely
Logic controller realizes identical function with the form of embedded microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure in hardware component can also be considered as to the device for realizing various functions included in it.Or
Even, it not only can be able to will be the software module of implementation method for realizing that the device of various functions is considered as but also can be Hardware Subdivision
Structure in part.
System, device, module or unit that above-described embodiment is illustrated, can specifically be realized by computer chip or entity,
Or realized by the product with certain function.
For convenience of description, it is divided into various units during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Understood based on such, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft
Part product can include some instructions to cause a computer equipment (can be personal computer, server, or network
Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can
To be stored in internal memory, internal memory potentially includes the volatile memory in computer-readable medium, random access memory
And/or the form, such as read-only storage (ROM) or flash memory (flash RAM) such as Nonvolatile memory (RAM).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes 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 instruction, 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 transition internal memory (PRAM), static RAM
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only storage
(ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, the information that can be accessed by a computing device available for storage.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.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.It is real especially for system
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part 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.For example：Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
Positioned at including in the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many distortion and
Change is without departing from spirit herein, it is desirable to which appended claim includes these distortion and changed without departing from the application's
Spirit.
Claims (12)
1. a kind of method for judging arch structure authenticity, it is characterised in that including：
The time-domain seismic profile data in purpose work area are obtained, carrying out layer thickness division to the purpose work area obtains multiple layers of position,
Based on the time-domain seismic profile data, the time-domain horizon texture cross-sectional data of the multiple layer position is determined；
The time-domain horizon texture cross-sectional data of the multiple layer position is converted into Depth Domain horizon texture cross-sectional data；
Based on the Depth Domain horizon texture cross-sectional data, it is determined that with the multiple layer position respectively it is corresponding distortion part and with
The multiple layer position corresponding part of not distorting respectively；Determine the distortion area of the distortion part, and the portion that do not distort
Divide the distance away from preset reference line；
According to it is described distortion area and it is described do not distort part the distance away from preset reference line, determine a plurality of piecewise fitting straight line with
And with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；A piecewise fitting in a plurality of piecewise fitting straight line
Straight line is corresponding with a construction interval in the purpose work area；One construction interval includes at least two layers of position；
Based on default criterion and the relevant parameter, judge and the corresponding techonosphere of a plurality of piecewise fitting straight line difference
The arch structure authenticity of section.
2. it is according to claim 1 it is a kind of judge arch structure authenticity method, it is characterised in that it is described will be described many
The time-domain horizon texture cross-sectional data of individual layer position is converted to Depth Domain horizon texture cross-sectional data, including：
According to default interval velocity, the time-domain horizon texture cross-sectional data of the multiple layer position is converted into the Depth Domain layer
Bit architecture cross-sectional data.
3. a kind of method for judging arch structure authenticity according to claim 2, it is characterised in that the default layer speed
Degree field is the conforming layer velocity field that speed is preset value.
4. a kind of method for judging arch structure authenticity according to claim 1, it is characterised in that described in the determination
Do not distort the distance away from preset reference line of part, including：
According to the depth of stratum of the part of not distorting, and the depth of stratum residing for the preset reference line, determine described abnormal
Become distance of the part away from preset reference line.
5. a kind of method for judging arch structure authenticity according to claim 4, it is characterised in that described in the basis
Do not distort the depth of stratum of part, and the depth of stratum residing for the preset reference line, it is determined that the part of not distorting is away from pre-
If the distance of reference line, including：
When part of not distorting corresponding with a layer position is with plane-parallel, by the ground of part of not distorting corresponding with this layer of position
Layer depth subtracts the absolute value after the depth of stratum residing for the preset reference line, not abnormal as layer bit line corresponding with this layer of position
Become distance of the part away from preset reference line；Or,
When part of not distorting corresponding with layer position is not parallel with horizontal plane, corresponding with this layer of position will not distort part and
The average value of intersecting vertical range of all intersection points away from the preset reference line in distortion part, as it is corresponding with this layer of position not
Distance of the distortion part away from preset reference line.
6. a kind of method for judging arch structure authenticity according to claim 1, it is characterised in that described in the determination
The distortion area of distortion part, including：
According to the depth of stratum of the distortion part and the depth of stratum of the part of not distorting, the abnormal of the distortion part is determined
Variable area.
7. a kind of method for judging arch structure authenticity according to claim 1, it is characterised in that described in the basis
Distortion area and it is described do not distort part the distance away from preset reference line, determine a plurality of piecewise fitting straight line and with it is described a plurality of
Piecewise fitting straight line distinguishes corresponding relevant parameter, including：
Linear fit processing is carried out to the distortion area and the distance of the part away from preset reference line that do not distort, it is determined that described
A plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；The relevant parameter bag
Include：Slope and linear fit coefficient.
8. a kind of method for judging arch structure authenticity according to claim 7, it is characterised in that the default differentiation
Condition includes：
When being more than or equal to 0.3 with the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences,
Arch structure with the construction interval of the piecewise fitting line correspondences is real；Or,
When linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences absolute value be less than 0.3 and slope it is absolute
When value is less than 0.01, the arch structure with the construction interval of the piecewise fitting line correspondences is false.
9. a kind of method for judging arch structure authenticity according to claim 7, it is characterised in that the default differentiation
Condition also includes：
When more than or equal to 0.3 and oblique with the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences
When rate is more than 0, the arch structure type with the construction interval of the piecewise fitting line correspondences is slippage arch structure；Or,
When more than or equal to 0.3 and oblique with the absolute value of linearly dependent coefficient in the relevant parameter of a piecewise fitting line correspondences
When rate is less than 0, the arch structure type with the construction interval of the piecewise fitting line correspondences is growth arch structure.
10. a kind of device for judging arch structure authenticity, it is characterised in that described device includes：Time-domain layer position section number
Sentence according to determining module, time and depth transfer module, distance and area determining module, fitting a straight line determining module and arch structure authenticity
Disconnected module；Wherein,
The time-domain layer position cross-sectional data determining module, the time-domain seismic profile data for obtaining purpose work area, to institute
State purpose work area progress layer thickness division and obtain multiple layers of position, based on the time-domain seismic profile data, determine the multiple layer
The time-domain horizon texture cross-sectional data of position；
The time and depth transfer module, for the time-domain horizon texture cross-sectional data of the multiple layer position to be converted into Depth Domain layer
Bit architecture cross-sectional data；
The distance and area determining module, for based on the Depth Domain horizon texture cross-sectional data, it is determined that with it is the multiple
Layer position distinguish corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively；Determine the distortion part
Distortion area, and it is described do not distort part the distance away from preset reference line；
The fitting a straight line determining module, for according to it is described distortion area and it is described do not distort part away from preset reference line away from
From, determine a plurality of piecewise fitting straight line and with the corresponding relevant parameter of a plurality of piecewise fitting straight line difference；It is described a plurality of
A piecewise fitting straight line is corresponding with a construction interval in the purpose work area in piecewise fitting straight line；One construction interval
Including at least two layers of position；
The arch structure authenticity determining module, for based on default criterion and the relevant parameter, judge with it is described
A plurality of piecewise fitting straight line distinguishes the arch structure authenticity of corresponding construction interval.
11. a kind of device for judging arch structure authenticity according to claim 10, it is characterised in that the time-domain
Layer position cross-sectional data determining module includes：Seismic profile data acquisition module, layer thickness division module and horizon texture cross-sectional data
Determining module；Wherein,
The seismic profile data acquisition module, the time-domain seismic profile data for obtaining purpose work area；
The layer thickness division module, multiple layers of position are obtained for carrying out layer thickness division to the purpose work area；
The horizon texture cross-sectional data determining module, for based on the time-domain seismic profile data, determining the multiple
The time-domain horizon texture cross-sectional data of layer position.
12. it is according to claim 10 it is a kind of judge arch structure authenticity device, it is characterised in that the distance with
Area determining module includes：Layer position distortion determining module, distortion and do not distort apart from determining module at area determining module；Wherein,
The layer position distortion determining module, for based on the Depth Domain horizon texture cross-sectional data, it is determined that dividing with the multiple layer position
Not corresponding distortion part and with the multiple layer position corresponding part of not distorting respectively；
The distortion area determining module, the distortion area for determining the distortion part；Specifically, do not distorted according to described
Partial depth of stratum, and the depth of stratum residing for the preset reference line, it is determined that the part of not distorting is away from preset reference
The distance of line；
It is described not distort apart from determining module, for determining the distance of the part away from preset reference line that do not distort；Specifically, root
According to the depth of stratum and the depth of stratum of the part of not distorting of the distortion part, the distortion face of the distortion part is determined
Product.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710407665.8A CN107238863B (en) | 2017-06-02 | 2017-06-02 | A kind of method and device judging arch structure authenticity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710407665.8A CN107238863B (en) | 2017-06-02 | 2017-06-02 | A kind of method and device judging arch structure authenticity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107238863A true CN107238863A (en) | 2017-10-10 |
| CN107238863B CN107238863B (en) | 2019-01-18 |
Family
ID=59986002
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710407665.8A Active CN107238863B (en) | 2017-06-02 | 2017-06-02 | A kind of method and device judging arch structure authenticity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107238863B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2005239730A1 (en) * | 2005-12-06 | 2007-07-19 | Leven, James Howard Dr | Determination of the geometry and compartmentalization of underground hydrocarbon reservoirs using differential wide aperture seismic surveying. |
| US8780132B1 (en) * | 2009-03-21 | 2014-07-15 | Charles Saron Knobloch | Enhanced assimilation of orientation-dependent data in a multidimensional data volume |
| CN104074512A (en) * | 2013-03-26 | 2014-10-01 | 中国石油大学(北京) | Method for measuring accumulation probability of oil-gas anticlinal reservoir |
| CN105507905A (en) * | 2016-02-03 | 2016-04-20 | 长沙有色冶金设计研究院有限公司 | Internal discharge method during anticline passing period of open-cast mining |
-
2017
- 2017-06-02 CN CN201710407665.8A patent/CN107238863B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2005239730A1 (en) * | 2005-12-06 | 2007-07-19 | Leven, James Howard Dr | Determination of the geometry and compartmentalization of underground hydrocarbon reservoirs using differential wide aperture seismic surveying. |
| US8780132B1 (en) * | 2009-03-21 | 2014-07-15 | Charles Saron Knobloch | Enhanced assimilation of orientation-dependent data in a multidimensional data volume |
| CN104074512A (en) * | 2013-03-26 | 2014-10-01 | 中国石油大学(北京) | Method for measuring accumulation probability of oil-gas anticlinal reservoir |
| CN105507905A (en) * | 2016-02-03 | 2016-04-20 | 长沙有色冶金设计研究院有限公司 | Internal discharge method during anticline passing period of open-cast mining |
Non-Patent Citations (1)
| Title |
|---|
| 梁顺军 等: "高陡构造地震勘探圈闭要素畸变分析", 《石油物探》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107238863B (en) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107092031B (en) | A kind of sound wave measuring well curve bearing calibration and device | |
| CN106324676B (en) | A kind of method and device of determining fault sealing property | |
| CN107765298B (en) | A kind of method and device of determining near-surface velocity model | |
| CN106353818B (en) | A kind of method and system of seismic data TEC time error correction | |
| CN106054251B (en) | A kind of first break picking method and device | |
| US11906696B2 (en) | Point-vector based modeling of petroleum reservoir properties for a gridless reservoir simulation model | |
| CN107783185B (en) | A kind of processing method and processing device of tomographic statics | |
| CN106896406A (en) | The method and apparatus that High-quality Reservoir is predicted based on impedance Domain Properties | |
| CN107783186A (en) | A kind of near-surface velocity modeling method and device | |
| CN107179546B (en) | A kind of method and device of determining velocity field | |
| CN107829727A (en) | A kind of method and device for determining crustal stress | |
| CN107246900A (en) | A kind of method and device for determining oil-water interface position | |
| Berlin et al. | Steepened channels upstream of knickpoints: Controls on relict landscape response | |
| CN106970423A (en) | A kind of method and device for determining shale gas stratigraphic structure model | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| CN107329174A (en) | The method and device of well logging information is shown on seismic profile | |
| CN107144878A (en) | A kind of fault recognizing method and device | |
| CN104459774B (en) | A kind of geologic lithology difference identification method and system | |
| CN103543478A (en) | Geologic morphological interpolation KM (Kriging and Multiple-point geostatistics) method | |
| CN106908842B (en) | A kind of method and device of determining individual well sedimentary micro curve | |
| CN107238863B (en) | A kind of method and device judging arch structure authenticity | |
| CN106772596B (en) | A kind of method and device of determining pre-stack time migration velocity field | |
| CN106896419B (en) | A kind of method and device of determining petroleum resources amount | |
| CN106932839B (en) | A kind of method and device of pair of rock density and rate curve adjustment | |
| CN105467452B (en) | A kind of method and device that static correction is carried out to transitional region |
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 |