CN111948709A - Fracture zone structure identification method, device and system - Google Patents
Fracture zone structure identification method, device and system Download PDFInfo
- Publication number
- CN111948709A CN111948709A CN202010749194.0A CN202010749194A CN111948709A CN 111948709 A CN111948709 A CN 111948709A CN 202010749194 A CN202010749194 A CN 202010749194A CN 111948709 A CN111948709 A CN 111948709A
- Authority
- CN
- China
- Prior art keywords
- variance
- fracture zone
- information
- feature information
- identified
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 27
- 238000010586 diagram Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 16
- 239000011435 rock Substances 0.000 claims description 17
- 230000015654 memory Effects 0.000 claims description 8
- 230000007704 transition Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 238000005755 formation reaction Methods 0.000 description 23
- 238000005070 sampling Methods 0.000 description 16
- 238000011161 development Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 238000013467 fragmentation Methods 0.000 description 3
- 238000006062 fragmentation reaction Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The specification provides a method, a device and a system for identifying a fracture zone structure. The method comprises the steps of obtaining a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation; matching the variance characteristic graph with preset structure variance characteristic information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map; and determining the structure of the fracture zone to be identified according to the matching result. The accuracy of the fracture zone structure identification can be improved by the aid of the method and the device.
Description
Technical Field
The application relates to the technical field of tectonic geology and oil gas exploration and development, in particular to a fracture zone structure identification method, device and system.
Background
The fault fracture zone is a geologic body consisting of derived cracks near the fracture surface and fracture surface fillers, and is a complex three-dimensional fracture zone structure. Accurate identification and description of a fractured band structural pattern play an important role in horizontal well deployment through faults.
At present, the fault fracture band structure is mainly identified by means of a synthetic wave model, natural change analysis of earthquake activity rate, earthquake position, earthquake velocity distribution and the like. However, these methods indirectly depict the fracture band structure, which may cause inaccurate identification of the fracture band structure.
Therefore, there is a need for a solution to the above technical problems.
Disclosure of Invention
The embodiment of the specification provides a method, a device and a system for identifying a broken belt structure, which can improve the accuracy of identifying the broken belt structure.
The method, the device and the system for identifying the structure of the fracture zone are realized in the following modes.
A method of fracture zone structure identification, comprising: acquiring a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation; matching the variance characteristic graph with preset structure variance characteristic information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map; and determining the structure of the fracture zone to be identified according to the matching result.
A fracture zone structure identification device comprising: the acquisition module is used for acquiring a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation; the matching module is used for matching the variance feature map with preset structure variance feature information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map; and the determining module is used for determining the structure of the fracture zone to be identified according to the matching result.
A fracture zone structure identification system comprising at least one processor and a memory storing computer-executable instructions that when executed by the processor perform the steps of the method of any one of the method embodiments of the present specification.
The specification provides a method, a device and a system for identifying a fracture zone structure. In some embodiments, the accuracy of subsequent fracture zone structure identification can be effectively improved by depicting the fracture zone to be identified by using the seismic variance attribute. The internal structure of the fracture zone and the deformation pattern of the fracture zone are depicted in advance by utilizing the seismic variance attribute, so that the accuracy of fracture zone structure identification can be improved. By pre-storing the corresponding relation between the fracture band structure and the variance feature information, the method can be used for directly matching when obtaining the variance feature map of the fracture band to be identified, and determining the structure style corresponding to the fracture band to be identified, so that the identification efficiency is effectively improved. With the embodiments provided in the present specification, the accuracy of fracture zone structure identification can be improved.
Drawings
The accompanying drawings, which are included to provide a further understanding of the specification, are incorporated in and constitute a part of this specification, and are not intended to limit the specification. In the drawings:
FIG. 1 is a schematic flow chart diagram illustrating one embodiment of a method for identifying a fracture zone structure provided herein;
FIG. 2 is a schematic illustration of one type of acquisition of seismic variance attributes provided herein;
FIG. 3 is a plot of the in-plane variance attribute of a target formation as provided by the present description;
FIG. 4 is a variance feature map obtained after processing a planar variance attribute map of a target formation as provided herein;
FIG. 5 is a schematic diagram of an embodiment of a relationship between a fracture zone structure and a variance feature map provided in the present specification;
FIG. 6 is a schematic block diagram of an embodiment of a fracture zone structure identification apparatus provided in the present specification.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in the embodiments of the present specification will be clearly and completely described below with reference to the drawings in the embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments in the present specification, and not all of the embodiments. All other embodiments that can be obtained by a person skilled in the art on the basis of one or more embodiments of the present description without inventive step shall fall within the scope of protection of the embodiments of the present description.
At present, the fault fracture band structure is mainly identified by means of a synthetic wave model, natural change analysis of earthquake activity rate, earthquake position, earthquake velocity distribution and the like. For example, the analysis of natural variation of seismic activity rate obtains that the farther the seismic activity rate is from the fault in the direction of the vertical fault strike, the smaller the seismic activity rate by analyzing the relationship between the seismic activity rate and the fracture zone. However, these methods indirectly depict the fracture band structure, which may cause inaccurate identification of the fracture band structure.
With the development of three-dimensional seismic technology, information such as geological structures and formation characteristics in the underground can be visualized and analyzed in a manner similar to the earth surface through geophysical data. Seismic attributes, which are components of the seismic data, may be measured, calculated, or otherwise extracted from the seismic data. Seismic geometry was first proposed by Taner et al, which classifies seismic attributes into two broad categories, geometric and physical. Geometric seismic attributes are understood to be attributes related to the geometry of the subsurface formations, and physical seismic attributes are understood to be attributes related to dynamics and kinematics.
Coherence properties, one of the geometric seismic properties, can be used primarily to characterize continuity between subsurface formations, and are most useful for identifying faults and fractures in subsurface formations. The coherence value is obtained by calculating the similarity degree of the waveforms of adjacent seismic traces through a certain algorithm, the stratum where the fault or fracture develops is discontinuous, and the coherence coefficient is low. The variance attribute and the coherence attribute are similar, the basic principle is that in a uniform and continuous stratum, the waveforms of reflected waves of adjacent seismic channels are similar, and the seismic waveforms have difference when the stratum is not uniform and discontinuous due to the conditions of fault, crack development or lithologic mutation, and the like, and the detection of the difference can achieve the purpose of extracting fault and other special structure development information.
The following describes an embodiment of the present disclosure with a specific application scenario as an example. Specifically, fig. 1 is a schematic flow chart of an embodiment of a method for identifying a fracture zone structure provided in the present specification. Although the present specification provides the method steps or apparatus structures as shown in the following examples or figures, more or less steps or modules may be included in the method or apparatus structures based on conventional or non-inventive efforts. In one embodiment, as shown in fig. 1, the present specification provides a method for identifying a fracture zone structure, which may include the following steps.
S0: acquiring a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation.
In embodiments of the present disclosure, the fracture zone to be identified may be a fracture zone, a composite fracture zone, or the like. The fault zone may also be referred to as a fault fracture zone, which consists of fault nuclei and fracture zones around the fault nuclei. Fault nuclei are typically narrow, localized slip bands containing high strain products such as fault mud, cracked rock, glutenite, and possibly wider regions containing multiple parallel or intersecting fault nuclei. The fracture zone is generally a relatively developed rock with fractures that exhibit significantly lower deformation strength than the fault nuclei. The width of the fault fracture zone is the sum of the widths of the fault nucleus and the fracture zones around the fault nucleus. It is understood that a fault-fractured zone may include a single fault nucleus or may include multiple fault nuclei. If the fault fracture zone comprises a single fault nucleus, the fault fracture zone is a single-nucleus fault fracture zone; if the fault fracture zone includes a plurality of fault nuclei, the fault fracture zone is a multi-nuclear fault fracture zone.
In an embodiment of the present description, the variance feature map may be obtained based on a quantification of a planar variance attribute map of the target formation. The plane variance attribute map may be obtained based on the seismic variance attribute.
In general, variance can be used to measure the deviation of a random variable from its averageThe case of a mean, which can be used to characterize the degree of dispersion of a set of data. A set of data is set: x is the number of1，x2，…，xi，…，xnN represents the number of data, and the average value is
The variance D of the group of data2Can be obtained by the formula (1). As can be seen from equation (1), the more discrete (the greater the difference) the set of data, the greater the variance value, and vice versa.
In the embodiments of the present description, the seismic variance attribute is primarily a technique developed based on variance for detecting formation discontinuities. Generally, if abnormal geological structures such as faults, cracks, karst caves and the like develop on the reflection interface, the waveforms between adjacent seismic channels theoretically have relatively large differences. The seismic variance attribute may be understood as the degree of difference from trace to trace by calculating the variance between the seismic trace waveforms. The development condition of the underground fault can be detected through the seismic variance attribute.
In some implementations, the process of obtaining the seismic variance attribute may include: acquiring seismic trace data of a target area, and selecting a proper time window length and a proper sampling interval (for example, the sampling interval is 1ms, and the time window length is 30 ms); taking the current sampling point as the center, and taking the sampling point with the length of one half time window from top to bottom according to the sampling interval, namely the current sampling point is positioned at the central position of one time window length on the section (as shown in figure 2); calculating the average amplitude value of sampling points in each seismic channel within the length of each time window; calculating the sum of the variances of each seismic trace within the selected time window length based on the amplitude mean; and obtaining the variance value of a single sampling point after weighting normalization processing.
In some implementations, the weighting function is typically a trigonometric function. The variance value of the seismic variance attribute is quantitative. The weighted normalization process can mainly control the variance value between 0 and 1. In the fault development area, the stronger the structural deformation is, the larger the waveform difference between adjacent seismic channels is, the higher the variance is, and the closer the value is to 1.
In some implementation scenarios, the variance of each seismic trace sampling point and the variance of each seismic trace sampling point after the weighting normalization processing can be calculated according to the following formula:
wherein,
is the variance of the sampling points and is,
is the weighted normalized sampling point variance uiIs the amplitude value of the ith seismic trace,
is the average value of the amplitudes of all seismic channels, L is the length of a time window, T is the number of adjacent seismic channels involved in calculating the variance, T represents a certain sampling point, j represents all the sampling points in a time window range, w is a weighting function, w is the average value of the amplitudes of all seismic channelsj-tA triangular weighting function of a certain sampling point in a certain time window, wherein w is more than or equal to 0j-t≤1。
Based on the above description, the variance value of each sampling point of the target formation can be calculated according to the formulas (2) to (4), so as to obtain the plane variance attribute map of the target formation. The plane variance attribute map may exhibit distribution morphology and variation trend of faults.
In some implementation scenarios, since the seismic trace data volume is large and the seismic variance attribute calculation is complex, the seismic variance attribute of the target formation may be calculated by computer software, so as to obtain a plane variance attribute map of the target formation. For example, in some implementation scenarios, the petrel software may be used to import the seismic trace data of the target formation into the petrel software, and then the plane variance attribute map of the target formation may be obtained according to the output result of the petrel software. The petrel software is developed by Schlumberger company, is an exploration and development integrated platform taking a three-dimensional geological model as a center, and belongs to geophysical professional software.
As shown in fig. 3, fig. 3 is a plane variance attribute graph of a target formation provided in this specification, where broken lines corresponding to F1, F2, F3, F4, F5, and F6 in an upper graph represent faults, a lower graph is an enlarged graph corresponding to a broken-line box in the upper graph, solid lines corresponding to F3-1, F3-2b, F4, and F4-b represent survey lines arranged when fracture zones of the faults F3 and F4 are depicted, and black dots represent well positions.
In an embodiment of the present description, after obtaining the planar variance attribute map of the target formation, a variance feature map may be obtained based on the planar variance attribute map. The variance characteristic diagram can be understood as a result obtained by depicting a fracture zone developed on a target stratum on the basis of acquiring a plane seismic variance attribute diagram of the target stratum.
In some embodiments of the present description, the obtaining the variance feature map based on the planar variance attribute map processing of the target formation may include: acquiring a plane variance attribute map of a target stratum; taking the position of a preset distance from a fracture zone in the plane variance attribute graph as a baseline position; the base line is parallel to the extension direction of the fracture zone; making a plurality of measuring lines perpendicular to the fracture zone by taking the base line as a starting point; each measuring line comprises a plurality of measuring points with equal spacing distance; and obtaining a variance characteristic diagram according to the relation between the variance value of the measuring point on each measuring line and the distance between the measuring point and the base line.
Fig. 4 is a variance feature map obtained after processing a planar variance attribute map of a target formation provided by the present specification. The abscissa represents the distance from the measurement point to the base line, the unit is meter (m), the ordinate represents the variance, the solid lines corresponding to L1, L2, and L3 represent the measurement lines, and the upper dotted line represents the overall trend of all the measurement lines, which may also be referred to as an envelope. The lower dotted line is a baseline indicating variance values of the surrounding rock portion.
Based on the above-described embodiments, the variance feature map of the fracture zone to be identified may be obtained, so as to identify the fracture zone according to the variance feature map subsequently.
In the embodiment of the specification, the variance attribute can well reflect the change of the fracture zone in a three-dimensional space, and the variance attribute value is obtained by calculation from seismic data, so that the fracture zone to be identified is depicted by using the seismic variance attribute, and the accuracy of subsequent fracture zone structure identification can be effectively improved.
S2: matching the variance characteristic graph with preset structure variance characteristic information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map.
In the embodiment of the present specification, after obtaining the variance feature map of the fractured zone to be identified, the variance feature map may be matched with preset structure variance feature information, so as to determine the structure of the fractured zone to be identified.
In some embodiments of the present disclosure, the preset structure variance feature information may include a corresponding relationship between a fracture band structure and variance feature information, where the variance feature information is obtained according to a variance feature map.
In some implementations, the fracture zone structure can include a pulse-like structure, a ribbon-like structure, a box-like structure, a dome-like structure.
In the embodiments of the present specification, the variance feature information may be understood as feature information related to a fracture zone structure extracted from a variance feature map, for example, the presence or absence of a peak, the presence of several peaks, the shape of the whole representation, and the like.
In some implementation scenarios, the variance feature information may be divided into variance information corresponding to different regions according to the variance feature map. For example, in some implementation scenarios, the variance characteristic information may include variance information corresponding to a fault nucleus region, variance information corresponding to a fracture zone region, and variance information corresponding to a surrounding rock region. The fault nucleus region can also be called a variance high value region, the broken belt region can also be called a variance median region, and the surrounding rock region can also be called a variance low value region.
In some implementation scenarios, the variance corresponding to the fault nucleus region is greater than a first threshold, the variance corresponding to the broken zone region is greater than a second threshold and less than the first threshold, and the variance corresponding to the surrounding rock region is greater than a third threshold and less than the second threshold. For example, in some implementations, the first threshold may be 0.4, the second threshold may be 0.1, and the third threshold may be 0.
In some embodiments of the present description, the variance characteristic information corresponding to the pulse-like structure may include a peak value, where the peak value is located in a fault nucleus region, and a variance value on both sides of the peak value exponentially decays to a surrounding rock region; the variance characteristic information corresponding to the strip-shaped structure can comprise at least two peak values, each peak value is located in different fault nucleus areas, and the different fault nucleus areas are separated by a crushing belt area; the variance characteristic information corresponding to the box-shaped structure can include that a variance value in a preset range is located in a fault nucleus area, and the variance information of the fault nucleus area is box-shaped; the variance characteristic information corresponding to the dome-shaped structure is in a dome shape as a whole, and the variance information of the crushing zone area and the variance information of the surrounding rock area are in smooth transition.
In this embodiment of the present description, a corresponding relationship between a fracture zone structure and a variance feature map may be obtained first, and then variance feature information corresponding to the fracture zone structure is extracted from the variance feature map, so as to obtain preset structure variance feature information.
As shown in fig. 5, fig. 5 is a schematic diagram of an embodiment of a corresponding relationship between a fracture zone structure and a variance feature map provided in the present specification. In fig. 5, F3-1, F3-2B, F4, and F4-B indicate cross-sectional positions arranged in 4 sections across the faults F3 and F4, a graph corresponding to a mark a indicates a variance feature map corresponding to a pulse-like structure, a graph corresponding to a mark B indicates a pulse-like structure, a graph corresponding to a mark C indicates a variance feature map corresponding to a band-like structure, a graph corresponding to a mark D indicates a band-like structure, a graph corresponding to a mark E indicates a variance feature map corresponding to a box-like structure, a graph corresponding to a mark F indicates a box-like structure, a graph corresponding to a mark G indicates a variance feature map corresponding to a dome-like structure, and a graph corresponding to a mark H indicates a dome-like structure. Wherein, can include broken zone in the fracture band structure, can include main fault nuclear secondary fault in the broken zone. The solid lines corresponding to L1, L2, L3, etc. represent the lines, and the upper dashed line represents the overall trend of all lines, which may also be referred to as an envelope. The lower dotted line is a baseline indicating variance values of the surrounding rock portion.
Further, variance feature information corresponding to the fracture zone structure can be extracted from the variance feature map, so that preset structure variance feature information is obtained.
For example, since the variance characteristic map corresponding to the pulse-like structure has a large variance value (>0.4, fault nucleus) in the middle fault nucleus region and two intermediate variance value bands (0.1-0.4, fragmentation band) beside the middle fault nucleus region, which exponentially attenuate to the matrix (0-0.1, surrounding rock region), the variance characteristic information corresponding to the pulse-like structure can be extracted to include a peak value, which is located in the fault nucleus region, and the variance values on both sides of the peak value exponentially attenuate to the surrounding rock region. Since a plurality of pulse-shaped high-intensity strain regions (fault nuclei) develop in the variance feature map corresponding to the strip-shaped structure and are separated by the corresponding fragmentation zones, it can be extracted that the variance feature information corresponding to the strip-shaped structure may include at least two peaks, each peak being located in a different fault nucleus region, and the different fault nucleus regions being separated by the fragmentation zone regions. The variance characteristic graph corresponding to the box-shaped structure shows that a series of tighter fault nuclei share a box-shaped crushing area, no obvious pulse or cusp exists on the variance attribute value distribution graph, the actual condition of the pattern is that the distance between the fault nuclei is short or the fault nuclei are intersected and merged and the secondary fault nuclei are not too developed, the overall deformation degree of the fault crushing area is strong, the variance characteristic information corresponding to the box-shaped structure can be extracted, the variance characteristic information can include that the variance value in a preset range is located in the fault nucleus area, and the variance information of the fault nucleus area is in the box shape. The variance characteristic map corresponding to the dome-shaped structure shows that the high-strength deformation area of the fractured zone is in a smooth transition to the surrounding rock, the whole variance characteristic map is in a dome shape and is related to the development of the secondary fault nucleus around the main fault nucleus, the variance characteristic information corresponding to the dome-shaped structure can be extracted, the whole variance characteristic information is in a dome shape, and the variance information of the fractured zone is in a smooth transition with the variance information of the surrounding rock area.
In this embodiment, after the variance feature information is extracted from the variance feature map, the corresponding relationship between the fracture zone structure and the variance feature information may be stored, so that the fracture zone structure may be correspondingly identified after the variance feature map or the variance feature information is obtained.
Based on the above-described embodiment, the preset structure variance feature information can be correspondingly obtained, and then the variance feature map of the fracture zone to be identified is matched with the preset structure variance feature information, so as to obtain the structure of the fracture to be identified.
In some embodiments of the present specification, the matching the variance feature map with preset structural feature information may include: extracting variance feature information of the fracture zone to be identified from the variance feature map to obtain variance feature information to be identified; and matching the variance characteristic information to be identified with variance characteristic information in preset structure variance characteristic information.
In the embodiment of the specification, the internal structure of the fracture zone and the deformation pattern of the fracture zone are described in advance by using the seismic variance attribute, so that the accuracy of fracture zone structure identification can be improved.
S4: and determining the structure of the fracture zone to be identified according to the matching result.
In the embodiment of the present specification, after the variance feature map is matched with the preset structure variance feature information, the structure of the fracture to be identified may be determined according to the matching result.
In some implementation scenarios, if the preset structure variance feature information does not include a fracture band structure corresponding to the variance feature map, feature information related to the fracture band structure may be extracted from the variance feature map, so as to obtain variance feature information, and define a fracture band structure based on the variance feature information, and finally, the corresponding relationship may be saved, so that the next identification is more accurate.
According to the fracture zone structure identification method provided by the specification, the fracture zone to be identified is depicted by utilizing the seismic variance attribute, so that the accuracy of subsequent fracture zone structure identification can be effectively improved. The internal structure of the fracture zone and the deformation pattern of the fracture zone are depicted in advance by utilizing the seismic variance attribute, so that the accuracy of fracture zone structure identification can be improved. By pre-storing the corresponding relation between the fracture band structure and the variance feature information, the method can be used for directly matching when obtaining the variance feature map of the fracture band to be identified, and determining the structure style corresponding to the fracture band to be identified, so that the identification efficiency is effectively improved.
In the present specification, each embodiment of the method is described in a progressive manner, and the same and similar parts in each embodiment may be joined together, and each embodiment focuses on the differences from the other embodiments. Reference is made to the description of the method embodiments.
Based on the method for identifying the structure of the broken belt, one or more embodiments of the present specification further provide a device for identifying the structure of the broken belt. The apparatus may include systems (including distributed systems), software (applications), modules, components, servers, clients, etc. that use the methods described in the embodiments of the present specification in conjunction with any necessary apparatus to implement the hardware. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
Specifically, fig. 6 is a schematic block diagram of an embodiment of a fracture zone structure identification device provided in the present specification, and as shown in fig. 6, the fracture zone structure identification device provided in the present specification may include: an acquisition module 120, a matching module 122, and a determination module 124.
The obtaining module 120 may be configured to obtain a variance feature map of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation;
a matching module 122, configured to match the variance feature map with preset structure variance feature information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map;
the determining module 124 may be configured to determine the structure of the fracture zone to be identified according to the matching result.
It should be noted that the above-mentioned description of the apparatus according to the method embodiment may also include other embodiments, and specific implementation manners may refer to the description of the related method embodiment, which is not described herein again.
The present specification also provides an embodiment of a fracture zone structure identification system, comprising at least one processor and a memory storing computer-executable instructions, which when executed by the processor, implement the steps of the method described in any one or more of the above embodiments, for example, comprising: acquiring a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation; matching the variance characteristic graph with preset structure variance characteristic information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map; and determining the structure of the fracture zone to be identified according to the matching result. The system may be a single server, or may include a server cluster, a system (including a distributed system), software (applications), an actual operating device, a logic gate device, a quantum computer, etc. using one or more of the methods or one or more of the example devices of the present specification, in combination with a terminal device implementing hardware as necessary.
The method or apparatus provided by the present specification and described in the foregoing embodiments may implement service logic through a computer program and record the service logic on a storage medium, where the storage medium may be read and executed by a computer, so as to implement the effect of the solution described in the embodiments of the present specification. The storage medium may include a physical device for storing information, and typically, the information is digitized and then stored using an electrical, magnetic, or optical media. The storage medium may include: devices that store information using electrical energy, such as various types of memory, e.g., RAM, ROM, etc.; devices that store information using magnetic energy, such as hard disks, floppy disks, tapes, core memories, bubble memories, and usb disks; devices that store information optically, such as CDs or DVDs. Of course, there are other ways of storing media that can be read, such as quantum memory, graphene memory, and so forth.
The embodiments of the method or the apparatus for identifying a fractured band structure provided in this specification may be implemented in a computer by a processor executing corresponding program instructions, for example, implemented in a PC end using a c + + language of a windows operating system, implemented in a linux system, or implemented in an intelligent terminal using android, iOS system programming languages, implemented in processing logic based on a quantum computer, or the like.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that the implementation can be by computer program instructions which can be provided to a processor of a general purpose computer, special purpose computer, embedded processor or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified.
The above description is merely exemplary of one or more embodiments of the present disclosure and is not intended to limit the scope of one or more embodiments of the present disclosure. Various modifications and alterations to one or more embodiments described herein will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims.
Claims (10)
1. A method for identifying a fracture zone structure, comprising:
acquiring a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation;
matching the variance characteristic graph with preset structure variance characteristic information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map;
and determining the structure of the fracture zone to be identified according to the matching result.
2. The method of claim 1, wherein the obtaining a variance profile based on a planar variance attribute map process for the target formation comprises:
acquiring a plane variance attribute map of a target stratum;
taking the position of a preset distance from a fracture zone in the plane variance attribute graph as a baseline position; the base line is parallel to the extension direction of the fracture zone;
making a plurality of measuring lines perpendicular to the fracture zone by taking the base line as a starting point; each measuring line comprises a plurality of measuring points with equal spacing distance;
and obtaining a variance characteristic diagram according to the relation between the variance value of the measuring point on each measuring line and the distance between the measuring point and the base line.
3. The method of claim 1, wherein the fractured band structure comprises a pulse-like structure, a ribbon-like structure, a box-like structure, a dome-like structure.
4. The method according to claim 1, wherein the variance feature information comprises variance information corresponding to a fault nucleus region, variance information corresponding to a fracture zone region and variance information corresponding to a surrounding rock region.
5. The method of claim 4, wherein the variance corresponding to the fault nucleus region is greater than a first threshold, the variance corresponding to the fracture zone region is greater than a second threshold and less than the first threshold, and the variance corresponding to the wall rock region is greater than a third threshold and less than the second threshold.
6. The method of claim 5, wherein the first threshold is 0.4, the second threshold is 0.1, and the third threshold is 0.
7. The method of claim 3, wherein the variance characterizing information corresponding to the pulse-like structure comprises a peak, the peak is located in a fault nucleus region, and the variance value on both sides of the peak exponentially decays to a surrounding rock region;
the variance characteristic information corresponding to the strip-shaped structure comprises at least two peak values, each peak value is located in different fault nucleus areas, and the different fault nucleus areas are separated through a crushing belt area;
the variance characteristic information corresponding to the box-shaped structure comprises a variance value in a preset range and is located in a fault nucleus area, and the variance information of the fault nucleus area is box-shaped;
the variance characteristic information corresponding to the dome-shaped structure is in a dome shape as a whole, and the variance information of the crushing zone area and the variance information of the surrounding rock area are in smooth transition.
8. The method of claim 1, wherein the matching the variance feature map with preset structural feature information comprises:
extracting variance feature information of the fracture zone to be identified from the variance feature map to obtain variance feature information to be identified;
and matching the variance characteristic information to be identified with variance characteristic information in preset structure variance characteristic information.
9. A fracture zone structure identification device, comprising:
the acquisition module is used for acquiring a variance characteristic diagram of a fracture zone to be identified; wherein the variance feature map is obtained based on a planar variance attribute map processing of the target formation;
the matching module is used for matching the variance feature map with preset structure variance feature information; the preset structure variance feature information comprises a corresponding relation between a fracture zone structure and variance feature information, and the variance feature information is obtained according to a variance feature map;
and the determining module is used for determining the structure of the fracture zone to be identified according to the matching result.
10. A fracture zone structure identification system comprising at least one processor and memory storing computer-executable instructions that when executed by the processor implement the steps of the method of any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010749194.0A CN111948709B (en) | 2020-07-30 | 2020-07-30 | Fracture zone structure identification method, device and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010749194.0A CN111948709B (en) | 2020-07-30 | 2020-07-30 | Fracture zone structure identification method, device and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111948709A true CN111948709A (en) | 2020-11-17 |
| CN111948709B CN111948709B (en) | 2021-08-06 |
Family
ID=73337999
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010749194.0A Active CN111948709B (en) | 2020-07-30 | 2020-07-30 | Fracture zone structure identification method, device and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111948709B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112198552A (en) * | 2020-09-25 | 2021-01-08 | 中国石油大学(北京) | Width determination method and device for pleat breaking system |
| CN113676945A (en) * | 2021-08-12 | 2021-11-19 | 锐捷网络股份有限公司 | Method, system, access controller and storage medium for identifying poor quality terminal |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103076623A (en) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | Crack detection method based on prestack coherence |
| CN105022092A (en) * | 2014-09-15 | 2015-11-04 | 中国矿业大学(北京) | Seismic data interpretation method and equipment used for determining fault distribution pattern |
| US10139508B1 (en) * | 2016-03-24 | 2018-11-27 | EMC IP Holding Company LLC | Methods and apparatus for automatic identification of faults on noisy seismic data |
| CN110488353A (en) * | 2019-09-04 | 2019-11-22 | 中海石油(中国)有限公司 | It is a kind of based on the fault interpretation method cuing open flat interaction and being combined with tectonic style guidance |
| CN110749924A (en) * | 2018-07-24 | 2020-02-04 | 中国石油化工股份有限公司 | Fracture zone identification method |
-
2020
- 2020-07-30 CN CN202010749194.0A patent/CN111948709B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103076623A (en) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | Crack detection method based on prestack coherence |
| CN105022092A (en) * | 2014-09-15 | 2015-11-04 | 中国矿业大学(北京) | Seismic data interpretation method and equipment used for determining fault distribution pattern |
| US10139508B1 (en) * | 2016-03-24 | 2018-11-27 | EMC IP Holding Company LLC | Methods and apparatus for automatic identification of faults on noisy seismic data |
| CN110749924A (en) * | 2018-07-24 | 2020-02-04 | 中国石油化工股份有限公司 | Fracture zone identification method |
| CN110488353A (en) * | 2019-09-04 | 2019-11-22 | 中海石油(中国)有限公司 | It is a kind of based on the fault interpretation method cuing open flat interaction and being combined with tectonic style guidance |
Non-Patent Citations (1)
| Title |
|---|
| 牛成民: "渤海海域隐性走滑断层形成机理、识别方法与控藏作用", 《中国海上油气》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112198552A (en) * | 2020-09-25 | 2021-01-08 | 中国石油大学(北京) | Width determination method and device for pleat breaking system |
| CN112198552B (en) * | 2020-09-25 | 2022-02-18 | 中国石油大学(北京) | Width determination method and device for pleat breaking system |
| CN113676945A (en) * | 2021-08-12 | 2021-11-19 | 锐捷网络股份有限公司 | Method, system, access controller and storage medium for identifying poor quality terminal |
| CN113676945B (en) * | 2021-08-12 | 2024-02-23 | 锐捷网络股份有限公司 | Method, system, access controller and storage medium for identifying bad terminal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111948709B (en) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5691958A (en) | Method for determining formation properties from seismic attributes | |
| EP2697734B1 (en) | Systems and methods for hydraulic fracture characterization using microseismic event data | |
| EP2502093B1 (en) | System and method for microseismic analysis | |
| EP3358376A1 (en) | Method and apparatus for unambiguously estimating seismic anisotropy parameters | |
| CN111948709B (en) | Fracture zone structure identification method, device and system | |
| EP2818896A2 (en) | Method and apparatus for estimating a parameter of a subsurface volume | |
| CN110320576A (en) | Method for forecasting advanced geology, device and storage medium | |
| US20150355354A1 (en) | Method of analyzing seismic data | |
| US20220120933A1 (en) | Method of detection of hydrocarbon horizontal slippage passages | |
| KR101324285B1 (en) | Method of modelling well log | |
| US20040220744A1 (en) | Method of measuring local similarities between several seismic trace cubes | |
| CN113970785A (en) | Method and system for predicting development of underground river crack, storage medium and electronic equipment | |
| CN116027454A (en) | Method for imaging contour map of sandstone thickness and intelligent drawing device | |
| AU709540B2 (en) | Determination of formation properties from seismic attributes | |
| CN112444871A (en) | Method and equipment for quantitatively analyzing crack spacing based on seismic response characteristics of scattered waves | |
| US11905809B2 (en) | Determining reservoir heterogeneity for optimized drilling location | |
| CN118091786B (en) | Method and device for correcting linear density of cracks based on underground logging data | |
| CN111856572B (en) | Method and device for determining width of fault fracture belt | |
| CN112147698B (en) | Crack development zone identification and feature determination method and system | |
| RU2201606C1 (en) | Method of typification and correlation of oil and gas productive rocks by borehole spectral-time parameters | |
| CN116482762A (en) | Method for calculating sub-seismic fault length | |
| Song et al. | Deciphering Cross-Well Strain Data from Single-Use Fiber with Slippage Effects During Fracturing | |
| CN116299664A (en) | Method, device and equipment for determining fracture-controlled fracture-cavity type oil reservoir body | |
| CN118501977A (en) | Method and device for detecting positions and sizes of goaf of strip mine | |
| CN117270050A (en) | Microseism positioning method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |