Crack extraction method and crack extraction device

Abstract

The invention provides a crack extraction method and a crack extraction device, and the method comprises the steps: carrying out the line-by-line peak scanning of a plane space of seismic attribute data in at least two directions, so as to obtain peak points, marking the peak points on the plane space of the seismic attribute data, and determining a graph formed by the peak points as a first crack trace; enhancing the continuity of the first crack trace by using an expansion-erosion method to form a second crack trace, determining a part with a length greater than or equal to a crack length preset value in the second crack trace as a third crack trace, performing scale division on a crack corresponding to the crack trace, and aiming at different scales, performing scale division on the third crack trace; and performing thickness and / or color characterization processing on each part of the third crack trace to form a fourth crack trace, and outputting the fourth crack trace. According to the crack extraction method provided by the invention, the seismic attribute data is used as input, and the fine crack trace is extracted by using a multi-direction peak scanning technology. And multi-scale fracture extraction is realized according to the seismic attribute intensity of the position of the fracture trace.

Description

Technical Field

[0001] The invention relates to the technical field of oil and gas geological exploration and development, and in particular to a fracture extraction method and a fracture extraction device. Background Art

[0002] Fractures are important seepage channels and storage spaces in tight reservoirs. Improving the ability to predict and characterize fractures is of great significance to oil and gas reservoir evaluation, well site deployment, and development plan formulation. Among them, reservoir fracture extraction based on seismic attributes is one of the important fracture prediction methods. Its general process is to identify fractures based on seismic attributes, represent the fracture results in binary form, skeletonize the binary fracture identification results, obtain fracture traces, and complete fracture extraction.

[0003] When extracting cracks from seismic attributes using the above method, there are two problems: First, weak seismic attributes are difficult to effectively use in crack identification, resulting in too many small and medium-sized cracks being lost during crack extraction. When predicting cracks through seismic attributes, the strength of the attribute generally indicates the scale and intensity of crack development, but when the attribute intensity is weak, its characteristics will be suppressed by the high-intensity attributes nearby, making it difficult to effectively identify; Second, when extracting cracks from the seismic attribute crack identification results, the two existing extraction methods, the morphological-based corrosion algorithm and the centerline trajectory based on the maximum inscribed circle, are based on binary images. There are large deviations in the extraction results, and the crack traces are distorted and deformed, making it difficult to accurately extract cracks at different scales. Summary of the invention

[0004] The present invention aims to solve at least one of the technical problems existing in the prior art. To this end, the present invention provides a crack extraction method and a crack extraction device, which aims to solve the problem that the crack extraction method of the related art loses too many small and medium-sized cracks when extracting cracks, and the extraction result has a large deviation, making it difficult to accurately extract the cracks by scale.

[0005] The present invention provides a crack extraction method, comprising: Obtain seismic attribute data and preset values ​​of fracture trace length; Based on the seismic attribute data, performing row-by-row peak scanning along at least two directions on the plane space of the seismic attribute data to obtain peak points, and marking the peak points on the plane space of the seismic attribute data, and determining that a graph composed of the peak points is a first fracture trace; Based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as a second crack trace; Based on the second crack trace and the preset value of the crack trace length, determining a portion of the second crack trace having a length greater than or equal to the preset value of the crack trace length as a third crack trace; Based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the third fracture trace is scaled, wherein the larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale is; Based on the third crack trace after scale division, thickness and / or color characterization processing is performed on each part of the third crack trace according to different scales to form a fourth crack trace and output it.

[0006] According to the fracture extraction method provided by the present invention, the seismic attribute data is plane seismic attribute data or slice data in three-dimensional seismic attribute data.

[0007] According to the fracture extraction method provided by the present invention, performing row-by-row peak scanning along at least two directions in the plane space of the seismic attribute data to obtain peak points includes: Determine a first scanning direction, and scan the plane space of the seismic attribute data line by line along a direction perpendicular to the first scanning direction, each scan is along the first scanning direction, and each scan extracts the seismic attribute data corresponding to the scanned line, and records the seismic attribute data as V i , where V i represents the earthquake attribute data of the i-th row; Determine a second scanning direction, and scan the plane space of the seismic attribute data line by line along a direction perpendicular to the second scanning direction, each scan is along the second scanning direction, and each scan extracts the seismic attribute data corresponding to the scanned line, and the seismic attribute data is V j , where V j represents the earthquake attribute data of the jth row; For all V i and V j Perform a peak scan and identify all peak points.

[0008] According to the crack extraction method provided by the present invention, the first scanning direction is orthogonal to the second scanning direction.

[0009] According to the fracture extraction method provided by the present invention, based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the fracture corresponding to the fracture trace is scaled, including: The seismic attribute intensity values ​​are graded, wherein the higher the seismic attribute intensity value is, the higher the corresponding grade is; Superimposing the third fracture trace with the seismic intensity attribute strength value, and extracting the seismic attribute strength value at each position on the third fracture trace; Based on the division level of the seismic attribute intensity value and the extracted seismic attribute intensity value of each position on the third fracture trace, each part of the third fracture trace is scaled. The higher the level of the seismic attribute intensity value corresponding to each position of the third fracture trace, the larger the scale of the corresponding position of the third fracture trace.

[0010] According to the fracture extraction method provided by the present invention, the third fracture trace is divided into at least large-scale fracture traces, medium-scale fracture traces and small-scale fracture traces.

[0011] The present invention also provides a crack extraction device, comprising: A data acquisition module, used to acquire seismic attribute data and a preset value of fracture trace length; A first crack trace extraction module is used to perform row-by-row peak scanning along at least two directions on the plane space of the seismic attribute data based on the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and extract and determine a graph composed of the peak points as a first crack trace; A second crack trace generation module, configured to enhance the continuity of the first crack trace by using a dilation-erosion method based on the first crack trace, and generate the first crack trace after the enhanced continuity as a second crack trace; A third crack trace generating module, configured to generate, based on the second crack trace and the preset value of the crack trace length, a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length as a third crack trace; A scale division module, for performing scale division on the cracks corresponding to the third crack trace based on the third crack trace and the seismic attribute intensity value in the seismic attribute data, wherein the larger the seismic attribute intensity value corresponding to the third crack trace is, the larger the scale of the corresponding part is; The fourth crack trace output module is used to perform thickness and / or color characterization processing on each part of the third crack trace according to different scales based on the third crack trace after scale division, to form and output a fourth crack trace.

[0012] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the above-mentioned crack extraction method when executing the program.

[0013] The present invention also provides a computer program product, comprising a computer program stored on a non-transitory computer-readable storage medium, wherein the computer program comprises program instructions, and the program instructions are used to execute the steps of the above-mentioned crack extraction method through a computer.

[0014] The present invention also provides a non-transitory computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps of the crack extraction method described above are implemented.

[0015] The present invention adopts the above technical solution, which has the following advantages: The crack extraction method provided by the present invention comprises the following steps: obtaining seismic attribute data and a preset value of the length of a crack trace; based on the seismic attribute data, performing a line-by-line peak scan along at least two directions on the plane space of the seismic attribute data to obtain peak points, marking the peak points on the plane space of the seismic attribute data, and determining a graph composed of the peak points as a first crack trace; based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as a second crack trace; based on the second crack trace and the preset value of the length of the crack trace, determining a portion of the second crack trace whose length is greater than or equal to the preset value of the length of the crack trace as a third crack trace; based on the third crack trace and the seismic attribute intensity value in the seismic attribute data, performing scale division on the cracks corresponding to the crack trace, wherein the larger the seismic attribute intensity value corresponding to the third crack trace is, the larger the scale of the corresponding portion is; based on the third crack trace after the scale division, performing thickness and / or color characterization processing on each portion of the third crack trace for different scales, forming a fourth crack trace and outputting it. The fracture extraction method provided by the present invention uses seismic attribute data as input and multi-directional peak scanning technology to extract fine fracture traces. According to the seismic attribute intensity of the position where the fracture trace is located, multi-scale fracture extraction is achieved. The fracture extraction method provided by the present invention can identify weak fracture response information in seismic attributes, accurately locate fractures and realize multi-scale fracture extraction. The operation is simple, the extracted fracture position is more accurate, the connection relationship is more natural, and the result is reliable. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 is a schematic flow chart of a crack extraction method provided by an embodiment of the present invention; Figure 2 is a plane diagram of earthquake attribute data provided by an embodiment of the present invention; Figure 3 is a schematic diagram of a 90° scanning line provided by an embodiment of the present invention; Figure 4 is a peak point schematic diagram provided by an embodiment of the present invention; Figure 5 is a schematic diagram of a fourth crack trace provided by an embodiment of the present invention; Figure 6 is a schematic diagram of crack traces extracted by the existing method provided by the present invention; Figure 7 is a schematic diagram of crack traces extracted by the crack extraction method provided by the present invention; Figure 8 It is a structural schematic diagram of the electronic device provided by the present invention.

[0018] Reference numerals: 810: processor; 820: communication interface; 830: memory; 840: communication bus. DETAILED DESCRIPTION

[0019] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0020] In the description of the present invention, it is to be understood that the terms “center”, “longitudinal”, “lateral”, “length”, “width”, “thickness”, “up”, “down”, “front”, “back”, “left”, “right”, “vertical”, “horizontal”, “top”, “bottom”, “inside”, “outside”, “clockwise”, “counterclockwise”, “axial”, “radial”, “circumferential”, etc., indicating orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the referred device or element must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as limiting the present invention.

[0021] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. In the description of the present invention, the meaning of "plurality" is more than two, unless otherwise clearly and specifically defined.

[0022] In the present invention, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", "fixed" and the like should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium, it can be the internal connection of two elements or the interaction relationship between two elements. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0023] In the present invention, unless otherwise clearly specified and limited, a first feature being "above" or "below" a second feature may mean that the first and second features are in direct contact, or the first and second features are in indirect contact through an intermediate medium. Moreover, a first feature being "above", "above" or "above" a second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is higher in level than the second feature. A first feature being "below", "below" or "below" a second feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature is lower in level than the second feature.

[0024] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0025] The crack extraction method provided by the present invention comprises the following steps: obtaining seismic attribute data and a preset value of the length of a crack trace; based on the seismic attribute data, performing a line-by-line peak scan on the plane space of the seismic attribute data in at least two directions to obtain peak points, marking the peak points on the plane space of the seismic attribute data, and determining a graph composed of the peak points as a first crack trace; based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as a second crack trace; based on the second crack trace and a preset value of the length of the crack trace, determining a portion of the second crack trace whose length is greater than or equal to the preset value of the length of the crack trace as a third crack trace; based on the third crack trace and the seismic attribute intensity value in the seismic attribute data, performing scale division on the cracks corresponding to the crack trace, wherein the larger the seismic attribute intensity value corresponding to the third crack trace is, the larger the scale of the cracks in the corresponding portion is; based on the third crack trace after the scale division, performing thickness and / or color characterization processing on each portion of the third crack trace for different scales, forming a fourth crack trace and outputting it. The fracture extraction method provided by the present invention uses seismic attribute data as input and multi-directional peak scanning technology to extract fine fracture traces. According to the seismic attribute intensity of the position where the fracture trace is located, multi-scale fracture extraction is achieved. The fracture extraction method provided by the present invention can identify weak fracture response information in seismic attributes, accurately locate fractures and realize multi-scale fracture extraction. The operation is simple, the extracted fracture position is more accurate, the connection relationship is more natural, and the result is reliable.

[0026] Combine the following Figures 1 to 7 The fracture extraction method of the present invention is described.

[0027] An embodiment of the present invention provides a crack extraction method, which comprises the following steps: Step S100, obtaining seismic attribute data and a preset value of the fracture trace length; Step S200: Based on the seismic attribute data, perform row-by-row peak scanning in at least two directions on the plane space of the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and determine that a graph composed of the peak points is a first fracture trace; Step S300: Based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as the second crack trace; Step S400: Based on the second crack trace and the preset value of the crack trace length, a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length is determined as a third crack trace; Step S500: based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the fracture corresponding to the fracture trace is scaled, wherein the larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale of the corresponding part is; Step S600: Based on the scale-divided third crack trace, for different scales, perform thickness and / or color characterization processing on each part of the third crack trace to form a fourth crack trace and output it.

[0028] Specifically, according to Figure 1 As shown in the process, first obtain the seismic attribute data and the preset value of the fracture trace length. The seismic attribute data is displayed as follows Figure 2 Then, on the plane space of the seismic attribute data, peak scans are performed in at least two directions to obtain peak points, and then the obtained peak points are marked on the plane space of the seismic attribute data, and the graph composed of the peak points is determined to be the first fracture trace.

[0029] Taking peak scanning in two directions as an example, the specific operation process is as follows: Select a scan line direction in the plane space of seismic attribute data, such as the 90° direction. Figure 3 As shown; Seismic attribute data are extracted row by row from the top of the plane space of seismic attribute data from top to bottom. The scanning direction of each time is 90°. The extracted seismic attribute data is recorded as v i , v i Represents the i-th row of data; For the seismic attribute data v extracted above i Perform peak scan and determine all peak points, such as Figure 4 , and mark all peak point locations in the plane space of seismic attribute data; For all earthquake attribute data v i After completing the peak scan, select another scan line direction that is best orthogonal to the last selected scan line direction, such as a 180° direction; Seismic attribute data are extracted row by row from left to right from the left side of the plane space of seismic attribute data. The scanning direction of each time is 180°. The extracted seismic attribute data is recorded as v j , v j Indicates the jth row of data; For the seismic attribute data v extracted above j Perform peak scan, determine all peak points, and mark all peak point positions in the plane space of seismic attribute data; After completing the peak scan in two directions, the graph composed of the peak point positions marked in the plane space of the seismic attribute data is determined as the first crack trace. For example, the position corresponding to the first crack trace can be assigned a value of 1, and the unmarked position can be assigned a value of 0, which is convenient for subsequent calculation and processing.

[0030] Of course, the scanning direction may be increased and the above steps may be repeated to enhance the display effect of the first crack trace. Generally, scanning in two orthogonal directions may obtain the first crack trace with a better display effect.

[0031] Then, the continuity of the first crack trace is enhanced by a morphological dilation-erosion method, and the first crack trace after the continuity enhancement is determined as the second crack trace.

[0032] Then, based on the second crack trace and the preset value of the crack trace length, the part whose length is less than the preset value of the crack trace length is eliminated, as it is more likely to be a noise point, and the part of the second crack trace whose length is greater than or equal to the preset value of the crack trace length is determined as the third crack trace.

[0033] After the third fracture trace is determined, the fracture corresponding to the third fracture trace is divided into multiple scales based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data. The larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale of the fracture in the corresponding part is.

[0034] Specifically, the specific operation of dividing the crack scale is as follows: The third fracture trace is superimposed on the seismic attribute data to extract the seismic attribute intensity value at the location of the third fracture trace.

[0035] The seismic attribute intensity values ​​are graded. For example, the seismic attribute intensity values ​​can be divided into three levels. Each level corresponds to a threshold. The seismic attribute intensity values ​​are divided into three levels: weak, medium, and strong. The level division can be given according to the requirements of the visual fracture research.

[0036] The third crack trace is divided according to the level of the seismic attribute intensity value at the location of the third crack trace. The higher the level of the seismic attribute intensity value, the larger the scale of the corresponding crack.

[0037] For example, in this example, the cracks are divided into three levels: large-scale cracks, medium-scale cracks, and small-scale cracks according to the seismic attribute intensity values ​​from large to small. Figure 5 .

[0038] Then, the third crack trace after scale division is post-processed, specifically, each part of the third crack trace is characterized by thickness and / or color, and the thickness is used to represent cracks of different scales. The thicker the crack, the larger the scale. At the same time, cracks of different scales are displayed in different colors. Thus, the final fourth crack trace is obtained and output. The result is as follows: Figure 7 shown.

[0039] In some embodiments, the above-mentioned seismic attribute data may be plane seismic attribute data or slice data in three-dimensional seismic attribute data.

[0040] The crack extraction device provided by the present invention is described below. The crack extraction device described below and the crack extraction method described above can be referred to each other.

[0041] The crack extraction device provided by the embodiment of the present invention includes a data acquisition module, a first crack trace extraction module, a second crack trace generation module, a third crack trace generation module, a scale division module and a fourth crack trace output module.

[0042] Among them, the data acquisition module is used to obtain seismic attribute data and preset values ​​of fracture trace length.

[0043] The first crack trace extraction module is used to perform line-by-line peak scanning along at least two directions in the plane space of the seismic attribute data based on the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and generate a graph composed of the peak points as a first crack trace; The second crack trace generation module is used to enhance the continuity of the first crack trace by using a dilation-erosion method based on the first crack trace, and generate the first crack trace after the continuity enhancement as the second crack trace; The third crack trace generation module is used to generate a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length as a third crack trace based on the second crack trace and the preset value of the crack trace length; The scale division module is used to divide the scale of the cracks corresponding to the third crack trace based on the third crack trace and the seismic attribute intensity value in the seismic attribute data, wherein the larger the seismic attribute intensity value corresponding to the third crack trace is, the larger the scale of the corresponding part is; The fourth crack trace output module is used to perform thickness and / or color characterization processing on each part of the third crack trace based on the scale-divided third crack trace for different scales, to form and output a fourth crack trace.

[0044] Figure 8 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 8As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication interface 820 and the memory 830 communicate with each other through the communication bus 840. The processor 810 may call the logic instructions in the memory 830 to execute the crack extraction method, which includes: Step S100, obtaining seismic attribute data and a preset value of the fracture trace length; Step S200: Based on the seismic attribute data, perform row-by-row peak scanning in at least two directions on the plane space of the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and determine that a graph composed of the peak points is a first fracture trace; Step S300: Based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as the second crack trace; Step S400: Based on the second crack trace and the preset value of the crack trace length, a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length is determined as a third crack trace; Step S500: based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the fracture corresponding to the fracture trace is scaled, wherein the larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale of the corresponding part is; Step S600: Based on the scale-divided third crack trace, for different scales, perform thickness and / or color characterization processing on each part of the third crack trace to form a fourth crack trace and output it.

[0045] In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0046] On the other hand, the present invention further provides a computer program product, the computer program product comprising a computer program stored on a non-transitory computer-readable storage medium, the computer program comprising program instructions, when the program instructions are executed by a computer, the computer can execute the crack extraction method provided by the above methods, the method comprising: Step S100, obtaining seismic attribute data and a preset value of the fracture trace length; Step S200: Based on the seismic attribute data, perform row-by-row peak scanning in at least two directions on the plane space of the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and determine that a graph composed of the peak points is a first fracture trace; Step S300: Based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as the second crack trace; Step S400: Based on the second crack trace and the preset value of the crack trace length, a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length is determined as a third crack trace; Step S500: based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the fracture corresponding to the fracture trace is scaled, wherein the larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale of the corresponding part is; Step S600: Based on the scale-divided third crack trace, for different scales, perform thickness and / or color characterization processing on each part of the third crack trace to form a fourth crack trace and output it.

[0047] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the method for extracting cracks provided above is implemented, and the method comprises: Step S100, obtaining seismic attribute data and a preset value of the fracture trace length; Step S200: Based on the seismic attribute data, perform row-by-row peak scanning in at least two directions on the plane space of the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and determine that a graph composed of the peak points is a first fracture trace; Step S300: Based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as the second crack trace; Step S400: Based on the second crack trace and the preset value of the crack trace length, a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length is determined as a third crack trace; Step S500: based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the fracture corresponding to the fracture trace is scaled, wherein the larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale of the corresponding part is; Step S600: Based on the scale-divided third crack trace, for different scales, perform thickness and / or color characterization processing on each part of the third crack trace to form a fourth crack trace and output it.

[0048] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0049] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0050] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them; although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present application, and they should all be included in the scope of the claims and specification of the present application. In particular, as long as there is no structural conflict, the various technical features mentioned in the various embodiments can be combined in any way. The present application is not limited to the specific embodiments disclosed herein, but includes all technical solutions that fall within the scope of the claims.

Claims

1. A crack extraction method, characterized in that: include: Obtain seismic attribute data and preset values ​​of fracture trace length; Based on the seismic attribute data, performing row-by-row peak scanning along at least two directions on the plane space of the seismic attribute data to obtain peak points, and marking the peak points on the plane space of the seismic attribute data, and determining that a graph composed of the peak points is a first fracture trace; Based on the first crack trace, using a dilation-erosion method to enhance the continuity of the first crack trace, and determining the first crack trace after the enhanced continuity as a second crack trace; Based on the second crack trace and the preset value of the crack trace length, determining a portion of the second crack trace having a length greater than or equal to the preset value of the crack trace length as a third crack trace; Based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the third fracture trace is scaled, wherein the larger the seismic attribute intensity value corresponding to the third fracture trace is, the larger the scale is; Based on the third crack trace after scale division, thickness and / or color characterization processing is performed on each part of the third crack trace according to different scales to form a fourth crack trace and output it.

2. The crack extraction method according to claim 1, characterized in that: The seismic attribute data are plane seismic attribute data or slice data in three-dimensional seismic attribute data.

3. The crack extraction method according to claim 1, characterized in that: Performing line-by-line peak scanning along at least two directions in the plane space of the seismic attribute data to obtain peak points includes: Determine a first scanning direction, and scan the plane space of the seismic attribute data line by line along a direction perpendicular to the first scanning direction, each scan is along the first scanning direction, and each scan extracts the seismic attribute data corresponding to the scanned line, and records the seismic attribute data as V i , where V i represents the earthquake attribute data of the i-th row; Determine a second scanning direction, and scan the plane space of the seismic attribute data line by line along a direction perpendicular to the second scanning direction, each scan is along the second scanning direction, and each scan extracts the seismic attribute data corresponding to the scanned line, and the seismic attribute data is V j , where V j represents the earthquake attribute data of the jth row; For all V i and V j Perform a peak scan and identify all peak points.

4. The crack extraction method according to claim 3, characterized in that: The first scanning direction is orthogonal to the second scanning direction.

5. The crack extraction method according to claim 1, characterized in that: Based on the third fracture trace and the seismic attribute intensity value in the seismic attribute data, the fracture corresponding to the fracture trace is scaled, including: The seismic attribute intensity values ​​are graded, wherein the higher the seismic attribute intensity value is, the higher the corresponding grade is; Superimposing the third fracture trace with the seismic intensity attribute strength value, and extracting the seismic attribute strength value at each position on the third fracture trace; Based on the division level of the seismic attribute intensity value and the extracted seismic attribute intensity value of each position on the third fracture trace, each part of the third fracture trace is scaled. The higher the level of the seismic attribute intensity value corresponding to each position of the third fracture trace, the larger the scale of the corresponding position of the third fracture trace.

6. The crack extraction method according to claim 5, characterized in that: The third fracture trace is divided into at least a large-scale fracture trace, a medium-scale fracture trace and a small-scale fracture trace.

7. A crack extraction device, characterized in that: include: A data acquisition module, used to acquire seismic attribute data and a preset value of fracture trace length; A first crack trace extraction module is used to perform row-by-row peak scanning along at least two directions on the plane space of the seismic attribute data based on the seismic attribute data to obtain peak points, mark the peak points on the plane space of the seismic attribute data, and extract and determine a graph composed of the peak points as a first crack trace; A second crack trace generation module, configured to enhance the continuity of the first crack trace by using a dilation-erosion method based on the first crack trace, and generate the first crack trace after the enhanced continuity as a second crack trace; A third crack trace generating module, configured to generate, based on the second crack trace and the preset value of the crack trace length, a portion of the second crack trace whose length is greater than or equal to the preset value of the crack trace length as a third crack trace; A scale division module, for performing scale division on the cracks corresponding to the third crack trace based on the third crack trace and the seismic attribute intensity value in the seismic attribute data, wherein the larger the seismic attribute intensity value corresponding to the third crack trace is, the larger the scale of the corresponding part is; The fourth crack trace output module is used to perform thickness and / or color characterization processing on each part of the third crack trace according to different scales based on the third crack trace after scale division, to form and output a fourth crack trace.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the crack extraction method according to any one of claims 1 to 6 are implemented.

9. A computer program product, comprising a computer program stored on a non-transitory computer-readable storage medium, the computer program comprising program instructions, characterized in that: The program instructions are used to execute the steps of the crack extraction method according to any one of claims 1 to 6 through a computer.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the crack extraction method according to any one of claims 1 to 6 are implemented.

Images