Carbonate rock fracture parameter calculation method and device, electronic equipment and medium

Abstract

The application discloses a carbonate rock fracture parameter calculation method and device, electronic equipment and medium. The method can comprise: determining a fracture development section according to core or imaging data, and then determining a plurality of fracture indexes; correcting the baseline of the curve of the fracture index through core or imaging interpretation of the fracture, to obtain a corrected fracture index; obtaining a comprehensive fracture index according to the corrected fracture index; and calculating a fracture parameter according to the comprehensive fracture index. The application obtains a comprehensive fracture index curve through the fracture index and the baseline drift method of the fracture index, and combines the difference method of the deep resistivity and the reconstructed matrix resistivity to obtain parameters such as fracture porosity, permeability and width, which is of great significance for improving the qualitative identification and quantitative interpretation accuracy of the fracture and the matching degree of the two.

Description

Technical Field

[0001] This invention relates to the field of well logging methods and technologies, and more specifically, to a method, apparatus, electronic device, and medium for calculating fracture parameters in carbonate rocks. Background Technology

[0002] In oil and gas field development, fractures are a crucial factor affecting the implementation of reservoir development plans. Fracture logging interpretation provides important small-scale fracture parameters for geological modeling and reservoir numerical simulation.

[0003] In existing technologies, fracture qualitative identification typically employs methods such as well core drilling and imaging logging to obtain fracture parameter data. However, these methods suffer from limited core and imaging logging data due to cost and other factors. In the study area, conventional logging curves lack shallow resistivity, making fracture identification using conventional logging data challenging; traditional qualitative and quantitative fracture identification methods, such as artificial neural networks, have shown poor performance.

[0004] The weighted method or comprehensive probability method of fracture indices has a certain effect on fracture identification in carbonate rocks. Commonly used fracture indices in this method include well diameter change, negative spontaneous potential anomaly, sonic transit rate, three-porosity ratio, pore structure index, Young's modulus, and deep-shallow resistivity difference ratio. However, density anomaly is not used as a fracture index, which may bring uncertainty to fracture identification in carbonate reservoirs.

[0005] Commonly used crack indices have undergone normal transformation and probability function calculation, but lack consistency correction processing based on crack response, resulting in unsatisfactory application effects.

[0006] In addition, the method of difference in resistivity between shallow and deep areas is widely used in the quantitative interpretation of parameters such as crack porosity, but it does not take the probability of cracks as a prerequisite, which often leads to inconsistencies between the quantitative interpretation of cracks and the actual occurrence of cracks.

[0007] Therefore, it is necessary to develop a method, device, electronic equipment, and medium for identifying fractures in carbonate rocks based on conventional well logging.

[0008] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention, and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0009] This invention proposes a method, apparatus, electronic device, and medium for calculating fracture parameters in carbonate rocks. It can obtain a comprehensive fracture index curve through fracture indices and the fracture index baseline drift method, and combine it with the deep resistivity and reconstructed matrix resistivity difference method to obtain parameters such as fracture porosity, permeability, and width. This is of great significance for improving the accuracy of qualitative identification and quantitative interpretation of fractures and the degree of matching between the two.

[0010] In a first aspect, embodiments of this disclosure provide a method for calculating fracture parameters in carbonate rocks, including:

[0011] Based on core or imaging data, determine the fracture development section, and then determine multiple fracture indices;

[0012] By interpreting the fractures through core samples or imaging, the curve of the fracture index is corrected for baseline offset to obtain the corrected fracture index.

[0013] Based on the corrected crack index, the comprehensive crack index is obtained;

[0014] Calculate the crack parameters based on the comprehensive crack index.

[0015] Preferably, the crack index includes a density anomaly ratio, which is calculated using formula (1):

[0016] ROHBbu_1=(RHOBbu-RHOB) / RHOBbu(1)

[0017] Wherein, ROHBbu_1 is the density anomaly ratio, RHOBbu is the density after filtering, and RHOB is the density before filtering.

[0018] Preferably, the comprehensive crack index is:

[0019] FRAC_ZH=∑a i FRAC_i (2)

[0020] Wherein, FRAC_ZH is the comprehensive crack index, FRAC_i is the i-th corrected crack index, and a i Let be the weight corresponding to the i-th corrected crack index.

[0021] Preferably, the crack parameters include porosity, opening degree, and permeability.

[0022] Preferably, the porosity is calculated based on the difference between deep resistivity and matrix resistivity, combined with the comprehensive crack index:

[0023]

[0024] Where, Φ f2 R represents porosity.mf R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 denoted as matrix resistivity, FRAC_ZH as comprehensive crack index, and c as a calculation coefficient.

[0025] Preferably, the opening of the low-angle seam or oblique seam is calculated using formula (4):

[0026]

[0027] Calculate the opening of high-angle seams or vertical seams using formula (5):

[0028]

[0029] Where ε is the opening degree, R m The resistivity of the mud cake.

[0030] Preferably, the permeability is calculated using formula (6):

[0031] K f =bε 2 Φ f2 (6)

[0032] Among them, K f denoted as ρ, and b is a calculation parameter.

[0033] As one specific implementation of this disclosure,

[0034] Secondly, this disclosure also provides a device for calculating fracture parameters in carbonate rocks, comprising:

[0035] The fracture index determination module identifies fracture development sections based on core or imaging data, and then determines multiple fracture indices.

[0036] The correction module interprets the cracks through core samples or imaging, performs baseline offset correction on the curve of the crack index, and obtains the corrected crack index.

[0037] The weighting module obtains the comprehensive crack index based on the corrected crack index.

[0038] The calculation module calculates the crack parameters based on the comprehensive crack index.

[0039] Preferably, the crack index includes a density anomaly ratio, which is calculated using formula (1):

[0040] ROHBbu_1=(RHOBbu-RHOB) / RHOBbu(1)

[0041] Wherein, ROHBbu_1 is the density anomaly ratio, RHOBbu is the density after filtering, and RHOB is the density before filtering.

[0042] Preferably, the comprehensive crack index is:

[0043] FRAC_ZH=∑a i FRAC_i (2)

[0044] Wherein, FRAC_ZH is the comprehensive crack index, FRAC_i is the i-th corrected crack index, and a i Let be the weight corresponding to the i-th corrected crack index.

[0045] Preferably, the crack parameters include porosity, opening degree, and permeability.

[0046] Preferably, the porosity is calculated based on the difference between deep resistivity and matrix resistivity, combined with the comprehensive crack index:

[0047]

[0048] Where, Φ f2 R represents porosity. mf R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 denoted as matrix resistivity, FRAC_ZH as comprehensive crack index, and c as a calculation coefficient.

[0049] Preferably, the opening of the low-angle seam or oblique seam is calculated using formula (4):

[0050]

[0051] Calculate the opening of high-angle seams or vertical seams using formula (5):

[0052]

[0053] Where ε is the opening degree, R m The resistivity of the mud cake.

[0054] Preferably, the permeability is calculated using formula (6):

[0055] K f =bε 2 Φ f2 (6)

[0056] Among them, K f denoted as ρ, and b is a calculation parameter.

[0057] Thirdly, embodiments of this disclosure also provide an electronic device, the electronic device comprising:

[0058] Memory, which stores executable instructions;

[0059] A processor that executes the executable instructions in the memory to implement the method for calculating carbonate rock fracture parameters.

[0060] Fourthly, embodiments of this disclosure also provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for calculating carbonate rock fracture parameters.

[0061] The methods and apparatus of the present invention have other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0062] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same parts.

[0063] Figure 1 A flowchart illustrating the steps of a method for calculating fracture parameters in carbonate rocks according to an embodiment of the present invention is shown.

[0064] Figure 2 A schematic diagram of the neutron-resistivity-crack density intersection is shown according to an embodiment of the present invention.

[0065] Figure 3 A schematic diagram of the crack inclination-neutron-acoustic wave intersection according to an embodiment of the present invention is shown.

[0066] Figure 4 A schematic diagram of the imaging logging interpretation results according to an embodiment of the present invention is shown.

[0067] Figure 5 A schematic diagram of conventional well logging fracture response according to an embodiment of the present invention is shown.

[0068] Figure 6 A schematic diagram of fracture parameters for well 010A according to an embodiment of the present invention is shown.

[0069] Figure 7 A schematic diagram of fracture parameters for well 61A according to an embodiment of the present invention is shown.

[0070] Figure 8A block diagram of a carbonate rock fracture parameter calculation device according to an embodiment of the present invention is shown.

[0071] Explanation of reference numerals in the attached figures:

[0072] 201. Crack index determination module; 202. Correction module; 203. Weighting module; 204. Calculation module. Detailed Implementation

[0073] Preferred embodiments of the invention will now be described in more detail. While preferred embodiments of the invention are described below, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0074] This invention provides a method for calculating fracture parameters in carbonate rocks, comprising:

[0075] Based on core or imaging data, determine the fracture development section, and then determine multiple fracture indices;

[0076] By interpreting the fractures through core samples or imaging, the curves of the fracture index are corrected for baseline shift to obtain the corrected fracture index.

[0077] Based on the corrected crack index, the comprehensive crack index is obtained;

[0078] Calculate crack parameters based on the comprehensive crack index.

[0079] In one example, the crack index includes the density anomaly ratio, which is calculated using formula (1):

[0080] ROHBbu_1=(RHOBbu-RHOB) / RHOBbu (1)

[0081] Wherein, ROHBbu_1 is the density anomaly ratio, RHOBbu is the density after filtering, and RHOB is the density before filtering.

[0082] In one example, the comprehensive crack index is:

[0083] FRAC_ZH=∑a i FRAC_i (2)

[0084] Wherein, FRAC_ZH is the comprehensive crack index, FRAC_i is the i-th corrected crack index, and a i Let be the weight corresponding to the i-th corrected crack index.

[0085] In one example, crack parameters include porosity, aperture, and permeability.

[0086] In one example, porosity is calculated based on the difference between deep resistivity and matrix resistivity, combined with a comprehensive crack index:

[0087]

[0088] Where, Φ f2 R represents porosity. mf R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 denoted as matrix resistivity, FRAC_ZH as comprehensive crack index, and c as a calculation coefficient.

[0089] In one example, the opening of a low-angle seam or oblique seam is calculated using formula (4):

[0090]

[0091] Calculate the opening of high-angle seams or vertical seams using formula (5):

[0092]

[0093] Where ε is the opening degree, R m The resistivity of the mud cake.

[0094] In one example, the permeability is calculated using formula (6):

[0095] K f =bε 2 Φ f2 (6)

[0096] Among them, K f denoted as ρ, and b is a calculation parameter.

[0097] Specifically, based on core or imaging data, the fracture development section is determined, and conventional curves or curve combinations that are more sensitive to fracture response are selected as fracture indicators, including density anomaly ratio, resistivity intrusion ratio, pore structure index, three-porosity ratio, and saturation ratio: the resistivity intrusion ratio is RD_2=(RD-MLL) / RD; the pore structure index is m_2=(A-BlogF) / (1+Blog Φ); the three-porosity ratio is PHI_C=(Φt-Φs) / Φt; the saturation ratio is SW_X=SXO / SW; and the density anomaly ratio is given by formula (1).

[0098] Interpreting fractures through core samples or imaging, the fracture index curve is corrected for baseline offset using the following formula to ensure that there are no fractures below the baseline, thus obtaining the corrected fracture index:

[0099] FRAC_1=RD_2+a

[0100] FRAC_2=m_2+b

[0101] FRAC_3=PHI_C+c

[0102] FRAC_4=SW_X+e

[0103] FRAC_5=ROHBbu_1+f

[0104] In the formula, FRAC_1 to FRAC_5 are the curves after correction of the crack index baseline offset, and a to f are constants.

[0105] By weighting the corrected crack indices, a comprehensive crack index curve is obtained, and crack development segments are identified. The comprehensive crack index is as follows:

[0106] FRAC_ZH=a1*FRAC_1+a2*FRAC_2+a3*FRAC_3+a4*FRAC_4+a5*FRAC_5

[0107] Reconstructing matrix resistivity through deep lateral intersection with compensating neutrons in both cracked and non-cracked areas:

[0108] R 基质 =3461,77*e (-27.0038*NPHI)

[0109] In a reservoir segment, when the matrix resistivity is significantly greater than the deep lateral resistivity, it indicates the presence of fractures in that segment; when the matrix resistivity is equal to or very close to the deep lateral resistivity, it indicates the absence of fractures in that reservoir.

[0110] By intersecting acoustic waves with compensated neutrons at low-angle and high-angle fractures, the acoustic transit time of the matrix is ​​reconstructed, and the fracture orientation is analyzed.

[0111] DT 基质 =46.0838 + 112.156 * NPHI

[0112] When the matrix acoustic transit time is significantly less than the acoustic transit time, low-angle or oblique-angle fractures exist in the reservoir. When the matrix acoustic transit time is equal to or very close to the acoustic transit time, high-angle fractures or no fractures exist in the reservoir.

[0113] Based on the difference between deep resistivity and matrix resistivity, and combined with comprehensive crack indices, crack porosity is calculated.

[0114] For fracture-porosity reservoirs, the fracture porosity model is derived based on the traditional deep-shallow resistivity difference:

[0115]

[0116] Considering the significant impact of comprehensive crack index on cracks, this invention improves the above formula and calculates porosity using formula (3).

[0117] The crack width is calculated using the finite element method simulation results, and the crack permeability is calculated based on the core calibration. The opening of low-angle or oblique cracks is calculated using formula (4), the opening of high-angle or vertical cracks is calculated using formula (5), and the permeability is calculated using formula (6).

[0118] The present invention also provides a device for calculating fracture parameters in carbonate rocks, comprising:

[0119] The fracture index determination module identifies fracture development sections based on core or imaging data, and then determines multiple fracture indices.

[0120] The correction module interprets the fractures through core samples or imaging, corrects the baseline offset of the fracture index curve, and obtains the corrected fracture index.

[0121] The weighted module obtains the comprehensive crack index based on the corrected crack index.

[0122] The calculation module calculates crack parameters based on the comprehensive crack index.

[0123] In one example, the crack index includes the density anomaly ratio, which is calculated using formula (1):

[0124] ROHBbu_1=(RHOBbu-RHOB) / RHOBbu (1)

[0125] Wherein, ROHBbu_1 is the density anomaly ratio, RHOBbu is the density after filtering, and RHOB is the density before filtering.

[0126] In one example, the comprehensive crack index is:

[0127] FRAC_ZH=∑a i FRAC_i (2)

[0128] Wherein, FRAC_ZH is the comprehensive crack index, FRAC_i is the i-th corrected crack index, and a i Let be the weight corresponding to the i-th corrected crack index.

[0129] In one example, crack parameters include porosity, aperture, and permeability.

[0130] In one example, porosity is calculated based on the difference between deep resistivity and matrix resistivity, combined with a comprehensive crack index:

[0131]

[0132] Where, Φ f2 R represents porosity. mf R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 denoted as matrix resistivity, FRAC_ZH as comprehensive crack index, and c as a calculation coefficient.

[0133] In one example, the opening of a low-angle seam or oblique seam is calculated using formula (4):

[0134]

[0135] Calculate the opening of high-angle seams or vertical seams using formula (5):

[0136]

[0137] Where ε is the opening degree, R m The resistivity of the mud cake.

[0138] In one example, the permeability is calculated using formula (6):

[0139] K f =bε 2 Φ f2 (6)

[0140] Among them, K f denoted as ρ, and b is a calculation parameter.

[0141] Specifically, based on core or imaging data, the fracture development section is determined, and conventional curves or curve combinations that are more sensitive to fracture response are selected as fracture indicators, including density anomaly ratio, resistivity intrusion ratio, pore structure index, three-porosity ratio, and saturation ratio: the resistivity intrusion ratio is RD_2=(RD-MLL) / RD; the pore structure index is m_2=(A-BlogF) / (1+Blog Φ); the three-porosity ratio is PHI_C=(Φt-Φs) / Φt; the saturation ratio is SW_X=SXO / SW; and the density anomaly ratio is given by formula (1).

[0142] Interpreting fractures through core samples or imaging, the fracture index curve is corrected for baseline offset using the following formula to ensure that there are no fractures below the baseline, thus obtaining the corrected fracture index:

[0143] FRAC_1=RD_2+a

[0144] FRAC_2=m_2+b

[0145] FRAC_3=PHI_C+c

[0146] FRAC_4=SW_X+e

[0147] FRAC_5=ROHBbu_1+f

[0148] In the formula, FRAC_1 to FRAC_5 are the curves after correction of the crack index baseline offset, and a to f are constants.

[0149] By weighting the corrected crack indices, a comprehensive crack index curve is obtained, and crack development segments are identified. The comprehensive crack index is as follows:

[0150] FRAC_ZH=a1*FRAC_1+a2*FRAC_2+a3*FRAC_3+a4*FRAC_4+a5*FRAC_5

[0151] Reconstructing matrix resistivity through deep lateral intersection with compensating neutrons in both cracked and non-cracked areas:

[0152] R 基质 =3461,77*e (-27.0038*NPHI)

[0153] In a reservoir segment, when the matrix resistivity is significantly greater than the deep lateral resistivity, it indicates the presence of fractures in that segment; when the matrix resistivity is equal to or very close to the deep lateral resistivity, it indicates the absence of fractures in that reservoir.

[0154] By intersecting acoustic waves with compensated neutrons at low-angle and high-angle fractures, the acoustic transit time of the matrix is ​​reconstructed, and the fracture orientation is analyzed.

[0155] DT 基质 =46.0838 + 112.156 * NPHI

[0156] When the matrix acoustic transit time is significantly less than the acoustic transit time, low-angle or oblique-angle fractures exist in the reservoir. When the matrix acoustic transit time is equal to or very close to the acoustic transit time, high-angle fractures or no fractures exist in the reservoir.

[0157] Based on the difference between deep resistivity and matrix resistivity, and combined with comprehensive crack indices, crack porosity is calculated.

[0158] For fracture-porosity reservoirs, the fracture porosity model is derived based on the traditional deep-shallow resistivity difference:

[0159]

[0160] Considering the significant impact of comprehensive crack index on cracks, this invention improves the above formula and calculates porosity using formula (3).

[0161] The crack width is calculated using the finite element method simulation results, and the crack permeability is calculated based on the core calibration. The opening of low-angle or oblique cracks is calculated using formula (4), the opening of high-angle or vertical cracks is calculated using formula (5), and the permeability is calculated using formula (6).

[0162] The present invention also provides an electronic device, comprising: a memory storing executable instructions; and a processor that executes the executable instructions in the memory to implement the above-described method for calculating carbonate rock fracture parameters.

[0163] The present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for calculating carbonate rock fracture parameters.

[0164] To facilitate understanding of the solutions and effects of the embodiments of the present invention, four specific application examples are given below. Those skilled in the art should understand that these examples are merely for the purpose of understanding the present invention, and any specific details therein are not intended to limit the present invention in any way.

[0165] Example 1

[0166] Figure 1 A flowchart illustrating the steps of the method for calculating carbonate rock fracture parameters according to the present invention is shown.

[0167] like Figure 1 As shown, the method for calculating fracture parameters in carbonate rocks includes: Step 101, determining the fracture development section based on core or imaging data, and then determining multiple fracture indices; Step 102, interpreting the fractures through core or imaging data, and performing baseline offset correction on the fracture index curves to obtain corrected fracture indices; Step 103, obtaining a comprehensive fracture index based on the corrected fracture index; Step 104, calculating fracture parameters based on the comprehensive fracture index.

[0168] Figure 2 A schematic diagram of neutron-resistivity-crack density intersection according to an embodiment of the present invention is shown. The horizontal axis represents neutron porosity, and the vertical axis represents deep lateral resistivity. The matrix resistivity is reconstructed by intersecting the deep lateral resistivity with compensated neutrons. The matrix data is represented within the polygons in the diagram. When the matrix resistivity is significantly greater than the deep lateral resistivity, it indicates the presence of cracks; when the matrix resistivity is close to the deep lateral resistivity, it indicates the absence of cracks.

[0169] Figure 3A schematic diagram of fracture dip angle-neutron-acoustic wave intersection according to an embodiment of the present invention is shown. The horizontal axis represents neutron porosity, and the vertical axis represents acoustic wave transit time. The matrix acoustic wave is reconstructed by the intersection of acoustic waves and compensating neutrons. The polygons in the diagram represent matrix data. When the matrix acoustic wave transit time is significantly less than the acoustic wave transit time, low-angle fractures or oblique fractures exist in the reservoir; when the matrix acoustic wave transit time is equal to or very close to the acoustic wave transit time, high-angle fractures or no fractures exist in the reservoir.

[0170] Figure 4 A schematic diagram of the imaging logging interpretation results according to an embodiment of the present invention is shown. From left to right, the images are depth, static image, fracture dip and azimuth tadpole, and dynamic image. The static image reflects the overall changes of fractures throughout the well section, while the dynamic image reflects the changes of fractures in local well sections.

[0171] Figure 5 A schematic diagram of conventional well logging fracture response according to an embodiment of the present invention is shown, from left to right: gamma ray, depth, shallow and deep resistivity, three-porosity curves (density, neutron, acoustic wave), fracture dip angle, and azimuth tadpole. The conventional well logging fracture response is characterized by increased resistivity invasion, abnormal density, and differences in the three porosities.

[0172] Figure 6 A schematic diagram of fracture parameters in well 010A according to an embodiment of the present invention is shown. FRAC_ZH is the comprehensive fracture index curve, AllFr_Den is the fracture density interpreted by imaging logging, and the last two channels are the fracture strike rose diagram and fracture tadpole diagram (dip and azimuth). It can be seen that the comprehensive fracture index curve is basically consistent with the imaging fracture density and fracture tadpole diagram. width_f is the fracture width, width is the fracture width interpreted by imaging, PHI_T is the total porosity, PHIE_f is the fracture porosity, and perm_f is the fracture permeability. It can be seen that the fracture width is basically consistent with the fracture width interpreted by imaging. The calculated fracture parameters are consistent with the qualitative fracture identification results.

[0173] Figure 7 A schematic diagram of fracture parameters in well 61A according to an embodiment of the present invention is shown. The comprehensive fracture index curve is basically consistent with the imaging fracture density and fracture tadpole. Perf2 represents the perforation section, and QO_2004 represents the production profile. It can be seen that fractures have a promoting effect on reservoir production. The calculated fracture parameters are consistent with the qualitative fracture identification results.

[0174] Example 2

[0175] Figure 8 A block diagram of a carbonate rock fracture parameter calculation device according to an embodiment of the present invention is shown.

[0176] like Figure 8 As shown, the carbonate rock fracture parameter calculation device includes:

[0177] The fracture index determination module 201 determines the fracture development section based on core or imaging data, and then determines multiple fracture indices.

[0178] The correction module 202 interprets the fractures through core samples or imaging, corrects the baseline offset of the fracture index curve, and obtains the corrected fracture index.

[0179] The weighting module 203 obtains the comprehensive crack index based on the corrected crack index;

[0180] Calculation module 204 calculates crack parameters based on the comprehensive crack index.

[0181] As an optional option, the crack index includes the density anomaly ratio, which is calculated using formula (1):

[0182] ROHBbu_1=(RHOBbu-RHOB) / RHOBbu (1)

[0183] Wherein, ROHBbu_1 is the density anomaly ratio, RHOBbu is the density after filtering, and RHOB is the density before filtering.

[0184] As an optional solution, the comprehensive crack index is:

[0185] FRAC_ZH=∑a i FRAC_i (2)

[0186] Wherein, FRAC_ZH is the comprehensive crack index, FRAC_i is the i-th corrected crack index, and a i Let be the weight corresponding to the i-th corrected crack index.

[0187] As an optional option, crack parameters include porosity, aperture, and permeability.

[0188] As an alternative, porosity can be calculated based on the difference between deep resistivity and matrix resistivity, combined with comprehensive crack indices.

[0189]

[0190] Where, Φ f2 R represents porosity. mf R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 denoted as matrix resistivity, FRAC_ZH as comprehensive crack index, and c as a calculation coefficient.

[0191] As an alternative, the opening of the low-angle seam or oblique seam can be calculated using formula (4):

[0192]

[0193] Calculate the opening of high-angle seams or vertical seams using formula (5):

[0194]

[0195] Where ε is the opening degree, R m The resistivity of the mud cake.

[0196] As an alternative, the permeability can be calculated using formula (6):

[0197] K f =bε 2 Φ f2 (6)

[0198] Among them, K f denoted as ρ, and b is a calculation parameter.

[0199] Example 3

[0200] This disclosure provides an electronic device comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the aforementioned method for calculating carbonate rock fracture parameters.

[0201] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0202] This memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include 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) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0203] The processor may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

[0204] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this disclosure.

[0205] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0206] Example 4

[0207] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for calculating carbonate rock fracture parameters.

[0208] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0209] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0210] Those skilled in the art should understand that the above description of the embodiments of the present invention is only intended to illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any of the examples given.

[0211] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for calculating fracture parameters in carbonate rocks, characterized in that, include: Based on core or imaging data, determine the fracture development section, and then determine multiple fracture indices; By interpreting the fractures through core samples or imaging, the curve of the fracture index is corrected for baseline offset to obtain the corrected fracture index. Based on the corrected crack index, the comprehensive crack index is obtained; Calculate the crack parameters based on the comprehensive crack index; The crack parameters include porosity, opening degree, and permeability; Specifically, the porosity is calculated based on the difference between deep resistivity and matrix resistivity, combined with the comprehensive crack index: （3） in, Porosity R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 FRAC_ZH is the matrix resistivity, FRAC_ZH is the comprehensive crack index, and c is the calculation coefficient; The opening of low-angle or oblique seams is calculated using formula (4): （4） Calculate the opening of high-angle seams or vertical seams using formula (5): （5） Where ε is the opening degree, The resistivity of the mud cake; The permeability is calculated using formula (6): （6） Among them, K f denoted as ρ, and b is a calculation parameter.

2. The method for calculating fracture parameters in carbonate rocks according to claim 1, wherein, The crack index includes the density anomaly ratio, which is calculated using formula (1): ROHBbu_1=(RHOBbu-RHOB) / RHOBbu (1) Wherein, ROHBbu_1 is the density anomaly ratio, RHOBbu is the density after filtering, and RHOB is the density before filtering.

3. The method for calculating fracture parameters in carbonate rocks according to claim 1, wherein, The comprehensive index of the crack is: （2） Wherein, FRAC_ZH is the comprehensive crack index, FRAC_i is the i-th corrected crack index, and a i Let be the weight corresponding to the i-th corrected crack index.

4. A device for calculating fracture parameters in carbonate rocks, characterized in that, include: The fracture index determination module identifies fracture development sections based on core or imaging data, and then determines multiple fracture indices. The correction module interprets the cracks through core samples or imaging, performs baseline offset correction on the curve of the crack index, and obtains the corrected crack index. The weighting module obtains the comprehensive crack index based on the corrected crack index. The calculation module calculates the crack parameters based on the comprehensive crack index. The crack parameters include porosity, opening degree, and permeability; Specifically, the porosity is calculated based on the difference between deep resistivity and matrix resistivity, combined with the comprehensive crack index: （3） in, Porosity R is the resistivity of the mud filtrate. LLD For deep lateral resistivity, R 基质 FRAC_ZH is the matrix resistivity, FRAC_ZH is the comprehensive crack index, and c is the calculation coefficient; The opening of low-angle or oblique seams is calculated using formula (4): （4） Calculate the opening of high-angle seams or vertical seams using formula (5): （5） Where ε is the opening degree, The resistivity of the mud cake; The permeability is calculated using formula (6): （6） Among them, K f denoted as ρ, and b is a calculation parameter.

5. 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 carbonate rock fracture parameters according to any one of claims 1-3.

6. 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 carbonate rock fracture parameters as described in any one of claims 1-3.

Images