CN113267815A - Filtering method and device for repeated broken edge data - Google Patents
Filtering method and device for repeated broken edge data Download PDFInfo
- Publication number
- CN113267815A CN113267815A CN202110769204.1A CN202110769204A CN113267815A CN 113267815 A CN113267815 A CN 113267815A CN 202110769204 A CN202110769204 A CN 202110769204A CN 113267815 A CN113267815 A CN 113267815A
- Authority
- CN
- China
- Prior art keywords
- broken
- broken edges
- coordinate axis
- data
- edges
- 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
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Abstract
The embodiment of the invention discloses a method and a device for filtering repeated broken edge data, wherein the method comprises the following steps: determining a main survey line direction coordinate axis and a contact survey line direction coordinate axis according to the seismic data; the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis are mutually vertical; extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data; aiming at any line number in a plurality of main measuring line numbers and a plurality of contact measuring line numbers, acquiring any two broken edges on the line number, and calculating the space distance between the two broken edges; judging whether the space distance meets a repetition condition; if yes, determining the two broken edges as repeated broken edges, and deleting any broken edge. Any two broken edge space distances on the main survey line number and the contact survey line number are calculated, repeated broken edges can be accurately determined according to the space distances for filtering, the processing efficiency is improved, and the geological structure characteristics can be more accurately determined.
Description
Technical Field
The embodiment of the invention relates to the technical field of seismic exploration, in particular to a method and a device for filtering repeated broken edge data.
Background
The broken edge is a line segment which reflects fault characteristics and is comprehensively interpreted by an interpreter according to broken points and section characteristics on a seismic section along an Inline (main survey line) or a Crossline (Crossline) by combining with a geological background of a research area structure in the seismic data interpretation process. The multiple broken edges of the same fault are combined in space to construct a plane, and then the plane forms a cross section, so the broken edges are basic units for forming the cross section.
The broken edge explanation is an important link in the construction explanation, and the result directly influences the subsequent geological modeling of the oil and gas reservoir and the describing and developing effects of the oil and gas reservoir. In the primary stage of exploration, the influence of the imaging effect of the existing seismic data and the geological recognition difference of a seismic interpreter on a research area causes the multi-solution of seismic data interpretation. The same block structure interpretation by different interpreters may have the problem of repetitive work, i.e. repetitive broken edges. Repeated broken edges can cause errors when the broken edges are combined into the broken surfaces during geological modeling, great inconvenience is caused to geological modeling and structural description, and therefore repeated broken edge data need to be filtered.
Disclosure of Invention
In view of the above problems, embodiments of the present invention are provided to provide a method and an apparatus for filtering repeated broken edge data, which overcome or at least partially solve the above problems.
According to an aspect of the embodiments of the present invention, there is provided a method for filtering repeated broken edge data, the method including:
determining a main survey line direction coordinate axis and a contact survey line direction coordinate axis according to the seismic data; the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis are mutually vertical;
extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data;
aiming at any line number in a plurality of main measuring line numbers and a plurality of contact measuring line numbers, acquiring any two broken edges on the line number, and calculating the space distance between the two broken edges;
judging whether the space distance meets a repetition condition;
if yes, determining the two broken edges as repeated broken edges, and deleting any broken edge.
According to another aspect of the embodiments of the present invention, there is provided a device for filtering repeated broken edge data, including:
the determining module is suitable for determining a main measuring line direction coordinate axis and an interconnection measuring line direction coordinate axis according to the seismic data; the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis are mutually vertical;
the extraction module is suitable for extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data;
the calculation module is suitable for acquiring any two broken edges on the line number according to any one of the main survey line numbers and the contact survey line numbers and calculating the space distance between the two broken edges;
the judging module is suitable for judging whether the space distance meets the repetition condition or not; if yes, determining the two broken edges as repeated broken edges, and deleting any broken edge.
According to still another aspect of an embodiment of the present invention, there is provided a computing device including: the system comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete mutual communication through the communication bus;
the memory is used for storing at least one executable instruction, and the executable instruction enables the processor to execute the operation corresponding to the filtering method of the repeated broken edge data.
According to another aspect of the embodiments of the present invention, there is provided a computer storage medium, where at least one executable instruction is stored, and the executable instruction causes a processor to perform an operation corresponding to the above method for filtering repeated broken edge data.
According to the filtering method and device for repeated broken edge data provided by the embodiment of the invention, any two broken edge space distances on the main measuring line number and the contact measuring line number are calculated, and the repeated broken edges can be accurately determined according to the space distances for filtering, so that the processing efficiency is improved, and the method and device are also beneficial to more finely determining the geological structure characteristics.
The foregoing description is only an overview of the technical solutions of the embodiments of the present invention, and the embodiments of the present invention can be implemented according to the content of the description in order to make the technical means of the embodiments of the present invention more clearly understood, and the detailed description of the embodiments of the present invention is provided below in order to make the foregoing and other objects, features, and advantages of the embodiments of the present invention more clearly understandable.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the embodiments of the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a flow chart illustrating a method for filtering repeat broken edge data according to one embodiment of the present invention;
FIG. 2a is a schematic diagram showing the raw data of a fault edge break;
FIG. 2b is a data diagram illustrating sampling statistics of broken edge raw data of a fault;
FIG. 2c shows a schematic of fault break data after processing;
FIG. 2d is a data schematic illustrating sample statistics of processed fault broken edge data;
FIG. 3a shows a schematic diagram of the broken edge raw data of a fault 2;
FIG. 3b is a data diagram illustrating sampling statistics of the broken edge raw data of the fault 2;
FIG. 3c shows a schematic diagram of the presence of repeated fault data between fault 2 and fault 1;
FIG. 3d shows a schematic of two fault break data after processing;
FIG. 3e shows a schematic view of the edge break data of the processed fault 2;
FIG. 3f is a data diagram illustrating sampling statistics of processed fault 2 edge break data;
fig. 4 is a schematic structural diagram of a filtering apparatus for repeated broken edge data according to an embodiment of the present invention;
FIG. 5 illustrates a schematic structural diagram of a computing device, according to an embodiment of the invention.
Detailed Description
Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention can be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Fig. 1 is a flowchart illustrating a method for filtering repeated broken edge data according to an embodiment of the present invention, and as shown in fig. 1, the method includes the following steps:
and S101, determining a main survey line direction coordinate axis and a contact survey line direction coordinate axis according to the seismic data.
The seismic data are various data obtained by processing and interpreting seismic data. According to the seismic data, the coordinates can be determined firstly, a space three-dimensional coordinate system is constructed, and data of each sampling point in the broken edge can be marked conveniently. Specifically, an active measurement line direction coordinate axis and a contact measurement line direction coordinate axis are determined, the active measurement line direction coordinate axis and the contact measurement line direction coordinate axis are perpendicular to each other, the active measurement line direction coordinate axis corresponds to an x coordinate, and the contact measurement line direction coordinate axis corresponds to a y coordinate. Further, a z coordinate axis can be established according to the direction corresponding to the time data or the depth data, the z coordinate axis is a coordinate axis perpendicular to both the main measurement line direction coordinate axis and the contact measurement line direction coordinate axis, a space stereo coordinate system is obtained, and the x coordinate, the y coordinate and the z coordinate are arranged in the space stereo coordinate system corresponding to each sampling point data in the broken edge.
The main measuring line direction coordinate axis and the contact measuring line direction coordinate axis form a plane, and the plane is a horizontal plane.
Preferably, when the main survey line direction coordinate axis and the crossline direction coordinate axis are established, the main survey line direction coordinate axis and the crossline direction coordinate axis may be created in a true north-south direction, or may be created in a proper direction according to an implementation situation, which is not limited herein.
And S102, extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data.
According to the seismic data, a plurality of main survey line numbers and a plurality of tie survey line numbers are directly extracted based on the established main survey line direction coordinate axes and tie survey line direction coordinate axes.
Step S103, for any one of the main survey line numbers and the contact survey line numbers, acquiring any two broken edges on the line number, and calculating the spatial distance between the two broken edges.
In the prior art, repeated broken edges are manually and sequentially removed in a three-dimensional space, when a work area is large and data points are numerous, considerable time and energy are consumed for manually filtering data, and reasonable broken edges are difficult to keep.
In view of the above, in the present embodiment, the distance between the broken edges is calculated, and when two broken edges are repeated broken edges, the distance is 0, so that the repeated broken edges can be detected. In order to further improve the efficiency, the distance between every two broken edges is not required to be calculated, the main line number and the interconnection line number are used for detecting line numbers one by one, and any two broken edges are selected for calculating the space distance of all the broken edges on each line number, so that the repeated broken edges on one line number can be quickly and accurately detected.
Specifically, all main survey line numbers and contact survey line numbers are detected one by one, any two broken edges on the main survey line numbers are obtained for any line number, for example, for any main survey line number, and whether the projections of any two broken edges on the horizontal plane are parallel to the coordinate axis of the main survey line direction or not is judged, or whether the projections of any two broken edges on the horizontal plane are overlapped or not is judged; if the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction, the x coordinates of the two broken edges on the same main measuring line number are the same data, processing is not needed, and the first spatial distance can be calculated by using a first formula according to the y coordinates of various sample points in any two broken edges. If the projections of any two broken edges on the horizontal plane are superposed, the projections of the two broken edges on the horizontal plane are the same point or the same line, and the first spatial distance is calculated by using a first formula according to the y coordinates of various points in any two broken edges. In addition to the above situation, it is necessary to calculate both the x coordinate and the y coordinate of each sample point in two broken edges, calculate a first spatial distance by using a first formula according to the y coordinate of each sample point in any two broken edges, and calculate a second spatial distance by using a second formula according to the x coordinate of each sample point in any two broken edges. For any one tie survey line number, acquiring any two broken edges on the tie survey line number; judging whether the projections of any two broken edges on the horizontal plane are parallel to the coordinate axis of the direction of the contact measuring line or not, or judging whether the projections of any two broken edges on the horizontal plane are overlapped or not; if the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the direction of the contact survey line, the y coordinates of the two broken edges on the same contact survey line number are the same data, processing is not needed, and the second spatial distance can be calculated by using a second formula according to the x coordinates of various sample points in any two broken edges. If the two broken edges are overlapped, the projection of the two broken edges on the horizontal plane is the same point or the same line, and the second space distance is calculated by using a second formula according to the x coordinates of various points in any two broken edges. In addition to the above, it is necessary to calculate both the x coordinate and the y coordinate of each sample point in two broken edges, calculate a first spatial distance by using a first formula according to the y coordinate of each sample point in any two broken edges, and calculate a second spatial distance by using a second formula according to the x coordinate of each sample point in any two broken edges.
Wherein the first formula is:
the second formula is:
in the formula s1,y(t) and s2,y(t) y coordinates corresponding to the t data points of any two broken edges on the main measuring line number respectively; s1,x(t) and s2,x(t) x coordinates corresponding to any two broken edges on the line number of the contact measuring line at the t data point respectively; the t data points are time data points or depth data points; t is t1Is the minimum time data point or depth data point, t, of two broken edge samplesNThe maximum time data point or the depth data point of the sampling points of the two broken edges. The abs function is an absolute value function of the data.
Seismic data interpretation is typically performed in the time domain, so the t data points are typically time of day data points prior to the time-depth transform, and the t data points are depth data points when the seismic data is depth domain data. In the following description, the t data point is taken as a time data point as an example, and the processing of the depth data point is the same as the processing of the time data point. Acquiring the y coordinate, s of the two broken edges 1 and 2 on the main measuring line number at the same t data point1,y(t) and s2，y(t) of (d). Calculating the absolute value of the distance between two y coordinates, t, using a first formula1Is the minimum time data point, t, of two broken edge samplesNAccumulating the maximum time data point of two broken edge sample points1To tNTwo seats corresponding to data points at each momentThe absolute value of the target distance is obtained to obtain a first space distance Disty. Since the calculated distance in the first formula is an absolute value of the coordinate difference, the finally obtained first spatial distance is 0 only when the coordinate difference of each t data point is 0. For the second formula, the x coordinate, s, of the same t data point of the two broken edges 1 and 2 on the line number of the contact measuring line is obtained1,x(t) and s2，x(t) of (d). Calculating the absolute value of the distance between two x coordinates, t, using a second formula1Is the minimum time data point, t, of two broken edge samplesNAccumulating the maximum time data point of two broken edge sample points1To tNThe absolute value of the distance between the two x coordinates corresponding to the data point at each moment is used for obtaining a second spatial distance Distx. Since the calculated distance in the second formula is an absolute value of the coordinate difference, the finally obtained second spatial distance is 0 only when the coordinate difference of each t data point is 0.
And step S104, judging whether the space distance meets the repeated condition.
When the judgment is carried out, whether the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction or not is judged, or whether the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the contact measuring line direction or not is judged, or whether the projections of any two broken edges on the horizontal plane are superposed or not is judged, if the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction and the x coordinates are the same, whether the first space distance obtained by calculating the y coordinate is a specified value or not is judged, and if the first space distance is the specified value, the repetition condition is met. If the projections of the two broken edges on the horizontal plane are parallel to coordinate axes in the direction of the contact measuring line and the y coordinates are the same, judging whether the second space distance obtained by the calculation of the x coordinate is a specified value or not, and if the second space distance is the specified value, meeting the repeated condition; or, if the projections of the two broken edges on the horizontal plane are overlapped and the projections are the same point or the same line, judging whether the first distance is a specified value or not, or judging whether the second distance is the specified value or not, and judging that the first distance meets the repetition condition by selecting either one of the first distance and the second distance.
Except the above situation, it is necessary to simultaneously determine whether the first spatial distance and the second spatial distance are both specified values, and only when the first spatial distance and the second spatial distance are both specified values, the repetition condition can be satisfied, and the two broken edges are repeated data.
Wherein the specified value is 0, i.e. the spatial distance is 0.
And step S105, determining that the two broken edges are repeated broken edges, and deleting any broken edge.
When the spatial data meet the repeated condition, two broken edges can be determined as repeated broken edges, and at the moment, any broken edge can be directly deleted to filter the repeated data.
Steps S103 to S105 are executed in a cyclic manner, and after the spatial distance between any two broken edges on one main measurement line number is determined, the spatial distances between other broken edges are determined until the spatial distances between all broken edges on each main measurement line number and the spatial distances between all broken edges on each contact measurement line number are determined, the detection of all broken edges is completed, and repeated data is filtered out. For the main survey line number or the cross survey line number, the spatial distance of any two broken edges on the main survey line number can be calculated firstly, the spatial distance of any two broken edges on the cross survey line number can also be calculated firstly, and the calculation sequence is not limited.
In specific implementation of this embodiment, the seismic data may include one or more faults, and the fault includes a plurality of broken edges.
When seismic data is processed, one fault can be processed, as shown in fig. 2a, 1 fault is included, and the broken edge raw data of the fault 1 is counted as shown in fig. 2b, and the statistical data shows that the fault includes 265 sampling points, wherein the abscissa in the graph is depth, the unit is meter, and the ordinate is proportion. However, the data are found to be inconsistent with the statistical data during manual processing, and repeated broken edge data exist. The above steps are executed by the present embodiment, and the processed edge breaking data of the fault 1 shown in fig. 2c and the statistical data of the processed edge breaking data shown in fig. 2d are obtained. Although the edge break in fig. 2c after the repeated data is filtered out is visually not different from the original edge break data in fig. 2a, it can be found from the data statistics result in fig. 2d that the number of sampling points is reduced from 265 to 53, and all the repeated edge break data is filtered out. The repeated data can be deleted manually for 1 hour, and the repeated broken edge data can be processed in only 1 minute by using the method, so that the processing efficiency is greatly improved.
When seismic data is processed, a plurality of faults can be processed. As shown in fig. 3a, in the broken edge original data of the fault 2, not only the repeated broken edge data exists in the fault 2 itself, but also the repeated data exists in the adjacent fault 1. As shown in fig. 3c, the thick line broken edge represents a fault 1, the dotted line represents a fault 2, and the broken edge where the thick line and the dotted line overlap is a portion where both lines overlap. After the fault 2 filters the repeated data, as shown in fig. 3e, compared with the broken edge original data of the fault 2 in fig. 3a, 11 broken edges are visually filtered, and 425 repeated data points are actually filtered. Fig. 3b is a diagram illustrating data statistics on the original data to obtain 508 samples, and 83 samples remain after filtering out the repeated data, as shown in fig. 3 f. The fault 2 is displayed together with the fault 1 after repeated data are filtered out, as shown in fig. 3d, all non-repeated data are retained, and the fault 1 and the fault 2 can be merged because the fault 1 and the fault 2 are similar in shape and small in closing difference. For the situation that the fault 2 and the fault 1 have repeated data, the repeated data can be deleted manually for 4 hours, the processing of repeated broken edge data can be completed in only 1 minute by using the method for processing, and the processing efficiency is greatly improved.
According to the filtering method for the repeated broken edge data, provided by the embodiment of the invention, any two broken edge space distances on the main measuring line number and the contact measuring line number are calculated, so that the repeated broken edge can be accurately determined according to the space distances for filtering, the processing efficiency is improved, and the method is also beneficial to more finely determining the geological structure characteristics.
Fig. 4 is a schematic structural diagram illustrating a device for filtering repeated broken edge data according to an embodiment of the present invention. As shown in fig. 4, the apparatus includes:
a determination module 410 adapted to determine a main line direction coordinate axis and a crossline direction coordinate axis from the seismic data; the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis are mutually vertical;
the extraction module 420 is suitable for extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data;
the calculating module 430 is adapted to obtain any two broken edges on a plurality of main survey line numbers and any one line number in a plurality of contact survey line numbers, and calculate the spatial distance between the two broken edges;
a judging module 440 adapted to judge whether the spatial distance satisfies a repetition condition; if yes, determining the two broken edges as repeated broken edges, and deleting any broken edge.
Optionally, the main survey line direction coordinate axis corresponds to an x coordinate; the coordinate axis of the direction of the contact line corresponds to the y coordinate.
Optionally, the calculation module 430 is further adapted to:
acquiring any two broken edges on any main measuring line number;
judging whether the projections of any two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction or not, or judging whether the projections of any two broken edges on the horizontal plane are overlapped or not;
if yes, calculating by using a first formula according to the y coordinates of various sample points in any two broken edges to obtain a first space distance;
if not, calculating by using a first formula according to the y coordinates of the various sample points in any two broken edges to obtain a first space distance; and according to the x coordinates of various points in any two broken edges, a second formula is used for calculating to obtain a second space distance.
Optionally, the calculation module 430 is further adapted to:
acquiring any two broken edges on any one contact measuring line number;
judging whether the projections of any two broken edges on the horizontal plane are parallel to the coordinate axis of the direction of the contact measuring line or not, or judging whether the projections of any two broken edges on the horizontal plane are overlapped or not;
if yes, calculating by using a second formula according to the x coordinates of all sample points in any two broken edges to obtain a second space distance;
if not, calculating by using a first formula according to the y coordinates of the various sample points in any two broken edges to obtain a first space distance; and according to the x coordinates of various points in any two broken edges, a second formula is used for calculating to obtain a second space distance.
Optionally, the first formula is:
the second formula is:
wherein s is1,y(t) and s2,y(t) y coordinates corresponding to the t data points of any two broken edges on the main measuring line number respectively; s1，x(t) and s2，x(t) x coordinates corresponding to any two broken edges on the line number of the contact measuring line at the t data point respectively; the t data points are time data points or depth data points; t is t1Is the minimum time data point or depth data point, t, of the sample points in two broken edgesNThe maximum time data point or the depth data point of the sampling points in the two broken edges.
Optionally, the determining module 440 is further adapted to:
judging whether the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction or not, or whether the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the contact measuring line direction or not, or judging whether the projections of any two broken edges on the horizontal plane are superposed or not;
if yes, judging whether the first space distance or the second space distance is a specified value;
if not, judging whether the first space distance and the second space distance are both specified values.
Optionally, the seismic data comprises one or more faults; the fault comprises a plurality of broken edges.
The descriptions of the modules refer to the corresponding descriptions in the method embodiments, and are not repeated herein.
The embodiment of the invention also provides a nonvolatile computer storage medium, wherein the computer storage medium stores at least one executable instruction, and the executable instruction can execute the filtering method of the repeated broken edge data in any method embodiment.
Fig. 5 is a schematic structural diagram of a computing device according to an embodiment of the present invention, and a specific embodiment of the present invention does not limit a specific implementation of the computing device.
As shown in fig. 5, the computing device may include: a processor (processor)502, a Communications Interface 504, a memory 506, and a communication bus 508.
The method is characterized in that:
the processor 502, communication interface 504, and memory 506 communicate with one another via a communication bus 508.
A communication interface 504 for communicating with network elements of other devices, such as clients or other servers.
The processor 502 is configured to execute the program 510, and may specifically execute the relevant steps in the above embodiment of the method for filtering repeated broken edge data.
In particular, program 510 may include program code that includes computer operating instructions.
The processor 502 may be a central processing unit CPU, or an Application Specific Integrated Circuit (ASIC), or one or more Integrated circuits configured to implement an embodiment of the present invention. The computing device includes one or more processors, which may be the same type of processor, such as one or more CPUs; or may be different types of processors such as one or more CPUs and one or more ASICs.
And a memory 506 for storing a program 510. The memory 506 may comprise high-speed RAM memory, and may also include non-volatile memory (non-volatile memory), such as at least one disk memory.
The program 510 may be specifically configured to enable the processor 502 to execute the method for filtering repeated broken edge data in any of the method embodiments described above. For specific implementation of each step in the program 510, reference may be made to corresponding steps and corresponding descriptions in units in the foregoing embodiments for filtering repeated broken edge data, which are not described herein again. It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described devices and modules may refer to the corresponding process descriptions in the foregoing method embodiments, and are not described herein again.
The algorithms or displays presented herein are not inherently related to any particular computer, virtual system, or other apparatus. Various general purpose systems may also be used with the teachings herein. The required structure for constructing such a system will be apparent from the description above. In addition, embodiments of the present invention are not directed to any particular programming language. It is appreciated that a variety of programming languages may be used to implement the teachings of embodiments of the present invention as described herein, and any descriptions of specific languages are provided above to disclose preferred embodiments of the invention.
In the description provided herein, numerous specific details are set forth. It is understood, however, that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Similarly, it should be appreciated that in the foregoing description of exemplary embodiments of the invention, various features of the embodiments of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. However, the disclosed method should not be interpreted as reflecting an intention that: that is, the claimed embodiments of the invention require more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this invention.
Those skilled in the art will appreciate that the modules in the device in an embodiment may be adaptively changed and disposed in one or more devices different from the embodiment. The modules or units or components of the embodiments may be combined into one module or unit or component, and furthermore they may be divided into a plurality of sub-modules or sub-units or sub-components. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or elements of any method or apparatus so disclosed, may be combined in any combination, except combinations where at least some of such features and/or processes or elements are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
Furthermore, those skilled in the art will appreciate that while some embodiments herein include some features included in other embodiments, rather than other features, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the following claims, any of the claimed embodiments may be used in any combination.
The various component embodiments of the invention may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that a microprocessor or Digital Signal Processor (DSP) may be used in practice to implement some or all of the functionality of some or all of the components according to embodiments of the present invention. Embodiments of the invention may also be implemented as apparatus or device programs (e.g., computer programs and computer program products) for performing a portion or all of the methods described herein. Such programs implementing embodiments of the present invention may be stored on computer-readable media or may be in the form of one or more signals. Such a signal may be downloaded from an internet website or provided on a carrier signal or in any other form.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. Embodiments of the invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The usage of the words first, second and third, etcetera do not indicate any ordering. These words may be interpreted as names. The steps in the above embodiments should not be construed as limiting the order of execution unless specified otherwise.
Claims (10)
1. A method for filtering repeated broken edge data is characterized by comprising the following steps:
determining a main survey line direction coordinate axis and a contact survey line direction coordinate axis according to the seismic data; the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis are mutually vertical;
extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data;
aiming at any line number in the main measuring line numbers and the contact measuring line numbers, acquiring any two broken edges on the line number, and calculating the space distance between the two broken edges;
judging whether the space distance meets a repetition condition;
if yes, determining the two broken edges as repeated broken edges, and deleting any broken edge.
2. The method of claim 1, wherein the inline direction coordinate axis corresponds to an x coordinate; the coordinate axis of the direction of the contact measuring line corresponds to a y coordinate; and a plane formed by the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis is a horizontal plane.
3. The method of claim 2, wherein the step of obtaining any two broken ribs on the line number for any one of the main line numbers and the cross line numbers further comprises:
acquiring any two broken edges on any main measuring line number;
judging whether the projections of any two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction or not, or judging whether the projections of any two broken edges on the horizontal plane are overlapped or not;
if yes, calculating by using a first formula according to the y coordinates of the various sample points in the any two broken edges to obtain a first space distance;
if not, calculating by using a first formula according to the y coordinates of the various sample points in the any two broken edges to obtain a first space distance; and according to the x coordinates of the various points in any two broken edges, calculating by using a second formula to obtain a second space distance.
4. The method of claim 2, wherein the step of obtaining any two broken ribs on the line number for any one of the main line numbers and the cross line numbers further comprises:
acquiring any two broken edges on any one contact measuring line number;
judging whether the projections of any two broken edges on the horizontal plane are parallel to the coordinate axis of the direction of the contact measuring line or not, or judging whether the projections of any two broken edges on the horizontal plane are overlapped or not;
if yes, calculating by using a second formula according to the x coordinates of the various sample points in the any two broken edges to obtain a second space distance;
if not, calculating by using a first formula according to the y coordinates of the various sample points in the any two broken edges to obtain a first space distance; and according to the x coordinates of the various points in any two broken edges, calculating by using a second formula to obtain a second space distance.
5. The method according to claim 3 or 4,
the first formula is:
the second formula is:
wherein s is1,y(t) and s2,y(t) y coordinates corresponding to the t data points of any two broken edges on the main measuring line number respectively; s1，x(t) and s2,x(t) x coordinates corresponding to any two broken edges on the line number of the contact measuring line at the t data point respectively; the t data points are time data points or depth data points; t is t1Is the minimum time data point or depth data point, t, of the sample points in two broken edgesNThe maximum time data point or the depth data point of the sampling points in the two broken edges.
6. The method of claim 3 or 4, wherein the determining whether the spatial distance satisfies a repetition condition further comprises:
judging whether the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the main measuring line direction, or whether the projections of the two broken edges on the horizontal plane are parallel to the coordinate axis of the contact measuring line direction, or judging whether the projections of any two broken edges on the horizontal plane are overlapped;
if yes, judging whether the first space distance or the second space distance is a specified value;
if not, judging whether the first space distance and the second space distance are both specified values.
7. The method of claim 1, wherein the seismic data comprises one or more faults; the fault comprises a plurality of broken edges.
8. A device for filtering repeated broken edge data, the device comprising:
the determining module is suitable for determining a main measuring line direction coordinate axis and an interconnection measuring line direction coordinate axis according to the seismic data; the main measuring line direction coordinate axis and the contact measuring line direction coordinate axis are mutually vertical;
the extraction module is suitable for extracting a plurality of main survey line numbers and a plurality of cross survey line numbers from the seismic data;
the calculation module is suitable for acquiring any two broken edges on the line number according to any one line number in the main measuring line numbers and the contact measuring line numbers and calculating the space distance between the two broken edges;
the judging module is suitable for judging whether the space distance meets a repetition condition or not; if yes, determining the two broken edges as repeated broken edges, and deleting any broken edge.
9. A computing device, comprising: the system comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete mutual communication through the communication bus;
the memory is used for storing at least one executable instruction, and the executable instruction causes the processor to execute the operation corresponding to the repeated broken edge data filtering method of any one of claims 1-7.
10. A computer storage medium having at least one executable instruction stored therein, the executable instruction causing a processor to perform operations corresponding to the method for filtering repeated broken edge data according to any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110769204.1A CN113267815B (en) | 2021-07-07 | 2021-07-07 | Method and device for filtering repeated broken edge data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110769204.1A CN113267815B (en) | 2021-07-07 | 2021-07-07 | Method and device for filtering repeated broken edge data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113267815A true CN113267815A (en) | 2021-08-17 |
| CN113267815B CN113267815B (en) | 2022-05-10 |
Family
ID=77236485
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110769204.1A Active CN113267815B (en) | 2021-07-07 | 2021-07-07 | Method and device for filtering repeated broken edge data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113267815B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166163A (en) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube |
| CN104199099A (en) * | 2014-09-19 | 2014-12-10 | 东北石油大学 | Efficient seismic interpretation method through multi-window continuous sections |
| CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
| CN106526675A (en) * | 2016-12-28 | 2017-03-22 | 中国石油化工股份有限公司 | Automatic extraction method for fault space parameter |
| CN108254784A (en) * | 2017-12-11 | 2018-07-06 | 中国石油天然气股份有限公司 | A kind of fault recognizing method, apparatus and system based on two-dimension earthquake data |
| CN110727025A (en) * | 2019-08-21 | 2020-01-24 | 中国石油化工股份有限公司 | Hidden fault recognition method |
| CN111913219A (en) * | 2020-08-13 | 2020-11-10 | 中海石油(中国)有限公司 | Method for determining fault control deposition |
| CN112305585A (en) * | 2019-07-26 | 2021-02-02 | 中国石油天然气集团有限公司 | Three-dimensional data volume acquisition method and device |
| CN113031066A (en) * | 2021-01-24 | 2021-06-25 | 东北石油大学 | Method for identifying hidden fault |
-
2021
- 2021-07-07 CN CN202110769204.1A patent/CN113267815B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166163A (en) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube |
| CN104199099A (en) * | 2014-09-19 | 2014-12-10 | 东北石油大学 | Efficient seismic interpretation method through multi-window continuous sections |
| CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
| CN106526675A (en) * | 2016-12-28 | 2017-03-22 | 中国石油化工股份有限公司 | Automatic extraction method for fault space parameter |
| CN108254784A (en) * | 2017-12-11 | 2018-07-06 | 中国石油天然气股份有限公司 | A kind of fault recognizing method, apparatus and system based on two-dimension earthquake data |
| CN112305585A (en) * | 2019-07-26 | 2021-02-02 | 中国石油天然气集团有限公司 | Three-dimensional data volume acquisition method and device |
| CN110727025A (en) * | 2019-08-21 | 2020-01-24 | 中国石油化工股份有限公司 | Hidden fault recognition method |
| CN111913219A (en) * | 2020-08-13 | 2020-11-10 | 中海石油(中国)有限公司 | Method for determining fault control deposition |
| CN113031066A (en) * | 2021-01-24 | 2021-06-25 | 东北石油大学 | Method for identifying hidden fault |
Non-Patent Citations (2)
| Title |
|---|
| 毛凤军等: "尼日尔Termit盆地三维地质构造建模研究与应用", 《地学前缘》 * |
| 田静: "基于聚类分析的断层自动拆分技术及应用", 《国外测井技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113267815B (en) | 2022-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103325085B (en) | The methods, devices and systems of automated graphics alignment | |
| EP3317852B1 (en) | Method in constructing a model of a scenery and device therefor | |
| Vuan et al. | Improving the detection of low‐magnitude seismicity preceding the Mw 6.3 L’Aquila earthquake: Development of a scalable code based on the cross correlation of template earthquakes | |
| CN106772592A (en) | The analysis method and device of diffracted wave focus energy | |
| CN116340072B (en) | Fault detection method and device for multi-core chip | |
| CN110969092A (en) | Pulse signal denoising method and device and terminal equipment | |
| CN112526470A (en) | Method and device for calibrating radar parameters, electronic equipment and storage medium | |
| CN113537023B (en) | Remote sensing image semantic change detection method | |
| CN113267815B (en) | Method and device for filtering repeated broken edge data | |
| CN111027601A (en) | Plane detection method and device based on laser sensor | |
| CN112380308B (en) | Geochemical anomaly delineating method and system based on data regularization | |
| US20130191071A1 (en) | System and method for automatic modal parameter extraction in structural dynamics analysis | |
| CN115248993A (en) | Method and device for detecting reality of simulation scene model and storage medium | |
| CN111929730B (en) | Small-scale geological anomaly detection method and device | |
| CN109387873A (en) | A kind of fracture and cave reservoir inversion method and system | |
| NO20181090A1 (en) | Visualizing attributes of multiple fault surfaces in real time | |
| CN112764102A (en) | Seismic data acquisition method and device | |
| CN109471161A (en) | The method for building up and device of geological structure attribute section | |
| CN115903043A (en) | Diffracted wave separation method and device | |
| CN111258894B (en) | Method and device for evaluating software risk, storage medium and electronic equipment | |
| CN115937652A (en) | Underground pipeline identification method, device, equipment and storage medium | |
| CN117853902A (en) | Slope stability analysis method | |
| CN107589444B (en) | Seismic data processing method and device | |
| CN117472491A (en) | Map rendering method, system, medium and electronic device for pipeline data | |
| CN117910883A (en) | Logging data quality detection method, device, computing equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |