CN105137485A - Automatic extraction and conversion method and device for interpretation information in seismic section bitmap - Google Patents
Automatic extraction and conversion method and device for interpretation information in seismic section bitmap Download PDFInfo
- Publication number
- CN105137485A CN105137485A CN201510546827.7A CN201510546827A CN105137485A CN 105137485 A CN105137485 A CN 105137485A CN 201510546827 A CN201510546827 A CN 201510546827A CN 105137485 A CN105137485 A CN 105137485A
- Authority
- CN
- China
- Prior art keywords
- bitmap
- seismic section
- value
- terrestrial coordinate
- point
- 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
Abstract
The invention relates to an automatic extraction and conversion method and device for interpretation information in a seismic section bitmap. According to the technical scheme, the method comprises the steps: determining two-dimensional geodetic coordinates and starting time of the head and tail ends of the seismic section bitmap after the preprocessing of the seismic section bitmap; automatically calculating and reading actual three-dimensional geodetic coordinates corresponding to an interpretation information point through a small gap; enabling the same type of interpretation results to be discriminated and merged into one text file, thereby bringing convenience to a user for importing the text file into seismic software so as to make spatial contrastive analysis and generate a precise map, achieving the quick digitalization of the interpretation information in the seismic section bitmap, and improving the precision and utilization rate of the interpretation result data in the seismic section bitmap.
Description
Technical field
The present invention relates to exploration engineering field, particularly the automatic Picking method for transformation of explain information and device in a kind of seismic section bitmap.
Background technology
In petroleum and natural gas industry, hydrocarbon-bearing pool is enriched among the ad-hoc location on cover stratum, underground.The exploration of modern petroleum and rock gas, exploitation need to utilize seismic section to determine the position of formation at target locations, buried depth, form to a great extent, and calculate with reserves for well location design.Thus by seismic section, sub-surface structural feature, tomography space characteristics are significant definitely rapidly and accurately.
But it is in seismic interpretation process, due to the limitation of data, low for original seismic cross-section sheet utilization factor; Especially, in overseas in exploration, the seismic section often obtained is all much the form of picture, is difficult to be formed on main geologic layer bit space to explanation personnel be familiar with intuitively, have impact on the actual effect of seismic prospecting.Section picture prints by mode in the past exactly, is sampled by artificial read point, and this mode efficiency is low, and error is large, can not be combined into figure with the interpretation results of existing geological data; Another kind of mode is that whole seismic section is changed into geological data, i.e. SegY file, and this mode requires high to picture quality, and debugging is slow for operation time.
In order to seismic interpretation personnel are freed from hard work, improve accuracy and efficiency, seismic section achievement is carried out high-density sampling by existing method in the mapping softwares such as Coreldraw, obtain rapidly and accurately and adopt some XYZ coordinate value spatially, and the interpretation results of multiple section can be merged into a volume coordinate file, accurate geologic structure diagram can be generated.But the weak point of this method is: the distortion of original seismic cross-section sheet large area, time fuzzy, performance data may exist certain error.
Summary of the invention
For solving the problem of prior art, the present invention proposes automatic Picking method for transformation and the device of explain information in a kind of seismic section bitmap, improves degree of accuracy and the utilization factor of interpretation results data in seismic section bitmap.
For achieving the above object, the invention provides the automatic Picking method for transformation of explain information in a kind of seismic section bitmap, comprising:
Carry out pre-service to seismic section bitmap, make the section horizontal line of seismic section bitmap be in horizontality, and the time shaft of seismic section bitmap is perpendicular to ground level, the border motif of seismic section bitmap is rectangle;
Pretreated seismic section bitmap is analyzed, obtain two-dimentional terrestrial coordinate value corresponding to the start time value of seismic section bitmap, the initial end points of seismic section and two-dimentional terrestrial coordinate value corresponding to seismic section distal point, the two-dimentional terrestrial coordinate value that the two-dimentional terrestrial coordinate value utilizing the initial end points of seismic section corresponding, seismic section distal point are corresponding obtains horizontal range coefficient, utilizes the start time value of seismic section bitmap to obtain vertical range coefficient;
Vectorized process is carried out to explain information in pretreated seismic section bitmap, and utilizes on the seismic section bitmap of the start time value of described seismic section bitmap after vector quantization and determine horizontal reference line; Wherein, the initial end points of described horizontal reference line is overlapping with the initial end points of seismic section, and the distal point of described horizontal reference line is overlapping with seismic section distal point;
With described horizontal reference line for benchmark, determine the vertical curve perpendicular to described horizontal reference line, obtain the bitmap vertical range of intersection point to described horizontal reference line of the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, described vertical curve and explain information;
Utilize the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, horizontal range coefficient, the intersection point of described vertical curve and explain information is to the bitmap vertical range of described horizontal reference line, the start time value of seismic section bitmap, vertical range coefficient, the intersection point of the two-dimentional terrestrial coordinate value that the initial end points of seismic section is corresponding and described vertical curve and explain information obtains the terrestrial coordinate X value of the intersection point of described vertical curve and explain information to the bitmap vertical range of described horizontal reference line, terrestrial coordinate Y value and time Z value,
The terrestrial coordinate X value of all vertical curves of identical layer position and tomography and the intersection point of explain information, terrestrial coordinate Y value and time Z value is utilized to be merged into figure.
Preferably, the step of the terrestrial coordinate X value of the intersection point of the described vertical curve of described acquisition and explain information, terrestrial coordinate Y value and time Z value comprises:
The bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, horizontal range coefficient, the terrestrial coordinate X of the initial end points of seismic section, Y value is utilized to obtain terrestrial coordinate X, the Y value of the intersection point of described vertical curve and explain information;
The intersection point of the start time value of seismic section bitmap, vertical range coefficient, described vertical curve and explain information is utilized to obtain the time Z value of the intersection point of described vertical curve and explain information to the bitmap vertical range of described horizontal reference line.
Preferably, the step of described acquisition horizontal range coefficient comprises:
The two-dimentional terrestrial coordinate value utilizing the initial end points of seismic section corresponding and the actual range between the initial end points of two-dimentional terrestrial coordinate value determination seismic section corresponding to seismic section distal point and distal point;
Measure the bitmap distance obtained between the initial end points of seismic section and distal point in seismic section bitmap;
Actual range between the initial end points of described seismic section and distal point obtains horizontal range coefficient divided by bitmap distance.
Preferably, the step of described vertical range coefficient comprises:
With start time value corresponding point for starting point, acquisition set time section end corresponding point and start time are worth the bitmap distance between corresponding point in seismic section bitmap;
Described set time section obtains vertical range coefficient divided by the bitmap distance between described set time section end corresponding point and start time value corresponding point.
Preferably, described the step that explain information in pretreated seismic section bitmap carries out vectorized process to be comprised:
Pretreated seismic cross-section is imported in vector software, use line of vector to carry out the layer position and the tomography that redraw seismic section;
Layer position, the tomography of different earthquake section are distinguished according to different colors.
For achieving the above object, present invention also offers the automatic Picking reforming unit of explain information in a kind of seismic section bitmap, comprising:
Pretreatment unit, for carrying out pre-service to seismic section bitmap, make the section horizontal line of seismic section bitmap be in horizontality, and the time shaft of seismic section bitmap is perpendicular to ground level, the border motif of seismic section bitmap is rectangle;
Scale-up factor acquiring unit, for analyzing pretreated seismic section bitmap, obtain two-dimentional terrestrial coordinate value corresponding to the start time value of seismic section bitmap, the initial end points of seismic section and two-dimentional terrestrial coordinate value corresponding to seismic section distal point, utilize the terrestrial coordinate value that the initial end points of seismic section is corresponding, the terrestrial coordinate value acquisition horizontal range coefficient that seismic section distal point is corresponding, utilize the start time value of seismic section bitmap to obtain vertical range coefficient;
Seismic section bitmap vector unit, for carrying out vectorized process to explain information in pretreated seismic section bitmap, and utilizes on the seismic section bitmap of the start time value of described seismic section bitmap after vector quantization and determines horizontal reference line; Wherein, the initial end points of described horizontal reference line is overlapping with the initial end points of seismic section, and the distal point of described horizontal reference line is overlapping with seismic section distal point;
Explain information dot picking first module, for with described horizontal reference line for benchmark, determine the vertical curve perpendicular to described horizontal reference line, obtain the bitmap vertical range of intersection point to described horizontal reference line of the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, described vertical curve and explain information;
Explain information dot picking second unit, the intersection point for the bitmap horizontal range of the initial end points to described vertical curve that utilize described horizontal reference line, horizontal range coefficient, described vertical curve and explain information obtains the terrestrial coordinate X value of the intersection point of described vertical curve and explain information, terrestrial coordinate Y value and time Z value to the bitmap vertical range of described horizontal reference line, the start time value of seismic section bitmap, vertical range coefficient, two-dimentional terrestrial coordinate value that the initial end points of seismic section is corresponding;
Explain information pickup unit, is merged into figure for the terrestrial coordinate X value of the intersection point of all vertical curves and explain information that utilize identical layer position and tomography, terrestrial coordinate Y value and time Z value.
Preferably, described explain information dot picking second unit comprises:
The two-dimentional terrestrial coordinate value determination module of explain information point, obtains terrestrial coordinate X, the Y value of the intersection point of described vertical curve and explain information for the terrestrial coordinate X in the two-dimentional terrestrial coordinate value that the bitmap horizontal range of the initial end points to described vertical curve that utilize described horizontal reference line, horizontal range coefficient, the initial end points of seismic section are corresponding, Y value;
The time Z value determination module of explain information point, the bitmap vertical range for the intersection point to described horizontal reference line that utilize the start time value of seismic section bitmap, vertical range coefficient, described vertical curve and explain information obtains the time Z value of the intersection point of described vertical curve and explain information.
Preferably, described scale-up factor acquiring unit comprises:
Actual range acquisition module, for the two-dimentional terrestrial coordinate value that utilizes the initial end points of seismic section corresponding and the actual range between the initial end points of two-dimentional terrestrial coordinate value determination seismic section corresponding to seismic section distal point and distal point;
Bitmap distance acquisition module, for measuring the bitmap distance obtained between the initial end points of seismic section and distal point in seismic section bitmap;
Horizontal range coefficient module, obtains horizontal range coefficient for the actual range between the initial end points of described seismic section and distal point divided by bitmap distance.
Preferably, described scale-up factor acquiring unit comprises:
Bitmap time gap acquisition module, for being worth corresponding point for starting point with the start time, acquisition set time section end corresponding point and start time are worth the bitmap distance between corresponding point in seismic section bitmap;
Vertical range coefficient module, obtains vertical range coefficient for described set time section divided by the bitmap distance between described set time section end corresponding point and start time value corresponding point in seismic section bitmap.
Preferably, described seismic section bitmap vector unit comprises:
Drafting module, for importing in vector software by pretreated seismic cross-section, uses line of vector to carry out the layer position and the tomography that redraw seismic section;
Staining module, distinguishes according to different colors for the layer position of different earthquake section, tomography.
Technique scheme has following beneficial effect: the technical program is by after the pre-service of seismic section bitmap, determine two-dimentional terrestrial coordinate value and the start time value at section bitmap head and the tail two ends, and automatically calculated by Small Distance and read actual three-dimensional terrestrial coordinate corresponding to explain information point, and the differentiation of same class interpretation results is merged into a text, be convenient to user and import to space comparative analysis and the accurate map of generation in seismic software, realize fast digitizing by explain information in seismic section bitmap, thus improve degree of accuracy and the utilization factor of interpretation results data in seismic section bitmap.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the automatic Picking method for transformation process flow diagram of explain information in a kind of seismic section bitmap of proposing of the present invention;
Fig. 2 is the schematic diagram before the seismic section bitmap pre-service of the present embodiment;
Fig. 3 is for the pretreated schematic diagram of the bitmap shown in Fig. 2;
Fig. 4 is that the terrestrial coordinate X value of explain information point in the present embodiment, terrestrial coordinate Y value and time Z value calculate schematic diagram;
Fig. 5 is the present embodiment seismic section line planimetric position distribution plan;
Fig. 6 is the pickup schematic diagram of the present embodiment explain information point;
Fig. 7 is the automatic Picking reforming unit block diagram of explain information in a kind of seismic section bitmap of proposing of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
As shown in Figure 1, be the automatic Picking method for transformation process flow diagram of explain information in a kind of seismic section bitmap of the present invention's proposition.Comprise:
Step 101): carry out pre-service to seismic section bitmap, make the section horizontal line of seismic section bitmap be in horizontality, and the time shaft of seismic section bitmap is perpendicular to ground level, the border motif of seismic section bitmap is rectangle;
In a step 101, utilize the software process seismic section bitmap pictures such as Photoshop, try to achieve real angle between horizontal line and section horizontal line in section.As shown in Figure 2, be the schematic diagram before the seismic section bitmap pre-service of the present embodiment.By rotating, section horizontal line is become horizontality; As picture exists certain distortion, also will reduce, make cross sectional boundary form remain rectangle, and the time shaft of seismic section bitmap be perpendicular to ground level, the border motif of seismic section bitmap is rectangle.As shown in Figure 3, for the pretreated schematic diagram of the bitmap shown in Fig. 2.
Step 102): pretreated seismic section bitmap is analyzed, obtain two-dimentional terrestrial coordinate value corresponding to the start time value of seismic section bitmap, the initial end points of seismic section and two-dimentional terrestrial coordinate value corresponding to seismic section distal point, the two-dimentional terrestrial coordinate value that the two-dimentional terrestrial coordinate value utilizing the initial end points of seismic section corresponding, seismic section distal point are corresponding obtains horizontal range coefficient, utilizes the start time value of seismic section bitmap to obtain vertical range coefficient;
As shown in Figure 5, be the present embodiment seismic section line planimetric position distribution plan.In a step 102, on the basis of step 101, pretreated seismic section bitmap and related data (as seismic line distribution plan) thereof are analyzed, obtains two-dimentional terrestrial coordinate value corresponding to seismic section head and the tail two-end-point and section start time value thereof; Wherein, two-dimentional terrestrial coordinate value comprises terrestrial coordinate X value and terrestrial coordinate Y value.
On in horizontal in earthquake three-dimensional coordinate system, by terrestrial coordinate X value, the actual range of terrestrial coordinate Y value acquisition from beginning to end between two-end-point that seismic section head and the tail two-end-point is corresponding, record the bitmap distance of head and the tail two-end-point in bitmap by assistant software, actual range is obtained the horizontal range coefficient of this seismic section bitmap divided by bitmap distance; In vertical direction in earthquake three-dimensional coordinate system, appoint and get 1000ms as Temporal-distance value, namely the time value of starting point is 200ms, the time value of distal point is 1200ms, and the bitmap distance measured between this starting point and distal point on seismic section bitmap, above-mentioned 1000ms is obtained vertical range coefficient divided by bitmap distance value.
Step 103): vectorized process is carried out to explain information in pretreated seismic section bitmap, and utilizes on the seismic section bitmap of the start time value of described seismic section bitmap after vector quantization and determine horizontal reference line; Wherein, the initial end points of described horizontal reference line is overlapping with the initial end points of seismic section, and the distal point of described horizontal reference line is overlapping with seismic section distal point;
In step 103, imported in vector software by seismic section bitmap pretreated in step 101, in the present embodiment, vector software is Coreldraw.Carry out by line of vector the layer position and the tomography that redraw seismic section bitmap, different layers position or tomography are distinguished according to color.In addition, draw one obtains start time value horizontal reference line through step 102, the head and the tail two-end-point of this horizontal reference line is strictly corresponding overlapping respectively with the seismic section head and the tail two-end-point that step 102 confirms.
Step 104): with described horizontal reference line for benchmark, determine the vertical curve perpendicular to described horizontal reference line, obtain the bitmap vertical range of intersection point to described horizontal reference line of the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, described vertical curve and explain information;
As shown in Figure 6, be the pickup schematic diagram of the present embodiment explain information point.In figure 6, a represents the bitmap horizontal range of initial end points to vertical curve of horizontal reference line, and b represents the bitmap vertical range of the intersection point of vertical curve and explain information to horizontal reference line.
Step 105): the bitmap horizontal range of initial end points to described vertical curve utilizing described horizontal reference line, horizontal range coefficient, the intersection point of described vertical curve and explain information is to the bitmap vertical range of described horizontal reference line, the start time value of seismic section bitmap, vertical range coefficient, the intersection point of the two-dimentional terrestrial coordinate value that the initial end points of seismic section is corresponding and described vertical curve and explain information obtains the terrestrial coordinate X value of the intersection point of described vertical curve and explain information to the bitmap vertical range of described horizontal reference line, terrestrial coordinate Y value and time Z value,
As shown in Figure 4, for the terrestrial coordinate X value of explain information point in the present embodiment, terrestrial coordinate Y value and time Z value calculate schematic diagram.As shown in Figure 4, for the terrestrial coordinate X value of explain information point in the present embodiment, terrestrial coordinate Y value and time Z value calculate schematic diagram.Composition graphs 6, the earth XY coordinate that certain point is actual and time Z value are calculated as follows formula 1, formula 2, formula 3.
Wherein,
the terrestrial coordinate X value of the impact point that expression will calculate,
represent the terrestrial coordinate X value of first point on seismic section,
represent the terrestrial coordinate X value of tail point on seismic section,
represent in tail point X-coordinate value corresponding to vector software Bitmap,
represent in first point X-coordinate value corresponding to vector software Bitmap,
represent the X-coordinate value that interpretation results is corresponding in vector software.
Wherein,
the terrestrial coordinate Y value of the impact point that expression will calculate,
represent the terrestrial coordinate Y value of first point on seismic section,
represent the terrestrial coordinate Y value of tail point on seismic section,
represent in tail point Y-coordinate value corresponding to vector software Bitmap,
represent in first point Y-coordinate value corresponding to vector software Bitmap,
represent the Y-coordinate value that interpretation results is corresponding in vector software.
For the calculating of Z value, adopt formula 3:
Wherein,
the actual time value of the impact point that expression will calculate,
represent seismic section start time value,
represent the time value that seismic section the preceding paragraph set time length ends is corresponding, (
) represent at bitmap distance value corresponding to one period of set time that vector software Bitmap is corresponding,
represent the bitmap distance value of the intersection point in vector software between vertical curve and interpretation results to horizontal reference line.
Vertical curve is constantly moved to section tail point direction, often move once, obtain the intersection point that a perpendicular line is corresponding with interpretation results, obtain the terrestrial coordinate X value of each intersection point, terrestrial coordinate Y value and time Z value in the manner described above, and be output into text with certain text formatting (the layer bit format as seismic interpretation software Landmark).
Step 106): utilize the terrestrial coordinate X value of all vertical curves of identical layer position and tomography and the intersection point of explain information, terrestrial coordinate Y value and time Z value to be merged into figure.
Accordingly, as shown in Figure 7, be the automatic Picking reforming unit block diagram of explain information in a kind of seismic section bitmap of the present invention's proposition.Comprise:
Pretreatment unit 701, for carrying out pre-service to seismic section bitmap, make the section horizontal line of seismic section bitmap be in horizontality, and the time shaft of seismic section bitmap is perpendicular to ground level, the border motif of seismic section bitmap is rectangle;
Scale-up factor acquiring unit 702, for analyzing pretreated seismic section bitmap, obtain two-dimentional terrestrial coordinate value corresponding to the start time value of seismic section bitmap, the initial end points of seismic section and two-dimentional terrestrial coordinate value corresponding to seismic section distal point, utilize the terrestrial coordinate value that the initial end points of seismic section is corresponding, the terrestrial coordinate value acquisition horizontal range coefficient that seismic section distal point is corresponding, utilize the start time value of seismic section bitmap to obtain vertical range coefficient;
Seismic section bitmap vector unit 703, for carrying out vectorized process to explain information in pretreated seismic section bitmap, and utilizes on the seismic section bitmap of the start time value of described seismic section bitmap after vector quantization and determines horizontal reference line; Wherein, the initial end points of described horizontal reference line is overlapping with the initial end points of seismic section, and the distal point of described horizontal reference line is overlapping with seismic section distal point;
Explain information dot picking first module 704, for with described horizontal reference line for benchmark, determine the vertical curve perpendicular to described horizontal reference line, obtain the bitmap vertical range of intersection point to described horizontal reference line of the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, described vertical curve and explain information;
Explain information dot picking second unit 705, the intersection point for the bitmap horizontal range of the initial end points to described vertical curve that utilize described horizontal reference line, horizontal range coefficient, described vertical curve and explain information obtains the terrestrial coordinate X value of the intersection point of described vertical curve and explain information, terrestrial coordinate Y value and time Z value to the bitmap vertical range of described horizontal reference line, the start time value of seismic section bitmap, vertical range coefficient, two-dimentional terrestrial coordinate value that the initial end points of seismic section is corresponding;
Explain information pickup unit 706, is merged into figure for the terrestrial coordinate X value of the intersection point of all vertical curves and explain information that utilize identical layer position and tomography, terrestrial coordinate Y value and time Z value.
In the present embodiment, start time value is 200ms, scale-up factor 9.8814, amounts to calculation 111 intersection points, and some distance is 2.5mm, only several seconds consuming time.Partial results is as table 1, and the text formatting of generation is Landmark interpretation horizon acquiescence introduced layer bit data form.In Table 1, data rows is from left to right respectively terrestrial coordinate X value, terrestrial coordinate Y value, time Z value.
Table 1
Seismic section picture directly utilizes by the technical program, does not need to print, and just can read X value, Y value, Z value that interpretation results point is corresponding fast.In addition, the technical program can read X value, Y value, the Z value in corresponding space in a large number exactly, and classifies to different interpretation results, generates different texts.Further, the interpretation results of different section picture is merged, generate the text of certain interpretation results on a large scale, and can be input in Landmark interpretation software and to be independently unified into figure to relevant seismic data volume.
Above-described embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only the specific embodiment of the present invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. the automatic Picking method for transformation of explain information in seismic section bitmap, is characterized in that, comprising:
Carry out pre-service to seismic section bitmap, make the section horizontal line of seismic section bitmap be in horizontality, and the time shaft of seismic section bitmap is perpendicular to ground level, the border motif of seismic section bitmap is rectangle;
Pretreated seismic section bitmap is analyzed, obtain two-dimentional terrestrial coordinate value corresponding to the start time value of seismic section bitmap, the initial end points of seismic section and two-dimentional terrestrial coordinate value corresponding to seismic section distal point, the two-dimentional terrestrial coordinate value that the two-dimentional terrestrial coordinate value utilizing the initial end points of seismic section corresponding, seismic section distal point are corresponding obtains horizontal range coefficient, utilizes the start time value of seismic section bitmap to obtain vertical range coefficient;
Vectorized process is carried out to explain information in pretreated seismic section bitmap, and utilizes on the seismic section bitmap of the start time value of described seismic section bitmap after vector quantization and determine horizontal reference line; Wherein, the initial end points of described horizontal reference line is overlapping with the initial end points of seismic section, and the distal point of described horizontal reference line is overlapping with seismic section distal point;
With described horizontal reference line for benchmark, determine the vertical curve perpendicular to described horizontal reference line, obtain the bitmap vertical range of intersection point to described horizontal reference line of the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, described vertical curve and explain information;
Utilize the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, horizontal range coefficient, the intersection point of described vertical curve and explain information is to the bitmap vertical range of described horizontal reference line, the start time value of seismic section bitmap, vertical range coefficient, the intersection point of the two-dimentional terrestrial coordinate value that the initial end points of seismic section is corresponding and described vertical curve and explain information obtains the terrestrial coordinate X value of the intersection point of described vertical curve and explain information to the bitmap vertical range of described horizontal reference line, terrestrial coordinate Y value and time Z value,
The terrestrial coordinate X value of all vertical curves of identical layer position and tomography and the intersection point of explain information, terrestrial coordinate Y value and time Z value is utilized to be merged into figure.
2. the method for claim 1, is characterized in that, the step of the terrestrial coordinate X value of the intersection point of the described vertical curve of described acquisition and explain information, terrestrial coordinate Y value and time Z value comprises:
The bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, horizontal range coefficient, the terrestrial coordinate X of the initial end points of seismic section, Y value is utilized to obtain terrestrial coordinate X, the Y value of the intersection point of described vertical curve and explain information;
The intersection point of the start time value of seismic section bitmap, vertical range coefficient, described vertical curve and explain information is utilized to obtain the time Z value of the intersection point of described vertical curve and explain information to the bitmap vertical range of described horizontal reference line.
3. method as claimed in claim 1 or 2, it is characterized in that, the step of described acquisition horizontal range coefficient comprises:
The two-dimentional terrestrial coordinate value utilizing the initial end points of seismic section corresponding and the actual range between the initial end points of two-dimentional terrestrial coordinate value determination seismic section corresponding to seismic section distal point and distal point;
Measure the bitmap distance obtained between the initial end points of seismic section and distal point in seismic section bitmap;
Actual range between the initial end points of described seismic section and distal point obtains horizontal range coefficient divided by bitmap distance.
4. method as claimed in claim 1 or 2, it is characterized in that, the step of described vertical range coefficient comprises:
With start time value corresponding point for starting point, acquisition set time section end corresponding point and start time are worth the bitmap distance between corresponding point in seismic section bitmap;
Described set time section obtains vertical range coefficient divided by the bitmap distance between described set time section end corresponding point and start time value corresponding point.
5. the method for claim 1, is characterized in that, describedly comprises the step that explain information in pretreated seismic section bitmap carries out vectorized process:
Pretreated seismic cross-section is imported in vector software, use line of vector to carry out the layer position and the tomography that redraw seismic section;
Layer position, the tomography of different earthquake section are distinguished according to different colors.
6. the automatic Picking reforming unit of explain information in seismic section bitmap, is characterized in that, comprising:
Pretreatment unit, for carrying out pre-service to seismic section bitmap, make the section horizontal line of seismic section bitmap be in horizontality, and the time shaft of seismic section bitmap is perpendicular to ground level, the border motif of seismic section bitmap is rectangle;
Scale-up factor acquiring unit, for analyzing pretreated seismic section bitmap, obtain two-dimentional terrestrial coordinate value corresponding to the start time value of seismic section bitmap, the initial end points of seismic section and two-dimentional terrestrial coordinate value corresponding to seismic section distal point, utilize the terrestrial coordinate value that the initial end points of seismic section is corresponding, the terrestrial coordinate value acquisition horizontal range coefficient that seismic section distal point is corresponding, utilize the start time value of seismic section bitmap to obtain vertical range coefficient;
Seismic section bitmap vector unit, for carrying out vectorized process to explain information in pretreated seismic section bitmap, and utilizes on the seismic section bitmap of the start time value of described seismic section bitmap after vector quantization and determines horizontal reference line; Wherein, the initial end points of described horizontal reference line is overlapping with the initial end points of seismic section, and the distal point of described horizontal reference line is overlapping with seismic section distal point;
Explain information dot picking first module, for with described horizontal reference line for benchmark, determine the vertical curve perpendicular to described horizontal reference line, obtain the bitmap vertical range of intersection point to described horizontal reference line of the bitmap horizontal range of initial end points to described vertical curve of described horizontal reference line, described vertical curve and explain information;
Explain information dot picking second unit, the intersection point for the bitmap horizontal range of the initial end points to described vertical curve that utilize described horizontal reference line, horizontal range coefficient, described vertical curve and explain information obtains the terrestrial coordinate X value of the intersection point of described vertical curve and explain information, terrestrial coordinate Y value and time Z value to the bitmap vertical range of described horizontal reference line, the start time value of seismic section bitmap, vertical range coefficient, two-dimentional terrestrial coordinate value that the initial end points of seismic section is corresponding;
Explain information pickup unit, is merged into figure for the terrestrial coordinate X value of the intersection point of all vertical curves and explain information that utilize identical layer position and tomography, terrestrial coordinate Y value and time Z value.
7. device as claimed in claim 6, it is characterized in that, described explain information dot picking second unit comprises:
The two-dimentional terrestrial coordinate value determination module of explain information point, obtains terrestrial coordinate X, the Y value of the intersection point of described vertical curve and explain information for the terrestrial coordinate X in the two-dimentional terrestrial coordinate value that the bitmap horizontal range of the initial end points to described vertical curve that utilize described horizontal reference line, horizontal range coefficient, the initial end points of seismic section are corresponding, Y value;
The time Z value determination module of explain information point, the bitmap vertical range for the intersection point to described horizontal reference line that utilize the start time value of seismic section bitmap, vertical range coefficient, described vertical curve and explain information obtains the time Z value of the intersection point of described vertical curve and explain information.
8. device as claimed in claim 6, it is characterized in that, described scale-up factor acquiring unit comprises:
Actual range acquisition module, for the two-dimentional terrestrial coordinate value that utilizes the initial end points of seismic section corresponding and the actual range between the initial end points of two-dimentional terrestrial coordinate value determination seismic section corresponding to seismic section distal point and distal point;
Bitmap distance acquisition module, for measuring the bitmap distance obtained between the initial end points of seismic section and distal point in seismic section bitmap;
Horizontal range coefficient module, obtains horizontal range coefficient for the actual range between the initial end points of described seismic section and distal point divided by bitmap distance.
9. device as claimed in claim 6, it is characterized in that, described scale-up factor acquiring unit comprises:
Bitmap time gap acquisition module, for being worth corresponding point for starting point with the start time, acquisition set time section end corresponding point and start time are worth the bitmap distance between corresponding point in seismic section bitmap;
Vertical range coefficient module, obtains vertical range coefficient for described set time section divided by the bitmap distance between described set time section end corresponding point and start time value corresponding point in seismic section bitmap.
10. device as claimed in claim 6, it is characterized in that, described seismic section bitmap vector unit comprises:
Drafting module, for importing in vector software by pretreated seismic cross-section, uses line of vector to carry out the layer position and the tomography that redraw seismic section;
Staining module, distinguishes according to different colors for the layer position of different earthquake section, tomography.
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 true CN105137485A (en) | 2015-12-09 |
| CN105137485B 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 |
|---|
| 李超玲 等: "基于CorelDraw的地震信息数据处理", 《中原工学院院报》 * |
| 王志胜 等: "coreIdraw在地震地质制图中的应用", 《科技咨询》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105137485B (en) | 2017-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bisdom et al. | Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling | |
| US7986319B2 (en) | Method and system for dynamic, three-dimensional geological interpretation and modeling | |
| CN106779417A (en) | The collection of engineering investigation information digitalization, management and integrated application method | |
| CN103810333B (en) | Engineering geology informatization work system | |
| CN104008252A (en) | Method for automatically generating plane-section drawing of underground cable duct lines | |
| Fey et al. | Deriving 3D displacement vectors from multi-temporal airborne laser scanning data for landslide activity analyses | |
| CN103699634A (en) | Database creation method and database creation system for one-stop pipeline detection | |
| CN104134002A (en) | Clastic rock reservoir modeling method and device based on a digital geological outcrop | |
| CN109270567B (en) | Method for making exploration groove gamma-irradiation rate equivalent map in uranium ore exploration | |
| CN103207419A (en) | Three-dimensional measurement method for tunnel rock formation attitude | |
| Bertotti et al. | Toward a quantitative definition of mechanical units: New techniques and results from an outcropping deep-water turbidite succession (Tanqua-Karoo Basin, South Africa) | |
| CN101788686B (en) | Point-by-point excitation well depth designing method based on uniform near-surface model | |
| CN114549774A (en) | Three-dimensional stratum modeling method based on drilling data | |
| CN105069162A (en) | Informatized exploratory adit catalog data acquisition system and method | |
| CN109872393B (en) | A kind of three-dimensional geological data processing method based on ground, subsurface geology information | |
| CN111812730A (en) | Resistivity data fusion three-dimensional imaging method and system for landslide detection | |
| CN103310061B (en) | Coal seismic prospecting drafting system | |
| CN104280764A (en) | Seismic exploration auxiliary data quality control method based on three-dimensional GIS | |
| CN108121853A (en) | A kind of System and method for that excavated-in t4 quantities is calculated based on AutoCAD | |
| CN202257688U (en) | Three-dimensional geologic body visual modeling and explaining system based on gravity field or magnetic field data | |
| CN105137485A (en) | Automatic extraction and conversion method and device for interpretation information in seismic section bitmap | |
| Lageson et al. | Applications of Google Earth Pro to fracture and fault studies of Laramide anticlines in the Rocky Mountain foreland | |
| Rana et al. | Identification of lineaments using Google tools | |
| CN106504319A (en) | Reservoir Three Dimensional Contrast map generalization method and device between well | |
| Zhang et al. | Semi-stochastic generation of rock discontinuity networks based on traces exposed on cavern roof |
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 |