Dynamic constraint fault control fracture-cavity type carbonate reservoir unit division method and device

Abstract

The invention relates to the technical field of seismic data interpretation, and particularly discloses a dynamically constrained fault control fracture-vuggy carbonate reservoir unit division method and device, and the method comprises the steps: carrying out the data processing of static seismic data, and obtaining the processed seismic data; based on the processed seismic data, strike-slip fracture depiction and fracture-cavity reservoir depiction are carried out respectively; dynamically constraining data of the strike-slip fracture depiction and the fracture-cavity reservoir depiction to obtain a dynamic constraint basis; performing oil reservoir unit division based on the dynamic constraint basis to obtain an oil reservoir unit division result; and guiding well pattern optimization deployment according to an oil reservoir unit division result. According to the method, strike-slip fracture and fault-control fracture-cavity reservoir feature parameters are extracted based on static data, fault-control fracture-cavity type carbonate reservoir unit boundary preliminary depiction is carried out, the reservoir unit boundary is further corrected in combination with dynamic data constraint on the basis, the reservoir unit boundary depiction precision is improved, well pattern optimization deployment is guided, and the reservoir unit boundary depiction efficiency is improved. And a technical support is provided for efficient development of the fault control fracture-vuggy carbonate reservoir.

Description

Technical Field

[0001] This invention relates to the field of seismic data interpretation technology, specifically to a method and apparatus for dynamically constrained fault-controlled fracture-vuggy carbonate reservoir unit division. Background Technology

[0002] Fault-controlled fracture-vuggy carbonate reservoirs are characterized by complex internal structures and strong heterogeneity, making the delineation of these reservoir units challenging. Current methods primarily rely on static 3D seismic data, using seismic attributes as the main basis. However, existing techniques suffer from low accuracy and poor compatibility with dynamic data, failing to effectively guide well network optimization in these reservoirs and exhibiting limited applicability. This innovative technology extracts strike-slip faults and characteristic parameters of fault-controlled fracture-vuggy reservoirs from static data to initially delineate the boundaries of fault-controlled fracture-vuggy carbonate reservoir units. Further refinement of these boundaries using dynamic data constraints significantly improves the accuracy of the delineation, closely matching actual development results and providing strong technical support for the efficient development of these reservoirs.

[0003] Based on this technical background, this invention studies a method and apparatus for dividing fractured-vuggy carbonate reservoir units under dynamic constraints. Summary of the Invention

[0004] To address the shortcomings of existing technologies, this invention provides a method and apparatus for dynamically constrained fault-controlled fracture-vuggy carbonate reservoir unit delineation. This method extracts characteristic parameters of strike-slip fractures and fault-controlled fracture-vuggy reservoirs based on static data, and conducts preliminary delineation of the boundaries of fault-controlled fracture-vuggy carbonate reservoir units. On this basis, it further combines dynamic data constraints to correct the reservoir unit boundaries, improves the accuracy of reservoir unit boundary delineation, guides the optimized deployment of well networks, and provides technical support for the efficient development of fault-controlled fracture-vuggy carbonate oil and gas reservoirs.

[0005] To achieve the above objectives, a first aspect of the present invention provides a dynamically constrained method for dividing fractured-vuggy carbonate reservoir units, comprising:

[0006] Processed seismic data is obtained by processing static seismic data.

[0007] Based on the processed seismic data, strike-slip faults and fracture-vuggy reservoirs are characterized respectively.

[0008] Dynamic constraints are applied to the data of strike-slip fracture characterization and fracture-vuggy reservoir characterization to obtain the dynamic constraint basis;

[0009] Based on the aforementioned dynamic constraints, reservoir unit division results are obtained by dividing the reservoir units.

[0010] The results of the reservoir unit division described above will guide the optimization and deployment of the well network.

[0011] A second aspect of the present invention provides a dynamically constrained fault-controlled fracture-vuggy carbonate reservoir unit delineation device, comprising:

[0012] The data processing module is used to process static seismic data to obtain processed seismic data.

[0013] The characterization module is used to characterize strike-slip faults and fracture-vuggy reservoirs based on the processed seismic data.

[0014] The dynamic constraint module is used to dynamically constrain the data of the strike-slip fracture characterization and the fractured-vuggy reservoir characterization to obtain the dynamic constraint basis.

[0015] The unit partitioning module is used to partition reservoir units based on the dynamic constraint criteria to obtain reservoir unit partitioning results.

[0016] The well network deployment module is used to guide the optimized deployment of the well network based on the reservoir unit division results.

[0017] A third aspect of the present invention provides an electronic device, the electronic device comprising:

[0018] Memory, which stores executable instructions;

[0019] A processor that executes the executable instructions in the memory to implement the dynamic constraint fault-controlled fracture-vuggy carbonate reservoir unit division method described in the first aspect.

[0020] A fourth aspect of the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the dynamic constraint fault-controlled fracture-vuggy carbonate reservoir unit division method described in the first aspect.

[0021] The beneficial effects of this invention include:

[0022] (1) The dynamic constraint method for dividing fault-controlled fracture-vuggy carbonate reservoir units proposed in this invention extracts the characteristic parameters of strike-slip faults and fault-controlled fracture-vuggy reservoirs based on static data, and conducts preliminary delineation of the boundaries of fault-controlled fracture-vuggy carbonate reservoir units. On this basis, the boundaries of reservoir units are further corrected by combining dynamic data constraints, thereby improving the accuracy of reservoir unit boundary delineation, guiding the optimized deployment of well networks, and providing technical support for the efficient development of fault-controlled fracture-vuggy carbonate oil and gas reservoirs.

[0023] (2) The dynamic constraint method for dividing fault-controlled fracture-vuggy carbonate reservoir units proposed in this invention takes into account the characteristics of fault-controlled fracture-vuggy carbonate rocks, starts from the characteristics of fracture and fracture-vuggy reservoirs, and obtains the static characterization basis of fault-controlled fracture-vuggy carbonate reservoir units. At the same time, dynamic data is used as the constraint condition of reservoir units, which improves the accuracy of reservoir unit division and has a high degree of consistency with the actual development effect. It can effectively guide the actual development of fault-controlled fracture-vuggy carbonate reservoirs.

[0024] Other features and advantages of the present invention will be described in detail in the following detailed description section. Attached Figure Description

[0025] 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.

[0026] Figure 1 This is a flowchart illustrating the dynamic constraint-based fault-controlled fracture-vuggy carbonate reservoir unit division method proposed in this invention.

[0027] Figure 2 This is a flowchart illustrating a specific implementation of the dynamic constraint-controlled fracture-vuggy carbonate reservoir unit division method proposed in this invention.

[0028] Figure 3 This is a schematic diagram illustrating the superposition of coherent top coherence and strike-slip fault characterization in a specific embodiment of the dynamic constraint fault-controlled fracture-vuggy carbonate reservoir unit division method proposed in this invention.

[0029] Figure 4 This is a schematic diagram of the porosity distribution in the target layer fractured-vuggy reservoir in a specific embodiment of the dynamic constraint fractured-vuggy carbonate reservoir unit division method proposed in this invention.

[0030] Figure 5 This is a schematic diagram of the oil pressure (top) and production (bottom) of a drilled well in a specific embodiment of the dynamic constraint fracture-vuggy carbonate reservoir unit division method proposed in this invention.

[0031] Figure 6 This is a schematic diagram of the carbonate reservoir unit division plan in a specific embodiment of the dynamic constraint fault-controlled fracture-vuggy carbonate reservoir unit division method proposed in this invention. Detailed Implementation

[0032] 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.

[0033] This invention provides a dynamically constrained method for dividing fractured-vuggy carbonate reservoir units, such as... Figure 1 As shown, it includes:

[0034] Processed seismic data is obtained by processing static seismic data.

[0035] Based on the processed seismic data, strike-slip faults and fracture-vuggy reservoirs were characterized respectively.

[0036] Dynamic constraints are obtained by applying dynamic constraints to the data characterizing strike-slip faults and fracture-vuggy reservoirs;

[0037] The reservoir unit division results are obtained by dividing the reservoir units based on dynamic constraints.

[0038] The results of reservoir unit division guide the optimization and deployment of well networks.

[0039] In this invention, based on static data, characteristic parameters of strike-slip fractures and fault-controlled fracture-vuggy reservoirs are extracted to conduct preliminary delineation of the boundaries of fault-controlled fracture-vuggy carbonate reservoir units. On this basis, dynamic data constraints are further combined to correct the reservoir unit boundaries, improve the accuracy of reservoir unit boundary delineation, guide the optimized deployment of well networks, and provide technical support for the efficient development of fault-controlled fracture-vuggy carbonate oil and gas reservoirs.

[0040] According to the present invention, data processing of static seismic data to obtain processed seismic data includes:

[0041] Using 3D seismic data and well logging data, a synthetic record calibration is performed to determine the seismic geological strata.

[0042] Processed seismic data are obtained by interpreting the target stratigraphic layers of the seismic geological horizon.

[0043] The accuracy of the layer interpretation is 2 lines × 2 channels.

[0044] According to the present invention, strike-slip fault characterization based on processed seismic data includes:

[0045] By processing seismic data, coherence volume properties are calculated, and coherence plane properties of the top surface of the target layer are extracted.

[0046] Based on the coherence plane properties of the top surface of the target layer, the distribution of strike-slip fault planes is characterized, and the distribution of fault planes is used as the basis for dividing reservoir units.

[0047] According to the present invention, characterizing fractured-vuggy reservoirs based on processed seismic data includes:

[0048] By processing seismic data, post-stack sparse pulse inversion is performed to obtain the wave impedance volume;

[0049] Using known porosity and wave impedance data from well logging, cross-plot analysis is performed to obtain the wave impedance porosity function;

[0050] Based on the wave impedance volume and the wave impedance porosity function, the porosity volume is obtained;

[0051] Based on the processed seismic data and porosity volume, extract the fracture-vuggy reservoir properties of the target layer;

[0052] Based on the fractured-vuggy reservoir properties of the target layer, porosity contour lines were drawn. According to the well logging interpretation results, the contour line with a porosity value of 1.8% was taken as the lower limit of the effective reservoir.

[0053] These contour lines are used as the reservoir basis for dividing reservoir units.

[0054] According to the present invention, the dynamic constraint basis for dynamically constraining the data of strike-slip fracture characterization and fracture-vuggy reservoir characterization includes:

[0055] Based on the drilling and testing data, the inter-well connectivity analysis was conducted on the stroke-slip fault characterization and fracture-vuggy reservoir characterization to obtain connected well groups and independent single wells. The results of the connectivity analysis were used as the dynamic constraint basis for reservoir unit division.

[0056] Preferably, the reservoir unit division based on dynamic constraints includes:

[0057] ① Dynamic constraints are based on the highest priority level, classifying dynamically connected well groups into the same reservoir unit, and conversely, dynamically independent single wells should be treated as independent reservoir units.

[0058] ② Based on the priority of reservoirs, and on the basis of satisfying the division criteria ①, the same porosity contour line boundary is divided into one reservoir unit, and different porosity contour line boundaries are divided into different reservoir units;

[0059] ③ The fracture criteria have the lowest priority. Based on satisfying criteria ① and ②, the same fracture fragment is classified as one reservoir unit, and different fracture fragments are classified as different reservoir units.

[0060] According to the present invention, the reservoir unit division result obtained based on dynamic constraint criteria includes:

[0061] Based on the criteria for reservoir unit division, the fault criteria for reservoir unit division, the reservoir criteria for reservoir unit division, and the dynamic constraint criteria for reservoir unit division, reservoir unit division results are obtained.

[0062] The principles of well network optimization deployment are:

[0063] ① For existing well-controlled reservoir units, according to the principle of 1 km well spacing in the plane and different strata in the vertical direction, the well network is densified in the plane and in the vertical three-dimensional density to realize the construction of a three-dimensional injection and production well network, timely replenish the energy of the reservoir unit, and improve the recovery rate of single wells and reservoir units;

[0064] ② For reservoir units without well control, a three-dimensional well network is constructed in one go according to the principle of 1-kilometer well spacing in the plane and different strata in the vertical direction, so as to realize the efficient utilization of the reservoir unit.

[0065] In this invention, based on the characteristics of fault-controlled fracture-cavity carbonate rocks, the static characterization basis of fault-controlled fracture-cavity carbonate rock reservoir units is obtained from the characteristics of fracture and fracture-cavity reservoirs. At the same time, dynamic data is used as the constraint condition of reservoir units, which improves the accuracy of reservoir unit division and has a high degree of consistency with the actual development effect. It can effectively guide the actual development of fault-controlled fracture-cavity carbonate rock reservoirs.

[0066] The present invention will be described in more detail below through embodiments.

[0067] Example 1:

[0068] like Figure 2 As shown in the figure, this embodiment proposes a dynamically constrained method for dividing fault-controlled fracture-vuggy carbonate reservoir units. First, it establishes a description of strike-slip faults and fault-controlled fracture-vuggy reservoir characteristics using 3D seismic data. Then, it uses characteristic parameters to achieve static characterization of fault-controlled fracture-vuggy carbonate reservoir units. Furthermore, it uses well dynamic data to constrain and correct the boundaries of fault-controlled fracture-vuggy carbonate reservoir units, achieving a refined characterization of these units, improving the accuracy of unit division, and guiding the efficient development of fault-controlled fracture-vuggy carbonate reservoirs.

[0069] The specific implementation steps of this method are as follows:

[0070] Step 1: Basic Data Processing

[0071] 1.1. Using 3D seismic data and well logging data, perform synthetic record calibration to determine seismic geological strata;

[0072] 1.2. Perform horizon interpretation of the target layer on the 3D seismic data, with an interpretation accuracy of 2 lines × 2 traces;

[0073] Step 2: Characterization of the strike-slip fracture system:

[0074] 2.1 Calculate coherence volume properties using 3D seismic data volume;

[0075] 2.2. Based on the stratigraphic interpretation results obtained in step 1.2, extract the coherence plane properties of the top surface of the target layer;

[0076] 2.3. Based on the coherence plane properties of the top surface of the target layer obtained in step 2.2, characterize the distribution of the strike-slip fracture plane, such as... Figure 3 As shown, the fracture plane distribution is used as the fracture basis for dividing reservoir units;

[0077] Step 3: Characterization of fractured-vuggy reservoirs:

[0078] 3.1 Using the 3D seismic data volume, post-stack sparse pulse inversion is performed to obtain the wave impedance volume;

[0079] 3.2 Using the known porosity and wave impedance data from well logging, cross-plot analysis is performed to obtain the wave impedance porosity function;

[0080] 3.3. Based on the wave impedance volume obtained in step 3.1 and the wave impedance porosity function obtained in step 3.2, the porosity volume is obtained;

[0081] 3.4. Based on the stratigraphic interpretation results obtained in step 1.2 and the porosity volume obtained in step 3.3, extract the fracture-vuggy reservoir properties of the target layer, such as... Figure 4 As shown;

[0082] 3.5. Draw porosity contour lines for the target layer fractured-vuggy reservoir attributes obtained in step 3.4. Based on the well logging interpretation results, take the contour line with a porosity value of 1.8% as the lower limit of the effective reservoir, and use this contour line as the reservoir basis for dividing the reservoir unit.

[0083] Step 4: Dynamic Data Analysis

[0084] 4.1 Based on the drilling and production test data, well connectivity analysis was conducted to identify connected well groups and independent single wells. The connectivity analysis results were used as the dynamic constraint basis for reservoir unit division. Since the oil pressure values ​​and trends of wells W1, W3, and H2 are consistent, these three wells constitute a connected well group and belong to the same reservoir unit. Well M1 has significantly different oil pressure and production compared to the aforementioned three wells, and is not connected to them, thus constituting an independent single well. Figure 5 As shown;

[0085] Step 5: Reservoir Unit Division:

[0086] 5.1 The criteria for dividing reservoir units are as follows: ① Dynamic constraints have the highest priority. Well groups that are dynamically connected are divided into the same reservoir unit, while single wells that are dynamically independent should be considered as independent reservoir units; ② Reservoir criteria have the next highest priority. On the basis of satisfying the criteria in ①, wells within the same porosity contour line boundary are divided into one reservoir unit, and wells with different porosity contour lines are divided into different reservoir units; ③ Fracture criteria have the lowest priority. On the basis of satisfying the criteria in ① and ②, the same fracture section should be divided into one reservoir unit, and different fracture sections should be divided into different reservoir units.

[0087] 5.2. Based on the reservoir unit division criteria in step 5.1, and according to the fault division criteria obtained in step 2.3, the reservoir division criteria in step 3.5, and the dynamic constraint criteria for reservoir unit division obtained in step 4.1, dynamic constraint fault-controlled fracture-vuggy carbonate reservoir unit division is performed. Based on the dynamic constraint that W1, W3, and H2 are connected well groups belonging to the same reservoir unit, and according to both reservoir and fault criteria, reservoir unit ① is determined in the south. Since M1 is an independent single well, reservoir unit ② is determined in the central area based on both reservoir and fault criteria. Since there is no well control in the north, reservoir unit ③ is determined in the north based on both reservoir and fault criteria. Figure 6 As shown;

[0088] Step Six: Well Pattern Optimization Deployment:

[0089] 6.1 Based on the reservoir unit division results obtained in step 5.2, guide the optimized deployment of the well network;

[0090] 6.2 The deployment principles are as follows: ① For existing well-controlled reservoir units, well networks can be densified in both planar and vertical directions according to the principle of 1 km well spacing in the plane and different strata in the vertical direction, so as to realize the construction of a three-dimensional injection-production well network, timely replenish the energy of the reservoir unit, and improve the recovery rate of single wells and reservoir units; 3 new wells can be deployed in the southern No. 1 reservoir unit and 2 new wells can be deployed in the central No. 2 reservoir unit; ② For reservoir units without well control, a three-dimensional well network can be constructed in one go according to the principle of 1 km well spacing in the plane and different strata in the vertical direction, so as to realize the efficient utilization of the reservoir unit; the northern No. 3 reservoir unit can form a well network in one go and deploy 3 new wells.

[0091] Example 2:

[0092] This embodiment provides a dynamically constrained method for dividing fractured-vuggy carbonate reservoir units, such as... Figure 1 As shown, it includes:

[0093] Processed seismic data is obtained by processing static seismic data.

[0094] Based on the processed seismic data, strike-slip faults and fracture-vuggy reservoirs were characterized respectively.

[0095] Dynamic constraints are obtained by applying dynamic constraints to the data characterizing strike-slip faults and fracture-vuggy reservoirs;

[0096] The reservoir unit division results are obtained by dividing the reservoir units based on dynamic constraints.

[0097] The results of reservoir unit division guide the optimization and deployment of well networks;

[0098] In this embodiment, the data processing of static seismic data to obtain processed seismic data includes:

[0099] Using 3D seismic data and well logging data, a synthetic record calibration is performed to determine the seismic geological strata.

[0100] Processed seismic data are obtained by interpreting the target stratigraphic layers of the seismic geological horizon.

[0101] The accuracy of the stratigraphic interpretation is 2 lines × 2 channels;

[0102] In this embodiment, the strike-slip fault characterization based on processed seismic data includes:

[0103] By processing seismic data, coherence volume properties are calculated, and coherence plane properties of the top surface of the target layer are extracted.

[0104] Based on the coherence plane properties of the top surface of the target layer, the distribution of strike-slip fault planes is characterized, and the distribution of fault planes is used as the basis for dividing reservoir units.

[0105] In this embodiment, characterizing fractured-vuggy reservoirs based on processed seismic data includes:

[0106] By processing seismic data, post-stack sparse pulse inversion is performed to obtain the wave impedance volume;

[0107] Using known porosity and wave impedance data from well logging, cross-plot analysis is performed to obtain the wave impedance porosity function;

[0108] Based on the wave impedance volume and the wave impedance porosity function, the porosity volume is obtained;

[0109] Based on the processed seismic data and porosity volume, extract the fracture-vuggy reservoir properties of the target layer;

[0110] Based on the fractured-vuggy reservoir properties of the target layer, porosity contour lines were drawn. According to the well logging interpretation results, the contour line with a porosity value of 1.8% was taken as the lower limit of the effective reservoir.

[0111] This contour line is used as the reservoir basis for dividing reservoir units;

[0112] In this embodiment, the dynamic constraint basis for dynamically constraining the data of strike-slip fracture characterization and fracture-vuggy reservoir characterization includes:

[0113] Based on the drilling and testing data, the inter-well connectivity analysis was conducted on the stroke-slip fault characterization and fracture-vuggy reservoir characterization to obtain connected well groups and independent single wells. The results of the connectivity analysis were used as the dynamic constraint basis for reservoir unit division.

[0114] In this embodiment, the reservoir unit division based on dynamic constraints includes:

[0115] ① Dynamic constraints are based on the highest priority level, classifying dynamically connected well groups into the same reservoir unit, and conversely, dynamically independent single wells should be treated as independent reservoir units.

[0116] ② Based on the priority of reservoirs, and on the basis of satisfying the division criteria ①, the same porosity contour line boundary is divided into one reservoir unit, and different porosity contour line boundaries are divided into different reservoir units;

[0117] ③ The fracture criteria have the lowest priority. Based on satisfying criteria ① and ②, the same fracture section is classified as one reservoir unit, and different fracture sections are classified as different reservoir units.

[0118] In this embodiment, the reservoir unit division result obtained based on dynamic constraints includes:

[0119] Based on the criteria for reservoir unit division, the fault criteria for reservoir unit division, the reservoir criteria for reservoir unit division, and the dynamic constraint criteria for reservoir unit division, reservoir unit division results are obtained.

[0120] The principles of well network optimization deployment are:

[0121] ① For existing well-controlled reservoir units, according to the principle of 1 km well spacing in the plane and different strata in the vertical direction, the well network is densified in the plane and in the vertical three-dimensional density to realize the construction of a three-dimensional injection and production well network, timely replenish the energy of the reservoir unit, and improve the recovery rate of single wells and reservoir units;

[0122] ② For reservoir units without well control, a three-dimensional well network is constructed in one go according to the principle of 1-kilometer well spacing in the plane and different strata in the vertical direction, so as to realize the efficient utilization of the reservoir unit.

[0123] Example 3:

[0124] This embodiment provides a dynamically constrained fault-controlled fracture-vuggy carbonate reservoir unit delineation device, including:

[0125] The data processing module is used to process static seismic data to obtain processed seismic data.

[0126] The characterization module is used to characterize strike-slip faults and fracture-vuggy reservoirs based on processed seismic data.

[0127] The dynamic constraint module is used to dynamically constrain the data of strike-slip fracture characterization and fracture-vuggy reservoir characterization to obtain the dynamic constraint basis;

[0128] The unit partitioning module is used to partition reservoir units based on dynamic constraints to obtain reservoir unit partitioning results.

[0129] The well network deployment module is used to guide the optimized deployment of well networks based on the reservoir unit division results.

[0130] In this embodiment, the data processing of static seismic data to obtain processed seismic data includes:

[0131] Using 3D seismic data and well logging data, a synthetic record calibration is performed to determine the seismic geological strata.

[0132] Processed seismic data are obtained by interpreting the target stratigraphic layers of the seismic geological horizon.

[0133] The accuracy of the stratigraphic interpretation is 2 lines × 2 channels;

[0134] In this embodiment, the strike-slip fault characterization based on processed seismic data includes:

[0135] By processing seismic data, coherence volume properties are calculated, and coherence plane properties of the top surface of the target layer are extracted.

[0136] Based on the coherence plane properties of the top surface of the target layer, the distribution of strike-slip fault planes is characterized, and the distribution of fault planes is used as the basis for dividing reservoir units.

[0137] In this embodiment, characterizing fractured-vuggy reservoirs based on processed seismic data includes:

[0138] By processing seismic data, post-stack sparse pulse inversion is performed to obtain the wave impedance volume;

[0139] Using known porosity and wave impedance data from well logging, cross-plot analysis is performed to obtain the wave impedance porosity function;

[0140] Based on the wave impedance volume and the wave impedance porosity function, the porosity volume is obtained;

[0141] Based on the processed seismic data and porosity volume, extract the fracture-vuggy reservoir properties of the target layer;

[0142] Based on the fractured-vuggy reservoir properties of the target layer, porosity contour lines were drawn. According to the well logging interpretation results, the contour line with a porosity value of 1.8% was taken as the lower limit of the effective reservoir.

[0143] This contour line is used as the reservoir basis for dividing reservoir units;

[0144] In this embodiment, the dynamic constraint basis for dynamically constraining the data of strike-slip fracture characterization and fracture-vuggy reservoir characterization includes:

[0145] Based on the drilling and testing data, the inter-well connectivity analysis was conducted on the stroke-slip fault characterization and fracture-vuggy reservoir characterization to obtain connected well groups and independent single wells. The results of the connectivity analysis were used as the dynamic constraint basis for reservoir unit division.

[0146] In this embodiment, the reservoir unit division based on dynamic constraints includes:

[0147] ① Dynamic constraints are based on the highest priority level, classifying dynamically connected well groups into the same reservoir unit, and conversely, dynamically independent single wells should be treated as independent reservoir units.

[0148] ② Based on the priority of reservoirs, and on the basis of satisfying the division criteria ①, the same porosity contour line boundary is divided into one reservoir unit, and different porosity contour line boundaries are divided into different reservoir units;

[0149] ③ The fracture criteria have the lowest priority. Based on satisfying criteria ① and ②, the same fracture section is classified as one reservoir unit, and different fracture sections are classified as different reservoir units.

[0150] In this embodiment, the reservoir unit division result obtained based on dynamic constraints includes:

[0151] Based on the criteria for reservoir unit division, the fault criteria for reservoir unit division, the reservoir criteria for reservoir unit division, and the dynamic constraint criteria for reservoir unit division, reservoir unit division results are obtained.

[0152] The principles of well network optimization deployment are:

[0153] ① For existing well-controlled reservoir units, according to the principle of 1 km well spacing in the plane and different strata in the vertical direction, the well network is densified in the plane and in the vertical three-dimensional density to realize the construction of a three-dimensional injection and production well network, timely replenish the energy of the reservoir unit, and improve the recovery rate of single wells and reservoir units;

[0154] ② For reservoir units without well control, a three-dimensional well network is constructed in one go according to the principle of 1-kilometer well spacing in the plane and different strata in the vertical direction, so as to realize the efficient utilization of the reservoir unit.

[0155] Example 4:

[0156] This invention provides an electronic device including a memory and a processor, comprising:

[0157] Memory, which stores executable instructions;

[0158] The processor runs executable instructions in memory to implement a method for dividing fractured-vuggy carbonate reservoir units.

[0159] 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.

[0160] 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 the invention, the processor is used to execute computer-readable instructions stored in the memory.

[0161] Those skilled in the art should 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 invention.

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

[0163] Example 5:

[0164] This invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements a method for dividing fractured-vuggy carbonate reservoir units.

[0165] A computer-readable storage medium according to embodiments of the present invention 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 invention are performed.

[0166] 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).

[0167] The embodiment of this invention proposes a dynamic constraint method for dividing fault-controlled fracture-vuggy carbonate reservoir units. Based on static data, it extracts characteristic parameters of strike-slip faults and fault-controlled fracture-vuggy reservoirs, and conducts preliminary delineation of the boundaries of fault-controlled fracture-vuggy carbonate reservoir units. On this basis, it further combines dynamic data constraints to correct the reservoir unit boundaries, improves the accuracy of reservoir unit boundary delineation, guides the optimized deployment of well networks, and provides technical support for the efficient development of fault-controlled fracture-vuggy carbonate oil and gas reservoirs.

[0168] 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 dividing fractured-vuggy carbonate reservoir units under dynamic constraints, characterized in that, include: Processed seismic data is obtained by processing static seismic data. Based on the processed seismic data, strike-slip faults and fracture-vuggy reservoirs are characterized respectively. Dynamic constraints are applied to the data of strike-slip fracture characterization and fracture-vuggy reservoir characterization to obtain the dynamic constraint basis; Based on the aforementioned dynamic constraints, reservoir unit division is performed to obtain reservoir unit division results; The results of the reservoir unit division described above will guide the optimization and deployment of the well network.

2. The method according to claim 1, characterized in that, Processing static seismic data yields processed seismic data, including: Using 3D seismic data and well logging data, a synthetic record calibration is performed to determine the seismic geological strata. Processed seismic data are obtained by interpreting the target stratigraphic layers of the seismic geological layers. The accuracy of the layer interpretation is 2 lines × 2 channels.

3. The method according to claim 2, characterized in that, Characterizing strike-slip faults based on the processed seismic data includes: Using the processed seismic data, coherence volume properties are calculated, and coherence plane properties of the top surface of the target layer are extracted; Based on the coherence plane properties of the top surface of the target layer, the distribution of strike-slip fracture planes is characterized, and the distribution of fracture planes is used as the basis for dividing reservoir units.

4. The method according to claim 3, characterized in that, Characterizing fracture-vuggy reservoirs based on the processed seismic data includes: Using the processed seismic data, post-stack sparse pulse inversion is performed to obtain the wave impedance volume; Using known porosity and wave impedance data from well logging, cross-plot analysis is performed to obtain the wave impedance porosity function; Based on the wave impedance volume and the wave impedance porosity function, the porosity volume is obtained; Based on the processed seismic data and porosity volume, extract the fracture-vuggy reservoir properties of the target layer; Based on the fractured-vuggy reservoir properties of the target layer, porosity contour lines were drawn. According to the well logging interpretation results, the contour line with a porosity value of 1.8% was taken as the lower limit of the effective reservoir. These contour lines are used as the reservoir basis for dividing reservoir units.

5. The method according to claim 4, characterized in that, The dynamic constraint basis for applying dynamic constraints to the strike-slip fracture characterization and fracture-vuggy reservoir characterization data includes: Based on the drilling and testing data, the inter-well connectivity analysis is performed on the strike-slip fracture characterization and fractured-vuggy reservoir characterization to obtain connected well groups and independent single wells. The results of the connectivity analysis are used as the dynamic constraint basis for reservoir unit division.

6. The method according to claim 5, characterized in that, The division of reservoir units based on the aforementioned dynamic constraints includes: ① Dynamic constraints are based on the highest priority level, classifying dynamically connected well groups into the same reservoir unit, and conversely, dynamically independent single wells should be treated as independent reservoir units. ② Based on the priority of reservoirs, and on the basis of satisfying the division criteria ①, the same porosity contour line boundary is divided into one reservoir unit, and different porosity contour line boundaries are divided into different reservoir units; ③ The fracture criteria have the lowest priority. Based on satisfying criteria ① and ②, the same fracture fragment is classified as one reservoir unit, and different fracture fragments are classified as different reservoir units.

7. The method according to claim 6, characterized in that, Based on the aforementioned dynamic constraints, the reservoir unit division results include: Based on the criteria for dividing reservoir units, the fault criteria for dividing reservoir units, the reservoir criteria for dividing reservoir units, and the dynamic constraint criteria for dividing reservoir units, reservoir units are divided to obtain reservoir unit division results. The principle of well network optimization deployment is as follows: ① For existing well-controlled reservoir units, according to the principle of 1 km well spacing in the plane and different strata in the vertical direction, the well network is densified in the plane and in the vertical three-dimensional density to realize the construction of a three-dimensional injection and production well network, timely replenish the energy of the reservoir unit, and improve the recovery rate of single wells and reservoir units; ② For reservoir units without well control, a three-dimensional well network is constructed in one go according to the principle of 1-kilometer well spacing in the plane and different strata in the vertical direction, so as to realize the efficient utilization of the reservoir unit.

8. A dynamically constrained fault-controlled fracture-vuggy carbonate reservoir unit delineation device, characterized in that, include: The data processing module is used to process static seismic data to obtain processed seismic data. The characterization module is used to characterize strike-slip faults and fracture-vuggy reservoirs based on the processed seismic data. The dynamic constraint module is used to dynamically constrain the data of the strike-slip fracture characterization and the fractured-vuggy reservoir characterization to obtain the dynamic constraint basis. The unit partitioning module is used to partition reservoir units based on the dynamic constraint criteria to obtain reservoir unit partitioning results. The well network deployment module is used to guide the optimized deployment of the well network based on the reservoir unit division results.

9. 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 dynamic constraint method for dividing fractured-vuggy carbonate reservoir units according to any one of claims 1-7.

10. 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 fault-controlled fracture-vuggy carbonate reservoir unit division method according to any one of claims 1-7 with dynamic constraints.

Images