CN109444953B - Automatic correction method and device for fault polygon drawing graph - Google Patents
Automatic correction method and device for fault polygon drawing graph Download PDFInfo
- Publication number
- CN109444953B CN109444953B CN201811234279.4A CN201811234279A CN109444953B CN 109444953 B CN109444953 B CN 109444953B CN 201811234279 A CN201811234279 A CN 201811234279A CN 109444953 B CN109444953 B CN 109444953B
- Authority
- CN
- China
- Prior art keywords
- sample point
- symmetric
- points
- distance
- pair
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000012937 correction Methods 0.000 title claims description 12
- 238000005070 sampling Methods 0.000 claims abstract description 97
- 238000013507 mapping Methods 0.000 claims abstract description 31
- 238000009499 grossing Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims description 21
- 238000009877 rendering Methods 0.000 claims description 15
- 238000012216 screening Methods 0.000 claims description 12
- 238000012952 Resampling Methods 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 8
- 239000011295 pitch Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000003247 decreasing effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/59—Other corrections
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Generation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a method and a device for automatically correcting a fault polygon drawing picture, wherein the method comprises the following steps: acquiring an original fault polygon drawing chart automatically identified and drawn by a machine; respectively acquiring coordinate data of all the symmetrical sample point pairs; according to the size relationship between the sample point ordinate data of the upper disc in the symmetrical sample point pair and the sample point ordinate data of the lower disc corresponding to the sample point ordinate data, cross points are screened and removed, and a first fault polygon drawing is obtained; determining the distance between each pair of symmetrical sampling points in the first fault polygon mapping chart, and acquiring the maximum distance, the first sampling points and the second sampling points; after the reference distance is reduced according to a second preset rule, the distance between the symmetrical sample point pairs except the central symmetrical sample point pair is used as the distance between the symmetrical sample point pairs except the central symmetrical sample point pair, the longitudinal coordinates of the symmetrical sample point pairs except the central symmetrical sample point pair are adjusted, and a second fault polygon drawing is obtained; and encrypting sampling points in the second fault polygon drawing and smoothing.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a method and a device for automatically correcting a fault polygon drawing.
Background
In the oil and mineral exploration work in the geological industry at present, the establishment of a tectonic graph is an important work content, and the method mainly utilizes the kinematic information of seismic wave variation such as reflection time, phasic, wave velocity and the like to research the distribution range and fluctuation form of a stratum interface and the fault development condition, converts the travel time in a seismic event section into the depth of the stratum interface, draws the geological tectonic graph and provides data for searching a tectonic oil and gas reservoir. Among them, planar drawing of fault polygons (polygons) is an important component, and the drawing precision directly affects the drawing precision of the structural diagram, and the precision of the structural diagram directly affects the accurate identification of traps and reservoirs, so that it is very important to accurately draw fault polygons. The current drawing methods mainly comprise the following two methods: one is manually finished, and manual drawing is performed according to the opening positions of faults and horizons, which is common in practical application; the second is to automatically connect by using the fracture of the horizon at each fault. The first mode adopts manual drawing, although the accuracy is high, the time and the labor are wasted, especially in an exploration block with a large area and multiple faults, if fault polygons is edited manually, each drawing needs a large amount of mechanical and repeated work, the time and the labor are wasted, and the work progress is greatly influenced; the second method, although improving the speed of drawing by automatic machine recognition, has the following problems:
1. the upper and lower disks of the Polygons have cross points;
the horizontal offset of the Polygons is not balanced enough and does not accord with geological rules;
2. polygons are not smooth and have more mutation points;
3. polygons are not closed at one end, and the closed end is not a sharp trend, but is rectangular.
Disclosure of Invention
The invention aims to provide a method and a device for automatically correcting a fault polygon drawing map, which are used for solving the problem of low precision of automatically drawing a geological structure map by using the existing machine, thereby improving the working efficiency and improving the precision of the geological structure map.
In order to achieve the above object, the technical solution of the present invention provides an automatic correction method for a tomographic polygon rendering map, including:
acquiring an original fault polygon drawing chart automatically identified and drawn by a machine;
respectively acquiring coordinate data of all symmetric sample point pairs from an original fault polygon drawing picture, wherein the symmetric sample point pairs are symmetric relations between sample points of an upper disc and sample points of a lower disc;
screening cross points in the original fault polygon drawing according to the size relationship between the sample point longitudinal coordinate data of the upper disc in the symmetrical sample point pair and the sample point longitudinal coordinate data of the lower disc symmetrical to the sample points of the upper disc;
removing the intersection points according to a first preset rule to obtain a first fault polygon drawing;
determining the distance between each pair of symmetrical sampling points in the first fault polygon mapping chart, and acquiring the maximum distance, and a first sampling point on the upper disc and a second sampling point on the lower disc corresponding to the maximum distance;
respectively taking the first sample point and the second sample point as central symmetric sample point pairs, taking the maximum distance as a reference distance, reducing the reference distance according to a second preset rule, and then adjusting the vertical coordinates of the symmetric sample point pairs except the central symmetric sample point pair as the distance between the symmetric sample point pairs except the central symmetric sample point pair to obtain a second fault polygon drawing, wherein the distance far away from the central symmetric sample point pair is gradually reduced until the distance is zero;
encrypting the number of sampling points in the second fault polygon mapping chart by adopting a resampling method;
and smoothing the encrypted sampling points to obtain a final fault polygon drawing.
The invention has the following advantages: according to the method, the cross point can be determined according to the size relationship between the vertical coordinate data of the sample point of the upper tray point and the vertical coordinate data of the sample point of the lower tray point which is symmetrical with the vertical coordinate data of the upper tray point in the symmetrical sample point pairs, after the cross point is removed according to a first preset rule, the vertical coordinate data between different symmetrical sample point pairs are adjusted according to the distance relationship between the symmetrical sample point pairs, so that the horizontal distance of a fault polygon is changed more uniformly, the two ends can be extinguished naturally, and the geological rule is met. Finally, sampling points in the second fault polygon drawing are encrypted by adopting a resampling method and are subjected to smoothing processing, so that the fault can be smooth and natural.
In order to achieve the above object, an embodiment of the present invention provides an automatic correction device for a tomographic polygon rendering map, including:
an acquisition unit configured to acquire an original tomographic polygon drawing chart drawn by automatic machine recognition;
respectively acquiring coordinate data of all symmetric sample point pairs from an original fault polygon drawing picture, wherein the symmetric sample point pairs are symmetric relations between sample points of an upper disc and sample points of a lower disc;
the screening unit is used for screening the cross points in the original fault polygon mapping chart according to the size relationship between the sample point longitudinal coordinate data of the upper disc in the symmetrical sample point pair and the sample point longitudinal coordinate data of the lower disc symmetrical to the sample points of the upper disc;
the processing unit is used for removing the intersection points according to a first preset rule to obtain a first fault polygon drawing;
determining the distance between each pair of symmetrical sampling points in the first fault polygon mapping chart, and acquiring the maximum distance, and a first sampling point on the upper disc and a second sampling point on the lower disc corresponding to the maximum distance;
respectively taking the first sample point and the second sample point as central symmetric sample point pairs, taking the maximum distance as a reference distance, reducing the reference distance according to a second preset rule, and then adjusting the vertical coordinates of the symmetric sample point pairs except the central symmetric sample point pair as the distance between the symmetric sample point pairs except the central symmetric sample point pair to obtain a second fault polygon drawing, wherein the distance far away from the central symmetric sample point pair is gradually reduced until the distance is zero;
the sampling unit is used for encrypting sampling points in the second fault polygon mapping chart by adopting a resampling method;
and the processing unit is also used for carrying out smoothing processing on the encrypted sampling points to obtain a final fault polygon drawing.
The invention has the following advantages: according to the method, the cross point can be determined according to the size relationship between the vertical coordinate data of the sample point of the upper tray point and the vertical coordinate data of the sample point of the lower tray point which is symmetrical with the vertical coordinate data of the upper tray point in the symmetrical sample point pairs, after the cross point is removed according to a first preset rule, the vertical coordinate data between different symmetrical sample point pairs are adjusted according to the distance relationship between the symmetrical sample point pairs, so that the horizontal distance of a fault polygon is changed more uniformly, the two ends can be extinguished naturally, and the geological rule is met. Finally, sampling points in the second fault polygon drawing are encrypted by adopting a resampling method and are subjected to smoothing processing, so that the fault can be smooth and natural.
Drawings
Fig. 1 is a schematic flow chart of an automatic correction method for a fault polygon rendering according to an embodiment of the present invention;
FIG. 2 is a schematic diagram comparing cross-points before and after removal;
FIG. 3 is a schematic diagram showing a comparison between horizontal fault distances of fault polygons before and after equalization processing;
FIG. 4 is a schematic diagram showing a comparison between before and after smoothing a fault;
fig. 5 is a schematic structural diagram of an automatic correction device for a tomographic polygon rendering according to an embodiment of the present invention.
Detailed Description
The following examples are intended to illustrate the invention but are not intended to limit the scope of the invention.
Example 1
The embodiment 1 of the invention provides a flow schematic diagram of an automatic correction method of a fault polygon rendering map. The method is mainly used for correcting the original fault polygon drawing which is drawn after automatic recognition by a traditional machine. The problems that the upper plate and the lower plate in the fault polygon drawing chart drawn by the traditional method have crossing points and uneven horizontal fault distance, do not accord with geological rules, the polygon is not smooth enough, the number of the mutation points is large, the two ends of the polygon are not killed and the like are solved.
The schematic flow diagram of the method is specifically shown in fig. 1, and the method may include:
and step 110, acquiring an original fault polygon drawing which is automatically identified and drawn by a machine.
And step 120, respectively acquiring coordinate data of all the symmetrical sample point pairs from the original fault polygon drawing graph.
Specifically, when the coordinate data in the original fault polygon drawing automatically identified and drawn by the machine is output, the sample point coordinate data of the upper disc is output first, and then the sample point coordinate data of the lower disc is output. And the sampling points of the upper disc and the sampling points of the lower disc have a symmetrical relation, so that the output sampling point coordinate data can be easily distinguished, wherein the first half part is the sampling point coordinate data of the upper disc, and the second half part is the sampling point coordinate data of the lower disc. And the first sampling point coordinate data of the upper disc and the first sampling point coordinate data of the lower disc are corresponding. By analogy, the coordinate data of all the symmetric sample point pairs can be easily known.
And step 130, screening the cross points in the original fault polygon drawing according to the size relationship between the sample point longitudinal coordinate data of the upper disc in the symmetrical sample point pair and the sample point longitudinal coordinate data of the lower disc symmetrical to the sample points of the upper disc.
Specifically, the ordinate data of the upper disk of the fault is generally larger than the ordinate data of the lower disk of the fault. Then, comparing the ordinate data of each pair of symmetric sample points, and determining that the upper and lower discs have intersection points at the positions which do not conform to the rule. That is, when it is determined that the sample point ordinate data of the upper disk in the second symmetric sample point pair in the original fault polygon mapping graph is less than or equal to the sample point ordinate data of the lower disk, it is determined that the second symmetric sample point pair is the intersection point. Wherein the second symmetric spline pair is any symmetric spline pair in all symmetric spline pairs.
And 140, removing the intersection points according to a first preset rule to obtain a first fault polygon drawing.
Specifically, in order to remove the intersection, the ordinate of the upper inventory may be adjusted, and the abscissa is not changed, so as to remove the intersection. And changing the vertical coordinate of the upper inventory is realized according to a first preset rule:
and adding a preset decimal to the ordinate data of the sample points of the lower tray of the second symmetric sample point pair to serve as a new ordinate of the upper tray of the second symmetric sample point pair, and forming a new sample point of the upper tray together with the abscissa of the upper tray of the second symmetric sample point pair.
And then, combining the sampling points of the new upper disc and the sampling points of the lower disc of the second symmetric sampling point pair into a new symmetric sampling point pair. Fig. 2 shows a comparative example before and after removing the intersection, as shown in fig. 2.
And the second symmetric sample point pair is any one of all symmetric sample point pairs in the original fault polygon rendering map. In fact, after all the intersections in the original tomographic polygon map are removed, the first tomographic polygon map, that is, a new tomographic polygon map, can be obtained. The first tomographic polygon map and the second tomographic polygon map to be described later are both new tomographic polygon maps obtained after a certain change, and are not the first and second tomographic polygons in the actual logical sense.
After removing the intersection, the following steps may be performed:
And 160, respectively taking the first sampling point and the second sampling point as center symmetric sampling point pairs, taking the maximum distance as a reference distance, reducing the reference distance according to a second preset rule, and then adjusting the vertical coordinates of the symmetric sampling point pairs except the center symmetric sampling point pair as the distance between the symmetric sampling point pairs except the center symmetric sampling point pair to obtain a second fault polygon drawing, wherein the distance far away from the center symmetric sampling point pair is gradually reduced until the distance is zero.
Specifically, based on many detailed observations (including observations of outcrops, deposits, seismic reflection data, etc.), the ideal pattern of "independent" individual fault pitches is that the fault pitch decreases in all directions towards the fault edge (end line) at a location near the center, and the fault pitch is zero at the end line location. Based on this feature, the spacing between each pair of symmetric sampling point pairs in the first tomographic polygon rendering map can be first determined, and then the maximum spacing is obtained. And the position of a pair of symmetric sample point pairs with the largest distance is taken as the central position of the fault, and the pair of symmetric sample point pairs can also be understood as central symmetric sample point pairs. For convenience of description, the sampling point of the upper disc in the centrosymmetric sampling point pair is defined as a first sampling point, and the sampling point of the lower disc in the centrosymmetric sampling point pair is defined as a second sampling point. The distance between the pairs of sample points is symmetrical about the center as the basic distance.
Because the distance between the symmetrical pairs of the sample points at the center in the fault is the largest, the distance between the symmetrical pairs of the sample points at the left and right sides is gradually reduced. Therefore, the coordinates of the symmetric pairs other than the central symmetric pair may be adjusted as the distance between the symmetric pairs other than the central symmetric pair after the basic distance is decreased according to the second preset rule. Thereby acquiring a second tomographic polygon map. Wherein, the distance between the symmetrical sample point pairs far away from the center is gradually reduced until the distance is zero.
Specifically, the maximum distance is taken as a reference distance, and after a value of n multiplied by i is reduced, the distance is taken as the distance between the nth symmetric sample point pair on the left side and the right side of the centrosymmetric sample point pair, n is a positive integer which is greater than or equal to 1 and less than or equal to one half of the total number of all the symmetric sample point pairs, the values of n are sequentially increased in a gradually increasing manner, and the increasing distance is 1; i is the ratio between the maximum pitch and n.
However, in order to ensure that the longitudinal center of each pair of symmetric sample points does not change, when adjusting the longitudinal coordinate of each pair of symmetric sample points, it is necessary to adjust the longitudinal coordinates of the original two pairs of symmetric sample points up and down, respectively, with reference to the centers of the original two pairs of symmetric sample points. For example, the maximum pitch is 10, n is 1-5, and i is 2. Then, when the vertical coordinate of each pair of symmetric sample point pairs is adjusted, the vertical coordinate of the sample point of the upper disc is a numerical value after the original vertical coordinate plus (the reference distance-nxi)/2; and the ordinate of the sample point of the lower plate is the value after the original ordinate- (the reference spacing-nxi)/2. Of course, the adjustment may be performed by sequentially decreasing the pitches before the symmetric sampling point pairs at equal pitches in the above-described manner, or may not be performed at equal pitches, and in short, the principle of gradually decreasing the pitches between the symmetric sampling point pairs on the left and right sides of the central symmetric sampling point pair is followed. Until the distance between the last symmetrical sampling point pairs on the left side and the right side is zero.
By the mode, the horizontal fault distances of fault polygons can be more balanced, the distance between the fault polygons is the largest from the center, and the distances between the fault polygons and the left side and the right side are gradually reduced until the fault polygons are extinguished, so that geological rules are met. Fig. 3 shows a comparison example before and after equalization processing is performed on the horizontal fault distance of the fault polygon, which is specifically shown in fig. 3.
And step 170, encrypting sampling points in the second fault polygon mapping chart by adopting a resampling method.
And 180, smoothing the encrypted sampling points to obtain a final fault polygon drawing.
Specifically, after the second tomographic polygon rendering is acquired, in order to further smooth the tomographic polygon rendering, the number of sampling points may be encrypted by using a resampling method. I.e. adding symmetric pairs of sample points. And then smoothing all the symmetric sample points, wherein a median filtering method or a mean filtering method can be adopted for smoothing optionally. To improve the smoothness of the polygon. Therefore, the curve in the finally obtained fault polygon drawing chart can be smoother and more natural. Specifically, as shown in fig. 4, fig. 4 is a front-back comparison example of smoothing processing of a fault. The curve in fig. 4 is significantly smoother.
According to the automatic correction method for the fault polygon drawing graph, the cross point can be determined according to the size relationship between the vertical coordinate data of the sample point of the upper tray point and the vertical coordinate data of the sample point of the lower tray point which is symmetrical to the upper tray point in the symmetrical sample point pairs, after the cross point is removed according to the first preset rule, the vertical coordinate data between different symmetrical sample point pairs are adjusted by using the distance relationship between the symmetrical sample point pairs, so that the horizontal distance of the fault polygon is changed more uniformly, and the two ends can be extinguished naturally, and accord with geological rules. Finally, sampling points in the second fault polygon drawing are encrypted by adopting a resampling method and are subjected to smoothing processing, so that the fault can be smooth and natural.
Example 2
Corresponding to embodiment 1, an embodiment of the present invention further provides an automatic correction apparatus for a tomographic polygon rendering map, specifically as shown in fig. 5, where the apparatus includes: an acquisition unit 501, a screening unit 502, a processing unit 503, and a sampling unit 504.
An acquisition unit 501 for acquiring an original tomographic polygon rendering map rendered by automatic machine recognition;
respectively acquiring coordinate data of all symmetric sample point pairs from an original fault polygon drawing picture, wherein the symmetric sample point pairs are symmetric relations between sample points of an upper disc and sample points of a lower disc;
a screening unit 502, configured to screen an intersection in the original fault polygon mapping chart according to a size relationship between sample point ordinate data of the upper tray in the symmetric sample point pair and sample point ordinate data of the lower tray symmetric to the sample point of the upper tray;
the processing unit 503 is configured to remove the intersection point according to a first preset rule, and obtain a first fault polygon rendering map;
determining the distance between each pair of symmetrical sampling points in the first fault polygon mapping chart, and acquiring the maximum distance, and a first sampling point on the upper disc and a second sampling point on the lower disc corresponding to the maximum distance;
respectively taking the first sample point and the second sample point as central symmetric sample point pairs, taking the maximum distance as a reference distance, reducing the reference distance according to a second preset rule, and then adjusting the vertical coordinates of the symmetric sample point pairs except the central symmetric sample point pair as the distance between the symmetric sample point pairs except the central symmetric sample point pair to obtain a second fault polygon drawing, wherein the distance far away from the central symmetric sample point pair is gradually reduced until the distance is zero;
a sampling unit 504, configured to encrypt sampling points in the second tomographic polygonal drawing by using a resampling method;
the processing unit 503 is further configured to perform smoothing processing on the encrypted sampling points to obtain a final tomographic polygonal mapping chart.
Optionally, the screening unit 502 is specifically configured to:
and when the vertical coordinate data of the sample points of the upper disc in the second symmetric sample point pair in the original fault polygon mapping chart is determined to be less than or equal to the vertical coordinate data of the sample points of the lower disc, determining the second symmetric sample point pair as an intersection point, wherein the second symmetric sample point pair is any one symmetric sample point pair in all the symmetric sample point pairs.
Optionally, the processing unit 503 is specifically configured to: adding a preset decimal to the ordinate data of the sample points of the lower tray of the second symmetric sample point pair, taking the ordinate data as the new ordinate of the upper tray of the second symmetric sample point pair, and forming a new sample point of the upper tray together with the abscissa of the upper tray of the second symmetric sample point pair;
and forming a new symmetric sample point pair by the sample point of the new upper disc and the sample point of the lower disc of the second symmetric sample point pair, and acquiring the first fault polygon drawing, wherein the second symmetric sample point pair is any one symmetric sample point pair in all symmetric sample point pairs in the original fault polygon drawing.
Optionally, the second preset rule is:
taking the maximum distance as a reference distance, reducing the value of n multiplied by i, and then taking the maximum distance as the distance between the nth symmetric sample point pair on the left side and the right side of the centrosymmetric sample point pair, wherein n is a positive integer which is more than or equal to 1 and less than or equal to one half of the total number of all the symmetric sample point pairs, the values of n are sequentially increased in an increasing way, and the increasing distance is 1; i is the ratio between the maximum pitch and n and is adjusted to ensure that each pair of symmetric samples is centered longitudinally unchanged.
Optionally, the processing unit 503 is specifically configured to: and smoothing the sampling points by using median filtering or mean filtering to obtain a final fault polygon mapping chart.
The functions performed by each component in the automatic correction device for a tomographic polygon rendering map according to an embodiment of the present invention are described in detail in embodiment 1 above, and therefore will not be described in detail here.
According to the automatic correction device for the fault polygon drawing graph, the cross point can be determined according to the size relationship between the sample point longitudinal coordinate data of the upper tray point and the sample point longitudinal coordinate data of the lower tray point which are symmetrical to the upper tray point in the symmetrical sample point pairs, after the cross point is removed according to the first preset rule, the longitudinal coordinate data between different symmetrical sample point pairs are adjusted according to the distance relationship between the symmetrical sample point pairs, so that the horizontal distance of the fault polygon is changed more uniformly, the two ends can be extinguished naturally, and the fault polygon drawing graph conforms to geological rules. Finally, sampling points in the second fault polygon drawing are encrypted by adopting a resampling method and are subjected to smoothing processing, so that the fault can be smooth and natural.
Example 3
Furthermore, an embodiment of the present invention also provides a computer-readable storage medium, on which computer program instructions are stored, which when executed by a processor implement the method steps of embodiment 1 described above.
Although the invention has been described in detail above with reference to a general description and specific examples, it will be apparent to one skilled in the art that modifications or improvements may be made thereto based on the invention. Accordingly, such modifications and improvements are intended to be within the scope of the invention as claimed.
Claims (4)
1. A method for automatically correcting a tomographic polygon rendering map, the method comprising:
acquiring an original fault polygon drawing chart automatically identified and drawn by a machine;
respectively acquiring coordinate data of all symmetric sample point pairs from the original fault polygon drawing graph, wherein the symmetric sample point pairs are symmetric relations between sample points of an upper disc and sample points of a lower disc;
screening the cross points in the original fault polygon mapping chart according to the size relationship between the sample point longitudinal coordinate data of the upper disc in the symmetrical sample point pair and the sample point longitudinal coordinate data of the lower disc symmetrical to the sample points of the upper disc;
removing the intersection points according to a first preset rule to obtain a first fault polygon drawing;
determining the distance between each pair of symmetrical sampling points in the first fault polygon mapping chart, and acquiring the maximum distance, and a first sampling point on the upper disc and a second sampling point on the lower disc corresponding to the maximum distance;
respectively taking the first sampling point and the second sampling point as central symmetrical sampling point pairs, taking the maximum distance as a reference distance, reducing the reference distance according to a second preset rule, and then adjusting the longitudinal coordinates of the symmetrical sampling point pairs except the central symmetrical sampling point pairs as the distance between the symmetrical sampling point pairs except the central symmetrical sampling point pairs to obtain a second fault polygon drawing, wherein the distance far away from the central symmetrical sampling point pairs is gradually reduced until the distance is zero;
encrypting sampling points in the second fault polygon mapping chart by adopting a resampling method;
smoothing the encrypted sampling points to obtain a final fault polygon drawing;
the screening of the intersection points in the original fault polygon mapping chart according to the size relationship between the sample point longitudinal coordinate data of the upper disc in the symmetrical sample point pair and the sample point longitudinal coordinate data of the lower disc symmetrical to the sample points of the upper disc specifically comprises the following steps:
when the vertical coordinate data of the upper disc sample points in a second symmetric sample point pair in the original fault polygon mapping chart is determined to be less than or equal to the vertical coordinate data of the lower disc sample points, determining the second symmetric sample point pair as an intersection point, wherein the second symmetric sample point pair is any one of all symmetric sample point pairs;
removing the intersection point according to a first preset rule, and acquiring a first fault polygon drawing, which specifically comprises the following steps:
adding a preset decimal to the ordinate data of the sample points of the lower tray of the second symmetric sample point pair, and using the ordinate data as a new ordinate of the upper tray of the second symmetric sample point pair to form a new sample point of the upper tray together with the abscissa of the upper tray of the second symmetric sample point pair;
forming a new symmetric sample point pair by the sample points of the new upper disc and the sample points of the lower disc of the second symmetric sample point pair, and acquiring a first fault polygon mapping chart, wherein the second symmetric sample point pair is any one symmetric sample point pair in all symmetric sample point pairs in the original fault polygon mapping chart;
the second preset rule is as follows:
taking the maximum distance as a reference distance, reducing the value of n multiplied by i to be used as the distance between the nth symmetric sample point pair on the left side and the right side of the centrosymmetric sample point pair, wherein n is a positive integer which is greater than or equal to 1 and less than or equal to one half of the total number of all the symmetric sample point pairs, the value of n is sequentially increased in a gradually increasing manner, and the increasing distance is 1; i is the ratio between the maximum pitch and n and is adjusted to ensure that each pair of symmetric samples is centered longitudinally unchanged.
2. The method according to claim 1, wherein the smoothing of the encrypted sampling points to obtain a final tomographic polygon mapping map specifically comprises:
and smoothing the sampling points by using median filtering or mean filtering to obtain a final fault polygon mapping chart.
3. An automatic correction device for a tomographic polygon rendering map, the device comprising:
an acquisition unit configured to acquire an original tomographic polygon drawing chart drawn by automatic machine recognition;
respectively acquiring coordinate data of all symmetric sample point pairs from the original fault polygon drawing picture, wherein the symmetric sample point pairs are symmetric relations between sample points of an upper disc and sample points of a lower disc;
the screening unit is used for screening the intersection points in the original fault polygon mapping chart according to the size relationship between the sampling point longitudinal coordinate data of the upper disc in the symmetrical sampling point pair and the sampling point longitudinal coordinate data of the lower disc symmetrical to the sampling points of the upper disc;
the processing unit is used for removing the intersection points according to a first preset rule to obtain a first fault polygon drawing;
determining the distance between each pair of symmetrical sampling points in the first fault polygon mapping chart, and acquiring the maximum distance, and a first sampling point on the upper disc and a second sampling point on the lower disc corresponding to the maximum distance;
respectively taking the first sampling point and the second sampling point as central symmetrical sampling point pairs, taking the maximum distance as a reference distance, reducing the reference distance according to a second preset rule, and then adjusting the longitudinal coordinates of the symmetrical sampling point pairs except the central symmetrical sampling point pairs as the distance between the symmetrical sampling point pairs except the central symmetrical sampling point pairs to obtain a second fault polygon drawing, wherein the distance far away from the central symmetrical sampling point pairs is gradually reduced until the distance is zero;
the sampling unit is used for encrypting sampling points in the second fault polygon mapping chart by adopting a resampling method;
the processing unit is further used for smoothing the encrypted sampling points to obtain a final fault polygon mapping chart;
the screening unit is specifically configured to:
when the vertical coordinate data of the upper disc sample points in a second symmetric sample point pair in the original fault polygon mapping chart is determined to be less than or equal to the vertical coordinate data of the lower disc sample points, determining the second symmetric sample point pair as an intersection point, wherein the second symmetric sample point pair is any one of all symmetric sample point pairs;
the processing unit is specifically configured to: adding a preset decimal to the ordinate data of the sample points of the lower tray of the second symmetric sample point pair, and using the ordinate data as a new ordinate of the upper tray of the second symmetric sample point pair to form a new sample point of the upper tray together with the abscissa of the upper tray of the second symmetric sample point pair;
forming a new symmetric sample point pair by the sample points of the new upper disc and the sample points of the lower disc of the second symmetric sample point pair, and acquiring a first fault polygon mapping chart, wherein the second symmetric sample point pair is any one symmetric sample point pair in all symmetric sample point pairs in the original fault polygon mapping chart;
the second preset rule is:
taking the maximum distance as a reference distance, reducing the value of n multiplied by i to be used as the distance between the nth symmetric sample point pair on the left side and the right side of the centrosymmetric sample point pair, wherein n is a positive integer which is greater than or equal to 1 and less than or equal to one half of the total number of all the symmetric sample point pairs, the value of n is sequentially increased in a gradually increasing manner, and the increasing distance is 1; i is the ratio between the maximum pitch and n and is adjusted to ensure that each pair of symmetric samples is centered longitudinally unchanged.
4. The apparatus according to claim 3, wherein the processing unit is specifically configured to: and smoothing the sampling points by using median filtering or mean filtering to obtain a final fault polygon mapping chart.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811234279.4A CN109444953B (en) | 2018-10-23 | 2018-10-23 | Automatic correction method and device for fault polygon drawing graph |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811234279.4A CN109444953B (en) | 2018-10-23 | 2018-10-23 | Automatic correction method and device for fault polygon drawing graph |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109444953A CN109444953A (en) | 2019-03-08 |
| CN109444953B true CN109444953B (en) | 2020-06-12 |
Family
ID=65548211
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811234279.4A Expired - Fee Related CN109444953B (en) | 2018-10-23 | 2018-10-23 | Automatic correction method and device for fault polygon drawing graph |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109444953B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
| CN106097446A (en) * | 2016-06-07 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of three dimensional geological curved surface method for reconstructing based on extensive geological data |
| WO2017213631A1 (en) * | 2016-06-07 | 2017-12-14 | Halliburton Energy Services, Inc. | Systems and methods for unfaulting point clouds |
-
2018
- 2018-10-23 CN CN201811234279.4A patent/CN109444953B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
| CN106097446A (en) * | 2016-06-07 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of three dimensional geological curved surface method for reconstructing based on extensive geological data |
| WO2017213631A1 (en) * | 2016-06-07 | 2017-12-14 | Halliburton Energy Services, Inc. | Systems and methods for unfaulting point clouds |
Non-Patent Citations (1)
| Title |
|---|
| 一种断层多边形自动平滑方法;李荟 等;《计算机系统应用》;20130831;第22卷(第8期);第140,159-161页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109444953A (en) | 2019-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2603979C1 (en) | Identification and selection of layers of fluid and fluid-saturated formations of one or more geological solid masses representing geological structure | |
| US9377546B2 (en) | Automatic extraction and characterization of fault and fracture populations | |
| CN102759748B (en) | The cross analysis Seismic Reservoir Prediction method that based target is analyzed | |
| CN110333543B (en) | Low-resistance body interpretation and imaging method and system based on reflection coefficient analysis | |
| CN111007064A (en) | Intelligent logging lithology identification method based on image identification | |
| CN105089642B (en) | The recognition methods of oil-gas reservoir exception producing well and device | |
| CN115113299B (en) | Intelligent marking stratum positioning and dividing method, device, equipment and storage medium | |
| CN104656140B (en) | Median filtering method and system based on vertical earthquake attenuation laws | |
| CN113570652B (en) | Sandstone reservoir mineral intercrystalline pore quantitative analysis method based on SEM image | |
| CN109444953B (en) | Automatic correction method and device for fault polygon drawing graph | |
| CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
| CN109325987B (en) | Automatic correction method and device for fault polygon drawing graph | |
| CN116201535B (en) | Automatic dividing method, device and equipment for oil and gas reservoir target well sign stratum | |
| CN102903149A (en) | Block forming method and device for geologic model | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| US20150047903A1 (en) | Analysis of geological objects | |
| CN116027454A (en) | Method for imaging contour map of sandstone thickness and intelligent drawing device | |
| CN106772596B (en) | A kind of method and device of determining pre-stack time migration velocity field | |
| CN104199107B (en) | Depth prediction approach and system before brill based on the many wave datum of vertical seismic | |
| CN112859172A (en) | Longitudinal and transverse wave micro-logging data processing method and device | |
| CN104462792B (en) | Log data lithological stratum numerical value reduction method | |
| CN112419493A (en) | Shale reservoir three-dimensional attribute model establishing method and device | |
| CN113640875B (en) | Deposition phase identification method, device and system | |
| CN111781661A (en) | Method and device for predicting sedimentary microphase plane spread in small-well area | |
| CN109079797B (en) | Distance measuring method and system of virtual clamp and mechanical arm control method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chu Weiwen Inventor after: Li Ang Inventor before: Li Ang |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200514 Address after: Room 302, unit 1, building 27, No.2, tiesan street, Nangang District, Harbin City, Heilongjiang Province 1500006 Applicant after: Chu Weiwen Address before: 110000 Wendalu 221-11, Huanggu District, Shenyang City, Liaoning Province, 1-22-1 Applicant before: Li Ang |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200612 |