A method and device for obtaining an interwell fracture network model of a fractured-porous carbonate reservoir

Abstract

The application discloses a kind of interwell fracture network model obtaining methods of fracture pore type carbonate reservoir, comprising: obtaining the initial form of the interwell fracture network of the fracture pore type carbonate reservoir;Correct the characteristic parameter in the interwell fracture network, obtain single well water cut curve objective function;According to the single well water cut curve objective function, obtain the interwell fracture network model.The present application considers the development characteristics of fracture and pore, uses the water cut curve in production process as fitting target, and based on the influence relationship between single well water cut curve and interwell fracture information, the shape characteristics of interwell fracture network, permeability and equivalent cross-sectional area and other parameters are quantified, the possibility of fracture related information multiple solution is greatly reduced, and the accuracy and efficiency of interwell fracture network modeling are improved.

Description

Technical Field

[0001] This invention relates to the field of petroleum extraction technology, and in particular to a method and apparatus for obtaining an inter-well fracture network model of fractured pore carbonate reservoirs. Background Technology

[0002] Fractured and porous carbonate reservoirs are characterized by a wide variety of fracture occurrences. The formation of fractures increases matrix permeability to some extent, thereby improving well production and positively impacting reservoir development. However, the presence of inter-well fracture networks exacerbates reservoir heterogeneity, causing injected water to surge along fractures during waterflooding, resulting in directional water channeling between oil and water wells. This reduces the waterflood sweep efficiency, leading to development challenges such as rapid water cut increases and significant differences in well productivity. Therefore, understanding fracture characteristics provides crucial guidance for efficient reservoir development, making the accurate quantification of fracture networks essential.

[0003] Currently, the commonly used methods for quantitative analysis of fracture networks mainly include the Catherine effect method, principal curvature method, stress field method, and neural network method. Each method has its own advantages and disadvantages and its own applicable scope. For example, (1) Catherine effect method: This method can obtain the maximum stress experienced by the rock in history, rather than the current stress. Although there is a difference between the two, since the formation of fractures is mainly related to the maximum stress in history, it is more accurate to use this method to study fractures. (2) Principal curvature method: This method can simulate the entire reservoir layer using mathematical and physical methods to obtain the distribution of tectonic stress field. The application of this method has limitations, namely, the reservoir thickness should not exceed 50m, and the internal structure of the block should not be too intense. If there is a fault, the block simulation should be carried out with the fault as the boundary. (3) Stress field method: The geostress field is mainly composed of the coupling of gravity stress, tectonic stress, pore pressure, thermal stress, etc. On the basis of completing the geostress field calculation, methods such as Coulomb discrimination criterion, Griffith discrimination criterion, and energy method discrimination criterion can be used to quantitatively predict fractures. (4) Neural network method: This method is a highly effective and feasible intelligent method. It does not require human intervention in the process of analyzing problems and calculations. It has advantages that other methods cannot match, such as high vertical resolution, large horizontal measurable range, low cost and high promotion value. The establishment of this method needs to be combined with well logging data. After obtaining the mathematical model of the neural network and analyzing the indication of fractures by well logging data, the neural network and well logging data can be used to identify reservoir fractures.

[0004] Meanwhile, fractures exhibit changes in their opening or closing state during production. Furthermore, the development of fractures and pores in the same space can alter the conductivity of fractures. Therefore, it is necessary to consider these characteristics and establish a simple and accurate inter-well fracture network modeling method based on dynamic data to replace the aforementioned quantification method, thereby meeting the requirements for the identification and inversion of inter-well fracture networks. Summary of the Invention

[0005] The purpose of this invention is to provide a method for obtaining an inter-well fracture network model of fractured porous carbonate reservoirs, so as to accurately quantify the inter-well fracture network of fractured porous carbonate reservoirs.

[0006] To achieve the above objectives, the present invention provides a method for obtaining an inter-well fracture network model of fractured porous carbonate reservoirs, comprising: The initial morphology of the inter-well fracture network in the fractured porous carbonate reservoir was obtained; By correcting the fracture characteristic parameters in the inter-well fracture network, the objective function for the water cut curve of a single well is obtained. The inter-well fracture network model is obtained based on the objective function of the single-well water cut curve.

[0007] Optionally, obtaining the initial morphology of the inter-well fracture network in the fractured porous carbonate reservoir includes: A numerical model of fractured and porous carbonate reservoirs was obtained. Initial morphology of the inter-well fracture network is added to the injection-production wells in the numerical model of the fractured-porous carbonate reservoir to obtain the initial morphology of the inter-well fracture network of the fractured-porous carbonate reservoir.

[0008] Optionally, the inter-well fracture network information refers to the fracture characteristic parameters in the inter-well fracture network.

[0009] Optionally, the fracture characteristic parameters in the inter-well fracture network include: fracture morphology characteristic parameters, fracture equivalent cross-sectional area, and fracture permeability.

[0010] Optionally, the fracture characteristic parameters in the inter-well fracture network are modified to obtain the objective function of the single-well water cut curve, including: By correcting the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area of ​​the fracture network between wells, the objective function of the single-well water cut curve is obtained.

[0011] Optionally, the objective function for the single-well water cut curve is obtained by modifying the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area in the inter-well fracture network, including: By correcting the fracture morphology characteristic parameters in the inter-well fracture network, an objective function for the upward trend of the water cut curve of a single well is obtained. By correcting the fracture permeability in the inter-well fracture network, an objective function is obtained to determine the rise time of the water cut curve for a single well. The objective function for obtaining the maximum value of the water cut curve of a single well is obtained by modifying the equivalent cross-sectional area of ​​the fractures in the inter-well fracture network. The objective function for the single-well water cut curve is obtained based on the objective function of the upward trend of the single-well water cut curve, the objective function of the upward time of the single-well water cut curve, and the objective function of the maximum value of the single-well water cut curve.

[0012] Optionally, the method further includes: performing dynamic index verification on the correctness of the obtained inter-well fracture network model.

[0013] This invention also provides a device for obtaining an inter-well fracture network model of fractured porous carbonate reservoirs, comprising: A morphology acquisition unit is used to obtain the initial morphology of the inter-well fracture network in the fractured porous carbonate reservoir. The function acquisition unit is used to correct the characteristic parameters in the inter-well fracture network and obtain the objective function of the single-well water cut curve. The model acquisition unit is used to obtain the inter-well fracture network model based on the objective function of the single-well water cut curve.

[0014] Optionally, obtaining the initial morphology of the inter-well fracture network in the fractured porous carbonate reservoir includes: A numerical model of fractured and porous carbonate reservoirs was obtained. Initial morphology of the inter-well fracture network is added to the injection-production wells in the numerical model of the fractured-porous carbonate reservoir to obtain the initial morphology of the inter-well fracture network of the fractured-porous carbonate reservoir.

[0015] Optionally, the inter-well fracture network information refers to the fracture characteristic parameters in the inter-well fracture network.

[0016] Optionally, the fracture characteristic parameters in the inter-well fracture network include: fracture morphology characteristic parameters, fracture equivalent cross-sectional area, and fracture permeability.

[0017] Optionally, the fracture characteristic parameters in the inter-well fracture network are modified to obtain the objective function of the single-well water cut curve, including: By correcting the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area of ​​the fracture network between wells, the objective function of the single-well water cut curve is obtained.

[0018] Optionally, the objective function for the single-well water cut curve is obtained by modifying the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area in the inter-well fracture network, including: By correcting the fracture morphology characteristic parameters in the inter-well fracture network, an objective function for the upward trend of the water cut curve of a single well is obtained. By correcting the fracture permeability in the inter-well fracture network, an objective function is obtained to determine the rise time of the water cut curve for a single well. The objective function for obtaining the maximum value of the water cut curve of a single well is obtained by modifying the equivalent cross-sectional area of ​​the fractures in the inter-well fracture network. The objective function for the single-well water cut curve is obtained based on the objective function of the upward trend of the single-well water cut curve, the objective function of the upward time of the single-well water cut curve, and the objective function of the maximum value of the single-well water cut curve.

[0019] The technical effects and advantages of this invention are as follows: This invention discloses a method for obtaining an inter-well fracture network model of fractured-porosity carbonate reservoirs, comprising: obtaining the initial morphology of the inter-well fracture network of the fractured-porosity carbonate reservoir; correcting the characteristic parameters in the inter-well fracture network to obtain a single-well water cut curve objective function; and obtaining the inter-well fracture network model based on the single-well water cut curve objective function. This invention considers the development characteristics of fractures and pores, uses the water cut curve during production as the fitting target, and quantifies the morphological characteristics, permeability, and equivalent cross-sectional area of ​​the inter-well fracture network based on the influence relationship between the single-well water cut curve and inter-well fracture information. This significantly reduces the possibility of multiple solutions related to fractures and improves the accuracy and efficiency of inter-well fracture network modeling.

[0020] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures pointed out in the description and the drawings. Attached Figure Description

[0021] Figure 1 Flowchart of the method for obtaining the inter-well fracture network model of fractured porous carbonate reservoirs; Figure 2 This is an initial morphological diagram of the inter-well fracture network in INJ-PRO. Figure 3 Diagram of the INJ-PRO inter-well fracture network model; Figure 4 This is a graph showing the water cut changes in PRO wells during the deep displacement phase. Figure 5 A structural diagram of the device was obtained for the inter-well fracture network model of fractured porous carbonate reservoirs. Detailed Implementation

[0022] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0023] To address the shortcomings of existing technologies, this invention discloses a method for obtaining an inter-well fracture network model of fractured porous carbonate reservoirs, such as... Figure 1 As shown. The method includes: obtaining the initial morphology of the inter-well fracture network in the fractured-porosity carbonate reservoir; correcting the characteristic parameters in the inter-well fracture network to obtain the objective function of the single-well water cut curve; and obtaining the inter-well fracture network model based on the objective function of the single-well water cut curve. This invention considers the development characteristics of fractures and pores, uses the water cut curve during the production process as the fitting target, and quantifies the morphological characteristics, permeability, and equivalent cross-sectional area of ​​the inter-well fracture network based on the influence relationship between the single-well water cut curve and the inter-well fracture information. This significantly reduces the possibility of multiple solutions related to fractures and improves the accuracy and efficiency of inter-well fracture network modeling.

[0024] To better understand this approach, the following section provides a detailed explanation of the specific steps involved in obtaining the inter-well fracture network model for fractured porous carbonate reservoirs.

[0025] 1. Based on the numerical model of fractured porous carbonate reservoirs, an initial-form inter-well fracture network is added between injection and production wells.

[0026] It should be noted that the inclusion of the initial well-to-well fracture network refers to adding fracture information from the well-to-well fracture network into the numerical model. This fracture information includes: fracture connectivity, equivalent cross-sectional area, and permeability. The fracture connectivity refers to the spatial network coordinates and depth at both ends of the fracture.

[0027] In fractured-porosity carbonate reservoirs, the co-development of fractures and pores in the same space leads to an increase in the cross-sectional area of ​​the fracture system. The equivalent cross-sectional area of ​​the fractures, proposed to describe this increase, is defined as follows: (1) In the formula, A This represents the equivalent cross-sectional area of ​​the crack. n The number of holes that develop in the same space as the crack; r The average radius of the pores that develop in the same space as the cracks; b The opening of the crack; h The height of the crack; l This refers to the distance between injection and production wells.

[0028] It should also be noted that the initial form of the inter-well fracture network refers to the initial value of fracture information in the inter-well fracture network. Since it is only the initial value of fracture information, researchers can determine it artificially based on their own understanding or geological research results, and this step does not need to pursue the accuracy of the initial form.

[0029] 2. Correct the morphological characteristic parameters in the inter-well fracture network and fit the upward trend of the water cut curve of a single well.

[0030] It should be noted that the inter-well fracture network morphology refers to the spatial structure of the inter-well fracture network at the locations where each fracture is connected. The characteristic parameters of the inter-well fracture network morphology refer to the locations where each fracture is connected.

[0031] In the process of correcting the above characteristic parameters, the connection positions of each fracture are adjusted to achieve changes in the shape of different fracture lengths and angles, thereby fitting the upward trend of the single-well water cut curve.

[0032] When fitting the upward trend of the water cut curve of a single well, the objective function for fitting is as follows: (2) In the formula, F 1 represents the objective function for fitting the upward trend segment of the water cut curve of a single well; f 1k , f wk The first k The actual moisture content at each moment and the moisture content obtained by numerical calculation; k Use time sequence numbers; it is recommended that one time be considered as one day. K 1 represents the number of moments in the upward trend segment of the water cut curve for a single well; it is recommended to use the number of days in the upward trend segment.

[0033] Based on fundamental theoretical research, this invention divides the single-well water cut curve into an initial horizontal segment, an upward trend segment, and a maximum horizontal segment, corresponding to three stages: the initial state of the reservoir, the start of water production in the inter-well fracture network, and the full water production of the inter-well fracture network. The upward trend segment of the single-well water cut curve refers to the stage in which the actual water cut curve of the single well rises from the initial horizontal segment to the maximum horizontal segment.

[0034] Considering the common data fluctuations in the actual water cut curves of single wells, the fitting accuracy requires the objective function to... F 1 is less than 0.05.

[0035] Previous studies have found that crack length is the main controlling factor affecting the upward trend of the moisture content curve. The recommended fitting criterion is: when the rate of increase of the moisture content calculated by numerical calculation is slower than that of the actual moisture content, adjust the crack connection position to make the length of each crack in the network smaller; otherwise, increase the crack length.

[0036] 3. Correct the fracture permeability in the inter-well fracture network and fit the rise time of the water cut curve of a single well.

[0037] It should be noted that the fracture permeability in the inter-well fracture network refers to the permeability of each fracture in the fracture network.

[0038] In the process of correcting fracture permeability, the rise time of the single-well water cut curve is fitted by adjusting the permeability of fractures at different locations. The rise time of the single-well water cut curve refers to the starting moment of the continuous rise phase of the single-well water cut curve. The continuous rise phase refers to the period of sustained water cut increase within 3-6 months.

[0039] When fitting the rise time of the water cut curve of a single well, the objective function for fitting is as follows: (3) In the formula, F 2 is the objective function for fitting the rise time of the water cut curve of a single well. The fitting accuracy requires the objective function to... F 2 is less than 25; f 1u , f wu These represent the rise time of the actual moisture content and the moisture content curve obtained from numerical calculation, respectively.

[0040] Previous studies have found that crack permeability is the main controlling factor affecting the rise time of the moisture content curve. The recommended fitting criterion is: when the rise time of the moisture content calculated by numerical calculation is later than that of the actual moisture content, increase the crack permeability; otherwise, decrease the crack permeability.

[0041] 4. Correct the equivalent cross-sectional area of ​​fractures in the inter-well fracture network and fit the maximum value of the water cut curve of a single well.

[0042] It should be noted that when correcting the equivalent cross-sectional area of ​​fractures in the inter-well fracture network, the maximum value of the single-well water cut curve is fitted by adjusting the cross-sectional area of ​​fractures at different locations. The maximum value of the single-well water cut curve refers to the water cut value corresponding to the end of the continuous rising phase of the single-well water cut curve.

[0043] When fitting the maximum value of the water cut curve of a single well, the objective function is as follows: (4) In the formula, F 3 represents the objective function for fitting the maximum value of the water cut curve of a single well, and the required fitting accuracy is equal to the objective function. F 3 is less than 0.04; f 1m , f wm These represent the actual moisture content and the maximum value of the moisture content curve obtained from numerical calculation, respectively.

[0044] Previous studies have found that the equivalent cross-sectional area of ​​cracks is the main controlling factor affecting the maximum value of the moisture content curve. The recommended fitting criterion is: when the maximum value of the moisture content calculated by numerical calculation is smaller than the actual moisture content, the equivalent cross-sectional area of ​​cracks should be increased; otherwise, the equivalent cross-sectional area of ​​cracks should be decreased.

[0045] 5. Repeat steps 2-4, fine-tuning the fracture location, permeability, and equivalent cross-sectional area to obtain the best single-well water cut fitting effect.

[0046] Previous studies have found that the upward trend, rise time, and maximum value of the single-well water cut curve are influenced not only by the main controlling factor but also by other fracture information. For example, the main controlling factor for the upward trend of the single-well water cut curve is the morphology of the inter-well fracture network, but it is also affected by fracture permeability and equivalent cross-sectional area. Therefore, steps 2-4 need to be repeated to fine-tune the fracture location, permeability, and equivalent cross-sectional area to obtain the best fitting effect for the single-well water cut. The objective function at this time is: (5) In the formula, F 4 represents the objective function for fitting the water cut curve of a single well, and the required fitting accuracy is equal to the objective function. F 1 is less than 4; c 1. c 2 and c 3 represents the weighting coefficients of the objective function. Previous studies have found that the fracture uplift trend and the corresponding inter-well fracture network morphology parameters have the greatest impact on the final fitting results. Therefore, the recommended weighting coefficients are as follows: c 1 = 0.7 c 2 = 0.1 c 3 = 0.2.

[0047] 6. Conduct dynamic index verification to confirm the correctness of the inter-well fracture network model.

[0048] After obtaining the inter-well fracture network model by fitting the water cut curve, dynamic index verification is required to determine the correctness of the model.

[0049] It should be noted that the dynamic index verification refers to calculating the dynamic index and comparing it with the actual value of the dynamic index. Based on previous research results, the dynamic index refers to the daily liquid production rate during the waterflooding stage or the water cut variation pattern during deep subsurface regulation.

[0050] To address the challenge of quantifying inter-well fracture networks in fractured and porous carbonate reservoirs, this invention, based on the influence relationship between single-well water cut curves and inter-well fracture information, summarizes the characteristics of the water cut curve as its upward trend, rise time, and maximum value, which correspond to the morphological characteristic parameters, permeability, and equivalent cross-sectional area of ​​the inter-well fracture network, respectively. This essentially achieves a one-to-one correspondence between curve characteristics and inversion information, significantly eliminating the ambiguity in fracture information determination and improving the accuracy and efficiency of inter-well fracture network modeling.

[0051] To better explain this solution, embodiments are also provided below.

[0052] The steps for establishing the INJ-PRO inter-well fracture network model for fractured porous carbonate reservoirs are as follows. Here, INJ represents injection wells, and PRO represents production wells.

[0053] First, based on the numerical model of fractured porous carbonate reservoirs, an initial-form inter-well fracture network is added between injection and production wells, such as... Figure 2 As shown in Table 1, the information on fractures in the inter-well fracture network is as follows.

[0054] Table 1. Information on fractures in the INJ-PRO inter-well fracture network.

[0055] Then, steps 2-5 of the above method are executed to correct the morphological characteristic parameters, permeability, and equivalent cross-sectional area of ​​the inter-well fracture network, and to fine-tune the parameters, thereby obtaining the inter-well fracture network model and corresponding parameters, such as... Figure 3 As shown in Table 2.

[0056] Table 2. Corresponding information of the established INJ-PRO inter-well fracture network

[0057] Finally, dynamic index verification was conducted to confirm the correctness of the inter-well fracture network model.

[0058] It should be noted that, comparing the water cut change pattern during the deep displacement stage of INJ-PRO, the simulation showed that the water cut of the PRO well decreased by 12.3 percentage points during the deep displacement stage, which is close to the actual water cut decrease of 10.5 percentage points. This verifies the correctness of the established inter-well fracture network. Figure 4 As shown.

[0059] This invention also provides a device for obtaining an inter-well fracture network model of fractured porous carbonate reservoirs, such as... Figure 5 As shown. The device includes: a morphology acquisition unit for acquiring the initial morphology of the inter-well fracture network in the fractured-pore carbonate reservoir; a function acquisition unit for correcting the characteristic parameters in the inter-well fracture network to obtain the objective function of the water cut curve of a single well; and a model acquisition unit for obtaining the model of the inter-well fracture network based on the objective function of the water cut curve of a single well.

[0060] Since the content protected by this device corresponds to the content protected by the method described above, it will not be described in detail here. Please refer to the explanation section of the method described above for details.

[0061] Finally, it should be noted that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for obtaining an interwell fracture network model of a fractured-porous carbonate reservoir well, characterized in that, include: The initial morphology of the inter-well fracture network in the fractured porous carbonate reservoir was obtained; By correcting the fracture characteristic parameters in the inter-well fracture network, the objective function for the water cut curve of a single well is obtained. Based on the objective function of the single-well water cut curve, the inter-well fracture network model is obtained; Correcting the fracture characteristic parameters in the inter-well fracture network to obtain the objective function of the single-well water cut curve includes: By correcting the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area of ​​the inter-well fracture network, the objective function of the single-well water cut curve is obtained. By correcting the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area of ​​the inter-well fracture network, the objective function for obtaining the water cut curve of the single well is obtained, including: By correcting the fracture morphology characteristic parameters in the inter-well fracture network, an objective function for the upward trend of the water cut curve of a single well is obtained. When fitting the upward trend of the water cut curve of a single well, the objective function for fitting is as follows: In the formula, F 1 is the objective function of fitting the rising trend section of the single well water cut curve; f 1k , f wk are the real water cut and the water cut obtained by numerical calculation at the i th moment, respectively; k is the moment number; k is the moment number; K 1 is the number of moments of the rising trend section of the single well water cut curve; By correcting the fracture permeability in the inter-well fracture network, an objective function is obtained to determine the rise time of the water cut curve for a single well. When fitting the rise time of the water cut curve of a single well, the objective function for fitting is as follows: In the formula, F 2 represents the objective function for fitting the rise time of the water cut curve of a single well; f 1u , f wu The rise time of the moisture content curves is shown for the actual moisture content and the moisture content curves obtained from numerical calculations, respectively. The objective function for obtaining the maximum value of the water cut curve of a single well is obtained by modifying the equivalent cross-sectional area of ​​the fractures in the inter-well fracture network. When fitting the maximum value of the water cut curve of a single well, the objective function is as follows: In the formula, F 3 represents the objective function for fitting the maximum value of the water cut curve of a single well; f 1m , f wm These represent the actual moisture content and the maximum value of the moisture content curve obtained from numerical calculation, respectively. The objective function of the single well water cut curve is obtained based on the objective function of the upward trend of the single well water cut curve, the objective function of the upward time of the single well water cut curve, and the objective function of the maximum value of the single well water cut curve. The objective function of the single-well water cut curve F 4 is represented as: In the formula, c 1. c 2 and c 3 represents the weight coefficients of the objective function.

2. The method according to claim 1, characterized in that, Obtaining the initial morphology of the inter-well fracture network in the fractured-pore carbonate reservoir includes: A numerical model of fractured and porous carbonate reservoirs was obtained. Initial morphology of the inter-well fracture network is added to the injection-production wells in the numerical model of the fractured-porous carbonate reservoir to obtain the initial morphology of the inter-well fracture network of the fractured-porous carbonate reservoir.

3. The method according to claim 2, characterized in that, The inter-well fracture network information refers to the fracture characteristic parameters in the inter-well fracture network.

4. The method according to claim 3, characterized in that, The fracture characteristic parameters in the inter-well fracture network include: fracture morphology characteristic parameters, fracture equivalent cross-sectional area, and fracture permeability.

5. The method according to claim 1, characterized in that, The method further includes: performing dynamic index verification on the correctness of the obtained inter-well fracture network model.

6. A device for obtaining an inter-well fracture network model of fractured porous carbonate oil reservoirs, characterized in that, include: A morphology acquisition unit is used to obtain the initial morphology of the inter-well fracture network in the fractured porous carbonate reservoir. The function acquisition unit is used to correct the characteristic parameters in the inter-well fracture network and obtain the objective function of the single-well water cut curve. The model acquisition unit is used to obtain the inter-well fracture network model based on the objective function of the single-well water cut curve; Correcting the fracture characteristic parameters in the inter-well fracture network to obtain the objective function of the single-well water cut curve includes: By correcting the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area of ​​the inter-well fracture network, the objective function of the single-well water cut curve is obtained. By correcting the fracture morphology parameters, fracture permeability, and equivalent cross-sectional area of ​​the inter-well fracture network, the objective function for obtaining the water cut curve of the single well is obtained, including: By correcting the fracture morphology characteristic parameters in the inter-well fracture network, an objective function for the upward trend of the water cut curve of a single well is obtained. When fitting the upward trend of the water cut curve of a single well, the objective function for fitting is as follows: In the formula, F 1 represents the objective function for fitting the upward trend segment of the water cut curve of a single well; f 1k , f wk The first k The actual moisture content at each moment and the moisture content obtained by numerical calculation; k For time sequence number; K 1 represents the number of moments in the upward trend segment of the water cut curve for a single well; By correcting the fracture permeability in the inter-well fracture network, an objective function is obtained to determine the rise time of the water cut curve for a single well. When fitting the rise time of the water cut curve of a single well, the objective function for fitting is as follows: In the formula, F 2 represents the objective function for fitting the rise time of the water cut curve of a single well; f 1u , f wu The rise time of the moisture content curves is shown for the actual moisture content and the moisture content curves obtained from numerical calculations, respectively. The objective function for obtaining the maximum value of the water cut curve of a single well is obtained by modifying the equivalent cross-sectional area of ​​the fractures in the inter-well fracture network. When fitting the maximum value of the water cut curve of a single well, the objective function is as follows: In the formula, F 3 represents the objective function for fitting the maximum value of the water cut curve of a single well; f 1m , f wm These represent the actual moisture content and the maximum value of the moisture content curve obtained from numerical calculation, respectively. The objective function of the single well water cut curve is obtained based on the objective function of the upward trend of the single well water cut curve, the objective function of the upward time of the single well water cut curve, and the objective function of the maximum value of the single well water cut curve. The objective function of the single-well water cut curve F 4 is represented as: In the formula, c 1. c 2 and c 3 represents the weight coefficients of the objective function.

7. The apparatus according to claim 6, characterized in that, Obtaining the initial morphology of the inter-well fracture network in the fractured-pore carbonate reservoir includes: A numerical model of fractured and porous carbonate reservoirs was obtained. Initial morphology of the inter-well fracture network is added to the injection-production wells in the numerical model of the fractured-porous carbonate reservoir to obtain the initial morphology of the inter-well fracture network of the fractured-porous carbonate reservoir.

8. The apparatus according to claim 7, characterized in that, The inter-well fracture network information refers to the fracture characteristic parameters in the inter-well fracture network.

9. The apparatus according to claim 8, characterized in that, The fracture characteristic parameters in the inter-well fracture network include: fracture morphology characteristic parameters, fracture equivalent cross-sectional area, and fracture permeability.

Images