CN111830558B - Fracture zone engraving method - Google Patents
Fracture zone engraving method Download PDFInfo
- Publication number
- CN111830558B CN111830558B CN201910308488.7A CN201910308488A CN111830558B CN 111830558 B CN111830558 B CN 111830558B CN 201910308488 A CN201910308488 A CN 201910308488A CN 111830558 B CN111830558 B CN 111830558B
- Authority
- CN
- China
- Prior art keywords
- fracture
- data
- analyzed
- cave
- determining
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000000605 extraction Methods 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000007476 Maximum Likelihood Methods 0.000 claims description 5
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 230000002265 prevention Effects 0.000 claims description 2
- 238000005553 drilling Methods 0.000 abstract description 19
- 230000010354 integration Effects 0.000 abstract 1
- 206010017076 Fracture Diseases 0.000 description 151
- 208000010392 Bone Fractures Diseases 0.000 description 146
- 239000011324 bead Substances 0.000 description 20
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 208000012322 Raynaud phenomenon Diseases 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method 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
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method of fracture zone engraving comprising: determining the fracture zone boundary of the region to be analyzed, and further obtaining the fracture zone profile of the region to be analyzed; determining cave data in an area to be analyzed through wave impedance inversion; extracting a fault plane in the area to be analyzed by utilizing a fault automatic extraction technology, and determining fracture data and crack data in the area to be analyzed; intersecting the cave data, the fracture data and the fracture data with the fracture zone profile and taking an intersection to respectively obtain a cave data body, a fracture data body and a fracture data body within the fracture zone profile range; and rescaling the cave data body, the fracture data body and the fracture data body in the fracture zone outline range according to respective value range ranges, and fusing the cave data body, the fracture data body and the fracture data body into a data body to obtain the fracture zone three-dimensional space carving of the zone to be analyzed. According to the method, the carving of the fracture zone is realized by optimizing the sensitive attribute and adopting a 'ternary integration' mode, and an effective and reliable basis can be provided for knowing the deployment of the drilling track.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration and development, in particular to a fracture zone engraving method.
Background
The oil field along the north of the Tarim basin is different from the surface karst fractured-cavern type oil reservoir in the Tahe area, and is a fracture-control fractured-cavern type oil reservoir. Through comprehensive comparison of drilling and recording well logging, earthquake, test, production dynamics and the like, the Ordovician carbonate rock reservoir matrix in the northward region has poor physical properties and underdeveloped pores, main reservoir spaces are 'cavity' -shaped caves and crack zones in sliding fracture zones, and the main reservoir type is a crack-cave type.
At present, good oil and gas results are revealed by trunk fracture drilling in the northward region, but the internal heterogeneity of a fracture zone is extremely strong, transverse and longitudinal segmented characteristics exist, and the reservoir space distribution is very complex. Due to the lack of large-scale reservoir section logging and coring data, the type and the distribution of the reservoir are controversial at present, and therefore an effective reservoir geological model is not established yet. Meanwhile, the obvious and effective seismic response characteristics of the reservoir are not established in the geophysical emptying of the high-quality reservoir at the leakage level.
When drilling and deploying, a horizontal well technology is often adopted to cross a fracture zone, and a drill bit is used for searching a high-quality reservoir stratum. The method can realize the search of a high-quality reservoir in the horizontal direction, but the fracture zone also has heterogeneity in the longitudinal direction, the fracture zone is not a good reservoir in the vertical direction, and the vertical target point cannot be calibrated. Meanwhile, in the drilling blow-down level reservoirs, the premium reservoirs are only a few meters in the fracture zone, which is hundreds of meters wide, even thousands of meters wide. Therefore, various sensitive attributes are required to be searched from seismic data to represent various reservoirs in the fracture, and the internal structure of the fracture is displayed more intuitively and vividly through body fusion and three-dimensional carving, so that drilling deployment and resource quantity calculation are guided.
Disclosure of Invention
In order to solve the above problems, the present invention provides a method for engraving a fracture zone, the method comprising:
firstly, identifying the boundary of a fracture zone of a region to be analyzed, determining the boundary of the fracture zone of the region to be analyzed, and further obtaining the profile of the fracture zone of the region to be analyzed;
step two, determining cave data in the area to be analyzed through wave impedance inversion;
extracting a fault plane in the region to be analyzed by utilizing a fault automatic extraction technology, and determining fracture data and crack data in the region to be analyzed according to the fault plane;
intersecting the cave data, the fracture data and the fracture data with the fracture zone profile and taking an intersection to respectively obtain a cave data body, a fracture data body and a fracture data body within the fracture zone profile range;
and fifthly, rescaling the cave data body, the fracture data body and the fracture data body in the fracture zone contour range according to respective value range, and fusing the cave data body, the fracture data body and the fracture data body into a data body to obtain the fracture zone three-dimensional space carving of the zone to be analyzed.
According to an embodiment of the invention, in the first step, the envelope curve of the maximum probability fracture-dense region of the region to be analyzed is determined through the maximum likelihood attribute, and the fracture zone boundary of the region to be analyzed is obtained.
According to an embodiment of the present invention, in the first step, the fracture zone boundary is smoothed in the line direction, the track direction and the time by using the structure-oriented filtering, so as to obtain the fracture zone profile of the region to be analyzed.
According to one embodiment of the present invention, in the second step,
and intersecting wave impedance data smaller than a preset wave impedance threshold with tensor data to determine cave data in the area to be analyzed.
According to an embodiment of the present invention, in the second step, the step of determining the preset wave impedance threshold includes:
determining the lower limit of the porosity of the cave reservoir of the area to be analyzed according to the logging data;
determining a corresponding lower threshold value of logging wave impedance according to the porosity lower limit of the cave reservoir;
and determining the corresponding seismic scale wave impedance threshold value through forward modeling analysis according to the lower threshold value of the logging wave impedance, so as to obtain the preset wave impedance threshold value.
According to one embodiment of the present invention, in the third step,
extracting a fault plane of which the AFE value is larger than a first preset AFE threshold value, determining the extracted fault plane as a fracture space area, and correspondingly obtaining fracture data;
and extracting a fault plane of which the AFE value is less than or equal to a first preset AFE threshold value and greater than a second preset AFE threshold value, determining the extracted fault plane as a fracture reservoir range, and correspondingly obtaining fracture data.
According to an embodiment of the present invention, determining an AFE value corresponding to the control leakage prevention of the region to be analyzed to obtain the first preset AFE threshold value; and/or the like, and/or,
and determining an AFE value corresponding to the small crack to obtain the second preset AFE threshold value.
According to an embodiment of the invention, in the fifth step, the value range of the fractured data volume is the highest, the value range of the cave data volume is the next to the value range of the fractured data volume, and the value range of the fractured data volume is the lowest.
According to an embodiment of the invention, in the fifth step, when the cave data volume, the broken data volume and the crack data volume with the rescaled value ranges are merged, the overlapping area is merged according to the preset priority.
According to one embodiment of the invention, the priority of the fracture data volume, the cavity data volume and the fracture data volume is sequentially decreased.
The method provided by the invention realizes fracture zone carving by optimizing the sensitive property and adopting a 'three-in-one' mode. Specifically, the method comprises the steps of enabling fracture zone profiles (unitary), cave reservoir spreading (binary) and fracture and crack spreading (ternary) proportions to be in different ranges, fusing the three data into a data file, and carving by using a three-dimensional visualization technology. Experiments show that the method can achieve 100% of success rate of trunk fracture drilling, and can provide effective and reliable basis for knowing drilling track deployment.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required in the description of the embodiments or the prior art:
FIG. 1 is a schematic illustration of "fracture + strong beads", "fracture + medium strong beads", "fracture + bead aggregate", "fracture + random weak reflections" according to one embodiment of the present invention;
FIG. 2 is a schematic flow diagram of a method for engraving a fractured zone according to one embodiment of the invention;
FIG. 3 is a schematic diagram of a fracture likelihood according to an embodiment of the invention;
FIG. 4 is a schematic view of a fracture zone profile of a region to be analyzed according to one embodiment of the present invention;
FIG. 5 is a schematic diagram of a cavern in accordance with an embodiment of the present disclosure;
FIG. 6 is a schematic diagram of seismic sections and AFE properties of three fracture patterns according to one embodiment of the present invention;
fig. 7 is a schematic diagram of the result of engraving a three-dimensional space according to an embodiment of the present invention.
Detailed Description
The following detailed description will be given with reference to the accompanying drawings and examples to explain how to apply the technical means to solve the technical problems and to achieve the technical effects. It should be noted that, as long as there is no conflict, the embodiments and the features in the embodiments of the present invention may be combined with each other, and the technical solutions formed are within the scope of the present invention.
In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced without some of these specific details or with other methods described herein.
Additionally, the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions, and while a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than here.
In recent years, carbonate rock in northern regions of towers mainly uses energy bodies or inversion bodies and ant properties to carve the fractured fluid, but the northern regions are different from the regions of towers and rivers, so that the phenomena of obvious corrosion expansion are not seen, the fracture activity is weak, and the reservoir is more complex. Meanwhile, the general burial depth of the northward region exceeds 7500m, the influence of desert earth surface and two-tier igneous rock is caused, the earthquake signal-to-noise ratio is low, the energy of beads on an earthquake section is weak and is not prominent, and the recognition accuracy of a fracture internal structure is lower. Therefore, a method for finely engraving the internal structure of the walk-slip fracture with higher precision is urgently needed.
In view of the above problems in the prior art, the present invention provides a new method for engraving fractured zones, which is particularly suitable for fine engraving of internal results of slippery fractures.
Currently, most wells are drilled in the northward region to reach the empty or distance lost circulation completion, and the bottom of the well basically corresponds to the position of a good-quality reservoir. The four categories of fracture + strong beads, fracture + medium-strength beads, fracture + bead aggregate, fracture + disordered weak reflection are basically established by finely calibrating the drilling.
As shown in figure 1, bead emptying loss of 'fracture + strong beads' and 'fracture + medium-strength beads' type drilling wells on fracture belts, bead emptying loss of 'fracture + bead aggregate' drilling wells on the edges of beads, emptying loss of 'fracture + disordered weak reflection' drilling wells on weak reflection areas can determine two types of geophysical response characteristics on the whole, fracture surfaces near the strong beads and in the weak reflection areas are high-quality reservoir concentrated areas, and blank weak reflection occurring on the top surface reflection of a room group can be determined as the position of a fracture surface basically.
Therefore, the internal structure of the northward fracture zone needs to be respectively carved by two types of bead phases and fracture phases, and analysis shows that the bead phases reflect karst cave, fracture opposite mapping fracture and fracture, the karst cave and fracture space correspond to a cave-type reservoir and form the fracture cave-type reservoir by combining with the small fractures.
Based on the analysis, the fracture zone carving method provided by the invention provides a 'three-in-one' fracture zone carving idea by optimizing the sensitive property, so that the space distribution of the three-dimensional carved fracture zone is realized.
Fig. 2 shows a schematic implementation flow diagram of the fracture zone engraving method provided by the present embodiment.
As shown in fig. 2, in the method for engraving a fracture zone according to this embodiment, a fracture zone boundary of a region to be analyzed is first identified in step S201, so as to determine a fracture zone boundary of the region to be analyzed. Subsequently, the method obtains a fracture zone profile of the region to be analyzed according to the fracture zone boundary of the region to be analyzed in step S202.
Specifically, in this embodiment, in step S201, the method preferably determines the envelope of the most probable fracture-dense region of the region to be analyzed through the maximum likelihood attribute, and further obtains the fracture zone boundary of the region to be analyzed.
The boundary identification of the fracture zone is an important reference factor for well position deployment and resource amount calculation. Analysis shows that whether the internal homogeneity degree of the fracture zone is high or low, the seismic section has difference, and weak reflection fault or local strong abnormity can be caused.
The method provided by the embodiment utilizes a maximum likelihood post-stack fracture detection technology, calculates the similarity between each sampling point in a specified inclination angle and trend range, obtains the probability of local fault and fracture development through similarity interception, and further generates a fault (fracture) trend, a fault (fracture) trend and a fault (fracture) maximum likelihood, so as to achieve the purpose of identifying and describing a fracture zone.
As shown in fig. 3, in this embodiment, an inclination angle body and an azimuth angle body may be obtained by calculation along a certain trend and inclination, and fracture likelihoods including Fault _ likelihood (light yellow) and thin _ Fault _ likelihood (green) may be obtained by optimization calculation, and have features of trend, inclination, and level. The three-dimensional fracture detection method is characterized in that the three-dimensional fracture detection method comprises the following steps of detecting the possible fracture information, and reflecting the outline of the fracture-dense area by the aid of the Fault _ likeliohood.
In response to the problem that the contour obtained through the above process is too irregular, the method provided in this embodiment preferably further applies structure-oriented filtering in step S202 to smooth the broken-zone boundary obtained in step S201 in the line direction, the track direction and the time, so as to obtain the broken-zone contour of the region to be analyzed. As shown in fig. 4, the profile of the fracture zone obtained at this time can better reflect the profile characteristics.
For the interior of the fracture zone, the present method preferably further subdivides the fracture-cavern reservoir into three types, namely: karst cave, fractures and cracks. As can be seen by combining the geophysical response characteristics of the typical well shown in figure 1, the three types of reservoirs are characterized by beads, staggered sections and disordered reflection on a seismic section, and the method uses the respective sensitive attributes to carve the spatial forms of the three types of reservoirs.
Specifically, as shown in fig. 2, in this embodiment, the method further determines cavern data in the area to be analyzed through wave impedance inversion in step S203.
Through research, when the volume is drawn by adopting wave impedance inversion, lower beads partially hidden in a strong interface are not easy to identify. If the beads are carved by the same threshold value, the strong interface is often carved out. In order to solve this problem, the method provided by this embodiment preferably identifies the small-scale cavern information in the strong interface by means of intersection of wave impedance and tensor attribute.
In this embodiment, the method preferably uses intersection of wave impedance data smaller than a preset wave impedance threshold and tensor attribute data, and a region satisfying both a strong tensor value and a low wave impedance (i.e., wave impedance is smaller than the preset wave impedance threshold and tensor value is greater than the preset tensor threshold) can be detected through intersection of the tensor attribute and the wave impedance data, so that small-scale karst cave information in the strong reflection interface can be effectively identified, as shown in fig. 5.
In this embodiment, the specific value of the preset tension value is preferably determined on the basis of enabling the slot-hole shape to be more prominent, and is preferably determined within the range of [30,50 ].
In this embodiment, in the process of determining the preset wave impedance threshold, the method preferably determines a cavern reservoir porosity lower limit of the region to be analyzed according to the logging data, then determines a corresponding logging wave impedance lower threshold according to the obtained cavern reservoir porosity lower limit, and then determines a corresponding seismic scale wave impedance lower threshold through forward model analysis according to the logging wave impedance lower threshold, so as to obtain the required preset wave impedance threshold.
For example, the cavern reservoir porosity threshold can be determined to be 5.2% from log data. The speed of pure limestone is 6100m/s, and the density is 2.7g/cm3From the time-averaged equation, the theoretical wave impedance corresponding to 5.2% of the porosity threshold was calculated to be 1.365 x 107m/s·kg/m3. However, since the frequency bandwidth of the seismic data (inversion data) is only about 0-80Hz, the wave impedance cannot be directly used as a threshold value during inversion carving, 80Hz high-cut filtering needs to be performed on the log curve, and the forward theoretical wave impedance of the obtained corresponding model is 1.365 x 107m/s·kg/m3. On a 10m cavern geologic model, the corresponding wave impedance on a seismic scale is 1.52 x 107m/s·kg/m3. Thus, 1.52 x 10 may be used7m/··kg/m3Threshold as wave impedance inversion (i.e. preset wave impedance threshold)
In this embodiment, the geostatistical inversion is adopted, so that the resolution is significantly improved, and the high-cut filtering range is also expanded (specifically, the high-cut filtering range needs to be adjusted according to the frequency bandwidth of the inversion wave impedance body).
Of course, in other embodiments of the present invention, the method may also use other reasonable ways to determine the cavern data of the area to be analyzed through wave impedance inversion, and the present invention is not limited thereto.
As shown in fig. 2 again, in this embodiment, the method further extracts a fault plane in the region to be analyzed by using an automatic fault extraction technique in the assistance S204, and determines fracture data and fracture data in the region to be analyzed according to the fault plane obtained in the step S204 in the step S205.
From the statistics of actual drilling, most wells in the northward region are emptied and lost on the main fracture surface in the fracture zone, but many main fracture surfaces are difficult to identify by naked eyes on the section, so that a geophysical mathematical operation method is required to be used for distinguishing. In this embodiment, the method preferably uses an Automatic Fault Extraction (AFE) method to process the three-dimensional seismic coherence volume data or the discontinuous attribute volume data (two steps of linear enhancement and fault enhancement) to automatically extract fault lines, so as to obtain a fault plane.
By adopting the automatic fault extraction technology, the fracture data and the Raynaud data obtained by the method can clearly reflect the fracture information, and compared with the coherent attribute, the transverse resolution is higher, and the method is better in fit with a drilling emptying loss section.
FIG. 6 shows seismic profiles and AFE properties of three fracture patterns, each labeled drilling loss location. It can be seen from the seismic profile shown in fig. 6 that the fracture characteristics are not obvious, the position of a main fracture surface with the width of a few meters is difficult to determine, and the centers of beads are not all favorable development zones of the reservoir. The left side and the right side are sections of the AFE attribute corresponding to the seismic section, and the position of a fracture surface depicted by the AFE attribute is well matched with the position of the atmospheric leakage when viewed from the section.
In this embodiment, in step S205, the method preferably extracts a fault plane having an AFE value greater than a first preset AFE threshold, and determines the extracted fault plane as a fracture space region, so as to obtain fracture data correspondingly. The method also extracts fault planes with AFE values smaller than or equal to a first preset AFE threshold value and larger than a second preset AFE threshold value, determines the extracted fault planes as a fracture reservoir range, and correspondingly obtains fracture data.
In this embodiment, the fracture threshold (i.e., the first predetermined AFE threshold) and the fracture threshold (i.e., the second predetermined AFE threshold) are obtained mainly by drilling calibration.
For example, counting 12 holes drilled in the No. 1 fracture zone and the No. 5 fracture zone in the northward direction at present, taking an AFE value larger than 201.2 as a fracture cavity area, matching 80% of drilling air loss, comparing fine fractures with AFE attributes through imaging logging and logging information, matching oil and gas development positions with AFE fractures well, selecting 28 as a threshold value through calibration, and selecting a fracture reservoir range between 201.2 and 28.
Of course, in other embodiments of the present invention, the method may also determine the fracture data and the crack data in the region to be analyzed in other reasonable manners according to actual needs, and the present invention is not limited thereto.
As shown in fig. 2, in this embodiment, after obtaining the fracture zone profile, the cave data, the fracture data, and the fracture data of the region to be analyzed, in step S206, the method intersects the cave data, the fracture data, and the fracture data with the fracture zone profile and takes an intersection, so as to obtain a cave data volume, a fracture data volume, and a fracture data volume within the fracture zone profile range, respectively.
Then, in step S207, the method rescales the cave data volume, the fracture data volume, and the fracture data volume within the fracture zone profile range according to respective value range ranges, and fuses the cave data volume, the fracture data volume, and the fracture data volume into one data volume, so as to obtain the fracture zone three-dimensional space carving of the region to be analyzed.
In this embodiment, preferably, the range of the fracture data volume is the highest, the range of the cave data volume is the second highest, and the range of the fracture data volume is the lowest. When the cave data body, the fracture data body and the crack data body after the value range is rescaled are merged, the data are likely to have overlapping areas, and for the overlapping areas, the method preferably merges the data according to preset priority.
Specifically, in this embodiment, for the data overlap region, the priorities of the fractured data volume, the cavern data volume, and the fractured data volume are sequentially decreased. That is, when the fractured data volume and the cave data volume are overlapped, the method takes the fractured data volume as the valid data volume; when the cave data body and the crack data body are overlapped, the method takes the cave data body as an effective data body.
It should be noted that, in this embodiment, according to actual needs, the method may further perform color coding according to the range of the different data in step S207, so as to obtain a colored three-dimensional space carving of the fracture zone, as shown in fig. 7.
Meanwhile, it should be noted that the present invention does not limit the order of determining the fracture zone profile, the cave data, the fracture data and the fracture data. In other embodiments of the present invention, according to actual needs, the specific order of the parameters may also be determined by using other reasonable orders.
As can be seen from the above description, the method provided by the present invention implements fracture zone engraving in a "three-in-one" manner by optimizing the sensitivity attribute. Specifically, the method comprises the steps of enabling fracture zone profiles (unitary), cave reservoir distribution (binary) and fracture and crack distribution (ternary) proportions to be in different ranges, fusing the three data into a data file, and carving by using a three-dimensional visualization technology. Experiments show that the success rate of trunk fracture drilling can be 100% by using the method, and effective and reliable bases can be provided for knowing the deployment of drilling tracks.
It is to be understood that the disclosed embodiments of this invention are not limited to the particular structures or process steps disclosed herein but are extended to equivalents thereof as would be understood by those ordinarily skilled in the relevant arts. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" or "an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
While the foregoing examples have been provided to illustrate the principles of the invention in one or more applications, it will be apparent to those skilled in the art that various changes in form, usage and details of implementation can be made without departing from the principles and concepts of the invention. Accordingly, the invention is defined by the appended claims.
Claims (8)
1. A method of fracture zone engraving, the method comprising:
firstly, identifying the boundary of a fracture zone of a region to be analyzed, determining the boundary of the fracture zone of the region to be analyzed, and further obtaining the contour of the fracture zone of the region to be analyzed; determining an envelope line of a maximum probability fracture concentration area of the area to be analyzed through a maximum likelihood attribute to obtain a fracture zone boundary of the area to be analyzed;
step two, determining cave data in the area to be analyzed through wave impedance inversion;
extracting a fault plane in the region to be analyzed by utilizing a fault automatic extraction technology, and determining fracture data and crack data in the region to be analyzed according to the fault plane;
intersecting the cave data, the fracture data and the fracture data with the fracture zone profile and taking an intersection to respectively obtain a cave data volume, a fracture data volume and a fracture data volume within the fracture zone profile range;
fifthly, rescaling the cave data body, the fracture data body and the fracture data body in the fracture zone profile range according to respective value range, and fusing the cave data body, the fracture data body and the fracture data body into a data body to obtain fracture zone three-dimensional space carving of the zone to be analyzed; the fracture data volume has the highest value range, the cave data volume has the next value range, and the fracture data volume has the lowest value range.
2. The method of claim 1, wherein in step one, the fracture zone boundary is smoothed in line direction, track direction and time by using texture oriented filtering to obtain the fracture zone profile of the region to be analyzed.
3. The method according to claim 1 or 2, wherein, in the second step,
and intersecting wave impedance data smaller than a preset wave impedance threshold with tensor data to determine cave data in the area to be analyzed.
4. The method according to claim 3, wherein in step two, the step of determining the preset wave impedance threshold comprises:
determining the lower limit of the porosity of the cave reservoir of the area to be analyzed according to the logging data;
determining a corresponding lower threshold value of logging wave impedance according to the lower limit of the porosity of the cavern reservoir;
and determining the corresponding seismic scale wave impedance threshold value through forward modeling analysis according to the lower threshold value of the logging wave impedance, so as to obtain the preset wave impedance threshold value.
5. The method of claim 1, wherein, in step three,
extracting a fault plane of which the AFE value is larger than a first preset AFE threshold value, determining the extracted fault plane as a fracture space region, and correspondingly obtaining fracture data;
and extracting fault planes with AFE values smaller than or equal to a first preset AFE threshold value and larger than a second preset AFE threshold value, determining the extracted fault planes as fracture reservoir ranges, and correspondingly obtaining fracture data.
6. The method of claim 5,
determining an AFE value corresponding to the control leakage prevention of the region to be analyzed to obtain the first preset AFE threshold value; and/or the like, and/or,
and determining an AFE value corresponding to the small crack to obtain the second preset AFE threshold value.
7. The method according to claim 1, wherein in the fifth step, when the cave data volume, the fractured data volume and the fractured data volume after the value range rescaling are merged, the overlapped area is merged according to the preset priority.
8. The method of claim 7, wherein the fracture data volume, cavern data volume, and fracture data volume are sequentially decreasing in priority.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910308488.7A CN111830558B (en) | 2019-04-17 | 2019-04-17 | Fracture zone engraving method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910308488.7A CN111830558B (en) | 2019-04-17 | 2019-04-17 | Fracture zone engraving method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111830558A CN111830558A (en) | 2020-10-27 |
| CN111830558B true CN111830558B (en) | 2022-07-22 |
Family
ID=72914877
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910308488.7A Active CN111830558B (en) | 2019-04-17 | 2019-04-17 | Fracture zone engraving method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111830558B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112698397B (en) * | 2020-11-20 | 2023-08-22 | 中国石油天然气股份有限公司 | Method for describing reservoir stratum with sliding fracture cavity in basin area |
| CN113156499B (en) * | 2021-02-26 | 2023-05-26 | 中国石油天然气股份有限公司 | Seismic data post-stack quantitative prediction method for fractured reservoir in basin area |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104818982A (en) * | 2015-04-20 | 2015-08-05 | 中国石油天然气股份有限公司 | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir |
| CN107390264A (en) * | 2017-07-20 | 2017-11-24 | 中国石油化工股份有限公司 | The characterizing method of the disconnected solution internal structure of carbonate rock |
| CN108241170A (en) * | 2016-12-27 | 2018-07-03 | 中国石油天然气股份有限公司 | Method and device for fitting carved volume and dynamic reserve of fractured-vuggy carbonate reservoir |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006056982A2 (en) * | 2004-11-24 | 2006-06-01 | Paradigm Geophysical Ltd. | System and method for fault identification |
| CN109583016B (en) * | 2018-10-22 | 2023-04-07 | 恒泰艾普集团股份有限公司 | Method for determining and quantitatively calculating space geometric form of fracture-cavity body |
-
2019
- 2019-04-17 CN CN201910308488.7A patent/CN111830558B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104818982A (en) * | 2015-04-20 | 2015-08-05 | 中国石油天然气股份有限公司 | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir |
| CN108241170A (en) * | 2016-12-27 | 2018-07-03 | 中国石油天然气股份有限公司 | Method and device for fitting carved volume and dynamic reserve of fractured-vuggy carbonate reservoir |
| CN107390264A (en) * | 2017-07-20 | 2017-11-24 | 中国石油化工股份有限公司 | The characterizing method of the disconnected solution internal structure of carbonate rock |
Non-Patent Citations (2)
| Title |
|---|
| 塔里木盆地走滑带碳酸盐岩断裂相特征及其与油气关系;邬光辉等;《地质学报》;20120229;第86卷(第2期);第219-227页 * |
| 多尺度缝洞型碳酸盐岩油藏不确定性建模方法;邓晓娟等;《石油学报》;20180930;第39卷(第9期);第1054-1058页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111830558A (en) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103529475B (en) | A kind of method identified and explain carbonate paleokarst reservoir three dimensional structure | |
| CN109116428B (en) | Fracture-cavity carbonate reservoir uncertainty modeling method and device | |
| CN106703779A (en) | Injection-production well pattern construction method suitable for carbonate fractured-cave reservoirs | |
| CN109541685B (en) | River channel sand body identification method | |
| CN111722277B (en) | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration | |
| CN109655896A (en) | The method and system that ancient stream channel identifies and portrays | |
| CN111611673B (en) | Modeling method for carbonate reservoir ancient underground river type reservoir | |
| CN112698392A (en) | Multilayer system reservoir three-dimensional drawing and oil-gas spatial distribution and constant volume method | |
| CN110984951B (en) | Shale oil development well pattern deployment method | |
| CN111830558B (en) | Fracture zone engraving method | |
| CN112946782B (en) | Earthquake fine depicting method for dense oil-gas storage seepage body | |
| CN112489208A (en) | Crack sheet extraction method based on ant algorithm and three-dimensional geological model construction method | |
| CN111190228B (en) | Submarine fan water channel reservoir configuration identification method | |
| CN112114358B (en) | Underground volcanic channel identification method based on three-dimensional seismic data representation | |
| CN110794477B (en) | Method for identifying overtaking stratum pinch-out line by combining well and seismic | |
| CN115857047B (en) | Comprehensive prediction method for earthquake reservoir | |
| CN107329168B (en) | Sedimentary body lithology identification method and system | |
| CN112433248B (en) | Method for detecting hidden reservoir stratum in carbonate rock deposition environment | |
| CN110764161B (en) | Comprehensive method for judging and identifying fractured and broken zone of carbonate rock through drilling data | |
| CN113514884A (en) | Compact sandstone reservoir prediction method | |
| CN114545498A (en) | Comprehensive depicting method for ancient karst fracture-cave structure | |
| CN115407407A (en) | Three-dimensional geological model construction method for carbonate rock ancient karst cave and filling thereof | |
| CN115387785A (en) | Sea-facies carbonate-cuttings limestone reservoir high-permeability strip identification method and device | |
| CN106353813A (en) | Method for identifying fluid properties based on array acoustic logging | |
| CN115032712A (en) | Method for rapidly judging development type of reservoir space under compact lithology background |
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 |