CN105137485B - The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap - Google Patents
The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap Download PDFInfo
- Publication number
- CN105137485B CN105137485B CN201510546827.7A CN201510546827A CN105137485B CN 105137485 B CN105137485 B CN 105137485B CN 201510546827 A CN201510546827 A CN 201510546827A CN 105137485 B CN105137485 B CN 105137485B
- Authority
- CN
- China
- Prior art keywords
- bitmap
- seismic profile
- value
- information
- geodetic coordinates
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to the automatic Picking method for transformation and device that information is explained in a kind of seismic profile bitmap, after the technical program pre-processes seismic profile bitmap, determine two-dimentional the geodetic coordinates value and time started value at section bitmap head and the tail both ends, and calculated and read automatically by small spacing and explain actual three-dimensional geodetic coordinates corresponding to information point, and same class interpretation results are distinguished and are merged into a text, it is easy to user to imported into hollow comparative analysis of seismic software and generates accurate map, the fast digitizing of information will be explained by realizing in seismic profile bitmap, so as to improve the accuracy and utilization rate of interpretation results data in seismic profile bitmap.
Description
Technical field
The present invention relates to exploration engineering field, explain that the automatic Picking of information turns in more particularly to a kind of seismic profile bitmap
Change method and device.
Background technology
In petroleum and natural gas industry, oil-gas reservoir is enriched among the ad-hoc location on underground set stratum.Modern petroleum
And exploration, the exploitation of natural gas largely need to determine the position of formation at target locations, buried depth, shape using seismic profile
State, and calculated for well location design with reserves.Stratigraphic structure form is thus rapidly and accurately definitely descended by seismic profile, is broken
Sheaf space feature is significant.
It is but low for original seismic cross-section piece utilization rate due to the limitation of data during seismic interpretation;Especially exist
In overseas exploration, the seismic profile frequently resulted in is all much the form of picture, it is difficult to forms main geologic layer to the personnel of explanation
Intuitively recognize on bit space, have impact on the actual effect of seismic prospecting.Mode in the past exactly prints section picture, passes through
Artificial read point is sampled, and this mode efficiency is low, and error is big, can not be combined into the interpretation results of existing geological data
Figure;Another way is that whole seismic profile is changed into geological data, i.e. SegY files, and this mode is to picture quality requirement
Height, debugging operation time are slow.
In order to which seismic interpretation personnel are freed from hard work, the degree of accuracy and efficiency, existing method are improved
Seismic profile achievement is subjected to high-density sampling in the mapping softwares such as Coreldraw, rapidly and accurately obtained using point in sky
Between on XYZ coordinate value, and the interpretation results of multiple sections can be merged into a space coordinates file, can generated accurate
Geologic structure diagram.But this method is disadvantageous in that:Original seismic cross-section piece large area distortion, when obscuring, achievement number
According to there may be certain error.
The content of the invention
To solve problem of the prior art, the automatic Picking that the present invention proposes to explain information in a kind of seismic profile bitmap turns
Change method and device, improve the accuracy and utilization rate of interpretation results data in seismic profile bitmap.
To achieve the above object, the invention provides the automatic Picking conversion side that information is explained in a kind of seismic profile bitmap
Method, including:
Seismic profile bitmap is pre-processed so that the section horizontal line of seismic profile bitmap is in horizontality, and
For the time shaft of seismic profile bitmap perpendicular to ground level, the border motif of seismic profile bitmap is rectangle;
Pretreated seismic profile bitmap is analyzed, between obtaining at the beginning of seismic profile bitmap value, earthquake cut open
Two-dimentional geodetic coordinates value corresponding to two-dimentional geodetic coordinates value corresponding to the starting endpoint of face and seismic profile distal point, utilizes earthquake
Two-dimentional geodetic coordinates value corresponding to section starting endpoint, two-dimentional geodetic coordinates value corresponding to seismic profile distal point obtain it is horizontal away from
From coefficient, using seismic profile bitmap at the beginning of between value obtain vertical range coefficient;
To explaining that information carries out vectorized process in pretreated seismic profile bitmap, and utilize the seismic profile position
Value determines horizontal reference line on the seismic profile bitmap after vector quantization between at the beginning of figure;Wherein, the horizontal reference line
Starting endpoint is overlapping with seismic profile starting endpoint, and the distal point of the horizontal reference line is overlapping with seismic profile distal point;
On the basis of the horizontal reference line, the vertical curve of the horizontal reference line is determined perpendicular to, obtains the level
The starting endpoint of reference line is to the bitmap horizontal range of the vertical curve, the vertical curve with explaining the intersection point of information to the water
The bitmap vertical range of flat reference line;
Using the horizontal reference line starting endpoint to the bitmap horizontal range of the vertical curve, horizontal range coefficient,
The vertical curve is with explaining the intersection point of information at the beginning of the bitmap vertical range of the horizontal reference line, seismic profile bitmap
Between value, two-dimentional geodetic coordinates value corresponding to vertical range coefficient, seismic profile starting endpoint and the vertical curve believe with explaining
The earth that the bitmap vertical range of the intersection point of breath to the horizontal reference line obtains intersection point of the vertical curve with explaining information is sat
Mark X values, geodetic coordinates Y value and time Z value;
Utilize geodetic coordinates X value of all vertical curves of identical layer position and tomography with the intersection point for explaining information, geodetic coordinates Y
Value and time Z value are merged into figure.
Preferably, the geodetic coordinates X values of the intersection point with explaining information that obtains the vertical curve, geodetic coordinates Y value and
The step of time Z value, includes:
Using the horizontal reference line starting endpoint to the bitmap horizontal range of the vertical curve, horizontal range coefficient,
The geodetic coordinates X of seismic profile starting endpoint, Y value obtain the geodetic coordinates X of intersection point of the vertical curve with explaining information, Y value;
The intersection point of value, vertical range coefficient, the vertical curve with explaining information between at the beginning of using seismic profile bitmap
Bitmap vertical range to the horizontal reference line obtains the time Z value of intersection point of the vertical curve with explaining information.
Preferably, the step of acquisition horizontal range coefficient includes:
Utilize two dimension corresponding to two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint and seismic profile distal point
Geodetic coordinates value determines the actual range between seismic profile starting endpoint and distal point;
Measurement obtains the bitmap distance in seismic profile bitmap between seismic profile starting endpoint and distal point;
Actual range divided by bitmap distance between the seismic profile starting endpoint and distal point obtain horizontal range system
Number.
Preferably, the step of vertical range coefficient includes:
Corresponding points are worth as starting point using the time started, obtain fixed time period end corresponding points and time started value corresponding points it
Between bitmap distance in seismic profile bitmap;
Position between the fixed time period divided by fixed time period end corresponding points and time started value corresponding points
Map distance obtains vertical range coefficient.
Preferably, it is described to explaining the step of information carries out vectorized process bag in pretreated seismic profile bitmap
Include:
Pretreated seismic cross-section is imported in vector software, using line of vector redraw the layer of seismic profile
Position and tomography;
The layer position of different earthquake section, tomography make a distinction according to different colors.
To achieve the above object, present invention also offers the automatic Picking conversion that information is explained in a kind of seismic profile bitmap
Device, including:
Pretreatment unit, for being pre-processed to seismic profile bitmap so that the section horizontal line of seismic profile bitmap
In horizontality, and the time shaft of seismic profile bitmap, perpendicular to ground level, the border motif of seismic profile bitmap is rectangle;
Proportionality coefficient acquiring unit, for analyzing pretreated seismic profile bitmap, obtain seismic profile position
It is worth between at the beginning of figure, two corresponding to two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint and seismic profile distal point
Geodetic coordinates value is tieed up, is greatly sat using corresponding to geodetic coordinates value corresponding to seismic profile starting endpoint, seismic profile distal point
Scale value obtain horizontal range coefficient, using seismic profile bitmap at the beginning of between value obtain vertical range coefficient;
Seismic profile bitmap vector unit, for explaining that information carries out vector in pretreated seismic profile bitmap
Change is handled, and horizontal ginseng is determined on the seismic profile bitmap after vector quantization using value between at the beginning of the seismic profile bitmap
Examine line;Wherein, the starting endpoint of the horizontal reference line is overlapping with seismic profile starting endpoint, the end of the horizontal reference line
Point is overlapping with seismic profile distal point;
Information dot picking first module is explained, on the basis of the horizontal reference line, being determined perpendicular to the level
The vertical curve of reference line, the starting endpoint of the horizontal reference line is obtained to the bitmap horizontal range of the vertical curve, described perpendicular
Straight line is with explaining the intersection point of information to the bitmap vertical range of the horizontal reference line;
Information dot picking second unit is explained, for the starting endpoint using the horizontal reference line to the vertical curve
The bitmap of intersection point to the horizontal reference line of bitmap horizontal range, horizontal range coefficient, the vertical curve with explaining information hangs down
Value between at the beginning of straight distance, seismic profile bitmap, two dimension the earth seat corresponding to vertical range coefficient, seismic profile starting endpoint
Scale value obtains geodetic coordinates X values, geodetic coordinates Y value and the time Z value of intersection point of the vertical curve with explaining information;
Information pickup unit is explained, for all vertical curves using identical layer position and tomography and the intersection point for explaining information
Geodetic coordinates X values, geodetic coordinates Y value and time Z value are merged into figure.
Preferably, the explanation information dot picking second unit includes:
The two-dimentional geodetic coordinates value determining module of information point is explained, for being arrived using the starting endpoint of the horizontal reference line
In the bitmap horizontal range of the vertical curve, horizontal range coefficient, two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint
Geodetic coordinates X, Y value obtain the geodetic coordinates X of the intersection point of the vertical curve with explaining information, Y value;
Explain the time Z value determining module of information point, for using be worth between at the beginning of seismic profile bitmap, vertical range
The bitmap vertical range of intersection point to the horizontal reference line of coefficient, the vertical curve with explaining information obtain the vertical curve and
Explain the time Z value of the intersection point of information.
Preferably, the proportionality coefficient acquiring unit includes:
Actual range acquisition module, for utilizing two-dimentional geodetic coordinates value and earthquake corresponding to seismic profile starting endpoint
Two-dimentional geodetic coordinates value determines the actual range between seismic profile starting endpoint and distal point corresponding to cross sections end point;
Bitmap is obtained between seismic profile starting endpoint and distal point in seismic profile position apart from acquisition module for measuring
Bitmap distance in figure;
Horizontal range coefficient module, for the actual range between the seismic profile starting endpoint and distal point divided by position
Map distance obtains horizontal range coefficient.
Preferably, the proportionality coefficient acquiring unit includes:
Bitmap time gap acquisition module, for being worth corresponding points as starting point using the time started, obtain fixed time period end
Bitmap distance between corresponding points and time started value corresponding points in seismic profile bitmap;
Vertical range coefficient module, for the fixed time period divided by fixed time period end corresponding points with starting
Bitmap distance between time value corresponding points in seismic profile bitmap obtains vertical range coefficient.
Preferably, the seismic profile bitmap vector unit includes:
Drafting module, for pretreated seismic cross-section to be imported in vector software, weight is carried out using line of vector
Paint layer position and the tomography of seismic profile;
Staining module, layer position, tomography for different earthquake section make a distinction according to different colors.
Above-mentioned technical proposal has the advantages that:After the technical program pre-processes seismic profile bitmap, it is determined that cuing open
The two-dimentional geodetic coordinates value at face bitmap head and the tail both ends and time started value, and calculated automatically by small spacing and read explanation information
Actual three-dimensional geodetic coordinates corresponding to point, and same class interpretation results are distinguished and are merged into a text, it is easy to user to lead
Enter to hollow comparative analysis of seismic software and generate accurate map, the quick number of information will be explained by realizing in seismic profile bitmap
Word, so as to improve the accuracy and utilization rate of interpretation results data in seismic profile bitmap.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the automatic Picking method for transformation flow chart that information is explained in a kind of seismic profile bitmap proposed by the present invention;
Fig. 2 is the schematic diagram before the seismic profile bitmap pretreatment of the present embodiment;
Fig. 3 is to the pretreated schematic diagram of bitmap shown in Fig. 2;
Fig. 4 is that the geodetic coordinates X values, geodetic coordinates Y value and time Z value that information point is explained in the present embodiment calculate signal
Figure;
Fig. 5 is the present embodiment seismic profile line plan-position distribution map;
Fig. 6 is the pickup schematic diagram that the present embodiment explains information point;
Fig. 7 is the automatic Picking reforming unit block diagram that information is explained in a kind of seismic profile bitmap proposed by the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other under the premise of creative work is not made
Embodiment, belong to the scope of protection of the invention.
As shown in figure 1, to explain the automatic Picking method for transformation of information in a kind of seismic profile bitmap proposed by the present invention
Flow chart.Including:
Step 101):Seismic profile bitmap is pre-processed so that the section horizontal line of seismic profile bitmap is in water
Level state, and the time shaft of seismic profile bitmap, perpendicular to ground level, the border motif of seismic profile bitmap is rectangle;
In a step 101, using the software processing seismic profile bitmap picture such as Photoshop, real water in section is tried to achieve
Angle between horizontal line and section horizontal line.As shown in Fig. 2 it is the signal before the seismic profile bitmap pretreatment of the present embodiment
Figure.Section horizontal line is changed into horizontality by rotating;As picture has certain distortion, also to be reduced so that cut open
Face border motif remains rectangle, and the time shaft of seismic profile bitmap is perpendicular to ground level, the border shape of seismic profile bitmap
State is rectangle.As shown in figure 3, for the pretreated schematic diagram of bitmap shown in Fig. 2.
Step 102):Pretreated seismic profile bitmap is analyzed, between obtaining at the beginning of seismic profile bitmap
Two-dimentional geodetic coordinates corresponding to two-dimentional geodetic coordinates value corresponding to value, seismic profile starting endpoint and seismic profile distal point
Value, utilize two-dimentional geodetic coordinates corresponding to two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint, seismic profile distal point
Value obtain horizontal range coefficient, using seismic profile bitmap at the beginning of between value obtain vertical range coefficient;
As shown in figure 5, it is the present embodiment seismic profile line plan-position distribution map.In a step 102, in step 101
On the basis of, pretreated seismic profile bitmap and its related data (such as seismic survey lines distribution map) are analyzed, obtain ground
Shake two-dimentional geodetic coordinates value corresponding to section head and the tail two-end-point and its section time started value;Wherein, two-dimentional geodetic coordinates value bag
Include geodetic coordinates X values and geodetic coordinates Y value.
On in terms of level in earthquake three-dimensional coordinate system, pass through geodetic coordinates X corresponding to seismic profile head and the tail two-end-point
Value, geodetic coordinates Y value obtain the actual range between head and the tail two-end-point, by aiding in software to measure head and the tail two-end-point in bitmap
Bitmap distance, actual range divided by bitmap distance are obtained into the horizontal range coefficient of the seismic profile bitmap;In earthquake three-dimensional
In vertical direction in coordinate system, appoint and take 1000ms as Temporal-distance value, be i.e. the time value of starting point is 200ms, distal point
Time value be 1200ms, and measure the bitmap distance between the starting point and distal point on seismic profile bitmap, will be above-mentioned
1000ms divided by bitmap distance be worth to vertical range coefficient.
Step 103):To explaining that information carries out vectorized process in pretreated seismic profile bitmap, and described in utilization
Value determines horizontal reference line on the seismic profile bitmap after vector quantization between at the beginning of seismic profile bitmap;Wherein, the water
The starting endpoint of flat reference line is overlapping with seismic profile starting endpoint, distal point and the seismic profile end of the horizontal reference line
Point is overlapping;
In step 103, pretreated seismic profile bitmap in step 101 is imported in vector software, in this reality
Apply in example, vector software Coreldraw.By line of vector redraw layer position and the tomography of seismic profile bitmap, it is different
Layer position or tomography make a distinction according to color.In addition, draw a horizontal reference that time started value is obtained through step 102
Line, the seismic profile head and the tail two-end-point that the head and the tail two-end-point of the horizontal reference line strictly confirms with step 102 correspond to overlapping respectively.
Step 104):On the basis of the horizontal reference line, the vertical curve of the horizontal reference line is determined perpendicular to, is obtained
Intersection point of the starting endpoint of the horizontal reference line to the bitmap horizontal range of the vertical curve, the vertical curve and explanation information
To the bitmap vertical range of the horizontal reference line;
As shown in fig. 6, the pickup schematic diagram of information point is explained for the present embodiment.In figure 6, a represents horizontal reference line
Starting endpoint is to the bitmap horizontal range of vertical curve, and b expression vertical curves are with explaining the intersection point of information to the bitmap of horizontal reference line
Vertical range.
Step 105):Using the horizontal reference line starting endpoint to the vertical curve bitmap horizontal range, level
Bitmap vertical range of the intersection point to the horizontal reference line, the seismic profile position of distance coefficient, the vertical curve with explaining information
Value between at the beginning of figure, two-dimentional geodetic coordinates value corresponding to vertical range coefficient, seismic profile starting endpoint and described vertical
The bitmap vertical range of intersection point to the horizontal reference line of the line with explaining information obtains the vertical curve with explaining the friendship of information
Geodetic coordinates X values, geodetic coordinates Y value and the time Z value of point;
As shown in figure 4, to explain geodetic coordinates X values, geodetic coordinates Y value and time Z the value meter of information point in the present embodiment
Calculate schematic diagram.As shown in figure 4, to explain geodetic coordinates X values, geodetic coordinates Y value and time Z the value meter of information point in the present embodiment
Calculate schematic diagram.With reference to Fig. 6, the actual the earth XY coordinates of certain point and time Z value are calculated as follows formula 1, formula 2, formula 3.
Wherein,The geodetic coordinates X values for the target point to be calculated are represented,Represent the earth of first point on seismic profile
Coordinate X values,The geodetic coordinates X values of tail point on seismic profile are represented,Represent in tail point corresponding to vector software Bitmap
X-coordinate value,Represent in first point X-coordinate value corresponding to vector software Bitmap,Represent the interpretation results in vector software
Corresponding X-coordinate value.
Wherein,The geodetic coordinates Y value for the target point to be calculated is represented,Represent the earth of first point on seismic profile
Coordinate Y value,The geodetic coordinates Y value of tail point on seismic profile is represented,Represent in tail point corresponding to vector software Bitmap
Y-coordinate value,Represent in first point Y-coordinate value corresponding to vector software Bitmap,Represent the interpretation results in vector software
Corresponding Y-coordinate value.
Calculating for Z values, using formula 3:
Wherein,The actual time value for the target point to be calculated is represented,Seismic profile time started value is represented,
Time value corresponding to seismic profile the preceding paragraph set time length ends is represented, () represent in vector software Bitmap
Bitmap distance value corresponding to corresponding one section of set time,Represent in vector software between vertical curve and interpretation results
Bitmap distance value of the intersection point to horizontal reference line.
Vertical curve is constantly moved to section tail point direction, often movement once, obtains a vertical line and interpretation results
Corresponding intersection point, geodetic coordinates X values, geodetic coordinates Y value and the time Z value of each intersection point are obtained in the manner described above, and with one
Fixed text formatting (such as seismic interpretation software Landmark layer bit format) is exported into text.
Step 106):Using all vertical curves of identical layer position and tomography with explain information intersection point geodetic coordinates X values,
Geodetic coordinates Y value and time Z value are merged into figure.
Accordingly, as shown in fig. 7, to explain that the automatic Picking of information turns in a kind of seismic profile bitmap proposed by the present invention
Change device block diagram.Including:
Pretreatment unit 701, for being pre-processed to seismic profile bitmap so that the section of seismic profile bitmap is horizontal
Line is in horizontality, and the time shaft of seismic profile bitmap, perpendicular to ground level, the border motif of seismic profile bitmap is square
Shape;
Proportionality coefficient acquiring unit 702, for analyzing pretreated seismic profile bitmap, obtain seismic profile
It is worth between at the beginning of bitmap, corresponding to two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint and seismic profile distal point
Two-dimentional geodetic coordinates value, utilize the earth corresponding to geodetic coordinates value corresponding to seismic profile starting endpoint, seismic profile distal point
Coordinate value obtain horizontal range coefficient, using seismic profile bitmap at the beginning of between value obtain vertical range coefficient;
Seismic profile bitmap vector unit 703, for explaining that information is carried out in pretreated seismic profile bitmap
Vectorized process, and using the seismic profile bitmap at the beginning of between value determine water on the seismic profile bitmap after vector quantization
Flat reference line;Wherein, the starting endpoint of the horizontal reference line is overlapping with seismic profile starting endpoint, the horizontal reference line
Distal point is overlapping with seismic profile distal point;
Information dot picking first module 704 is explained, on the basis of the horizontal reference line, being determined perpendicular to the water
The vertical curve of flat reference line, the starting endpoint of the horizontal reference line is obtained to the bitmap horizontal range of the vertical curve, described
Vertical curve is with explaining the intersection point of information to the bitmap vertical range of the horizontal reference line;
Explain information dot picking second unit 705, for the starting endpoint using the horizontal reference line to it is described vertically
The bitmap horizontal range of line, horizontal range coefficient, the vertical curve are with explaining the intersection point of information to the position of the horizontal reference line
Value between at the beginning of figure vertical range, seismic profile bitmap, two dimension is big corresponding to vertical range coefficient, seismic profile starting endpoint
Ground coordinate value obtains geodetic coordinates X values, geodetic coordinates Y value and the time Z value of intersection point of the vertical curve with explaining information;
Information pickup unit 706 is explained, for all vertical curves using identical layer position and tomography and the friendship for explaining information
Geodetic coordinates X values, geodetic coordinates Y value and the time Z value of point are merged into figure.
In the present embodiment, time started value is 200ms, proportionality coefficient 9.8814, calculates 111 intersection points altogether, puts distance
For 2.5mm, only several seconds are taken.Partial results such as table 1, the text formatting of generation give tacit consent to introduced layer position for Landmark interpretation horizons
Data format.In table 1, data row from left to right are respectively geodetic coordinates X values, geodetic coordinates Y value, time Z value.
Table 1
The technical program directly utilizes seismic profile picture, it is not necessary to prints, it is possible to which quick read is construed to
X values, Y value, Z values corresponding to fruit point.In addition, the technical program can largely be accurately read the X values, Y value, Z in corresponding space
Value, and different interpretation results are classified, generate different texts.Further, by the explanation of different section pictures
Achievement merges, and generates some interpretation results in a wide range of interior text, and it is soft to be input to Landmark explanations
It is independent in part to be unified into figure to related seismic data cube.
Above-described embodiment, the purpose of the present invention, technical scheme and beneficial effect are carried out further
Describe in detail, should be understood that the embodiment that the foregoing is only the present invention, be not intended to limit the present invention
Protection domain, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., all should include
Within protection scope of the present invention.
Claims (8)
1. the automatic Picking method for transformation of information is explained in a kind of seismic profile bitmap, it is characterised in that including:
Seismic profile bitmap is pre-processed so that the section horizontal line of seismic profile bitmap is in horizontality, and earthquake
For the time shaft of section bitmap perpendicular to ground level, the border motif of seismic profile bitmap is rectangle;
Pretreated seismic profile bitmap is analyzed, value, seismic profile rise between obtaining at the beginning of seismic profile bitmap
Two-dimentional geodetic coordinates value corresponding to two-dimentional geodetic coordinates value corresponding to top point and seismic profile distal point, utilizes seismic profile
Two-dimentional geodetic coordinates value corresponding to starting endpoint, two-dimentional geodetic coordinates value obtains horizontal range system corresponding to seismic profile distal point
Number, using seismic profile bitmap at the beginning of between value obtain vertical range coefficient;
To explaining that information carries out vectorized process in pretreated seismic profile bitmap, and utilize the seismic profile bitmap
Time started value determines horizontal reference line on the seismic profile bitmap after vector quantization;Wherein, the starting of the horizontal reference line
End points is overlapping with seismic profile starting endpoint, and the distal point of the horizontal reference line is overlapping with seismic profile distal point;
On the basis of the horizontal reference line, the vertical curve of the horizontal reference line is determined perpendicular to, obtains the horizontal reference
The starting endpoint of line is joined to the bitmap horizontal range of the vertical curve, the vertical curve and the intersection point of explanation information to the level
Examine the bitmap vertical range of line;
Using the starting endpoint of the horizontal reference line to the bitmap horizontal range of the vertical curve, horizontal range coefficient, described
Vertical curve and explain the intersection point of information between at the beginning of the bitmap vertical range of the horizontal reference line, seismic profile bitmap
Two-dimentional geodetic coordinates value corresponding to value, vertical range coefficient, seismic profile starting endpoint and the vertical curve are with explaining information
The bitmap vertical range of intersection point to the horizontal reference line obtain the geodetic coordinates X of the intersection point of the vertical curve with explaining information
Value, geodetic coordinates Y value and time Z value;
Using all vertical curves of identical layer position and tomography with explain the geodetic coordinates X values of intersection point of information, geodetic coordinates Y value and
Time Z value is merged into figure;
The step of acquisition horizontal range coefficient, includes:Using two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint with
And two-dimentional geodetic coordinates value corresponding to seismic profile distal point determine between seismic profile starting endpoint and distal point it is actual away from
From;Measurement obtains the bitmap distance in seismic profile bitmap between seismic profile starting endpoint and distal point;The earthquake is cutd open
Actual range divided by bitmap distance between face starting endpoint and distal point obtain horizontal range coefficient.
2. the method as described in claim 1, it is characterised in that the acquisition vertical curve is big with the intersection point of explanation information
The step of ground coordinate X values, geodetic coordinates Y value and time Z value, includes:
Using the starting endpoint of the horizontal reference line to the bitmap horizontal range of the vertical curve, horizontal range coefficient, earthquake
The geodetic coordinates X of section starting endpoint, Y value obtain the geodetic coordinates X of intersection point of the vertical curve with explaining information, Y value;
Value between at the beginning of using seismic profile bitmap, vertical range coefficient, the vertical curve are with explaining the intersection point of information to institute
The bitmap vertical range for stating horizontal reference line obtains the time Z value of intersection point of the vertical curve with explaining information.
3. method as claimed in claim 1 or 2, it is characterised in that include the step of the vertical range coefficient:
Using the time started be worth corresponding points as starting point, obtain fixed time period end corresponding points and the time started value corresponding points between
Bitmap distance in seismic profile bitmap;
The fixed time period divided by fixed time period end corresponding points and the time started value corresponding points between bitmap away from
From obtaining vertical range coefficient.
4. the method as described in claim 1, it is characterised in that described to explaining information in pretreated seismic profile bitmap
The step of carrying out vectorized process includes:
By pretreated seismic cross-section import vector software in, using line of vector redraw seismic profile layer position and
Tomography;
The layer position of different earthquake section, tomography make a distinction according to different colors.
5. the automatic Picking reforming unit of information is explained in a kind of seismic profile bitmap, it is characterised in that including:
Pretreatment unit, for being pre-processed to seismic profile bitmap so that the section horizontal line of seismic profile bitmap is in
Horizontality, and the time shaft of seismic profile bitmap, perpendicular to ground level, the border motif of seismic profile bitmap is rectangle;
Proportionality coefficient acquiring unit, for analyzing pretreated seismic profile bitmap, obtain seismic profile bitmap
Two dimension corresponding to two-dimentional geodetic coordinates value corresponding to time started value, seismic profile starting endpoint and seismic profile distal point is big
Ground coordinate value, utilize geodetic coordinates value corresponding to geodetic coordinates value corresponding to seismic profile starting endpoint, seismic profile distal point
Obtain horizontal range coefficient, using seismic profile bitmap at the beginning of between value obtain vertical range coefficient;
Seismic profile bitmap vector unit, for explaining that information is carried out at vector quantization in pretreated seismic profile bitmap
Reason, and using the seismic profile bitmap at the beginning of between value determine horizontal reference on the seismic profile bitmap after vector quantization
Line;Wherein, the starting endpoint of the horizontal reference line is overlapping with seismic profile starting endpoint, the distal point of the horizontal reference line
It is overlapping with seismic profile distal point;
Information dot picking first module is explained, on the basis of the horizontal reference line, being determined perpendicular to the horizontal reference
The vertical curve of line, the starting endpoint of the horizontal reference line is obtained to the bitmap horizontal range of the vertical curve, the vertical curve
With explaining the intersection point of information to the bitmap vertical range of the horizontal reference line;
Information dot picking second unit is explained, the bitmap for the starting endpoint using the horizontal reference line to the vertical curve
Horizontal range, horizontal range coefficient, the vertical curve with explain information intersection point to the horizontal reference line bitmap it is vertical away from
From two-dimentional geodetic coordinates value corresponding to value between at the beginning of, seismic profile bitmap, vertical range coefficient, seismic profile starting endpoint
Obtain geodetic coordinates X values, geodetic coordinates Y value and the time Z value of intersection point of the vertical curve with explaining information;
Information pickup unit is explained, the earth for all vertical curves using identical layer position and tomography and the intersection point for explaining information
Coordinate X values, geodetic coordinates Y value and time Z value are merged into figure;
The proportionality coefficient acquiring unit includes:Actual range acquisition module, for using corresponding to seismic profile starting endpoint
Two-dimentional geodetic coordinates value corresponding to two-dimentional geodetic coordinates value and seismic profile distal point determines seismic profile starting endpoint and end
Actual range between end points;Bitmap is obtained between seismic profile starting endpoint and distal point apart from acquisition module for measuring
Bitmap distance in seismic profile bitmap;Horizontal range coefficient module, for the seismic profile starting endpoint and distal point
Between actual range divided by bitmap distance obtain horizontal range coefficient.
6. device as claimed in claim 5, it is characterised in that the explanation information dot picking second unit includes:
The two-dimentional geodetic coordinates value determining module of information point is explained, for the starting endpoint using the horizontal reference line described in
It is big in the bitmap horizontal range of vertical curve, horizontal range coefficient, two-dimentional geodetic coordinates value corresponding to seismic profile starting endpoint
Ground coordinate X, Y value obtain the geodetic coordinates X of intersection point of the vertical curve with explaining information, Y value;
Explain the time Z value determining module of information point, for using be worth between at the beginning of seismic profile bitmap, vertical range system
The bitmap vertical range of intersection point to the horizontal reference line of several, the described vertical curve with explaining information obtains the vertical curve and solution
Release the time Z value of the intersection point of information.
7. device as claimed in claim 5, it is characterised in that the proportionality coefficient acquiring unit includes:
Bitmap time gap acquisition module, for being worth corresponding points as starting point using the time started, it is corresponding to obtain fixed time period end
Bitmap distance between point and time started value corresponding points in seismic profile bitmap;
Vertical range coefficient module, for the fixed time period divided by fixed time period end corresponding points and time started
Bitmap distance between value corresponding points in seismic profile bitmap obtains vertical range coefficient.
8. device as claimed in claim 5, it is characterised in that the seismic profile bitmap vector unit includes:
Drafting module, for pretreated seismic cross-section to be imported in vector software, carry out redrawing ground using line of vector
Shake layer position and the tomography of section;
Staining module, layer position, tomography for different earthquake section make a distinction according to different colors.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510546827.7A CN105137485B (en) | 2015-08-31 | 2015-08-31 | The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510546827.7A CN105137485B (en) | 2015-08-31 | 2015-08-31 | The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105137485A CN105137485A (en) | 2015-12-09 |
CN105137485B true CN105137485B (en) | 2017-12-05 |
Family
ID=54722891
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510546827.7A Active CN105137485B (en) | 2015-08-31 | 2015-08-31 | The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137485B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0570615A1 (en) * | 1992-05-15 | 1993-11-24 | Prakla-Seismos GmbH | Method for improvement of seismic sections |
JP4065946B2 (en) * | 2003-03-20 | 2008-03-26 | 独立行政法人産業技術総合研究所 | Initial motion detection method for multi-component AE waveform |
CN102509511A (en) * | 2011-11-25 | 2012-06-20 | 武汉大学 | Rapid map making system based on GIS (Geographic Information System) data |
CN103376463A (en) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | Inversion modeling method based on fault control |
CN104570098A (en) * | 2015-01-28 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | Method and device for acquiring common-azimuth-angle seismic data |
-
2015
- 2015-08-31 CN CN201510546827.7A patent/CN105137485B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0570615A1 (en) * | 1992-05-15 | 1993-11-24 | Prakla-Seismos GmbH | Method for improvement of seismic sections |
JP4065946B2 (en) * | 2003-03-20 | 2008-03-26 | 独立行政法人産業技術総合研究所 | Initial motion detection method for multi-component AE waveform |
CN102509511A (en) * | 2011-11-25 | 2012-06-20 | 武汉大学 | Rapid map making system based on GIS (Geographic Information System) data |
CN103376463A (en) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | Inversion modeling method based on fault control |
CN104570098A (en) * | 2015-01-28 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | Method and device for acquiring common-azimuth-angle seismic data |
Non-Patent Citations (2)
Title |
---|
coreIdraw在地震地质制图中的应用;王志胜 等;《科技咨询》;20121231(第29期);第7-8页 * |
基于CorelDraw的地震信息数据处理;李超玲 等;《中原工学院院报》;20071231;第18卷(第6期);第72-75页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105137485A (en) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102759745B (en) | Method for forecasting carbonate reservoir based on forward modeling of digital geological outcrop model | |
Bisdom et al. | Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling | |
CN104134002A (en) | Clastic rock reservoir modeling method and device based on a digital geological outcrop | |
CN104007464B (en) | A kind of TSP big gun hole measures and powder charge integrated apparatus and method | |
CN109884710B (en) | Micro-logging tomography method aiming at excitation well depth design | |
CN109270567B (en) | Method for making exploration groove gamma-irradiation rate equivalent map in uranium ore exploration | |
CN102096944A (en) | Cloud recognition method of three-dimensional laser scanning point in geological body structural plane | |
CN102901989A (en) | Gravity field or magnetic field data based geologic body three-dimensional visualized modeling and interpretation method | |
CN102879799A (en) | Multi-direction seismic energy gradient difference carbonate karst cave type reservoir identification method | |
CN109872393B (en) | A kind of three-dimensional geological data processing method based on ground, subsurface geology information | |
CN105911592B (en) | A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint | |
CN110764142B (en) | Method for assisting seismic data interpretation | |
US20140095078A1 (en) | Method and system for presenting seismic information | |
CN101788686A (en) | Point-by-point excitation well depth designing method based on uniform near-surface model | |
US20130327125A1 (en) | Method for geochemical grandient exploration | |
CN103310061B (en) | Coal seismic prospecting drafting system | |
CN106407500B (en) | A kind of preparation method and device of vector quantization section | |
CN105137485B (en) | The automatic Picking method for transformation and device of information are explained in a kind of seismic profile bitmap | |
CN108242069B (en) | Oil reservoir profile compiling and drawing method based on seismic profile and time-depth relation | |
CN202257688U (en) | Three-dimensional geologic body visual modeling and explaining system based on gravity field or magnetic field data | |
Lageson et al. | Applications of Google Earth Pro to fracture and fault studies of Laramide anticlines in the Rocky Mountain foreland | |
CN107917692A (en) | A kind of block identification method based on unmanned plane | |
CN106504319A (en) | Reservoir Three Dimensional Contrast map generalization method and device between well | |
CN107918160B (en) | Comprehensive scanning system and method for abandoned well | |
YU et al. | Control function of base level for facies and sand body distribution, Putaohua oil layer of X13 area, Songliao basin |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |