A method, device and electronic device for identifying cracks in different directions based on Gabor operator

Abstract

The present invention provides a method, device, computer-readable storage medium, and electronic device for azimuthally identifying fractures based on the Gabor operator. The method comprises: obtaining a coherent attribute data volume from a three-dimensional post-stack seismic data volume of a target area; determining a time-difference data volume in the directions of the main survey line and the connecting survey line based on the three-dimensional post-stack seismic data volume of the target area; performing azimuthally processing on the coherent attribute data volume using the Gabor operator based on the time-difference data volume in the directions of the main survey line and the connecting survey line to obtain fracture development data volumes in different azimuths; and determining fracture development characteristics in different azimuths based on the fracture development data volumes in different azimuths. The present invention can obtain information on fracture development and fracture development intensity and their azimuth in different azimuths to predict fracture and fracture development zones and their directions.

Description

Technical Field

[0001] The present invention relates to the field of geophysical exploration technology, and in particular to a Gabor operator-based azimuth fracture identification method, device, computer-readable storage medium and electronic equipment. Background Art

[0002] Fractured oil and gas reservoirs account for one-third of China's proven reserves. Fractures serve as both pathways for oil and gas migration and reservoirs. Accurately identifying and predicting fractures and their development zones, as well as accurately understanding the occurrence and distribution of fractures in these reservoirs, are crucial for the efficient development of these reservoirs. The region, orientation, and density of fracture development within these reservoirs are crucial for exploration. Currently, there are two main types of techniques for seismic fracture prediction. One is fracture detection based on pre-stack azimuthal anisotropy, using pre-stack partial stack data or OVT gather data. This type of technique can extract fracture intensity and direction, but it requires high pre-stack data quality, is computationally intensive, and has a long operation cycle. This type of technique is currently the main development direction for fracture prediction. The other is fracture characterization based on geometric attributes of post-stack seismic data, such as traditional coherence, curvature, and seismic tensor techniques. These techniques, based on post-stack seismic data, offer small data volumes and fast computational speeds, making them widely used in production. However, post-stack seismic data lacks a concept of orientation, and fracture prediction techniques based on post-stack seismic data generally cannot obtain information on fracture direction. Although some researchers have used the multi-directional characteristics of the curvelet transform to study fracture direction in post-stack seismic data and have achieved some success in combining fracture direction prediction with wellbore data, the curvelet transform principle is complex and computationally inefficient. Therefore, it is urgent to combine the idea of ​​multi-directional prediction of post-stack curved waves to carry out research on more computationally efficient fracture development intensity and direction prediction technology to improve the efficiency of fracture prediction. Summary of the Invention

[0003] To address the above problems, embodiments of the present invention provide a method, device, computer-readable storage medium, and electronic device for identifying cracks in different orientations based on a Gabor operator.

[0004] In a first aspect, an embodiment of the present invention provides a method for identifying cracks in different orientations based on a Gabor operator, comprising the following steps:

[0005] S100, acquiring a three-dimensional post-stack seismic data volume of a target area;

[0006] S200, obtaining a coherent attribute data volume from the three-dimensional post-stack seismic data volume:

[0007] S300, based on the 3D post-stack seismic data volume of the target area, determining the time difference data volume in two directions of the main survey line and the connecting survey line;

[0008] S400: Based on the time difference data volumes in the two directions of the main survey line and the connecting survey line, the coherent attribute data volume is processed by the Gabor operator to obtain the fault development data volume in different directions;

[0009] S500: Based on the fracture development data volumes at different orientations, determine the fracture development characteristics at different orientations.

[0010] According to an embodiment of the present invention, the above method further includes the following steps:

[0011] S600, based on the fracture development data volumes at different orientations, determining the maximum value of the fracture development data volumes coherent in different directions and the minimum value of the fracture development data volumes coherent in different directions;

[0012] S700: Determine a fracture development intensity data volume based on a maximum value of fracture development data volumes coherent in different directions and a minimum value of fracture development data volumes coherent in different directions.

[0013] According to an embodiment of the present invention, the above method further includes the following steps:

[0014] S800, determines the fracture development direction by comparing fracture development data volumes in different directions;

[0015] S900 , determining a data volume of a fracture development direction from the coherent attribute data volume to predict a fracture development zone and its direction.

[0016] According to an embodiment of the present invention, the above method further includes the following steps:

[0017] S1000 predicts fault and fracture development zones and their trends based on the fracture development intensity data volume and the fracture development direction data volume.

[0018] According to an embodiment of the present invention, in the above step S400, the coherent attribute data volume is processed by using the Gabor operator according to the following formula:

[0019]

[0020] Where x'=xcosθ+ysinθ; y'=-xsinθ+ycosθ; θ is the angle of the Gabor operator; γ is the aspect ratio of the Gabor operator; λ is the wavelength; σ is half of the filter element, π is a natural constant, x represents the position in the direction of the main survey line, y represents the position in the direction of the contact survey line, Coh θ It is a data body of fault development in different directions.

[0021] According to an embodiment of the present invention, in the above step S700, the crack development intensity data volume is calculated according to the following formula:

[0022] FraInt=(Coh Max -Coh Min ) / Coh Max

[0023] Among them, Coh Max is the maximum value of coherence in different directions, Coh Min is the minimum value of coherence in different directions, and FraInt is the fracture development intensity data volume.

[0024] According to an embodiment of the present invention, in the above step S800, the fracture development data volumes in different directions are compared, wherein the direction with the smallest fracture development data volume is the crack development direction.

[0025] According to a second aspect, the present invention further provides a device for implementing the method for identifying cracks in different orientations based on a Gabor operator as described in the first aspect, comprising:

[0026] A data acquisition module is used to acquire three-dimensional post-stack seismic data volumes in the target area;

[0027] A data extraction module is used to obtain a coherent attribute data volume from the three-dimensional post-stack seismic data volume:

[0028] A data determination module is used to determine the time difference data volume in two directions of the main survey line and the connecting survey line based on the three-dimensional post-stack seismic data volume of the target area;

[0029] Axis processing of a block is used to process the time difference data volume in the two directions of the main survey line and the connecting survey line. The Gabor operator is used to process the coherent attribute data volume in different directions to obtain the fault development data volume in different directions.

[0030] The feature analysis module is used to determine the fracture development characteristics in different orientations based on the fracture development data volumes in different orientations.

[0031] In a third aspect, an embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon. When the program is executed by a processor, the method for identifying cracks in different orientations based on the Gabor operator as described in the first aspect is implemented.

[0032] In a fourth aspect, an embodiment of the present invention provides an electronic device, comprising:

[0033] processor;

[0034] a memory for storing instructions executable by the processor;

[0035] The processor is configured to execute the instructions to implement the orientation-based crack identification method based on Gabor operator as described in the first aspect.

[0036] Compared with the prior art, the above technical solution of the present invention has the following beneficial effects:

[0037] An embodiment of the present invention proposes a method for identifying fractures in different orientations based on the Gabor operator. The method is based on post-stack seismic data and uses the Gabor operator in different orientations to extract fracture development characteristics in different orientations. Then, according to the fracture development characteristics in different orientations, faults and fracture development zones and their directions are predicted to improve the efficiency of fracture identification. An embodiment of the present invention proposes a method for identifying fractures in different orientations based on the Gabor operator. This method can analyze the fracture development characteristics in each orientation separately, or it can comprehensively analyze the fracture development characteristics in different directions to extract the fracture development intensity and main development direction to predict the fractures and fracture development zones and their directions. Compared with the traditional fracture prediction effect based on post-stack seismic data, the method for identifying fractures in different orientations based on the Gabor operator proposed in an embodiment of the present invention has higher accuracy, richer parameters, and higher computational efficiency compared with pre-stack data. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0039] Figure 1 The flowchart of the Gabor operator-based orientation crack identification method according to an embodiment of the present invention is shown;

[0040] Figure 2 (a) to (f) are schematic diagrams showing coherent slices at different orientations obtained by using the orientation-specific Gabor operator according to an embodiment of the present invention;

[0041] Figure 3 A schematic diagram showing the structure of an electronic device according to an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0042] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments 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 making creative efforts shall fall within the scope of protection of the present invention.

[0043] Example 1

[0044] The present invention provides an azimuth-based fracture identification method based on Gabor operator. The technology is based on post-stack seismic data and uses azimuth-based Gabor operator to extract fracture development characteristics in different orientations. According to the fracture development characteristics in different orientations, faults and fracture development zones and their directions are predicted.

[0045] like Figure 1 The specific implementation steps of the method include but are not limited to the following steps:

[0046] Step 1: Extract the coherent attribute volume Coh.

[0047] Input the 3D post-stack seismic data volume Cube and extract the coherent attribute data volume Coh of the 3D data volume.

[0048] Step 2: Extract the azimuth time difference data volumes DTInline and DTXline.

[0049] For the 3D post-stack seismic data volume Cube, calculate the time difference data volumes DTInline and DTXline in the two directions of the main survey line and the connecting survey line.

[0050] Step 3: Based on the time difference data volume DTInline and DTXline, the coherent attribute data volume Coh is processed by Gabor operator to obtain the fracture development data volume Coh at different directions. θ , according to the fracture development data body Coh in different directions θ It is possible to obtain the fault development characteristics in different directions.

[0051] Using the azimuth Gabor operator of formula (1), the coherent attribute data volume Coh is processed in azimuth by combining the time difference of DTInline and DTXline in two directions to obtain the fracture development data volume Coh in different azimuths. θ .

[0052]

[0053] Where x' = xcosθ + ysinθ; y' = -xsinθ + ycosθ; θ is the angle of the Gabor operator; γ is the aspect ratio of the Gabor operator; λ is the wavelength; σ is half the filter bin; π is a natural constant; x and y are the corresponding positions in the direction of the main survey line and the direction of the connecting survey line, respectively, and the values ​​of the corresponding positions are determined by DTInline and DTXline.

[0054] Step 4: Based on the fracture development data volumes in different directions, determine the maximum value of the fracture development data volumes with different directions of coherence and the minimum value of the fracture development data volumes with different directions of coherence, and then calculate the fracture development intensity data volume FraInt using formula (2):

[0055] FraInt=(Coh Max -Coh Min ) / Coh Max (2)

[0056] Among them, Coh Max is the maximum value of the fracture development data volume with different directions of coherence, Coh Min is the minimum value of fracture development data with coherence in different directions.

[0057] Step 5: Compare the fracture development data in different directions Coh Min To determine the direction of fracture development, the direction with the smallest fracture development data volume can be regarded as the main direction of fracture development. After determining the main direction of fracture development, the data volume FraOri of fracture development direction is found from the coherent attribute data volume Coh, and the fracture and fracture development zones and their directions are predicted based on this.

[0058] Example 2

[0059] The implementation process and technical effects of the technical solution of the present invention are described below with reference to an embodiment.

[0060] In this embodiment, a method for identifying cracks in different orientations based on Gabor operators can be implemented according to the following steps:

[0061] First, perform step 1:

[0062] The coherent attribute volume Coh is extracted from the 3D post-stack seismic data volume Cube in the target area.

[0063] In actual operation, the 3D post-stack seismic data volume Cube is input and the coherent attribute data volume Coh of the 3D data volume is extracted.

[0064] In this embodiment, coherent attribute extraction from post-stack seismic data is a traditional calculation method. First-generation three-channel cross-correlation-based coherence algorithms, second-generation similarity-based cross-correlation algorithms, and third-generation eigenvalue-based coherence algorithms can be used to extract coherent attributes. Coherent attribute extraction involves calculating bins and time windows. These three types of coherence algorithms are well-established technologies, and common commercial geophysical software includes coherent attribute extraction functionality.

[0065] Then, proceed to step 2:

[0066] Extract the azimuth time difference data volumes DTInline and DTXline from the 3D post-stack seismic data volume Cube in the target area.

[0067] In actual operation, for the 3D post-stack seismic data volume Cube, the time difference data volumes DTInline and DTXline in the two directions of the main survey line and the connecting survey line are calculated.

[0068] In this embodiment, a cross-correlation algorithm can be used to calculate two-way time difference data volumes along the main survey line and the connecting survey line. The cross-correlation time window and maximum scan length must be set. The cross-correlation algorithm is used to calculate the strike time difference of the event axis for adjacent traces. The time difference corresponding to the maximum similarity coefficient is the directional time difference, represented by the number of sample points. The time difference data in both directions will serve as input in step three to determine the extracted data for calculating the azimuth-specific fault characteristics.

[0069] Then, proceed to step three:

[0070] Extract the fracture development characteristics under different orientations Coh θ .

[0071] Using the azimuth Gabor operator of formula (1), the coherent attribute data volume is processed in azimuth by combining the time difference of DTInline and DTXline in two directions to obtain the fracture development data volume Coh in different directions. θ .

[0072]

[0073] Wherein, x'=xcosθ+ysinθ; y'=-xsinθ+ycosθ; θ is the angle of the Gabor operator; γ is the aspect ratio of the Gabor operator; λ is the wavelength; σ is half of the filter element, and π is a natural constant.

[0074] x and y represent the location of the seismic data. Their values ​​are determined by the DTInline and DTXline time difference data volumes. x represents the location along the main survey line, and y represents the location along the tie survey line. If the DTInline value at that location is 0, x takes the value at x; if the DTInline value at that location is 1, x takes the value at x+1; if the DTInline value at that location is -1, x takes the value at x-1. Similarly, tie survey lines are similar.

[0075] Then, proceed to step 4:

[0076] Based on the fracture development data volumes in different directions, the maximum value of the fracture development data volumes with different directions of coherence and the minimum value of the fracture development data volumes with different directions of coherence are determined, and then the fracture development intensity data volume FraInt is calculated using formula (2):

[0077] FraInt=(Coh Max -Coh Min ) / Coh Max (2)

[0078] Among them, Coh Max is the maximum value of coherence in different directions, Coh Min is the minimum value of coherence in different directions.

[0079] Then, execute step five:

[0080] By comparing the fracture development data in different directions, Coh Min To determine the direction of fracture development, the direction with the smallest fracture development data volume can be regarded as the main direction of fracture development. After determining the main direction of fracture development, the data volume FraOri of fracture development direction is found from the coherent attribute data volume Coh, and the fracture and fracture development zones and their directions are predicted based on this.

[0081] Application Examples

[0082] The beneficial effects of the embodiments of the present invention are described below through a specific example.

[0083] Figure 2 (a) to Figure 2 (f) shows the fracture development characteristics in different orientations extracted using the orientation-specific Gabor operator, where:

[0084] (a) is the extracted crack development characteristics at 0 degrees;

[0085] (b) is the extracted crack development characteristics at 30 degrees;

[0086] (c) is the extracted crack development characteristics at 60 degrees;

[0087] (d) is the extracted 90-degree crack development characteristics;

[0088] (e) is the extracted crack development characteristics at 120 degrees;

[0089] (f) is the extracted crack development characteristics at 150 degrees.

[0090] This comparison reveals the primary direction and corresponding intensity of fracture development at each location, allowing for the identification of additional secondary or sub-fractures. Fracture development directions extracted from seismic data are relative, requiring a comprehensive correction of the seismic predictions based on the fracture development direction in the wellbore to determine the fracture development direction for the entire work area.

[0091] In summary, the present invention provides an azimuth-based fracture identification method based on Gabor operator, which is based on post-stack seismic data and uses azimuth-based Gabor operator to extract fracture development characteristics in different orientations. That is, the fracture development characteristics in each orientation can be analyzed separately, or the fracture development intensity and main development direction can be extracted by comprehensively analyzing the fracture development characteristics in different directions to predict the fault and fracture development zone and its trend. Compared with the traditional fracture prediction method based on post-stack seismic data, the method has higher accuracy, richer parameters and higher calculation efficiency than pre-stack data.

[0092] Example 3

[0093] The following are embodiments of the apparatus of the present invention, which can be used to implement the method embodiments of the present invention. For details not disclosed in the apparatus embodiments of the present invention, please refer to the method embodiments of the present invention.

[0094] This embodiment provides a device for identifying cracks in different orientations based on a Gabor operator, which is characterized by comprising:

[0095] An acquisition module, used to acquire a three-dimensional seismic data volume of a target area;

[0096] a decomposition module for decomposing each seismic signal in the three-dimensional seismic data volume into a data set consisting of a plurality of IMF components, extracting the instantaneous amplitude and instantaneous frequency of each of the IMF components, calculating the instantaneous energy-frequency value of each of the IMF components based on the instantaneous amplitude and instantaneous frequency, and calculating the instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency value of each of the IMF components,

[0097] An analysis module is used to analyze the energy and frequency variation characteristics of the target area formation based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy-frequency value of each IMF component of each seismic signal, and to analyze the oil and gas content of the formation according to the energy and frequency variation characteristics.

[0098] Example 4

[0099] This embodiment provides a computer-readable medium having a computer program stored thereon. When the program is executed by a processor, the steps of the azimuth-based crack identification method based on the Gabor operator as described in the above embodiment are implemented.

[0100] It should be noted that the present invention can implement all or part of the processes in the above-mentioned embodiment method by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, it can implement the steps of the above-mentioned various method embodiments. The computer program includes computer program code, which can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal and software distribution medium. Of course, there are other types of readable storage media, such as quantum memory, graphene memory, etc. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

[0101] Example 5

[0102] Figure 3 FIG. 1 is a schematic diagram of the structure of an electronic device according to an embodiment of the present invention. Figure 3 As shown, at the hardware level, the electronic device includes a processor and, optionally, an internal bus, a network interface, and memory. The memory may include internal memory, such as high-speed random-access memory (RAM), and may also include non-volatile memory, such as at least one disk storage device. Of course, the electronic device may also include other hardware required for its services.

[0103] The processor, network interface, and memory can be interconnected via an internal bus, such as an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, or an EISA (Extended Industry Standard Architecture) bus. These buses can be categorized as address buses, data buses, and control buses. For ease of illustration, the diagram uses only line segments, but this does not imply that there is only one bus or only one type of bus.

[0104] The memory is configured to store a program. Specifically, the program may include program code, which includes computer operating instructions. The memory may include internal memory and non-volatile memory, and provides instructions and data to the processor. The processor reads the corresponding computer program from the non-volatile memory into the internal memory and then runs it. The processor executes the program stored in the memory to perform all steps of the aforementioned method for azimuth-based crack identification based on Gabor operators.

[0105] The communication bus mentioned in the above devices may be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. This communication bus can be divided into address buses, data buses, control buses, etc. For ease of illustration, only one thick line is used in the figure, but this does not mean that there is only one bus or only one type of bus. The communication interface is used for communication between the above electronic devices and other devices.

[0106] Bus comprises hardware, software or both, for above-mentioned parts are coupled together.For example, bus can comprise accelerated graphics port (AGP) or other graphics bus, enhanced industry standard architecture (EISA) bus, front side bus (FSB), hypertransport (HT) interconnection, industry standard architecture (ISA) bus, infinite bandwidth interconnection, low pin count (LPC) bus, memory bus, micro channel architecture (MCA) bus, peripheral component interconnection (PCI) bus, PCI-Express (PCI-X) bus, serial advanced technology attachment (SATA) bus, video electronics standard association local (VLB) bus or other suitable bus or two or more above these combinations.In suitable case, bus can comprise one or more buses.Although the embodiment of the present invention describes and shows specific bus, the present invention considers any suitable bus or interconnection.

[0107] The memory may include random access memory (RAM) or non-volatile memory (NVM), such as at least one disk storage. Alternatively, the memory may be at least one storage device located away from the processor.

[0108] The memory may include a large capacity memory for data or instructions. By way of example and not limitation, the memory may include a hard disk drive (HDD), a floppy disk drive, a flash memory, an optical disk, a magneto-optical disk, a magnetic tape, or a universal serial bus (USB) drive, or a combination of two or more of these. Where appropriate, the memory may include a removable or non-removable (or fixed) medium. In a specific embodiment, the memory is a non-volatile solid-state memory. In a specific embodiment, the memory includes a read-only memory (ROM). Where appropriate, the ROM may be a mask-programmed ROM, a programmable ROM (PROM), an erasable PROM (EPROM), an electrically erasable PROM (EEPROM), an electrically rewritable ROM (EAROM), or a flash memory, or a combination of two or more of these.

[0109] The above-mentioned processor can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, and discrete hardware components.

[0110] It should be noted that those skilled in the art can clearly understand that for the convenience and brevity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of the present invention. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, and will not be repeated here.

[0111] The devices, apparatuses, systems, modules, or units described in the above embodiments may be implemented by computer chips or entities, or by products having certain functions. A typical implementation device is a computer. Specifically, the computer may be, for example, a personal computer, a laptop computer, an in-vehicle human-computer interaction device, a cellular phone, a camera phone, a smartphone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

[0112] Although the present invention provides method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive means. The order of steps listed in the embodiments is only one way of executing the order of many steps and does not represent the only execution order. When an actual device or terminal product is executed, it can be executed in sequence or in parallel according to the method shown in the embodiments or the drawings (for example, in a parallel processor or multi-threaded processing environment, or even a distributed data processing environment).

[0113] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0114] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0115] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0116] It should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or apparatus comprising the element.

[0117] Each embodiment in this specification is described in a related manner. Similar portions between the various embodiments can be referenced to each other. Each embodiment focuses on the differences from other embodiments. In particular, the device, electronic device, and readable storage medium embodiments are generally similar to the method embodiments, so their descriptions are relatively simple. For related portions, reference can be made to the descriptions of the method embodiments.

[0118] The above description is only a preferred embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention are included in the scope of protection of the present invention.

Claims

1. A method for identifying cracks in different directions based on Gabor operator, characterized in that: The following steps are involved: S100, acquiring a three-dimensional post-stack seismic data volume of a target area; S200, obtaining a coherent attribute data volume from the three-dimensional post-stack seismic data volume: S300, based on the 3D post-stack seismic data volume of the target area, determining the time difference data volume in two directions of the main survey line and the connecting survey line; S400, based on the time difference data volumes in the two directions of the main survey line and the connecting survey line, the coherent attribute data volume is processed by using the Gabor operator according to the following formula to obtain the fault development data volume in different directions; Where x'=xcosθ+ysinθ; y'=-xsinθ+ycosθ; θ is the angle of the Gabor operator; γ is the aspect ratio of the Gabor operator; λ is the wavelength; σ is half of the filter bin, π is a natural constant, x and y represent the position of the seismic data, and their values ​​are determined by the time difference data volume in two directions. x represents the position in the direction of the main survey line, and y represents the position in the direction of the connecting survey line. Coh θ is the fracture development data volume under different orientations; x' is the position of the seismic data in the detection direction of the filter, and y' is the position of the seismic data in the orthogonal direction of the detection direction; S500: Based on the fracture development data volumes at different orientations, determine the fracture development characteristics at different orientations.

2. The method for identifying cracks in different orientations based on Gabor operators according to claim 1, characterized in that: The method further comprises: S600, based on the fracture development data volumes at different orientations, determining the maximum value of the fracture development data volumes coherent in different directions and the minimum value of the fracture development data volumes coherent in different directions; S700: Determine a fracture development intensity data volume based on a maximum value of fracture development data volumes coherent in different directions and a minimum value of fracture development data volumes coherent in different directions.

3. The method for identifying cracks in different orientations based on Gabor operators according to claim 2, characterized in that: The method further comprises: S800, determines the fracture development direction by comparing fracture development data volumes in different directions; S900: Determine a data volume of a crack development direction from the coherent attribute data volume.

4. The method for identifying cracks in different orientations based on Gabor operators according to claim 3, wherein: The method further comprises: S1000 predicts fault and fracture development zones and their trends based on the fracture development intensity data volume and the fracture development direction data volume.

5. The method for identifying cracks in different orientations based on Gabor operators according to claim 2, wherein: In step S700, the crack development intensity data volume is calculated according to the following formula: FraInt=(Coh Max - Coh Min ) / Co Max Among them, Coh Max is the maximum value of the fracture development data volume with different directions of coherence, Coh Min is the minimum value of the fracture development data body with coherence in different directions, and FraInt is the fracture development intensity data body.

6. The method for identifying cracks in different orientations based on Gabor operators according to claim 3, wherein: In step S800 , the fracture development data volumes in different directions are compared, wherein the direction with the smallest fracture development data volume is the fracture development direction.

7. A device for identifying cracks in different directions based on Gabor operator, characterized in that: include: A data acquisition module is used to acquire three-dimensional post-stack seismic data volumes in the target area; A data extraction module is used to obtain a coherent attribute data volume from the three-dimensional post-stack seismic data volume: A data determination module is used to determine the time difference data volume in two directions of the main survey line and the connecting survey line based on the three-dimensional post-stack seismic data volume of the target area; A block is processed in different directions. It is used for the time difference data volume based on the main survey line and the connecting survey line. According to the following formula, the coherent attribute data volume is processed in different directions using the Gabor operator to obtain the fault development data volume in different directions. Where x'=xcosθ+ysinθ; y'=-xsinθ+ycosθ; θ is the angle of the Gabor operator; γ is the aspect ratio of the Gabor operator; λ is the wavelength; σ is half of the filter bin, π is a natural constant, x and y represent the position of the seismic data, and their values ​​are determined by the time difference data volume in two directions. x represents the position in the direction of the main survey line, and y represents the position in the direction of the connecting survey line. Coh θ is the fracture development data volume under different orientations; x' is the position of the seismic data in the detection direction of the filter, and y' is the position of the seismic data in the orthogonal direction of the detection direction; The feature analysis module is used to determine the fracture development characteristics in different orientations based on the fracture development data volumes in different orientations.

8. A computer-readable storage medium, characterized in that A computer program is stored thereon, and when the program is executed by a processor, the method for identifying cracks in different orientations based on a Gabor operator as claimed in any one of claims 1 to 6 is implemented.

9. An electronic device comprising: processor; a memory for storing instructions executable by the processor; The processor is configured to execute the instructions to implement the orientation-based crack identification method based on Gabor operator according to any one of claims 1 to 6.

Images