Cave geology scale volume calculation method and device, electronic equipment and medium

Abstract

The application discloses a karst cave geological scale apparent volume calculation method and device, electronic equipment and medium. The method can include: performing karst cave plane boundary correction and description, calculating the maximum bottom projection area of the karst cave; determining the height of the karst cave based on wave impedance; and calculating the apparent volume of the karst cave according to the projection area and the height of the karst cave. The application can quickly and effectively calculate the geological scale apparent volume of the karst cave based on the karst cave column body approximate estimation method, and provides a basis for more accurate calculation of the static reserves of the fracture-cave reservoir.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of oil and gas geophysical exploration, and more particularly, to a cave body geological scale apparent volume calculation method and device, electronic equipment and medium. BACKGROUND

[0002] The Ordovician system of Tahe oilfield is a typical karst fracture-cave type reservoir, with various karst types, developed caves, pores and fracture-caves, deep burial (more than 5300 meters), strong heterogeneity, large scale variation and complex connectivity between fracture-cave units. The cave is the main reservoir space of the fracture-cave type reservoir, and the development scale, development morphology and combination mode of the cave are complex and diverse. Domestic scholars have done a lot of research on how to finely depict the cave, and have proposed that amplitude variation ratio, root mean square and other energy-based seismic attributes can effectively identify the development morphology of the cave, and the size and position of the cave can be determined based on forward modeling and multi-attribute fusion.

[0003] However, due to the different development scales, complex morphology and combination mode of the cave body, there is a difference between the seismic recognition scale and the actual geological scale. With the adjustment of the energy threshold value, the energy value decreases from large to small, the size of the cave recognized by the earthquake changes, and the number decreases. Uniform adjustment of the color threshold value of the energy attribute will lead to different scale variation rates of different scale caves, especially the small cave bodies appear from "yes" to "no". How to accurately determine the outline, scale and apparent volume of the cave is crucial for reserve calculation. Based on the seismic anomaly profile of the cave recognized by the earthquake, how to transform it into a geological profile close to the actual scale, that is, the "volume correction" of the cave, domestic scholars have done a lot of research. Basically, through forward modeling, an empirical relationship between the scale of the forward model and the scale of the seismic anomaly body is established, and height, width and volume correction coefficients are used to calculate the corrected volume. However, the correction coefficients are different for caves of different scales, making it difficult to land in actual production.

[0004] At present, a cave body geological scale apparent volume calculation method needs to be developed. According to the field requirements, relying on commonly used interpretation software, a fast and effective field practical process for calculating the volume of the cave body is established, which is of great significance for the reserve calculation and development plan formulation of the fracture-cave reservoir.

[0005] The information disclosed in the background section of this application is only intended to deepen the understanding of the general background of the application, and should not be regarded as acknowledging or implying in any form that the information constitutes prior art known to those skilled in the art. SUMMARY

[0006] The present application provides a cave geological scale apparent volume calculation method and device, electronic equipment and medium, which can quickly and effectively calculate the geological scale apparent volume of the cave body based on the cave body column approximation estimation method. Relying on the seismic interpretation software platform, the column approximation estimation method practical process suitable for field popularization and application is formed through the step-by-step correction idea innovation for different types of caves, which can efficiently complete the correction of the apparent volume of the cave and the accurate estimation of the static reserve.

[0007] In the first aspect, the present application provides a cave geological scale apparent volume calculation method, comprising:

[0008] correcting and depicting the cave plane boundary, and calculating the maximum bottom projection area of the cave body;

[0009] determining the cave height based on the wave impedance;

[0010] calculating the apparent volume of the cave body according to the projection area and the cave height.

[0011] Preferably, the step of correcting and depicting the cave plane boundary and calculating the maximum bottom projection area of the cave body comprises:

[0012] determining the advantage attribute of the cave shape and boundary identification for different geological targets of the cave identification;

[0013] determining the seismic anomaly shape and plane boundary of the cave based on the amplitude energy attribute and the coherent attribute, and depicting the seismic anomaly boundary on the plane;

[0014] based on the theoretical basis of forward simulation and multi-attribute superposition, adjusting the energy attribute value range of the cave identification, correcting the cave seismic anomaly boundary to the geological anomaly boundary, depicting the boundary contour, and calculating the maximum bottom projection area of the cave body.

[0015] Preferably, the step of adjusting the energy attribute value range of the cave identification comprises:

[0016] adjusting the energy attribute value range of the cave identification to 60% of the maximum value.

[0017] Preferably, the step of determining the cave height based on the wave impedance comprises:

[0018] identifying the cave shape based on the original seismic data and the amplitude energy attribute;

[0019] determining the cave seismic anomaly body based on the root mean square amplitude anomaly;

[0020] quantitatively correcting the height of the cave body by using the wave impedance inversion technology, combining the impedance threshold value of the typical cave reservoir body calibrated by well logging interpretation and drilling, determining the longitudinal height of the cave body, and then correcting the seismic anomaly height of the cave body to the geological height to determine the cave height.

[0021] Preferably, calculating the cave apparent volume according to the projection area and the cave height comprises:

[0022] Multiplying the projection area and the cave height, the cave apparent volume is calculated.

[0023] Preferably, it further comprises:

[0024] According to the cave apparent volume and the average porosity of the work area fracture-cave reservoir, the effective volume of the cave is determined, and then the static reserves of the cave geological body are calculated according to the reservoir parameters.

[0025] Preferably, it further comprises:

[0026] Based on the cave plane boundary correction and delineation of multi-attribute and forward simulation, the maximum bottom projection area of the cave before and after correction is calculated respectively;

[0027] According to the multiplication of the maximum bottom projection area before and after correction and the cave height, the apparent volume of the cave before and after correction is calculated;

[0028] According to the apparent volume of the cave before and after correction, the correction rate is calculated at different scales;

[0029] According to the correction rate, the uncorrected cave apparent volume of the corresponding scale is corrected.

[0030] Preferably, the correction rate is:

[0031] The correction rate=(1-corrected apparent volume / corrected apparent volume)×100%.

[0032] In the second aspect, the embodiments of the present disclosure also provide a cave body geological scale apparent volume calculation device, comprising:

[0033] The area calculation module corrects and delineates the cave plane boundary, and calculates the maximum bottom projection area of the cave body;

[0034] The height calculation module determines the cave height based on wave impedance;

[0035] The volume calculation module calculates the cave apparent volume according to the projection area and the cave height.

[0036] Preferably, the correction and delineation of the cave plane boundary and the calculation of the maximum bottom projection area of the cave body comprise:

[0037] The advantage attribute of cave shape and boundary identification is determined for different geological targets of cave identification;

[0038] The seismic anomaly shape and plane boundary of the cave are determined based on the amplitude energy attribute and the coherent attribute, and the seismic anomaly boundary is delineated on the plane.

[0039] The energy attribute value range for karst cave identification is adjusted based on the theoretical basis of forward modeling and multi-attribute stacking, the boundary of the karst cave seismic anomaly is corrected to the boundary of the geological anomaly, the boundary profile is depicted, and the maximum bottom projection area of the karst cave body is calculated.

[0040] Preferably, adjusting the energy attribute value range for karst cave identification comprises:

[0041] The energy attribute value range for karst cave identification is adjusted to 60% of the maximum value.

[0042] Preferably, determining the height of the karst cave based on wave impedance comprises:

[0043] Identifying the shape of the karst cave based on original seismic data and amplitude energy attributes;

[0044] Determining the karst cave seismic anomaly body based on the root mean square amplitude anomaly;

[0045] Using wave impedance inversion technology, the height of the karst cave body is quantitatively corrected, the impedance threshold value of the typical karst cave reservoir is determined by combining well interpretation and drilling calibration, the longitudinal height of the karst cave body is determined, and then the height of the karst cave body seismic anomaly is corrected to the geological height, and the height of the karst cave is determined.

[0046] Preferably, calculating the apparent volume of the karst cave body according to the projection area and the height of the karst cave comprises:

[0047] The projection area and the height of the karst cave are multiplied to calculate the apparent volume of the karst cave.

[0048] Preferably, it further comprises:

[0049] According to the apparent volume of the karst cave and the average value of the porosity of the fracture-cave reservoir in the work area, the effective volume of the karst cave is determined, and then the static reserves of the karst geological body are calculated according to the reservoir parameters.

[0050] Preferably, it further comprises:

[0051] Based on the plane boundary correction and depiction of the karst cave by multi-attribute and forward modeling, the maximum bottom projection area of the karst cave body before and after correction is calculated respectively;

[0052] According to the multiplication of the maximum bottom projection area before and after correction and the height of the karst cave, the apparent volume of the karst cave before and after correction is calculated;

[0053] According to the apparent volume of the karst cave before and after correction, the correction rate is calculated at different scales;

[0054] According to the correction rate, the uncorrected apparent volume of the karst cave body corresponding to the scale is corrected.

[0055] Preferably, the correction rate is:

[0056] The correction rate = (1 - the volume after correction / the volume before correction) * 100%.

[0057] In a third aspect, the embodiments of the present disclosure further provide an electronic device, which comprises:

[0058] A memory, which stores executable instructions;

[0059] A processor, which runs the executable instructions in the memory to realize the cavern geologic scale apparent volume calculation method.

[0060] In a fourth aspect, the embodiments of the present disclosure further provide a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the cavern geologic scale apparent volume calculation method.

[0061] Its beneficial effects are:

[0062] (1) Based on the preferred cavern recognition, the seismic anomaly boundary is determined by the seismic sensitive multi-attribute, the maximum energy is adjusted to fill the seismic anomaly boundary, then the maximum energy threshold is adjusted to 60% based on the forward simulation, the cavern geologic boundary is determined, and the theoretical basis of the fracture-cave body forward simulation width quantitative analysis is applied.

[0063] (2) For the description of single cavern plane profile, the energy change before and after correction is adjusted, and the explanation software polgon is manually carved to complete, that is, the maximum bottom area of the cavern body can be quickly obtained, and the unified correction threshold value is not affected (automatic carving), and the method has the advantages of practicality and quickness.

[0064] The method and device have other characteristics and advantages, which will be apparent or will be described in detail in the accompanying drawings and subsequent detailed description incorporated herein, which are collectively used to explain the specific principles of the present application. BRIEF DESCRIPTION OF DRAWINGS

[0065] The above and other objects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings, in which like reference characters refer to like parts throughout the figures, and wherein:

[0066] Figure 1 A flow chart showing the steps of the cavern geologic scale apparent volume calculation method according to one embodiment of the present application is shown.

[0067] Figure 2 A schematic diagram showing the technical flow of the cylinder approximation estimation method according to one embodiment of the present application is shown.

[0068] Figure 3 A schematic diagram of a cave plan profile correction process according to an embodiment of the present application is shown.

[0069] Figure 4 A schematic diagram of a single type cave volume distribution correction and calculation process according to an embodiment of the present application is shown.

[0070] Figure 5 A schematic diagram of a composite type cave volume correction and calculation process according to an embodiment of the present application is shown.

[0071] Figure 6 A block diagram of a cave volume geological scale apparent volume calculation device according to an embodiment of the present application is shown.

[0072] BRIEF DESCRIPTION OF DRAWINGS

[0073] 201, area calculation module; 202, height calculation module; 203, volume calculation module. DETAILED DESCRIPTION

[0074] The preferred embodiments of the present application will be described in more detail below. Although the preferred embodiments of the present application are described below, it is understood that the present application can be implemented in various forms and should not be limited by the embodiments set forth herein.

[0075] In order to understand the scheme and effects of the embodiments of the present application, four specific application examples are given below. Those skilled in the art should understand that the examples are only for the convenience of understanding the present application, and any specific details thereof are not intended to limit the present application in any way.

[0076] Example 1

[0077] Figure 1 A flowchart of the steps of a cave volume geological scale apparent volume calculation method according to an embodiment of the present application is shown.

[0078] As Figure 1 shown, the cave volume geological scale apparent volume calculation method includes: step 101, performing cave plan boundary correction and delineation, and calculating the maximum bottom surface projection area of the cave volume; step 102, determining the cave height based on wave impedance; and step 103, calculating the apparent volume of the cave volume according to the projection area and the cave height.

[0079] In one example, performing cave plan boundary correction and delineation, and calculating the maximum bottom surface projection area of the cave volume includes:

[0080] Determining the advantageous attributes of cave morphology and boundary identification for different geological targets identified by the cave;

[0081] determine the seismic anomaly shape and planar boundary of the cave based on the amplitude energy attribute and the coherent attribute, and delineate the seismic anomaly boundary of the cave on the planar surface;

[0082] based on the theoretical basis of forward modeling and the multi-attribute superposition, adjust the energy attribute value range of the cave identification, correct the cave seismic anomaly boundary to the geological anomaly boundary, delineate the boundary profile, and calculate the maximum bottom surface projection area of the cave body.

[0083] In one example, adjusting the energy attribute value range of the cave identification includes:

[0084] adjusting the energy attribute value range of the cave identification to 60% of the maximum value.

[0085] In one example, determining the cave height based on wave impedance includes:

[0086] based on the original seismic data and the amplitude energy attribute, identifying the cave shape;

[0087] based on the root mean square amplitude anomaly, determining the cave seismic anomaly body;

[0088] Using wave impedance inversion technology, the height of the cave body is quantitatively corrected, combined with the impedance threshold value of the typical cave reservoir body calibrated by well logging interpretation and drilling, the longitudinal height of the cave body is determined, and then the cave body seismic anomaly height is corrected to the geological height, and the cave height is determined.

[0089] In one example, according to the projection area and the cave height, the apparent volume of the cave body is calculated.

[0090] The projection area is multiplied by the cave height to calculate the apparent volume of the cave body.

[0091] In one example, it also includes:

[0092] According to the apparent volume of the cave body and the average value of the porosity of the fracture-cave reservoir in the work area, the effective volume of the cave is determined, and then the static reserves of the cave geological body are calculated according to the reservoir parameters.

[0093] In one example, it also includes:

[0094] Based on the cave planar boundary correction and delineation of multi-attribute and forward modeling, the maximum bottom surface projection area of the cave body before and after correction is calculated respectively;

[0095] According to the multiplication of the maximum bottom surface projection area before and after correction and the cave height, the apparent volume of the cave before and after correction is calculated;

[0096] According to the apparent volume of the cave before and after correction, the correction rate is calculated in different scales;

[0097] According to the correction rate, the uncorrected cave apparent volume of the corresponding scale is corrected.

[0098] In one example, the correction rate is:

[0099] The correction rate = (1 - corrected post-visual volume / corrected pre-visual volume) x 100%.

[0100] Specifically, the present application is based on the theoretical basis of fracture-vug width forward modeling, and relies on the interpretation software to realize the geological scale calculation of fracture-vug reservoirs by the columnar approximation estimation method of multiplying the bottom area by the height, thereby providing a fast and effective method for more accurate reserve calculation.

[0101] Based on multi-attribute and forward modeling, the cave plane boundary correction and depiction are carried out. For different geological targets of cave identification, the advantage attributes of cave shape and boundary identification are determined, the seismic anomaly shape and plane boundary of the cave are determined based on amplitude energy attribute and coherent attribute, and the seismic anomaly boundary is depicted on the plane based on the interpretation software; on this basis, based on the theoretical basis of forward modeling and multi-attribute superposition, the energy attribute value range of cave identification is adjusted to 60% of the maximum value, so as to correct the cave seismic anomaly boundary to the geological anomaly boundary, and the boundary profile is manually depicted, and the maximum bottom projection area of the cave body before and after correction is calculated. For cave identification, the amplitude change rate and root mean square energy attribute are better, and the coherent attribute is better for cave boundary identification.

[0102] Based on wave impedance, the height of the cave is determined. For the identification of the shape of the cave, the original seismic data and the amplitude energy attribute are mainly used, and the root mean square amplitude anomaly is determined as the seismic anomaly body of the cave, which is difficult to determine the actual height of the cave reservoir. In production, the height of the cave body is quantitatively corrected, and the height is mainly predicted by using the wave impedance inversion technology. Since the impedance inversion technology can eliminate the wavelet sidelobes and improve the longitudinal resolution, compared with other attributes, the impedance inversion is closer to the geological model height; combined with well logging interpretation and drilling calibration, the impedance threshold value of the typical cave reservoir is determined to determine the longitudinal height of the cave body, and then the seismic anomaly height of the cave body is corrected to the geological height.

[0103] Based on the columnar approximation estimation method, the apparent volume of the cave body is determined. Since the cave reservoir is approximately a column in three-dimensional space, the layer attribute of the energy attribute in the development time window of the cave body is extracted to determine the maximum bottom area of the cave, and then the corrected height of the cave body is obtained by the impedance body calibration, and then the approximate estimation method of multiplying the bottom area by the height can quickly calculate the geological apparent volume of the cave reservoir.

[0104] At the same time, the apparent volume of the cave reservoir before and after correction can be compared and the correction rate can be calculated. Since the correction rates of different scales of caves are different, according to the shape and scale of the cave in the research, the correction is classified and graded:

[0105] For single-hole, only different scale correction rates are needed, combined with forward simulation seismic identification results, it can be classified into greater than 100m, 50-100m, and less than 50m;

[0106] For complex hole, it is composed of two or more caves, before correction, seismic identification is a large energy group, and there is energy change in the local part, after correction, it can be obviously distinguished into two or several caves, if the scales of the caves constituting the complex hole are basically the same, the threshold value is uniformly adjusted and corrected according to the above method; if the scales of the caves constituting the complex hole are quite different, uniform adjustment of the energy threshold value will lead to disappearance of the small hole, and the small hole needs to be corrected separately, and then the correction rate is counted according to the scale.

[0107] Based on the advantage of the discovered well, the correction rate range of the dissolution hole of different types and different scales is counted, the geological abnormal volume of the hole is quickly estimated according to the classification of the dissolution hole identified by the seismic volume in the later period, and the basis for the accurate and efficient estimation of the reserve is provided.

[0108] According to the average value of the porosity of the fracture-cave reservoir in the working area, the effective volume of the corrected dissolution reservoir is obtained, and the static reserve of the dissolution geological body is estimated according to the reservoir parameters, so as to provide a practical calculation process for the reserve calculation of the fracture-cave reservoir in the working area.

[0109] The method first corrects the out-of-plane contour of the dissolution hole according to the seismic abnormal boundary of the dissolution hole identified by the seismic multi-attribute identification and the forward simulation, and manually describes the dissolution hole boundary relying on the seismic interpretation software; the height of the dissolution hole identified by the seismic attribute is corrected based on the wave impedance inversion, and on this basis, the geological apparent volume of the dissolution hole is estimated by the column approximation method, that is, the maximum bottom area of the corrected dissolution hole is multiplied by the height of the dissolution hole to obtain the geological apparent volume of the dissolution hole, and then the static reserve is calculated, and the effectiveness of the method is verified by comparing the dynamic reserve of the typical dissolution hole with the constant volume well. The method is simple to operate, cost-saving, fast and effective in serving production, and can effectively support the efficient development of the carbonate fracture-cave reservoir.

[0110] Example 2

[0111] The application also provides a dissolution hole geological scale apparent volume calculation device, which comprises:

[0112] The area calculation module corrects and describes the plane boundary of the dissolution hole, and calculates the maximum bottom projection area of the dissolution hole;

[0113] The height calculation module determines the height of the dissolution hole based on the wave impedance;

[0114] The volume calculation module calculates the apparent volume of the dissolution hole according to the projection area and the height of the dissolution hole.

[0115] In one example, the karst cave planar boundary correction and delineation is performed, and the calculation of the maximum bottom projection area of the karst cave body includes:

[0116] The advantage attribute of karst cave morphology and boundary identification is determined for different geological targets of karst cave identification;

[0117] The seismic anomaly morphology and planar boundary of the karst cave are determined based on the amplitude energy attribute and the coherence attribute, and the seismic anomaly boundary is delineated on the plane;

[0118] Based on the theoretical basis of forward simulation and multi-attribute superposition, the energy attribute value range of karst cave identification is adjusted, the karst cave seismic anomaly boundary is corrected to the geological anomaly boundary, the boundary profile is delineated, and the maximum bottom projection area of the karst cave body is calculated.

[0119] In one example, adjusting the energy attribute value range of karst cave identification includes:

[0120] The energy attribute value range of karst cave identification is adjusted to 60% of the maximum value.

[0121] In one example, determining the height of the karst cave based on wave impedance includes:

[0122] Based on the original seismic data and the amplitude energy attribute, the karst cave morphology is identified;

[0123] Based on the root mean square amplitude anomaly, the karst cave seismic anomaly body is determined;

[0124] Using wave impedance inversion technology, the height of the karst cave body is quantitatively corrected, combined with the impedance threshold value of the typical karst cave reservoir body determined by well logging interpretation and drilling calibration, the longitudinal height of the karst cave body is determined, and then the karst cave body seismic anomaly height is corrected to the geological height, and the height of the karst cave is determined.

[0125] In one example, according to the projection area and the height of the karst cave, the apparent volume of the karst cave body is calculated, which includes:

[0126] The projection area and the height of the karst cave are multiplied to calculate the apparent volume of the karst cave body.

[0127] In one example, it also includes:

[0128] According to the apparent volume of the karst cave and the average porosity of the work area fracture-cave reservoir, the effective volume of the karst cave is determined, and then the static reserves of the karst geological body are calculated according to the reservoir parameters.

[0129] In one example, it also includes:

[0130] Based on the karst cave planar boundary correction and delineation of multi-attribute and forward simulation, the maximum bottom projection area of the karst cave body before and after correction is calculated respectively;

[0131] According to the multiplication of the maximum bottom projection area before and after correction and the height of the karst cave, the apparent volume of the karst cave before and after correction is calculated;

[0132] According to the volume of the cave before and after the correction, the correction rate is calculated at different scales;

[0133] According to the correction rate, the uncorrected volume of the cave corresponding to the scale is corrected.

[0134] In one example, the correction rate is:

[0135] The correction rate = (1 - the volume after correction / the volume before correction) x 100%.

[0136] Specifically, the present application is based on the theoretical basis of the fracture-cave width forward simulation, and relies on the interpretation software to realize the geological scale calculation of the fracture-cave reservoir by the columnar approximation estimation method of the bottom area multiplied by the height, which provides a fast and effective method for more accurate reserve calculation.

[0137] Based on the multi-attribute and forward simulation, the cave plane boundary correction and description are carried out. For different geological targets of cave identification, the advantage attributes of cave shape and boundary identification are determined, the seismic anomaly shape and plane boundary of the cave are determined based on the amplitude energy attribute and the coherent attribute, and the seismic anomaly boundary of the cave is described on the plane based on the interpretation software; On this basis, based on the theoretical basis of forward simulation and multi-attribute superposition, the energy attribute value range of the cave identification is adjusted to 60% of the maximum value, so as to correct the cave seismic anomaly boundary to the geological anomaly boundary, and manually describe the boundary profile to calculate the maximum bottom projection area before and after the correction of the cave volume. For the identification of the cave, the amplitude change rate and the root mean square energy attribute are better, and the coherent attribute is good for the identification of the cave boundary.

[0138] Determine the height of the cave based on wave impedance. For the identification of the cave shape, it is mainly based on the original seismic data and the amplitude energy attribute, and the root mean square amplitude anomaly determined is the seismic anomaly body of the cave, which is difficult to determine the actual height of the cave reservoir. In production, the height of the cave is quantitatively corrected, mainly using wave impedance inversion technology to predict the height. Because the impedance inversion technology can eliminate the sidelobes of the wavelet and improve the longitudinal resolution, compared with other attributes, the impedance inversion is closer to the geological model height; Combined with well logging interpretation and drilling calibration, the impedance threshold value of the typical cave reservoir is determined to determine the longitudinal height of the cave, and then the seismic anomaly height of the cave is corrected to the geological height.

[0139] Determine the apparent volume of the cave based on the columnar approximation estimation method. Since the cave reservoir is approximately a column in three-dimensional space, the layer attribute of the energy attribute in the development time window of the cave is extracted to determine the maximum bottom area of the cave, and then the corrected height of the cave is obtained by the impedance body, and then the approximate estimation method of the bottom area multiplied by the height can quickly calculate the geological apparent volume of the cave reservoir.

[0140] At the same time, the volume contrast before and after the correction of the cave reservoir can be corrected and the correction rate can be calculated. Since the correction rates of different sizes of caves are different, the correction is classified and graded according to the shape and size of the cave in the study:

[0141] For a single cave, only the correction rate needs to be counted according to different sizes, and combined with the forward simulation seismic identification result, it can be classified into more than 100m, 50-100m, and less than 50m.

[0142] For the complex cave, which is composed of two or more caves, the seismic identification is a large energy group before correction, and there is energy change locally. After correction, it will be obviously distinguished as two or several caves. If the sizes of the caves constituting the complex cave are basically the same, the threshold value is uniformly adjusted and corrected according to the above method; if the sizes of the caves constituting the complex cave are very different, uniform adjustment of the energy threshold value will lead to the disappearance of the small cave, which needs to be corrected separately, and then the correction rate is counted according to the size.

[0143] Based on the advantage of the discovered well, the correction rate range of different types and different sizes of caves is counted, which can quickly estimate the geological abnormal volume of the cave according to the classified and graded seismic abnormal volume of the cave body in the later stage, and provide a basis for more accurate and efficient estimation of the reserve.

[0144] According to the average porosity of the fracture-cave reservoir in the working area, the effective volume of the corrected cave reservoir is obtained, and the static reserve of the cave geological body is estimated according to the reservoir parameters, which provides a practical calculation process for the reserve calculation of the fracture-cave reservoir in the working area.

[0145] Example 3

[0146] According to the method, the three-dimensional volume of a single cave is calculated.

[0147] Figure 2 A schematic diagram of a column approximation estimation method implementation technical process according to an embodiment of the present application is shown.

[0148] Figure 2 The seismic abnormal plane external contour correction, height correction, and geological scale three-dimensional volume calculation process of a typical cave are shown, and the column approximation estimation method is used to realize the rapid calculation of the cave geological scale volume.

[0149] Figure 3 A schematic diagram of a cave plane contour correction process according to an embodiment of the present application is shown.

[0150] Figure 3The preferred boundary attribute identifies the planar boundary profile of the cave, i.e., the seismic anomaly boundary, and then the boundary attribute is superimposed with the energy attribute. The energy value range is adjusted under the constraint of the boundary attribute to fill the seismic anomaly boundary, and then the energy is adjusted to 60% to determine the geological boundary profile of the cave. The seismic anomaly boundary and the corrected geological boundary of the cave are manually delineated, and the scale change of the diameter of the cave before and after correction can be measured.

[0151] Figure 4 A schematic diagram of the single-type cave volume distribution correction and calculation process according to an embodiment of the present application is shown.

[0152] Figure 4 The single-type cave planar external profile and height correction process is shown. For a typical independent single-type cave, the maximum bottom area correction on the plane, i.e., the external profile identified by seismic, is the seismic anomaly profile. How to correct it to a profile that is relatively close to the actual geological scale is based on the results of the fracture-cave forward modeling. First, the seismic anomaly planar profile is determined based on the boundary attribute (such as the coherence attribute), and the seismic anomaly boundary is manually carved. Then, based on the forward modeling theory, the energy attribute (such as the root mean square amplitude, total energy attribute, etc.) is adjusted to 60% to be closest to the actual size. Then, relying on the interpretation software, the manually carved planar boundary profile after correction is obtained, i.e., the maximum cross-sectional area of the fracture-cave.

[0153] According to the present method, the geological scale three-dimensional volume correction of the composite cave is carried out.

[0154] Figure 5 A schematic diagram of the composite-type cave volume correction and calculation process according to an embodiment of the present application is shown.

[0155] Figure 5 The planar profile correction process of the composite-type cave is shown. For a composite-type cave, the seismic recognition of the distribution pattern obtains the number and pattern of the composite cave. As shown in the figure, there are three caves. For No. 1 and No. 2 caves, the scales are roughly equivalent. According to the Figure 3 、 4 The planar profile correction process shown adjusts the threshold value uniformly, delineates the seismic anomaly boundary and the corrected geological profile boundary, and measures the scale change of the diameter before and after correction. The external profile of No. 1 and No. 2 caves before correction is in a large boundary, and after correction, the cave bodies are separated, and the geological boundary profiles of the two caves are delineated, and the bottom areas of the two caves are obtained. For No. 3 cave, since the scale is quite different from that of No. 1 and No. 2 caves, it needs to be corrected separately. The correction process and Figure 3 、 Figure 4The correction process is the same as that of the single cave, manually delineating the external profile of the plane before and after correction and calculating the bottom area. For the height correction of the composite cave, the height of each cave reservoir is calibrated based on impedance inversion; then the apparent volume of each cave after correction is obtained based on the cylinder approximation estimation method.

[0156] According to the present method, static reserves are calculated and compared with dynamic reserves of constant-volume wells.

[0157] Based on the apparent volume of the cave reservoir after correction, the effective volume of the cave reservoir after correction is calculated according to the porosity of the fracture-cave reservoir in the work area, and the static reserves of the cave geological body are estimated according to the reservoir parameters, and compared with the dynamic reserves of the typical cave constant-volume well. The typical cave well in the work area is, for example, Figure 4 The seismic anomaly delineation area of the medium cave well before correction is 25024m 2 , and the area after correction is 13920m 2 . The seismic anomaly height (such as RMS amplitude identification) of the cave reservoir before correction is 96m, and the cave reservoir height after correction is 39m. Based on the cylinder approximation estimation method, the apparent volume is obtained, the seismic anomaly apparent volume before correction is 975936m 3 , and the apparent volume after correction is 542880m 3 , the volume correction rate is 44%, and the geological reserves calculated based on the corrected volume are 100,000 tons. The dynamic reserves of the constant-volume well are 120,000 tons, the static reserves calculated after correction are close to the dynamic reserves of the constant-volume well, which shows the effectiveness of the cave reservoir geological scale volume estimation method.

[0158] Example 4

[0159] Figure 6 A block diagram of a cave body geological scale apparent volume calculation device according to an embodiment of the present application is shown.

[0160] As shown in Figure 6 , the cave body geological scale apparent volume calculation device includes:

[0161] The area calculation module 201 corrects and delineates the cave plane boundary, calculates the maximum bottom surface projection area of the cave body, and calculates the maximum bottom surface projection area of the cave body.

[0162] The height calculation module 202 determines the cave height based on wave impedance.

[0163] The volume calculation module 203 calculates the apparent volume of the cave body according to the projection area and the cave height.

[0164] As an optional solution, the cave plane boundary correction and delineation, and the calculation of the maximum bottom surface projection area of the cave body include:

[0165] Determine the advantageous attributes for cave morphology and boundary identification based on different geological targets in cave identification;

[0166] The seismic anomaly morphology and planar boundary of the karst cave are determined based on the amplitude energy attribute and coherence attribute, and the seismic anomaly boundary is delineated on the plane.

[0167] Based on the theoretical basis of forward modeling and the superposition of multiple attributes, the energy attribute value range of the karst cave identification is adjusted, the seismic anomaly boundary of the karst cave is corrected to the geological anomaly boundary, the boundary outline is delineated, and the maximum bottom surface projection area of ​​the karst cave is calculated.

[0168] As an optional solution, adjusting the energy attribute value range of the karst cave identification includes:

[0169] Adjust the energy attribute value range of the identified cave to 60% of the maximum value.

[0170] As an optional approach, determining the height of a cave based on wave impedance includes:

[0171] Based on raw seismic data and amplitude energy properties, the morphology of karst caves is identified;

[0172] Based on the root mean square amplitude anomaly, the seismic anomaly of the karst cave was identified;

[0173] By using wave impedance inversion technology, the height of the karst cave is quantitatively corrected. Combined with well logging interpretation and drilling calibration of the impedance threshold value of typical karst cave reservoirs, the longitudinal height of the karst cave is determined. Then, the seismic anomaly height of the karst cave is corrected to the geological height to determine the karst cave height.

[0174] As an optional approach, the apparent volume of the cave can be calculated based on the projected area and the cave's height, including:

[0175] Calculate the apparent volume of the cave by multiplying the projected area by the cave's height.

[0176] As an optional solution, the following are also included:

[0177] The effective volume of the karst cave is determined based on the apparent volume of the cave and the average porosity of the fractured-cavity reservoir in the work area. Then, the static reserves of the karst cave geological body are calculated based on the reservoir parameters.

[0178] As an optional solution, the following are also included:

[0179] Based on multi-attribute and forward modeling, the plane boundary of the cave is corrected and characterized, and the maximum bottom surface projection area of ​​the cave body before and after correction is calculated respectively.

[0180] The apparent volume of the cave before and after correction is calculated by multiplying the maximum bottom projection area before and after correction by the height of the cave.

[0181] The correction rate is calculated on different scales based on the apparent volumes of the karst cave before and after correction.

[0182] The uncorrected cavity volume is corrected according to the correction rate corresponding to the scale.

[0183] As an option, the correction rate is:

[0184] The correction rate = (1 - corrected volume / corrected volume before correction) x 100%.

[0185] Example 5

[0186] The embodiment provides an electronic device, which comprises a memory storing executable instructions, and a processor running the executable instructions in the memory to implement the above-mentioned cavity volume calculation method of geological scale.

[0187] The electronic device according to the embodiment of the present disclosure comprises a memory and a processor.

[0188] The memory is used to store non-transitory computer readable instructions. Specifically, the memory can comprise one or more computer program products, which can comprise various forms of computer readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM), cache memory, and the like. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, and the like.

[0189] The processor can be a central processing unit (CPU) or other forms of processing units with data processing capability and / or instruction execution capability, and can control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer readable instructions stored in the memory.

[0190] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience effect, the embodiment can also include well-known structures such as a communication bus, an interface, and the like, which should also be included in the protection scope of the present disclosure.

[0191] Detailed descriptions of the embodiment can refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0192] Example 6

[0193] The embodiment provides a computer readable storage medium storing a computer program, which is executed by a processor to implement the above-mentioned cavity volume calculation method of geological scale.

[0194] A computer readable storage medium according to embodiments of the present disclosure has non-transitory computer readable instructions stored thereon. When the non-transitory computer readable instructions are run by a processor, all or part of the steps of the method of the embodiments of the present disclosure described above are performed.

[0195] The computer readable storage medium described above includes, but is not limited to, an optical storage medium (for example, a CD-ROM and a DVD), a magneto-optical storage medium (for example, an MO), a magnetic storage medium (for example, a magnetic tape or a moving hard disk), a medium having a built-in rewritable nonvolatile memory (for example, a memory card), and a medium having a built-in ROM (for example, a ROM cartridge).

[0196] It will be understood by those skilled in the art that the above description of the embodiments of the present application is only for the purpose of exemplarily illustrating the beneficial effects of the embodiments of the present application, and is not intended to limit the embodiments of the present application to any of the examples given.

[0197] The embodiments of the present application have been described above, and the above description is exemplary and is not exhaustive, and is also not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for calculating the apparent volume of a karst cave at the geological scale, characterized in that, include: Correct and delineate the planar boundaries of the cave, and calculate the maximum projected area of ​​the bottom surface of the cave. Determining the height of a karst cave based on wave impedance; Calculate the apparent volume of the cave based on the projected area and the cave height; The process includes correcting and delineating the planar boundaries of the cave, and calculating the maximum projected area of ​​the cave's bottom surface. Determine the advantageous attributes for cave morphology and boundary identification based on different geological targets in cave identification; The seismic anomaly morphology and planar boundary of the karst cave are determined based on the amplitude energy attribute and coherence attribute, and the seismic anomaly boundary is delineated on the plane. Based on the theoretical basis of forward modeling and the superposition of multiple attributes, the energy attribute value range of the karst cave identification is adjusted, the seismic anomaly boundary of the karst cave is corrected to the geological anomaly boundary, the boundary outline is delineated, and the maximum bottom surface projection area of ​​the karst cave is calculated.

2. The method for calculating the apparent volume of a karst cave at geological scale according to claim 1, wherein, Adjusting the energy attribute value range for cave identification includes: Adjust the energy attribute value range of the cave identification to 60% of the maximum value.

3. The method for calculating the apparent volume of a karst cave at geological scale according to claim 1, wherein, Determining the height of a karst cave based on wave impedance includes: Based on raw seismic data and amplitude energy properties, the morphology of karst caves is identified; Based on the root mean square amplitude anomaly, the seismic anomaly of the karst cave was identified; Using wave impedance inversion technology, the height of the karst cave is quantitatively corrected. Combined with well logging interpretation and the impedance threshold value of typical karst cave reservoirs determined by drilling, the seismic anomaly height of the karst cave is corrected to the geological height, which is the karst cave height.

4. The method for calculating the apparent volume of a karst cave at geological scale according to claim 1, wherein, Also includes: The effective volume of the karst cave is determined based on the apparent volume of the cave and the average porosity of the fractured-cavity reservoir in the work area. Then, the static reserves of the karst cave geological body are calculated based on the reservoir parameters.

5. The method for calculating the apparent volume of a karst cave at geological scale according to claim 1, wherein, Also includes: Based on multi-attribute and forward modeling, the plane boundary of the cave is corrected and characterized, and the maximum bottom surface projection area of ​​the cave body before and after correction is calculated respectively. The apparent volume of the cave before and after correction is calculated by multiplying the maximum bottom surface projection area before and after correction by the height of the cave. Based on the apparent volume of the cavern before and after correction, the correction rate is calculated on different scales. The apparent volume of the uncorrected cavern at the corresponding scale is corrected according to the correction rate.

6. The method for calculating the apparent volume of a karst cave at geological scale according to claim 5, wherein, The correction rate is: Correction rate = (1 - corrected post-view volume / corrected pre-view volume) × 100%.

7. A device for calculating the apparent volume of a karst cave at the geological scale, characterized in that, include: The area calculation module performs plane boundary correction and delineation of the cave, and calculates the maximum bottom surface projection area of ​​the cave. The height calculation module determines the height of the cave based on wave impedance; The volume calculation module calculates the apparent volume of the cave based on the projected area and the height of the cave. The process includes correcting and delineating the planar boundaries of the cave, and calculating the maximum projected area of ​​the cave's bottom surface. Determine the advantageous attributes for cave morphology and boundary identification based on different geological targets in cave identification; The seismic anomaly morphology and planar boundary of the karst cave are determined based on the amplitude energy attribute and coherence attribute, and the seismic anomaly boundary is delineated on the plane. Based on the theoretical basis of forward modeling and the superposition of multiple attributes, the energy attribute value range of the karst cave identification is adjusted, the seismic anomaly boundary of the karst cave is corrected to the geological anomaly boundary, the boundary outline is delineated, and the maximum bottom surface projection area of ​​the karst cave is calculated.

8. An electronic device, characterized in that, The electronic device includes: Memory, which stores executable instructions; A processor that executes the executable instructions in the memory to implement the method for calculating the apparent volume of a karst cave at the geological scale according to any one of claims 1-7.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method for calculating the apparent volume of a karst cave at geological scale as described in any one of claims 1-7.

Images