Fractured reservoir type evaluation method and device, computer equipment and storage medium

Abstract

The invention provides a fractured reservoir type evaluation method and device, computer equipment and a storage medium, and the method comprises the steps: obtaining conventional logging information of a fractured reservoir, and solving the total porosity PorT, matrix porosity PorB and fracture porosity PorFra of the reservoir according to a conventional method; calculating a fracture reservoir development index RECP, and defining reservoir development levels (class I, class II and class III) through the RECP value; calculating a reservoir type main control factor index ZKYS, and defining fracture development characteristics (a high angle, a low angle and a medium-high angle) through the ZKYS; and establishing a reservoir type evaluation method through a cross plot technology. According to the method, quantitative evaluation of the fracture type is achieved by applying a conventional logging method, the coincidence rate is high compared with imaging logging, the logging interpretation evaluation level of the fractured reservoir is improved, and the application range and prospect are quite wide.

Description

Technical Field

[0001] The present invention relates to the technical field of fracture reservoir type judgment, and particularly relates to a method, device, computer device and storage medium for evaluating fracture reservoir types. Background Art

[0002] Fault structures are basic tectonic types widely developed in the earth's crust. The activities and distributions of faults control the development and evolution of basins and the formation of tectonic patterns, and play a decisive role in the generation, migration and accumulation of oil and gas resources. In the multi-stage and different-scale fractures formed by fault structures, cave-type and fracture network reservoir systems formed under the action of deep fluids or epigene karstification, or multi-stage strike-slip tectonic ruptures and material volume adjustment within the fault zone, are important influencing factors for the development of deep and ultra-deep reservoirs in the Tarim Basin and Mesozoic reservoirs on the southwestern margin of the Ordos Basin.

[0003] Accurately evaluating the fracture reservoir type based on logging data is of great significance for the selection of oil and gas exploration and development directions, the optimization of development methods, and the improvement of oil and gas well productivity. Usually, high-end logging technologies such as micro-resistivity imaging logging and dipole acoustic imaging logging are mainly used to evaluate reservoir types. The differences in imaging logging images are used to evaluate reservoir levels, the angles of holes / fractures, etc., to achieve quantitative analysis of the fracture-cave system. However, only the existence of fractures can be qualitatively evaluated through conventional dual laterolog, acoustic transit time, density, neutron and other curves, and it is difficult to achieve quantification by classifying the corresponding differences between dual laterolog, acoustic logging, density logging, neutron logging and the bedrock. Summary of the Invention

[0004] Based on this, it is necessary to provide a method, device, computer device and storage medium for evaluating fracture reservoir types in view of the above technical problems.

[0005] A method for evaluating fracture reservoir types includes:

[0006] Measuring and obtaining logging data of a fracture reservoir;

[0007] Based on the logging data of the fracture reservoir, calculating the total porosity PorT, matrix porosity PorB and fracture porosity PorFra of the fracture reservoir respectively;

[0008] Based on the logging data of the fracture reservoir and the calculation results of the total porosity PorT, matrix porosity PorB and fracture porosity PorFra, calculating the development index RECP;

[0009] Based on the measurement results of the logging data of the fracture reservoir, calculating the main control factor index ZKYS;

[0010] Intersect the calculated development index RECP with the main control factor index ZKYS to obtain a crossplot, and evaluate the development type of the fractured reservoir based on the crossplot.

[0011] In one embodiment, the step of measuring and obtaining the logging data of the fractured reservoir includes:

[0012] Measure and obtain the logging data of the fractured reservoir based on wireline logging.

[0013] In one embodiment, the logging data of the fractured reservoir includes logging curves of acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS.

[0014] In one embodiment, in the step of calculating the development index RECP based on the logging data of the fractured reservoir, as well as the calculation results of the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra, the development index RECP is calculated based on the following calculation formula;

[0015]

[0016] wherein, RD is the deep lateral logging data, PorT is the total porosity data, PorB is the matrix porosity data, PorFra is the fracture porosity data, and RS is the shallow lateral logging data.

[0017] In one embodiment, in the step of calculating the main control factor index ZKYS based on the measurement results of the logging data of the fractured reservoir, the main control factor index ZKYS is calculated based on the following calculation formula;

[0018]

[0019] wherein, AC is the acoustic logging data, RS is the shallow lateral logging data, and RD is the deep lateral logging data.

[0020] In one embodiment, the step of intersecting the calculated development index RECP with the main control factor index ZKYS to obtain a crossplot, and evaluating the development type of the fractured reservoir based on the crossplot includes:

[0021] Based on crossplot technology, establish an evaluation method for the development type of the fractured reservoir.

[0022] In one embodiment, after the step of intersecting the calculated development index RECP with the main control factor index ZKYS to obtain a crossplot, and evaluating the development type of the fractured reservoir based on the crossplot, it further includes:

[0023] The result graph is obtained based on the micro-resistivity imaging logging method, and the result of the fracture reservoir type evaluation method is verified and interpreted by using the result graph.

[0024] A fracture reservoir type evaluation device, comprising:

[0025] A logging data acquisition module, configured to measure and acquire logging data of a fracture reservoir;

[0026] A porosity calculation module, configured to calculate the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra of the fracture reservoir respectively based on the logging data of the fracture reservoir;

[0027] A development index calculation module, configured to calculate a development index RECP based on the logging data of the fracture reservoir and the calculation results of the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra;

[0028] A main control factor index calculation module, configured to calculate a main control factor index ZKYS based on the measurement result of the logging data of the fracture reservoir;

[0029] An evaluation result obtaining module, configured to cross the calculated development index RECP and the main control factor index ZKYS to obtain a cross plot, and evaluate the development type of the fracture reservoir based on the cross plot.

[0030] A computer device, comprising a memory and a processor, the memory stores a computer program, and is characterized in that when the processor executes the computer program, the following steps are implemented:

[0031] Measure and acquire logging data of a fracture reservoir;

[0032] Based on the logging data of the fracture reservoir, calculate the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra of the fracture reservoir respectively;

[0033] Based on the logging data of the fracture reservoir and the calculation results of the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra, calculate a development index RECP;

[0034] Based on the measurement result of the logging data of the fracture reservoir, calculate a main control factor index ZKYS;

[0035] Cross the calculated development index RECP and the main control factor index ZKYS to obtain a cross plot, and evaluate the development type of the fracture reservoir based on the cross plot.

[0036] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the following steps:

[0037] Measure and obtain logging data of fractured reservoirs;

[0038] Based on the logging data of the fractured reservoir, the total porosity PorT, the matrix porosity PorB and the fracture porosity PorFra of the fractured reservoir are calculated respectively;

[0039] Calculating the development index RECP based on the well logging data of the fractured reservoir and the calculation results of the total porosity PorT, the matrix porosity PorB and the fracture porosity PorFra;

[0040] Based on the measurement results of the well logging data of the fractured reservoir, a main controlling factor index ZKYS is calculated;

[0041] The calculated development index RECP is intersected with the main controlling factor index ZKYS to obtain an intersection diagram, and the development type of the fractured reservoir is evaluated based on the intersection diagram.

[0042] The above-mentioned fracture reservoir type evaluation method, device, computer equipment and storage medium have carried out in-depth exploration of conventional logging and reservoir parameters. Without adding new technical means, it can realize the classification evaluation of fracture reservoirs and improve the level of interpretation and application of logging data.

[0043] The present invention provides a method for quantitatively evaluating the type of fractured reservoir based on resistivity and acoustic logging, the method comprising: obtaining conventional logging data of fractured reservoirs, and obtaining the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of fractured reservoirs according to conventional methods; calculating the development index RECP of fractured reservoirs, and defining the development level (class I, class II, class III) of fractured reservoirs by the numerical value of the development index RECP; calculating the main controlling factor index ZKYS of fractured reservoirs, and defining the fracture development characteristics (high angle, low angle, medium-high angle) by the numerical value of the main controlling factor index ZKYS; and establishing a reservoir type evaluation method by using the intersection diagram technology. The method realizes the quantitative evaluation of fracture types by using conventional logging methods, has a high conformity rate compared with imaging logging, improves the logging interpretation and evaluation level of fractured reservoirs, and has a very broad application scope and prospect. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 A schematic diagram of a process of evaluating a fracture reservoir type in one embodiment;

[0045] Figure 2 A schematic diagram of the measurement results of Well 1 in one embodiment;

[0046] Figure 3 Schematic diagram of the measurement results of Well 2 in one embodiment;

[0047] Figure 4 Schematic diagram of the crossplot results of the development index RECP and the main control factor index ZKYS in one embodiment;

[0048] Figure 5 Schematic diagram of the results obtained from microresistivity imaging logging in one embodiment;

[0049] Figure 6 Structural block diagram of a fracture reservoir type evaluation device in one embodiment;

[0050] Figure 7 Internal structure diagram of a computer device in one embodiment. Detailed implementation manners

[0051] In order to make the objectives, technical solutions and advantages of the present application clearer and more understandable, the present application will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0052] Embodiment 1

[0053] In this embodiment, as Figure 1 shown, a fracture reservoir type evaluation method is provided, which includes:

[0054] Step 110, measuring and obtaining logging data of a fracture reservoir.

[0055] In this embodiment, the logging data of the fracture reservoir includes logging curves of acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS. Among them, the acoustic logging AC is an acoustic travel time curve of three porosities; the density logging DEN is a compensated density curve of three porosities; the deep lateral logging RD is a resistivity curve of the deep lateral; the shallow lateral logging RS is a resistivity curve of the shallow lateral.

[0056] In one embodiment, the logging data of the fracture reservoir is measured based on the cable logging method. Specifically, cable logging is a method of detecting the structure and properties of underground rock formations using a cable. Its principle is to install a logging tool with a cable in the wellbore, extend the cable from the wellhead to the bottom of the well, and then infer the physical properties and structure of the underground rock formations based on the resistance of the cable. The logging tool for cable logging includes a probe with a cable, a cable, and a cable output device. The probe is usually equipped with a sealed annular detector. According to the measurement requirements and the depth of the well, logging tools with different lengths and precisions can be selected.

[0057] In this embodiment, the steps of wireline logging include:

[0058] 1. First, check whether the logging tool and its related equipment are operating properly and confirm that the measurement parameter settings are reasonable.

[0059] 2. Then lower the logging tool into the wellbore and pay attention to the movement of the logging tool during the process to prevent damage.

[0060] 3. When the logging tool reaches the predetermined position, start data acquisition and recording. At this time, check whether the data of the logging equipment is stable to ensure the accuracy of the data.

[0061] 4. Conduct vertical logging or inclined logging, continuously acquire data as needed, and record the actual depth of the logging.

[0062] 5. After completing the logging, slowly raise the logging tool out of the well and observe whether the logging tool is damaged or deformed during the raising process.

[0063] Step 120: Based on the logging data of the fractured reservoir, calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of the fractured reservoir respectively.

[0064] In this embodiment, according to the logging curve data such as the acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS of the fractured reservoir obtained in step 110, substitute them into the relevant calculation formulas to calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of the fractured reservoir.

[0065] Specifically, use the density logging DEN to determine the total porosity PorT of the reservoir. Based on the following calculation formula, calculate the total porosity PorT:

[0066]

[0067] In the above calculation formula:

[0068] ρ ma represents the reservoir matrix density, with the unit of grams per cubic centimeter (g / cm3);

[0069] ρ b represents the reservoir bulk density, with the unit of grams per cubic centimeter (g / cm3);

[0070] ρ f represents the reservoir fluid density, with the unit of grams per cubic centimeter (g / cm3);

[0071] Use the acoustic transit time logging AC to determine the reservoir matrix porosity PorB. Based on the following calculation formula, calculate the total porosity PorB:

[0072]

[0073] In the above calculation formula:

[0074] Δt represents the time difference of acoustic logging, with the unit of microseconds per meter (μs / m);

[0075] Δt ma represents the acoustic wave time difference of the reservoir rock skeleton, with the unit of microseconds per meter (μs / m);

[0076] Δt f represents the acoustic wave time difference of the reservoir fluid, with the unit of microseconds per meter (μs / m);

[0077] C P represents the compaction correction coefficient of the reservoir. Select 3.7 for old formations and 1 for new formations.

[0078] Use the deep lateral RD and shallow lateral data RS to calculate the fracture porosity PorFra; first, judge the fracture state. Since the dual laterolog shows different curve characteristics for different fracture states, to solve for the fracture porosity PorFra, use the following calculation formula to judge the fracture state:

[0079]

[0080] In the above calculation formula:

[0081] Y represents the discrimination index, dimensionless, and is the calculation result;

[0082] RD represents the deep lateral resistivity logging value, with the unit of ohm-meter (Ω·m);

[0083] RS represents the shallow lateral resistivity logging value, with the unit of ohm-meter (Ω·m).

[0084] During the whole well section processing, first calculate the Y value, and then judge the fracture state. When Y > 0.1, it is a high-angle fracture; when 0.1 ≥ Y > 0, it is an inclined fracture; when Y < 0, it is a low-angle fracture. After the fracture state is judged, use the following interpretation model to calculate the fracture porosity:

[0085]

[0086] In the above calculation formula:

[0087] PORFra represents the fracture porosity, with the unit of decimal;

[0088] Rmf represents the resistivity of the mud filtrate, with the unit of ohm-meter (Ω·m);

[0089] A1, A2, and A3 are constants, and their values vary depending on the crack state. Based on the value of the above discrimination index Y, the values of A1, A2, and A3 are determined, as shown in Table 1.

[0090] Table 1 Constant Value Table for the Crack Porosity Interpretation Model

[0091] Crack state Y A1 A2 A3 Low-angle crack Y＜0 -0.992417 1.97247 0.000318291 Inclined crack 0≤Y≤0.1 -17.6332 20.36451 0.00093177 High-angle crack Y＞0.1 8.522532 -8.242788 0.00071236

[0092] Step 130: Calculate the development index RECP based on the well logging data of the fractured reservoir, as well as the calculation results of the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra.

[0093] In this embodiment, the development index RECP is calculated based on the following formula:

[0094]

[0095] In the formula, RD is the deep lateral logging data, PorT is the total porosity data, PorB is the matrix porosity data, PorFra is the fracture porosity data, and RS is the shallow lateral logging data.

[0096] Define the reservoir development level through the value of the development index RECP.

[0097] When RECP > 1, the reservoir development level is defined as Class I.

[0098] When 0.1 < RECP < 1, the reservoir development level is defined as Class II.

[0099] When 0.02 < RECP < 0.1, the reservoir development level is defined as Class III.

[0100] Step 140: Calculate the main control factor index ZKYS based on the measurement results of the well logging data of the fractured reservoir.

[0101] In this embodiment, the main control factor index ZKYS is also called the type main control factor index ZKYS.

[0102] In this embodiment, the main control factor index ZKYS is calculated based on the following formula:

[0103]

[0104] In the formula, AC is the acoustic logging data, RS is the shallow lateral logging data, and RD is the deep lateral logging data.

[0105] Define the fracture development characteristics through the value of the main control factor index ZKYS.

[0106] When ZKYS>0.04, the reservoir is defined as having high-angle fracture development characteristics;

[0107] When ZKYS<-0.04, the reservoir is defined as having low-angle fracture development characteristics;

[0108] When -0.04<ZKYS<0.04, the reservoir is defined as having fracture development characteristics of medium and high angles.

[0109] Step 150: Intersect the calculated development index RECP with the main controlling factor index ZKYS to obtain an intersection diagram, and evaluate the development type of the fracture reservoir based on the intersection diagram.

[0110] In this embodiment, a method for evaluating the development type of fractured reservoirs is established based on the crossplot technology. Specifically, the crossplot technology refers to a mapping interpretation technology for well logging data. It intersects two types of well logging data on a plane map and determines the value or range of the desired parameter based on the coordinates of the intersection point.

[0111] In one embodiment, after step 150, the method further includes:

[0112] Step 160, obtaining a result map based on microresistivity imaging logging, and using the result map to verify and interpret the result of the fracture reservoir type evaluation method.

[0113] Specifically, microresistivity imaging logging refers to the use of multiple rows of button-shaped small electrodes on a multi-electrode plate to emit current to the well wall formation. Due to the different rock composition, structure and fluid contained in the electrodes, the current changes, which reflects the change in rock resistivity at various locations on the well wall, and can display resistivity imaging of the well wall. A result map is obtained according to the microresistivity imaging logging method, and the result map can be used to verify whether the development type of the fracture reservoir evaluated based on the intersection map is correct.

[0114] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0115] Example Two

[0116] In this example, a method for evaluating the type of fractured reservoir is provided, which includes:

[0117] Select two wells in a certain well area, denoted as Well 1 and Well 2 respectively;

[0118] Step 1: Measure and obtain the logging data of Well 1 and Well 2 respectively.

[0119] In this example, the logging data of Well 1 and Well 2 include logging curves of acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS. Among them, the acoustic logging AC is the acoustic travel time curve of three porosities; the density logging DEN is the compensated density curve of three porosities; the deep lateral logging RD is the resistivity curve of the deep lateral; the shallow lateral logging RS is the resistivity curve of the shallow lateral.

[0120] In one embodiment, the logging data of the fractured reservoir is measured based on the wireline logging method. Specifically, wireline logging is a method of detecting the structure and properties of underground rock formations using a cable. Its principle is to install a logging tool with a cable in the wellbore, extend the cable from the wellhead to the bottom of the well, and then infer the physical properties and structure of the underground rock formations according to the resistance of the cable. The logging tools for wireline logging include a probe with a cable, a cable, and a cable output device. The probe is usually equipped with a sealed annular detector. According to the measurement requirements and the depth of the well, logging tools with different lengths and precisions can be selected.

[0121] Based on the logging data of Well 2, calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of Well 2 respectively;

[0122] In this example, according to the logging curve data such as acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS of Well 1 obtained in Step 1, calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of Well 1; according to the logging curve data such as acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS of Well 2 obtained in Step 1, calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of Well 2.

[0123] Specifically, use the density logging DEN to determine the total porosity PorT of the reservoir, and calculate the total porosity PorT based on the following calculation formula:

[0124]

[0125] In the above calculation formula:

[0126] ρ marepresents the reservoir matrix density, with the unit of grams per cubic centimeter (g / cm3);

[0127] ρ b represents the reservoir bulk density, with the unit of grams per cubic centimeter (g / cm3);

[0128] ρ f represents the reservoir fluid density, with the unit of grams per cubic centimeter (g / cm3);

[0129] The acoustic travel time log AC is used to determine the reservoir matrix porosity PorB. Based on the following calculation formula, the total porosity PorB is calculated:

[0130]

[0131] In the above calculation formula:

[0132] Δt represents the acoustic travel time difference of the acoustic log, with the unit of microseconds per meter (μs / m);

[0133] Δt ma represents the acoustic travel time of the reservoir rock matrix, with the unit of microseconds per meter (μs / m);

[0134] Δt f represents the acoustic travel time of the reservoir fluid, with the unit of microseconds per meter (μs / m);

[0135] C P represents the compaction correction coefficient of the reservoir. Select 3.7 for old formations and 1 for new formations.

[0136] The deep lateral RD and shallow lateral data RS are used to calculate the fracture porosity PorFra; first, the fracture state is discriminated. Since the dual laterolog shows different curve characteristics for different fracture states, to solve for the fracture porosity PorFra, the following calculation formula is used to discriminate the fracture state:

[0137]

[0138] In the above calculation formula:

[0139] Y represents the discrimination index, dimensionless, and is the calculation result;

[0140] RD represents the deep lateral resistivity log value, with the unit of ohm-meter (Ω·m);

[0141] RS represents the shallow lateral resistivity log value, with the unit of ohm-meter (Ω·m).

[0142] During the process of full-well section treatment, the Y value is first calculated, and then the state of the fracture is discriminated. When Y > 0.1, it is a high-angle fracture; when 0.1 ≥ Y > 0, it is an inclined fracture; when Y < 0, it is a low-angle fracture. After the fracture state is discriminated, the following interpretation model is used to calculate the fracture porosity:

[0143]

[0144] In the above calculation formula:

[0145] PORFra represents the fracture porosity, and the unit is a decimal;

[0146] Rmf represents the resistivity of the mud filtrate, and the unit is ohm-meter (Ω·m);

[0147] A1, A2, and A3 are constants, and their values vary according to different fracture states. Based on the value of the above discrimination index Y, the values of A1, A2, and A3 are taken, as shown in Table 1.

[0148] Step 3: Based on the logging data of Well 1, as well as the calculation results of the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra of Well 1, calculate the development index RECP of Well 1;

[0149] Based on the logging data of Well 2, as well as the calculation results of the total porosity PorT, the matrix porosity PorB, and the fracture porosity PorFra of Well 2, calculate the development index RECP of Well 2.

[0150] In this embodiment, based on the following calculation formula, calculate the development index RECP of Well 1 and Well 2 respectively;

[0151]

[0152] In the formula, RD is the deep lateral logging data, PorT is the total porosity data, PorB is the matrix porosity data, PorFra is the fracture porosity data, and RS is the shallow lateral logging data.

[0153] Define the reservoir development levels of Well 1 and Well 2 respectively through the values of the development index RECP of Well 1 and Well 2,

[0154] When RECP > 1, the development level of the reservoir is defined as Class I;

[0155] When 0.1 < RECP < 1, the development level of the reservoir is defined as Class II;

[0156] When 0.02 < RECP < 0.1, the development level of the reservoir is defined as Class III.

[0157] Step 4: Calculate the main control factor index ZKYS of Well 1 based on the measurement results of the logging data of Well 1.

[0158] Calculate the main control factor index ZKYS of Well 2 based on the measurement results of the logging data of Well 2.

[0159] In this embodiment, calculate the main control factor index ZKYS of Well 1 and Well 2 based on the following calculation formula;

[0160]

[0161] In the formula, AC is the acoustic logging data, RS is the shallow lateral logging data, and RD is the deep lateral logging data.

[0162] Define the fracture development characteristics through the value of the main control factor index ZKYS.

[0163] When ZKYS > 0.04, the fracture development characteristics of the reservoir are defined as high angle;

[0164] When ZKYS < -0.04, the fracture development characteristics of the reservoir are defined as low angle;

[0165] When -0.04 < ZKYS < 0.04, the fracture development characteristics of the reservoir are defined as medium-high angle.

[0166] After the above four steps, the logging data, total porosity PorT, matrix porosity PorB, fracture porosity PorFra, development index RECP, and main control factor index ZKYS of Well 1 are as Figure 2 and shown in Table 2;

[0167] The logging data, total porosity PorT, matrix porosity PorB, fracture porosity PorFra, development index RECP, and main control factor index ZKYS of Well 2 are as Figure 3 and shown in Table 2.

[0168] Table 2 Logging data table of Well 1 and Well 2

[0169]

[0170]

[0171] Step 5: As Figure 4 shown, cross the calculated development index RECP of Well 1 with the main control factor index ZKYS of Well 1; cross the calculated development index RECP of Well 2 with the main control factor index ZKYS of Well 2; obtain two cross plots respectively, and evaluate the development types of Well 1 and Well 2 based on the two cross plots.

[0172] In this step, based on the crossplot technology, an evaluation method for the development type of the fractured reservoir is established. Specifically, the crossplot technology refers to a mapping interpretation technology for well logging data. It intersects two types of well logging data on a plane map and determines the value or range of the required parameter based on the coordinates of the intersection point.

[0173] In one embodiment, Figure 5 As shown, after step five, the method further includes:

[0174] Step six, obtaining a result map based on microresistivity imaging logging, and using the result map to verify and interpret the result of evaluating the development type of the reservoirs of Well 1 and Well 2 based on the intersection map.

[0175] Specifically, microresistivity imaging logging refers to the use of multiple rows of button-shaped small electrodes on a multi-electrode plate to emit current to the well wall formation. Due to the different rock composition, structure and fluid contained in the electrodes, the current changes, which reflects the change in rock resistivity at various locations on the well wall, and can display resistivity imaging of the well wall. A result map is obtained according to the microresistivity imaging logging method, and the result map can be used to verify whether the development type of the fracture reservoir evaluated based on the intersection map is correct.

[0176] Embodiment 3

[0177] In this embodiment, Figure 6 As shown, a fracture reservoir type evaluation device is provided, comprising:

[0178] Well logging data acquisition module 210, used to measure and acquire well logging data of fractured reservoirs;

[0179] A porosity calculation module 220, for respectively calculating the total porosity PorT, matrix porosity PorB and fracture porosity PorFra of the fractured reservoir based on the well logging data of the fractured reservoir;

[0180] A development index calculation module 230 is used to calculate a development index RECP based on the well logging data of the fractured reservoir and the calculation results of the total porosity PorT, the matrix porosity PorB and the fracture porosity PorFra;

[0181] A main controlling factor index calculation module 240 is used to calculate a main controlling factor index ZKYS based on the measurement results of the well logging data of the fractured reservoir;

[0182] The evaluation result obtaining module 250 is used to intersect the calculated development index RECP with the main controlling factor index ZKYS to obtain an intersection diagram, and evaluate the development type of the fracture reservoir based on the intersection diagram.

[0183] In one embodiment, the logging data acquisition module includes:

[0184] An acoustic logging unit for measuring the logging curve of acoustic logging AC based on wireline logging;

[0185] A density logging unit for measuring the logging curve of density logging DEN based on wireline logging;

[0186] A deep lateral logging unit for measuring the logging curve of deep lateral logging RD based on wireline logging;

[0187] A shallow lateral logging unit for measuring the logging curve of shallow lateral logging RS based on wireline logging.

[0188] In one embodiment, the porosity calculation module includes:

[0189] A total porosity calculation unit for calculating the total porosity PorT of the fractured reservoir based on the logging data of the fractured reservoir;

[0190] A matrix porosity calculation unit for calculating the matrix porosity PorB of the fractured reservoir based on the logging data of the fractured reservoir;

[0191] A fracture porosity calculation unit for calculating the fracture porosity PorFra of the fractured reservoir based on the logging data of the fractured reservoir.

[0192] In one embodiment, the development index calculation module is used to calculate the development index RECP based on the following calculation formula;

[0193]

[0194] Wherein, RD is the deep lateral logging data, PorT is the total porosity data, PorB is the matrix porosity data, PorFra is the fracture porosity data, and RS is the shallow lateral logging data.

[0195] In one embodiment, the main control factor index calculation module is used to calculate the main control factor index ZKYS based on the following calculation formula;

[0196]

[0197] Wherein, AC is the acoustic logging data, RS is the shallow lateral logging data, and RD is the deep lateral logging data.

[0198] In one embodiment, the evaluation result acquisition module includes:

[0199] The intersection unit is used to perform an intersection of the calculated development index RECP and the main control factor index ZKYS based on crossplot technology to obtain a crossplot.

[0200] The development type evaluation unit is used to evaluate the development type of the fractured reservoir based on the crossplot.

[0201] In one embodiment, the fractured reservoir type evaluation device further includes a microresistivity imaging module.

[0202] The microresistivity imaging module is used to obtain a result map based on microresistivity imaging logging, and use the result map to verify and interpret the result of the fractured reservoir type evaluation method.

[0203] For the specific limitations of the fractured reservoir type evaluation device, reference may be made to the limitations of the fractured reservoir type evaluation method in the above text, which will not be elaborated here. Each unit in the above-mentioned fractured reservoir type evaluation device can be implemented in whole or in part by software, hardware, and their combination. The above-mentioned units can be embedded in the processor of the computer device in hardware form or be independent of it, or can be stored in the memory of the computer device in software form, so as to facilitate the processor to call and execute the operations corresponding to the above units.

[0204] Embodiment 4

[0205] In this embodiment, a computer device is provided. Its internal structure diagram can be as Figure 7 shown. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program, and a database is deployed on the non-volatile storage medium. The database is used to store logging data, total porosity PorT, matrix porosity PorB, fracture porosity PorFra, development index RECP, and main control factor index ZKYS. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with other computer devices on which application software is deployed. When the computer program is executed by the processor, it realizes a method for processing goods inbound and outbound data. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen. The input device of the computer device can be a touch layer covered on the display screen, or a button, a trackball, or a touchpad set on the shell of the computer device, or an external keyboard, a touchpad, or a mouse, etc.

[0206] Those skilled in the art can understand, Figure 7The structure shown is only a block diagram of some structures related to the solution of this application, and does not constitute a limitation on the computer device to which the solution of this application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have different component arrangements.

[0207] In one embodiment, a computer device is provided, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the following steps are implemented:

[0208] Measure and obtain well logging data of a fractured reservoir;

[0209] Based on the well logging data of the fractured reservoir, calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of the fractured reservoir respectively;

[0210] Based on the well logging data of the fractured reservoir and the calculation results of the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra, calculate the development index RECP;

[0211] Based on the measurement results of the well logging data of the fractured reservoir, calculate the main control factor index ZKYS;

[0212] Intersect the calculated development index RECP with the main control factor index ZKYS to obtain an intersection diagram, and evaluate the development type of the fractured reservoir based on the intersection diagram.

[0213] In one embodiment, when the processor executes the computer program, the following steps are also implemented:

[0214] Measure and obtain the well logging data of the fractured reservoir based on the wireline logging method.

[0215] In one embodiment, the well logging data of the fractured reservoir includes well logging curves of acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS.

[0216] In one embodiment, when the processor executes the computer program, the following steps are also implemented:

[0217] Based on the following calculation formula, calculate the development index RECP;

[0218]

[0219] In the formula, RD is the deep lateral logging data, PorT is the total porosity data, PorB is the matrix porosity data, PorFra is the fracture porosity data, and RS is the shallow lateral logging data.

[0220] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0221] Based on the following calculation formula, calculate the main control factor index ZKYS;

[0222]

[0223] In the formula, AC is the acoustic logging data, RS is the shallow lateral logging data, and RD is the deep lateral logging data.

[0224] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0225] Based on the crossplot technology, establish an evaluation method for the development type of the fractured reservoir.

[0226] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0227] Obtain the result map based on the micro-resistivity imaging logging method, and use the result map to verify and interpret the result of the evaluation method for the fractured reservoir type.

[0228] Embodiment Five

[0229] In this embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented:

[0230] Measure and obtain the logging data of the fractured reservoir;

[0231] Based on the logging data of the fractured reservoir, calculate the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of the fractured reservoir respectively;

[0232] Based on the logging data of the fractured reservoir and the calculation results of the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra, calculate the development index RECP;

[0233] Based on the measurement results of the logging data of the fractured reservoir, calculate the main control factor index ZKYS;

[0234] Cross the calculated development index RECP with the main control factor index ZKYS to obtain a crossplot, and evaluate the development type of the fractured reservoir based on the crossplot.

[0235] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0236] Measure and obtain the logging data of the fractured reservoir based on the wireline logging method.

[0237] In one embodiment, the well logging data of the fractured reservoir includes well logging curves of acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS.

[0238] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0239] Based on the following calculation formula, the development index RECP is calculated;

[0240]

[0241] In the formula, RD is deep lateral logging data, PorT is total porosity data, PorB is matrix porosity data, PorFra is fracture porosity data, and RS is shallow lateral logging data.

[0242] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0243] Based on the following calculation formula, the main control factor index ZKYS is calculated;

[0244]

[0245] In the formula, AC is acoustic logging data, RS is shallow lateral logging data, and RD is deep lateral logging data.

[0246] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0247] Based on crossplot technology, an evaluation method for the development type of the fractured reservoir is established.

[0248] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0249] Based on the micro-resistivity imaging logging method, a result map is obtained, and the result of the evaluation method for the fractured reservoir type is verified and interpreted by using the result map.

[0250] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, storage, database, or other medium used in the various embodiments provided in the present application can include non-volatile and / or volatile memories. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0251] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope described in this specification.

[0252] The above-described embodiments only represent several implementation manners of the present application. The description is relatively specific and detailed, but it should not be construed as a limitation on the scope of the invention patent. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can still be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application should be subject to the appended claims.

Claims

1. A method for evaluating fracture reservoir types, characterized in that including: Measuring and obtaining well logging data of a fractured reservoir; Based on the well logging data of the fractured reservoir, calculating the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of the fractured reservoir respectively; Based on the well logging data of the fractured reservoir and the calculation results of the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra, calculating the development index RECP; Based on the measurement results of the well logging data of the fractured reservoir, calculating the main controlling factor index ZKYS; Intersecting the calculated development index RECP with the main controlling factor index ZKYS to obtain a crossplot, and evaluating the development type of the fractured reservoir based on the crossplot.

2. The method according to claim 1, characterized in that, The step of measuring and obtaining well logging data of the fractured reservoir includes: Measuring and obtaining the well logging data of the fractured reservoir based on the cable well logging method.

3. The method according to claim 1, wherein The well logging data of the fractured reservoir includes well logging curves of acoustic logging AC, density logging DEN, deep lateral logging RD, and shallow lateral logging RS.

4. The method according to claim 3, characterized in that, In the step of calculating the development index RECP based on the well logging data of the fractured reservoir and the calculation results of the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra, the development index RECP is calculated based on the following calculation formula; In the formula, RD is the deep lateral logging data, PorT is the total porosity data, PorB is the matrix porosity data, PorFra is the fracture porosity data, and RS is the shallow lateral logging data.

5. The method according to claim 3, characterized in that In the step of calculating the main controlling factor index ZKYS based on the measurement results of the well logging data of the fractured reservoir, the main controlling factor index ZKYS is calculated based on the following calculation formula; In the formula, AC is the acoustic logging data, RS is the shallow lateral logging data, and RD is the deep lateral logging data.

6. The method according to claim 1, characterized in that, The step of intersecting the calculated development index RECP with the main controlling factor index ZKYS to obtain a crossplot, and evaluating the development type of the fractured reservoir based on the crossplot includes: Based on the crossplot technology, establishing an evaluation method for the development type of the fractured reservoir.

7. The method according to claim 1, wherein After the step of intersecting the calculated development index RECP with the main controlling factor index ZKYS to obtain a crossplot, and evaluating the development type of the fractured reservoir based on the crossplot, it further includes: Obtaining a result map based on the micro-resistivity imaging logging method, and using the result map to verify and interpret the result of the evaluation method for the fractured reservoir type.

8. An evaluation device for fracture reservoir types, characterized in that, including: A well logging data acquisition module for measuring and obtaining well logging data of a fractured reservoir; A porosity calculation module for calculating the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra of the fractured reservoir respectively based on the well logging data of the fractured reservoir; A development index calculation module for calculating the development index RECP based on the well logging data of the fractured reservoir and the calculation results of the total porosity PorT, matrix porosity PorB, and fracture porosity PorFra; The master control factor index calculation module is used to calculate the master control factor index ZKYS based on the measurement results of the logging data of the fractured reservoir; The evaluation result acquisition module is used to cross the calculated development index RECP with the master control factor index ZKYS to obtain a cross plot, and evaluate the development type of the fractured reservoir based on the cross plot.

9. A computer device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images