A five-dimensional trace gather discontinuity body carving method and device and electronic equipment

Abstract

Embodiments of the present application provide a five-dimensional fracture and fissure body carving method and device and electronic equipment. The five-dimensional fracture and fissure body carving method comprises: obtaining fracture attribute and fissure plane distribution attribute of a target work area from five-dimensional fracture and fissure data of the target work area according to the maximum principal stress direction of the target work area; using the correlation between seismic attribute and physical property of the target work area, calculating the fracture and fissure attribute volume of the target layer according to the seismic attribute of the target layer extracted from full stack data of the target work area; fusing the fracture attribute, the fissure plane distribution attribute and the fracture and fissure attribute volume to obtain a data volume for reflecting the spatial distribution of the fracture and fissure body; and obtaining the spatial distribution of the fracture and fissure body of the target layer according to the data volume.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of geophysical exploration, and in particular to a five-dimensional gather-based fault-fracture body carving method and device and electronic equipment. BACKGROUND

[0002] Fault-fracture body is a new type of trap related to faults and affected by the development degree of fractures in reservoirs. Carbonate rocks are prone to tectonic deformation and generate fractures and dissolution pores under tectonic stress, especially near faults.

[0003] In related technologies, a single attribute method is used to depict faults and fractures, for example, a fault and fracture can be depicted based on a full stack data volume. However, due to the anisotropy of fault-fracture bodies, it is difficult to depict the fault-fracture body and the distribution of effective reservoir space by using a single attribute method.

[0004] Therefore, how to accurately depict the fault-fracture body and the distribution of effective reservoir space is a technical problem to be solved. SUMMARY

[0005] The embodiments of the present application provide a five-dimensional gather-based fault-fracture body carving method, device, electronic equipment and storage medium, which are used to accurately depict the fault-fracture body and the distribution of effective reservoir space.

[0006] One of the embodiments of the present application provides a five-dimensional gather-based fault-fracture body carving method, which comprises the following steps: determining the wellbore peripheral stress field direction of a drilled well in a target work area; determining the maximum principal stress direction of the target work area according to the coherent plane map of a target layer in the target work area and the wellbore peripheral stress field direction of the drilled well in the target work area; wherein the coherent plane map of the target layer is obtained by extracting seismic coherent attributes from full stack seismic data of the target work area; obtaining the fault attribute and the fracture plane distribution attribute of the target work area from five-dimensional gather data of the target work area according to the maximum principal stress direction of the target work area; determining the correlation between seismic attributes and physical properties of the target work area according to the physical property information of the drilled well in the target work area and the full stack seismic data of the target work area; wherein the seismic attributes at least include the amplitude of seismic waves, and the physical properties at least include fracture data and dissolution pore data; calculating the fault-fracture attribute volume of the target layer according to the seismic attributes of the target layer extracted from the full stack data of the target work area by using the correlation between the seismic attributes and the physical properties of the target work area; wherein the fault-fracture attribute volume is used to reflect the physical properties of the fault-fracture body of the target layer; fusing the fault attribute, the fracture plane distribution attribute and the fault-fracture attribute volume to obtain a data volume for reflecting the spatial distribution of the fault-fracture body; and obtaining the spatial distribution of the fault-fracture body of the target layer according to the data volume.

[0007] In some embodiments, the determining the direction of the wellbore peripheral stress field of the drilled well in the target work area comprises:

[0008] The azimuth of the induced fracture of the drilled well is determined by using the imaging logging data of the drilled well in the target work area; and the direction of the wellbore peripheral stress field of the drilled well in the target work area is determined according to the azimuth of the induced fracture.

[0009] In some embodiments, the coherent plane map of the target layer is obtained by: determining the seismic reflection interface of the target layer by using the composite record calibration result of the drilled well in the target work area; obtaining the horizon data of the target layer from the full stack seismic data of the target work area according to the seismic reflection interface of the target layer; and extracting the coherent attribute of the target layer from the full stack seismic data according to the horizon data to obtain the coherent plane map of the target layer.

[0010] In some embodiments, the fracture attribute and the fracture plane distribution attribute of the target work area are obtained from the five-dimensional gather data of the target work area according to the direction of the maximum principal stress of the target work area, comprising: stacking the five-dimensional gather data according to the division azimuth to obtain the data volume of the angle division stacking; selecting the data volume perpendicular to the direction of the maximum principal stress of the target work area from the data volume of the angle division stacking to obtain the data volume for the fracture and fracture body attribute analysis; and extracting the coherent attribute from the data volume for the fracture and fracture body attribute analysis to obtain the fracture attribute and the fracture plane distribution attribute.

[0011] One embodiment of this application provides a five-dimensional gather-based fracture carving device. The device includes: a first determining module for determining the stress field direction around the drilled wellbore in the target area; a second determining module for determining the direction of the maximum principal stress in the target area based on the coherence plane map of the target layer and the stress field direction around the drilled wellbore; wherein the coherence plane map of the target layer is obtained by extracting seismic coherence attributes from the full-stack seismic data of the target area; a first acquiring module for obtaining the fracture attributes and fracture plane distribution attributes of the target area from the five-dimensional gather data of the target area based on the direction of the maximum principal stress; and a second acquiring module for obtaining the fracture attributes and fracture plane distribution attributes of the target area based on the physical property information of the drilled well and the stress field direction around the drilled wellbore in the target area. The following steps are taken: A first module overlays seismic data to determine the correlation between seismic attributes and physical properties of the target work area; wherein the seismic attributes include at least the amplitude of seismic waves, and the physical properties include at least fracture data and dissolution pore data; a second module uses the correlation between seismic attributes and physical properties of the target work area to calculate the fracture attribute volume of the target layer based on the seismic attributes of the target layer extracted from the full overlay data of the target work area; wherein the fracture attribute volume reflects the fracture physical properties of the target layer; a third module fuses the fracture attributes, the plane distribution attributes of the fractures, and the fracture attribute volume to obtain a data volume reflecting the spatial distribution of the fractures; a fourth module obtains the spatial distribution of the fractures in the target layer based on the data volume.

[0012] This application provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the method described above when running the program.

[0013] This application provides a storage medium for storing a computer-readable program, which, when run, performs the method described above.

[0014] The technical solutions provided in this application have at least the following advantages compared with the prior art:

[0015] In the implementation provided in this application, based on the direction of the maximum principal stress in the target area, the fracture properties and fracture planar distribution properties of the target area are obtained from the five-dimensional gather data of the target area. Utilizing the correlation between the seismic properties and physical properties of the target area, the fracture property volume of the target layer is calculated based on the seismic properties of the target layer extracted from the full stack data of the target area. The fracture properties, fracture planar distribution properties, and fracture property volume are then fused to obtain a data volume reflecting the spatial distribution of the fracture body. Based on this data volume, the spatial distribution of the fracture body in the target layer is obtained. This allows for accurate characterization of the fracture body and the effective reservoir spatial distribution. Attached Figure Description

[0016] The present application will be further described in the way of example embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not limiting, in these embodiments, the same numbers denote the same structures, wherein:

[0017] Figure 1 is an example flow chart of a five-dimensional gather fracture body sculpting method according to some embodiments of the present application;

[0018] Figure 2 is an example schematic diagram of a A3 well logging induced fracture map according to some embodiments of the present application;

[0019] Figure 3 is an example schematic diagram of a coherence flat of a target zone according to some embodiments of the present application;

[0020] Figure 4 is an example schematic diagram of a scheme of azimuthal stack according to some embodiments of the present application;

[0021] Figure 5 is an example schematic diagram of a root mean square amplitude attribute map extracted based on full stack data according to some embodiments of the present application;

[0022] Figure 6 is an example schematic diagram of a scatter plot of root mean square amplitude attribute versus porosity of dissolution pores according to some embodiments of the present application;

[0023] Figure 7 is an example schematic diagram of a fracture body plan map according to some embodiments of the present application;

[0024] Figure 8 is an example schematic diagram of a fracture body map on a seismic profile according to some embodiments of the present application;

[0025] Figure 9 is an example schematic diagram of a spatial distribution map of a fracture body according to some embodiments of the present application;

[0026] Figure 10 is an example schematic diagram of a five-dimensional gather fracture body sculpting device according to some embodiments of the present application;

[0027] Figure 11 is an example structural schematic diagram of an electronic device according to some embodiments of the present application. DETAILED DESCRIPTION

[0028] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some examples or embodiments of the present application, and for those skilled in the art, the present application can also be applied to other similar scenarios without creative labor on the basis of these drawings. Unless it is obvious from the language environment or otherwise stated, the same reference numbers in the drawings represent the same structure or operation.

[0029] It should be understood that the "system", "device", "unit" and / or "module" used herein is a method for distinguishing different components, elements, parts, sections or assemblies at different levels. However, if other words can achieve the same purpose, the words can be replaced by other expressions.

[0030] As shown in the present application and claims, unless the context clearly indicates otherwise, the words "one", "a", "an", and / or "the" do not refer to the singular, but can also include the plural. Generally speaking, the terms "comprise" and "include" only indicate the inclusion of the steps and elements explicitly identified, and these steps and elements do not constitute an exclusive list, and the method or device can also include other steps or elements.

[0031] Flowcharts are used in the present application to illustrate the operations performed by the system according to the embodiments of the present application. It should be understood that the preceding or subsequent operations are not necessarily performed in sequence. On the contrary, each step can be processed in reverse order or simultaneously. At the same time, other operations can be added to these processes, or one or more steps of operation can be removed from these processes.

[0032] In order to facilitate understanding, the technical solutions of the present application will be introduced below in combination with the drawings and embodiments.

[0033] Figure 1 is an exemplary flowchart of a five-dimensional fracture set discontinuity body carving method according to some embodiments of the present application. As shown in Figure 1 The five-dimensional fracture set discontinuity body carving method includes the following steps:

[0034] In step S110, the wellbore peripheral stress field direction of the drilled well of the target work area is determined.

[0035] In the specific implementation process, the azimuth of the induced fracture of the drilled well can be determined by using the imaging logging data of the drilled well of the target work area; and the wellbore peripheral stress field direction of the target work area is determined according to the azimuth of the induced fracture.

[0036] For example only, as Figure 2As shown, the induced fracture trend of the A3 well in the target work area is mainly northwest, and according to the relationship between the induced fracture and the wellbore peripheral principal stress, it can be determined that the stress field direction around the A3 well is northeast, so the wellbore peripheral stress field direction of the drilled well in the target work area is northeast.

[0037] In step S120, the maximum principal stress direction of the target work area is determined according to the coherence panel of the target layer in the target work area and the wellbore peripheral stress field direction of the drilled well; wherein the coherence panel of the target layer is obtained from the seismic coherence attribute extracted from the full stack seismic data of the target work area.

[0038] In the specific implementation process, the coherence panel of the target layer can be obtained by the following method:

[0039] The seismic reflection interface of the target layer is determined by using the synthetic record calibration result of the drilled well in the target work area; the horizon data of the target layer is interpreted from the full stack seismic data of the target work area according to the seismic reflection interface of the target layer; the coherence attribute of the target layer is extracted from the full stack seismic data according to the horizon data, and the coherence panel of the target layer is obtained, as shown in the coherence panel of the target layer. Figure 3

[0040] According to the Figure 3 It can be seen that the main fracture and the fracture trend of the target work area are northwest, that is, the target work area is affected by the northeast stress field, which is consistent with the wellbore peripheral stress field direction determined in step S120, so it can be determined that the maximum principal stress direction of the target work area is northeast.

[0041] In step S130, the fracture attribute and the fracture plane distribution attribute of the target work area are obtained from the five-dimensional gather data of the target work area according to the maximum principal stress direction of the target work area.

[0042] In some embodiments, the five-dimensional gather data is stacked according to the division azimuth angle to obtain an angle-stacked data volume; the data volume perpendicular to the maximum principal stress direction of the target work area is selected from the angle-stacked data volume to obtain a data volume for fracture and fracture body attribute analysis; the coherence attribute is extracted from the data volume for fracture and fracture body attribute analysis to obtain the fracture attribute and the fracture plane distribution attribute.

[0043] In the specific implementation process, the angle-stacked data volume can be obtained by stacking according to the signal-to-noise ratio of the five-dimensional gather data of the target work area, the fracture azimuth trend distribution, etc.

[0044] For example only, according to the maximum principal stress direction of the target work area being northeast and the fracture and fracture trend being mainly northwest, it is determined that Figure 4 ​The azimuthal superposition scheme is shown: the five-dimensional gathers of the target work area are superposed according to the mapping relationship of -22.5°~22.5°, 22.5°~67.5°, -67.5~-22.5°, 67.5~112.5° corresponding to 0°, 45°, -45°, 90° respectively, to obtain the azimuthal superposition data volume. The data volume perpendicular to the maximum principal stress direction of the target work area (0° azimuth data) is selected from the azimuthal superposition data volume as the data volume for fracture and crack attribute analysis; the coherence volume attribute is extracted from the data volume for fracture and crack attribute analysis to obtain the fracture attribute and crack plane distribution attribute of the target work area.

[0045] In step S140, the correlation between the seismic attribute and the physical property of the target work area is determined according to the physical property information of the drilled well in the target work area and the full stack seismic data of the target work area; wherein the seismic attribute at least includes the amplitude of the seismic wave, and the physical property at least includes the fracture data and the dissolution pore data.

[0046] For example only, on the basis of the calibration of the synthetic record of the drilled well, the response characteristics of the fracture and dissolution pore at the target layer of the drilled well in the full stack seismic data are analyzed, and it is found that several wells with good oil and gas shows in the target layer have medium-low frequency and medium-strong amplitude characteristics in the full stack seismic data, and the corresponding physical property is also relatively good. The root mean square amplitude extracted from the full stack seismic data has a good corresponding relationship with the physical property, and the root mean square amplitude attribute value corresponding to the well with good physical property is also high, such as Figure 5 As shown, the two are positively correlated.

[0047] According to the above analysis results, the root mean square amplitude value of the seismic data corresponding to the target layer of the drilled well in the work area and the porosity of the dissolution pore are counted respectively, and the correlation between the two is established as shown in Table 1 and Figure 6 As shown, the correlation between the two is established as shown in Table 1 and

[0048] y=7*e -6x +0.7338 (1)

[0049] Wherein, x is the root mean square amplitude value, y is the porosity size of the dissolution pore, unit %.

[0050] Table 1

[0051]

[0052] In step S150, the fracture attribute volume of the target layer is calculated according to the seismic attribute of the target layer extracted from the full stack data of the target work area by using the correlation between the seismic attribute and the physical property of the target work area; wherein the fracture attribute volume is used to reflect the physical property of the fracture volume of the target layer.

[0053] For example, the root-mean-square amplitude value of the seismic data corresponding to the target layer can be extracted from the full stack data of the target work area, and the root-mean-square amplitude value is substituted into formula (1) to calculate the fracture attribute volume representing the size of the porosity of the dissolution pore of the target layer.

[0054] In step S160, the fracture attribute, the fracture plane distribution attribute and the fracture attribute volume are fused to obtain a data volume for reflecting the spatial distribution of the fracture volume.

[0055] In the implementation process, the fracture attribute and the fracture plane distribution attribute of the target work area obtained in step S130 can be fused with the fracture attribute volume of the target layer reflecting the reservoir property obtained in step S150 to obtain a data volume for reflecting the spatial distribution of the fracture volume.

[0056] In step S170, the spatial distribution of the fracture volume of the target layer is obtained according to the data volume for reflecting the spatial distribution of the fracture volume.

[0057] In the implementation process, the spatial distribution of the fracture volume with double medium of the carbonate reservoir can be obtained according to the data volume for reflecting the spatial distribution of the fracture volume obtained in step S160 (as shown in Figure 7 、 Figure 8 、 Figure 9 ).

[0058] In the embodiments provided in the present application, starting from the stress field, five-dimensional gather data reflecting the anisotropy of the reservoir are used to carry out research, and the spatial distribution of the fracture volume with double medium is obtained. The method provided in the embodiments of the present application is applied in the Chaixi Sag, effectively predicts the spatial distribution of the fracture volume, and two newly drilled wells are consistent with the predicted spatial distribution of the fracture volume; and high yield is obtained according to the predicted spatial distribution of the fracture volume, which shows that the prediction of the spatial distribution of the fracture volume based on the method provided in the embodiments of the present application is accurate and feasible. The method provided in the embodiments of the present application can be used to depict faults, and can be further used to depict fractures and pores, and is suitable for various stages of oilfield development.

[0059] Figure 10 is an example schematic diagram of a five-dimensional gather fracture volume carving device based on some embodiments of the present application.

[0060] As shown in Figure 10 , the five-dimensional gather fracture volume carving device based on the five-dimensional gather fracture volume carving device includes a first determination module 1010, a second determination module 1020, a first acquisition module 1030, a second acquisition module 1040, a third acquisition module 1050, a fourth acquisition module 1060 and a fifth acquisition module 1070.

[0061] The first determination module 1010 is configured to determine a wellbore peripheral stress field direction of a drilled well in a target work area.

[0062] The second determination module 1020 is configured to determine a maximum principal stress direction of the target work area according to a coherent plane map of a target layer in the target work area and the wellbore peripheral stress field direction of the drilled well in the target work area; wherein the coherent plane map of the target layer is obtained from a seismic coherent attribute extracted from full stack seismic data of the target work area.

[0063] The first acquisition module 1030 is configured to obtain a fracture attribute and a fracture plane distribution attribute of the target work area from five-dimensional gather data of the target work area according to the maximum principal stress direction of the target work area.

[0064] The second acquisition module 1040 is configured to determine a correlation between seismic attributes and physical properties of the target work area according to physical property information of the drilled well in the target work area and full stack seismic data of the target work area; wherein the seismic attributes at least include amplitudes of seismic waves, and the physical properties at least include fracture data and dissolution pore data.

[0065] The third acquisition module 1050 is configured to calculate a fracture attribute volume of the target layer according to seismic attributes of the target layer extracted from the full stack data of the target work area by using the correlation between the seismic attributes and the physical properties of the target work area; wherein the fracture attribute volume is used to reflect physical properties of a fracture volume of the target layer.

[0066] The fourth acquisition module 1060 is configured to fuse the fracture attribute, the fracture plane distribution attribute, and the fracture attribute volume to obtain a data volume used to reflect spatial distribution of the fracture volume.

[0067] The fifth acquisition module 1070 is configured to obtain spatial distribution of the fracture volume of the target layer according to the data volume.

[0068] In some embodiments, the determination of the wellbore peripheral stress field direction of the drilled well in the target work area includes:

[0069] The azimuth of the induced fracture of the drilled well is determined by using imaging logging data of the drilled well in the target work area; and the wellbore peripheral stress field direction of the drilled well in the target work area is determined according to the azimuth of the induced fracture.

[0070] In some embodiments, the coherent plane map of the target layer is obtained by: determining the seismic reflection interface of the target layer by using the result of the synthetic record calibration of the drilled well of the target work area; obtaining the horizon data of the target layer from the full stack seismic data of the target work area according to the seismic reflection interface of the target layer; and obtaining the coherent plane map of the target layer by extracting the coherent attribute of the target layer from the full stack seismic data according to the horizon data.

[0071] In some embodiments, the fracture attribute and the fracture plane distribution attribute of the target work area are obtained from the five-dimensional gather data of the target work area according to the maximum principal stress direction of the target work area, including: performing stacking on the five-dimensional gather data according to the division azimuth to obtain an angle-stacked data volume; obtaining a fracture volume attribute analysis data volume from the angle-stacked data volume, which is perpendicular to the maximum principal stress direction of the target work area; and obtaining the fracture attribute and the fracture plane distribution attribute by extracting the coherent attribute from the fracture volume attribute analysis data volume.

[0072] In the above embodiment of the five-dimensional gather fracture volume carving device, the specific processing of each module and the technical effects brought by the specific processing can be referred to the related description in the corresponding method embodiment, which will not be repeated here.

[0073] Figure 11 is an exemplary structural schematic diagram of an electronic device according to some embodiments of the present application.

[0074] As shown in Figure 11 The electronic device includes at least one processor 1101, at least one communication interface 1102, at least one memory 1103, and at least one communication bus 1104. Optionally, the communication interface 1102 can be an interface of a communication module, such as an interface of a GSM module. The processor 1101 can be a processor CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement embodiments of the present application. The memory 1103 can include a high-speed RAM memory, and can also include a non-volatile memory, such as at least one disk memory. The memory 1103 stores a program, and the processor 1101 invokes the program stored in the memory 1103 to execute part or all of the above-mentioned method embodiments.

[0075] The present application relates to a storage medium for storing a computer readable program, which, when executed, performs part or all of the above-mentioned method embodiments.

[0076] Optionally, the storage medium can be a non-transitory computer-readable storage medium, for example, the non-transitory computer-readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, and an optical data storage device, etc.

[0077] Based on the same inventive concept, the embodiments of the present application further provide a computer program product comprising a computer program, which, when executed by a processor, implements some or all of the method embodiments described above.

[0078] The above has described the basic concept, and it is obvious that the above detailed disclosure is only taken as an example and does not limit the present application. Although it is not explicitly stated here, those skilled in the art can make various modifications, improvements and corrections to the present application. Such modifications, improvements and corrections are suggested in the present application, so such modifications, improvements and corrections still belong to the spirit and scope of the exemplary embodiments of the present application.

[0079] At the same time, specific words are used in the present application to describe the embodiments of the present application. For example, "one embodiment", "an embodiment", and / or "some embodiments" means a certain feature, structure or characteristic related to at least one embodiment of the present application. Therefore, it should be emphasized and noted that the "one embodiment" or "one embodiment" or "one alternative embodiment" mentioned in different places in the present application does not necessarily refer to the same embodiment. In addition, some features, structures or characteristics in one or more embodiments of the present application can be properly combined.

[0080] In addition, unless the claim explicitly states, the order of the processing elements and sequences described in the present application, the use of numerals and letters, or the use of other names, is not intended to limit the order of the processes and methods of the present application. Although some currently considered useful embodiments of the application are discussed in the above disclosure through various examples, it should be understood that such details are only for the purpose of illustration, and the additional claims are not limited to the disclosed embodiments, on the contrary, the claims are intended to cover all modifications and equivalent combinations that meet the spirit and scope of the embodiments of the present application. For example, although the system components described above can be realized by hardware devices, they can also be realized by software solutions only, such as installing the described system on existing servers or mobile devices.

[0081] Similarly, it should be noted that, in order to simplify the description of the present application and to help understand one or more embodiments of the present application, sometimes various features are combined into one embodiment, figure or description of the same. However, this disclosure method does not mean that the features required by the present application are more than the features mentioned in the claims. In fact, the features of the embodiments are less than all the features of the single embodiment disclosed above.

[0082] Some embodiments use numerical values to describe components, quantities of attributes, and the like. It should be understood that such numerical values used in the description of embodiments are, in some examples, modified by the adjectives "about," "approximately," or "substantially." Unless otherwise stated, "about," "approximately," or "substantially" indicate that the described dimension allows for ±20% variation. Accordingly, numerical values used in the specification and claims of some embodiments are approximations. Variations to these values can occur depending on the desired properties sought to be obtained by a particular embodiment. In some embodiments, the numerical values are those already described herein. In some embodiments, numerical values used in the specification and claims are one of the approximate values. In some embodiments, numerical values should be considered in the context and with the understanding that numbers of units in scientific or technical calculations typically contain certain numbers of significant figures and that rounding according to generally accepted significant figure rules can occur, as appropriate, to the extent that variations will occur in the feasible range.

[0083] Each patent, patent application, patent publication, and other material, such as articles, books, specifications, publications, documents, and the like, referenced herein are hereby incorporated by reference in their entirety for the teachings relevant to the sentence and / or paragraph in which such reference is made. Discrepancies between art history documents and the present disclosure, and between art history documents that limit the scope of the claims of the present application, are hereby excluded. It is specifically intended that the description, definitions, and / or terminology used herein be interpreted in accordance with the description, definitions, and / or terminology used in the present disclosure, and not a description, definitions, and / or terminology used in the art history documents.

[0084] Finally, it should be understood that the embodiments described herein are merely exemplary of the application. Other variations of the embodiments can be devised by those skilled in the art without departing from the scope of the present application. Accordingly, while the present application is depicted and described by reference to exemplary embodiments, it will be readily apparent to those of ordinary skill in the art that changes in form and detail can be made hereto without departing from the spirit and scope of the application. Accordingly, the disclosed embodiments are to be considered as illustrative and not restrictive, and the scope of the application is not to be determined by the description of the embodiments, but rather by the claims.

Claims

1. A method of carving based on five-dimensional gather discontinuity volume, characterized in that, The method comprises: determining a wellbore peripheral stress field direction of a drilled well in a target work area; determining a maximum principal stress direction of the target work area according to a coherent plane map of a target layer in the target work area and the wellbore peripheral stress field direction of the drilled well in the target work area, wherein the coherent plane map of the target layer is obtained by extracting a seismic coherent attribute from full stack seismic data of the target work area; obtaining a fracture attribute and a fracture plane distribution attribute of the target work area from five-dimensional gather data of the target work area according to the maximum principal stress direction of the target work area; determining a correlation between seismic attributes and physical properties of the target work area according to physical property information of the drilled well in the target work area and the full stack seismic data of the target work area, wherein the seismic attributes at least include amplitudes of seismic waves, and the physical properties at least include fracture data and dissolution pore data; calculating a fracture attribute volume of the target layer according to seismic attributes of the target layer extracted from the full stack data of the target work area by using the correlation between the seismic attributes and the physical properties of the target work area, wherein the fracture attribute volume is used to reflect physical properties of a fracture volume of the target layer; fusing the fracture attribute, the fracture plane distribution attribute and the fracture attribute volume to obtain a data volume used to reflect spatial distribution of the fracture volume; obtaining spatial distribution of the fracture volume of the target layer according to the data volume; obtaining the fracture attribute and the fracture plane distribution attribute of the target work area from the five-dimensional gather data of the target work area according to the maximum principal stress direction of the target work area, comprising: stacking the five-dimensional gather data according to a division azimuth to obtain a data volume of angle division stacking; selecting a data volume perpendicular to the maximum principal stress direction of the target work area from the data volume of angle division stacking to obtain a data volume for fracture attribute analysis; extracting a coherent attribute from the data volume for fracture attribute analysis to obtain the fracture attribute and the fracture plane distribution attribute.

2. The method of claim 1, wherein, The method comprises: determining an azimuth of an induced fracture of the drilled well in the target work area by using imaging logging data of the drilled well in the target work area; determining the wellbore peripheral stress field direction of the drilled well in the target work area according to the azimuth of the induced fracture.

3. The method of claim 2, wherein, The coherent plane map of the target layer is obtained by: determining a seismic reflection interface of the target layer by using a synthetic record calibration result of the drilled well in the target work area; obtaining horizon data of the target layer by interpreting the seismic reflection interface of the target layer from full stack seismic data of the target work area; extracting a coherent attribute of the target layer from the full stack seismic data according to the horizon data to obtain the coherent plane map of the target layer.

4. A five-dimensional gather discontinuity volume based sculpting apparatus, characterized by, The device comprises: a first determining module configured to determine a wellbore peripheral stress field direction of a drilled well in a target work area; a second determining module configured to determine a maximum principal stress direction of the target work area according to a coherent plane map of a target layer in the target work area and the wellbore peripheral stress field direction of the drilled well in the target work area, wherein the coherent plane map of the target layer is obtained by extracting a seismic coherent attribute from full stack seismic data of the target work area; The first obtaining module is configured to obtain fracture attributes and fracture plane distribution attributes of the target work area from five-dimensional gather data of the target work area according to a maximum principal stress direction of the target work area; The second obtaining module is configured to determine a correlation between seismic attributes and physical properties of the target work area according to physical property information of drilled wells of the target work area and full-stack seismic data of the target work area; the seismic attributes at least include amplitudes of seismic waves, and the physical properties at least include fracture data and dissolution pore data; The third obtaining module is configured to calculate fracture attribute volumes of the target layer according to seismic attributes of the target layer extracted from the full-stack data of the target work area by using the correlation between the seismic attributes and the physical properties of the target work area; the fracture attribute volumes are used to reflect physical properties of fracture bodies of the target layer; The fourth obtaining module is configured to fuse the fracture attributes, the fracture plane distribution attributes and the fracture attribute volumes to obtain a data volume reflecting spatial distribution of fracture bodies; The fifth obtaining module is configured to obtain spatial distribution of fracture bodies of the target layer according to the data volume. The fracture attributes and the fracture plane distribution attributes of the target work area are obtained from five-dimensional gather data of the target work area according to a maximum principal stress direction of the target work area, including: The five-dimensional gather data is stacked according to division azimuth to obtain a data volume of angle-stacked data; A data volume perpendicular to the maximum principal stress direction of the target work area is selected from the data volume of angle-stacked data to obtain a data volume for fracture body attribute analysis; Coherent attributes are extracted from the data volume for fracture body attribute analysis to obtain the fracture attributes and the fracture plane distribution attributes.

5. The apparatus of claim 4, wherein, The direction of the stress field of the wellbore periphery of the drilled well of the target work area is determined, including: The direction of induced fractures of the drilled well of the target work area is determined by using imaging logging data of the drilled well of the target work area; The direction of the stress field of the wellbore periphery of the drilled well of the target work area is determined according to the direction of the induced fractures.

6. The apparatus of claim 4, wherein, The coherent plane map of the target layer is obtained in the following manner: The seismic reflection interface of the target layer is determined by using a synthetic record calibration result of the drilled well of the target work area; The horizon data of the target layer are interpreted from full-stack seismic data of the target work area according to the seismic reflection interface of the target layer; The coherent attributes of the target layer are extracted from the full-stack seismic data according to the horizon data to obtain the coherent plane map of the target layer. 7.An electronic device, comprising a memory and a processor, the memory storing a computer program, and the processor executing the program to perform the method of any one of claims 1 to 3. 8.A storage medium for storing a computer readable program, the computer readable program being executed to perform the method of any one of claims 1 to 3.

Images