Interwell connectivity quantification and analysis method and device

Abstract

The invention discloses an inter-well connectivity quantification and analysis method and device, and relates to the technical field of water drive development, and the method comprises the steps: carrying out the interlayer division of an original model, and determining the distribution information of a water injection well and a producing well; an oil reservoir is divided into a plurality of grid units through a geometric segmentation method, and it is ensured that the water injection wells and the oil production wells are located in the grid units. And dividing a control area for each well by using a division rule, and determining the affiliation of the grid units through Euclidean distance comparison. And obtaining the conductivity between two adjacent grid units divided by the control boundary. And respectively evaluating the production and water injection capacities of the well by using the production index and the water injection index. And based on a finite volume method, constructing a mass conservation equation for the grid units, and carrying out iterative solution until convergence. And traversing the conductivity connectivity table to calculate a flow component, accumulating to obtain a total flow, deducing an inter-well connectivity coefficient according to the total flow, and realizing quantitative evaluation of the inter-well connectivity coefficient. The problems of how to optimize the injection-production well pattern design and how to improve the water drive development effect and the oil reservoir recovery efficiency are solved.

Description

Technical Field

[0001] This application relates to the technical field of water flooding development, and particularly to a method and device for quantifying and analyzing inter-well connectivity. Background Art

[0002] A key challenge in reservoir development is to efficiently manage the fluid flow in the underground reservoir to achieve optimal crude oil recovery.

[0003] The existing technologies for obtaining the layer system distribution information of a reservoir usually rely on geological exploration, logging data, and reservoir modeling techniques. These technologies provide a basic understanding of the reservoir structure, fluid properties, and distribution. However, when converting this information into the actual water flooding injection-production well pattern design, some challenges still exist. For example, there is a lack of sufficient fineness and dynamics in describing inter-well connectivity. The existing technologies often can only provide static and simplified descriptions of inter-well connectivity, and cannot accurately capture the actual fluid flow between wells. As a result, the optimization effect of the injection-production well pattern design is limited, which affects the water flooding development effect and the reservoir recovery rate.

[0004] Therefore, how to develop a more comprehensive and accurate method to analyze inter-well connectivity to optimize the injection-production well pattern design and improve the water flooding development effect and reservoir recovery rate has become an urgent problem to be solved. Summary of the Invention

[0005] In the embodiments of this application, by providing a method for quantifying and analyzing inter-well connectivity, the problem of how to develop a more comprehensive and accurate method to analyze inter-well connectivity to optimize the injection-production well pattern design and improve the water flooding development effect and reservoir recovery rate is solved.

[0006] In a first aspect, an embodiment of the present application provides a method for quantifying and analyzing inter-well connectivity, including: performing interlayer division on an original reservoir model to obtain distribution information of injection wells and production wells in a lateral layer system; processing the distribution information of injection wells and production wells in the lateral layer system by using a geometric segmentation method to divide the reservoir into multiple interconnected grid units, such that the injection wells and production wells are both located within the grid units; using a division rule to divide the control area for each injection well and production well; wherein, the control area of a well is a polygonal area surrounded by the intersection points of the control dividing lines of six pairs of injection wells or production wells around it; traversing all grid units, and determining whether the grid unit belongs to the control area of an injection well or a production well by comparing the Euclidean distance between the center point of each grid unit and all injection wells and production wells; obtaining the conductivity between two adjacent grid units divided by the control dividing line, and using a conductivity connection table to record the conductivity between all two adjacent grid units divided by the control dividing line; using a production index to represent the production capacity of a production well under a first pressure difference, and using a water injection index to represent the water injection capacity of an injection well under a second pressure difference; for all two adjacent grid units divided by the control dividing line, constructing a mass conservation equation according to the finite volume method, and iteratively solving the pressure vector and the total fluid saturation until the convergence condition is satisfied; traversing the conductivity connection table, and for all two adjacent grid units divided by the control dividing line, calculating the flow component between the injection well and the production well, accumulating all flow components, and calculating the total flow of the fluid between the injection well and the production well; obtaining an expression of the inter-well connectivity coefficient from the total flow calculation formula, and calculating the inter-well connectivity coefficient according to the expression of the inter-well connectivity coefficient.

[0007] In a possible implementation manner, the division rule includes: for any pair of injection well and production well, calculating the coordinates of their midpoint and the slope of the connecting line, and calculating the slope of the perpendicular bisector; obtaining the equation of the perpendicular bisector according to the slope of the perpendicular bisector, and using the perpendicular bisector as the control dividing line; the calculation formula for the coordinates of the midpoint is: ; wherein, is the coordinate of the midpoint , is the coordinate of the production well , is the coordinate of the injection well ; the calculation formula for the slope of the connecting line is: ; wherein, is the slope of the connecting line; the calculation formula for the slope of the perpendicular bisector is: ; wherein, is the slope of the perpendicular bisector; the equation of the perpendicular bisector is: ; wherein, and are variables in the equation of the perpendicular bisector, is the abscissa of any point in the horizontal layer system, is the ordinate of any point in the horizontal layer system.

[0008] In a possible implementation, traversing all grid cells and determining whether a grid cell belongs to the control area of an injection well or a production well by comparing the Euclidean distances between the center point of each grid cell and all injection wells and production wells includes: if the Euclidean distance from the center point of a grid cell to a certain injection well is less than its Euclidean distances to all production wells, it is determined that the grid cell belongs to the control area of the injection well; if the Euclidean distance from the center point of a grid cell to a certain production well is less than its Euclidean distances to all injection wells, it is determined that the grid cell belongs to the control area of the production well.

[0009] In a possible implementation, obtaining the conductivity between two adjacent grid cells divided by the control dividing line includes: the calculation formula for the conductivity between two adjacent grid cells divided by the control dividing line is: , , ; where is the conductivity between two adjacent grid cells divided by the control dividing line, and are intermediate variables, is the permeability of grid cell , is the permeability of grid cell , is the area of grid cell , is the area of grid cell , is the distance from the center point of grid cell to the control dividing line, is the distance from the center point of grid cell to the control dividing line, is the unit normal vector of the control dividing line pointing to grid cell , is the unit normal vector of the control dividing line pointing to grid cell , is the unit direction vector of the control dividing line from its center point to grid cell , is the unit direction vector of the control dividing line from its center point to grid cell .

[0010] In a possible implementation, using the production index to represent the production capacity of an oil production well under a first pressure difference and using the water injection index to represent the water injection capacity of a water injection well under a second pressure difference includes: The first pressure difference is: ; where is the first pressure difference, is the reservoir pressure, is the bottom hole pressure; The second pressure difference is: ; where is the second pressure difference; The calculation formula for the production index or the water injection index is: ; where is the production index or the water injection index , is the unit conversion factor, is the effective angle around the wellbore, is the reservoir capacity, is the pressure equivalent radius of the grid cell, is the actual radius of the wellbore, is the skin effect factor, is the logarithmic term, representing the flow resistance from the grid cell to the wellbore.

[0011] In a possible implementation, for all adjacent two grid cells divided by the controlled dividing line, constructing a mass conservation equation according to the finite volume method and iteratively solving the pressure vector and the total fluid saturation until the convergence condition is satisfied, includes: The mass conservation equation is: ; where is the fluid residual, is the grid cell 's porosity, is the grid cell 's fluid density, is the grid cell 's fluid 's saturation, is the time step 's change amount within, is the grid cell 's volume, is the mass flow rate flowing into the grid cell from the grid cell , is the set of grid cells adjacent to the grid cell , is the source-sink term; When the grid cell is an oil production well, the source-sink term is: ; where is the fluid 's relative permeability in the grid cell , is a fluid in the grid cell viscosity, is the production index ; when the grid cell is an injection well, the source-sink term is: , is the injection index ; from the grid cell flows into the adjacent grid cell mass flow rate The expression is: ; where is the conductivity between two adjacent grid cells separated by the controlled dividing line, is in the grid cell and the grid cell the average value of the mobility parameters between, is a fluid density, is a fluid relative permeability, is a fluid viscosity, represents the gravity term, is the acceleration due to gravity, is the height difference between adjacent grid cells, is the pressure difference between adjacent grid cells; ; where is the fluid compressibility, is the fluid density corresponding to the reference pressure under, is the pressure of the current grid cell, is the reference pressure, is the base of the natural logarithm; ; where is the fluid mobility parameter in the grid cell , is the fluid mobility parameter in the grid cell ; The pressure vector includes the pressures of all grid cells, and the total fluid saturation includes the saturations of the fluids in all grid cells.

[0012] In a possible implementation, traversing the conductivity connection table, for all adjacent two grid cells separated by the controlled dividing line, calculating the flow component between the injection well and the production well, accumulating all flow components, and calculating the total flow of the fluid between the injection well and the production well, including: The calculation formula for the flow component is: ; where is the flow component between two adjacent grid cells between the injection well and the production well; The calculation formula for the total flow is: ; wherein, is the total flow rate of the fluid between the injection well and the production well, is the control area of the injection well ; is the control area of the production well ; is the inter-well connectivity coefficient, is the average pressure within the control area of the injection well ; is the average pressure within the control area of the production well ;

[0013] In a possible implementation, the expression of the inter-well connectivity coefficient is obtained from the calculation formula of the total flow rate, and the inter-well connectivity coefficient is calculated according to the expression of the inter-well connectivity coefficient, including: the expression of the inter-well connectivity coefficient is: ; wherein, , , is the pressure of the grid cell ; is the pressure of the grid cell ; is the volume of the grid cell ; is the volume of the grid cell ; is the total volume of the control area of the injection well ; is the total volume of the control area of the production well ;

[0014] In a second aspect, an inter-well connectivity quantification and analysis device according to an embodiment of the present application includes: an acquisition distribution information module, configured to perform interlayer division on an original reservoir model to acquire the distribution information of injection wells and production wells in a lateral layer series; a grid cell division module, configured to process the distribution information of injection wells and production wells in the lateral layer series by using a geometric segmentation method, divide the reservoir into a plurality of interconnected grid cells, and enable both the injection wells and the production wells to be located within the grid cells; a control area division module, configured to use division rules to divide the control area for each injection well and production well; wherein, the control area of a well is a polygonal area surrounded by the intersection points of the control demarcation lines of six pairs of injection wells or production wells around it; a determination module, configured to traverse all grid cells, and determine whether the grid cell belongs to the control area of an injection well or the control area of a production well by comparing the Euclidean distances between the center point of each grid cell and all injection wells and production wells; an acquisition conductivity module, configured to acquire the conductivity between two adjacent grid cells divided by the control demarcation line, and record the conductivity between all two adjacent grid cells divided by the control demarcation line by using a conductivity connection table; a representation module, configured to represent the production capacity of a production well under a first pressure difference by using a production index, and represent the injection capacity of an injection well under a second pressure difference by using an injection index; an iterative solution module, configured to construct a mass conservation equation according to the finite volume method for all two adjacent grid cells divided by the control demarcation line, and iteratively solve the pressure vector and the total fluid saturation until the convergence condition is satisfied; a total flow calculation module, configured to traverse the conductivity connection table, calculate the flow components between an injection well and a production well for all two adjacent grid cells divided by the control demarcation line, accumulate all the flow components, and calculate the total flow of the fluid between the injection well and the production well; a well-to-well connectivity coefficient calculation module, configured to obtain an expression of the well-to-well connectivity coefficient from the total flow calculation formula, and calculate the well-to-well connectivity coefficient according to the expression of the well-to-well connectivity coefficient.

[0015] One or more technical solutions provided in the embodiments of the present application have at least the following technical effects: The embodiment of the present application provides a method for quantifying and analyzing inter-well connectivity. By performing inter-layer division on the original reservoir model, the distribution information of injection wells and production wells in the lateral layer system can be obtained more accurately. Using the geometric segmentation method, the lateral layer system of the reservoir is divided into multiple interconnected grid units to ensure that injection wells and production wells are both located within the grid units. By dividing the control area for each injection well and production well, the distribution of fluid flow in the reservoir can be understood more intuitively. In determining the attribution of grid units, a method of comparing the Euclidean distances between the center points of each grid unit and all injection wells and production wells is adopted. This method is both simple and effective, and can accurately classify the grid units into the control areas of the corresponding injection wells or production wells. The present application also obtains the conductivity between two adjacent grid units divided by the control dividing line and records it using a conductivity connection table. This step is crucial for subsequent calculation of the flow components between injection wells and production wells. In constructing the mass conservation equation, the present application adopts the finite volume method and iteratively solves the pressure vector and total fluid saturation until the convergence condition is met. This method can more accurately simulate the dynamic change process of fluids in the reservoir. By traversing the conductivity connection table, calculating the flow components between two adjacent grid units divided by the control dividing line between injection wells and production wells, and accumulating to obtain the total flow, the expression of the inter-well connectivity coefficient is obtained, and the inter-well connectivity coefficient is calculated accordingly. This result is of great significance for evaluating the inter-well connectivity and optimizing the production strategy. It solves the problem of how to develop a more comprehensive and accurate method to analyze the inter-well connectivity to optimize the injection-production well pattern design and improve the effect of water flooding development and the reservoir recovery factor. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings required for describing the embodiments of the present application or the prior art. Obviously, the drawings in the following description are some embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0017] Figure 1 It is a flowchart of the method for quantifying and analyzing inter-well connectivity provided by the embodiment of the present application; Figure 2 It is a schematic diagram of the original reservoir model provided by the embodiment of the present application; Figure 3 It is a schematic diagram of the divided control area provided by the embodiment of the present application; Figure 4 It is a schematic diagram of the conductivity between two adjacent grid units divided by the control dividing line provided by the embodiment of the present application; Figure 5The calculation data diagram of the inter-well connectivity coefficient provided by the embodiment of the present application; Figure 6 The comparison result diagram of reservoir productivity prediction between the method based on numerical simulation and the method based on material balance principle provided by the embodiment of the present application; Figure 7 The schematic diagram of the device for quantifying and analyzing inter-well connectivity provided by the embodiment of the present application; Figure 8 The schematic diagram of the server for quantifying and analyzing inter-well connectivity provided by the embodiment of the present application. Detailed implementation manners

[0018] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present application.

[0019] The following explanations are made for some technologies involved in the embodiments of the present application to facilitate understanding. It should be considered that they are only exemplary. Therefore, those of ordinary skill in the art should recognize that various changes and modifications can be made to the embodiments described here without departing from the scope and spirit of the present application. Similarly, for the sake of clarity and conciseness, the descriptions of some well-known functions and structures are omitted in the following descriptions.

[0020] The embodiment of the present application provides a method for quantifying and analyzing inter-well connectivity. As Figure 1 shown, the method includes steps S101 to S109. Among them, Figure 1 This is only an execution order shown in the embodiment of the present application and does not represent the only execution order of a method for quantifying and analyzing inter-well connectivity. Under the condition that the final result can be achieved, Figure 1 The steps shown can be executed in parallel or reversed.

[0021] S101: Perform interlayer division on the original reservoir model to obtain the distribution information of injection wells and production wells in the horizontal layer series.

[0022] Specifically, this application combines the layered information of the original reservoir model with the actual perforation conditions of injection and production wells. By carefully classifying the perforation grids of injection and production wells within each layer, a unique injection-production well network system for each layer is constructed. Specifically, this application first accurately extracts the structural characteristics of each layer based on the layered data of the original reservoir model, and simultaneously collects the detailed perforation data of injection and production wells. Subsequently, each geological layer (zone) is traversed to accurately locate the perforation positions of each well in each layer, and then the perforation grid distribution of these wells in different layers is determined. If there is a perforation grid of an injection-production well in a certain layer, it means that the well has injection-production functions in that layer, so it is officially included in the composition of the injection-production well network of that layer. This process ensures that the injection-production well network of each layer can be accurately depicted, realizing the all-round division of the injection-production well network from the top layer to the bottom layer and from the vertical direction to the horizontal direction of the reservoir. It should be noted that injection-production wells refer to injection wells and production wells, and the injection-production well network refers to multiple pairs of injection wells and production wells.

[0023] Figure 2 It is a schematic diagram of the original reservoir model provided by the embodiment of this application. Figure 2 In (a), after the interlayer division of the original reservoir model, the distribution of a pair of injection well and production well in a certain longitudinal layer system at the vertical well position is shown. Figure 2 In (b), after the interlayer division of the original reservoir model, the distribution of the injection well and production well in a certain horizontal layer system at the vertical well position is obtained. Figure 2 It shows the formation structure in the reservoir and the spatial distribution of geological bodies, and can intuitively reflect the overall situation of reservoir geology. Figure 2 In it, the injection well is marked in blue, the production well is marked in red, and the arrow indicates the direction of fluid flow, vividly depicting the fluid exchange process between injection and production wells. The distribution of the injection well and production well in a certain horizontal layer system at the vertical well position can refine the specific distribution of the injection-production well network in different layers.

[0024] S102: Use the geometric segmentation method to process the distribution information of injection wells and production wells in the horizontal layer system, divide the reservoir into multiple interconnected grid units, and make the injection wells and production wells all located within the grid units.

[0025] S103: Use the division rules to divide the control area for each injection well and production well. Among them, the control area of a well is a polygonal area surrounded by the intersection points of the control dividing lines of six pairs of surrounding injection wells or production wells.

[0026] The division rules include: for any pair of injection well and production well, calculate the coordinates of their midpoint and the slope of the connection line, and calculate the slope of the perpendicular bisector. Obtain the equation of the perpendicular bisector according to the slope of the perpendicular bisector, and use the perpendicular bisector as the control dividing line.

[0027] The calculation formula for the coordinates of the midpoint is: . Among them, is the coordinate of the midpoint , is the coordinate of the oil production well , is the coordinate of the water injection well I.

[0028] The calculation formula for the slope of the connecting line is: . Among them, is the slope of the connecting line.

[0029] The calculation formula for the slope of the perpendicular bisector is: . Among them, is the slope of the perpendicular bisector.

[0030] The equation of the perpendicular bisector is: . Among them, and are variables in the equation of the perpendicular bisector, is the abscissa of any point in the horizontal layer system, is the ordinate of any point in the horizontal layer system.

[0031] Furthermore, the division rule of the present application not only considers the distribution information of water injection wells and oil production wells in the horizontal layer system, but also deeply combines the actual geological conditions of the oil reservoir to ensure the accuracy and practicability of the division.

[0032] Specifically, the obtained distribution information is analyzed in detail, and then according to an innovative space segmentation strategy, the oil reservoir area in the horizontal layer system is carefully divided into multiple interconnected control areas. The division of these control areas is not arbitrary, but strictly follows the geometric rules mentioned above. For each pair of water injection wells and oil production wells, the midpoint coordinates and the slope of the connecting line are calculated, and then the slope of the perpendicular bisector is obtained. Using these geometric parameters, the equation of the perpendicular bisector is derived, and the perpendicular bisector is used as the control boundary. It should be noted that when dividing the control areas in the present application, special attention is paid to the adaptability of the shape of each control area to the actual geological conditions of the oil reservoir. This means that the boundaries of the control areas are not simple straight lines or regular polygons, but are finely adjusted according to the geological characteristics of the oil reservoir, the well location distribution, and the fluid flow characteristics. Such a design ensures that each control area can accurately reflect the complex situation inside the oil reservoir, and at the same time ensures that the connection relationship between the water injection wells and the oil production wells can be captured and calculated.

[0033] Figure 3 This is a schematic diagram of the divided control areas provided by the embodiment of the present application. Figure 3 In (a) shows multiple divided control areas, Figure 3In (b), the control area of a specific water injection well and production well is magnified. Figure 3 The square grid in (b) in it is multiple interconnected grid cells into which the reservoir is divided in S102. Figure 3 Well I in it is a water injection well, well J is a production well, and the other wells are marked in blue and red respectively, with blue marking water injection wells and red marking production wells. After determining the control areas of each well, the connectivity relationship between the water injection well and the production well is established using the pressure field distribution information in the reservoir. Specifically, after injecting fluid into the water injection well and the production well, the diffusion path of the fluid in the reservoir and the pumping situation of the production well are analyzed. Figure 3 In (b) in it, the pressure field distribution between the injection and production wells is shown, where the arrows clearly indicate the direction of fluid displacement. It should be noted that the fluid in this application is the aqueous phase and the oil phase.

[0034] S104: Traverse all grid cells, and by comparing the Euclidean distance between the center point of each grid cell and all water injection wells and production wells, determine whether the grid cell belongs to the control area of a water injection well or a production well.

[0035] Traverse all grid cells, and by comparing the Euclidean distance between the center point of each grid cell and all water injection wells and production wells, determine whether the grid cell belongs to the control area of a water injection well or a production well, including: If the Euclidean distance from the center point of the grid cell to a certain water injection well is less than its Euclidean distance to all production wells, it is determined that the grid cell belongs to the control area of this water injection well. If the Euclidean distance from the center point of the grid cell to a certain production well is less than its Euclidean distance to all water injection wells, it is determined that the grid cell belongs to the control area of this production well.

[0036] The formula for calculating the Euclidean distance from the center point of the grid cell to the production well is: . Among them, is the Euclidean distance from the center point of the grid cell to the production well, is the coordinate of the center point of the grid cell.

[0037] The formula for calculating the Euclidean distance from the center point of the grid cell to the water injection well is: . Among them, is the Euclidean distance from the center point of the grid cell to the water injection well.

[0038] S105: Obtain the conductivity between two adjacent grid cells divided by the control dividing line, and record the conductivity between all two adjacent grid cells divided by the control dividing line using the conductivity connection table.

[0039] Obtain the conductivity between two adjacent grid cells divided by the controlled dividing line, including: The calculation formula for the conductivity between two adjacent grid cells divided by the controlled dividing line is: , , . Among them, is the conductivity between two adjacent grid cells divided by the controlled dividing line, and are intermediate variables, is the permeability of grid cell , is the permeability of grid cell , is the area of grid cell , is the area of grid cell , is the distance from the center point of grid cell to the controlled dividing line, is the distance from the center point of grid cell to the controlled dividing line, is the unit normal vector of the controlled dividing line pointing to grid cell , is the unit normal vector of the controlled dividing line pointing to grid cell , is the unit direction vector of the center point of the controlled dividing line pointing to grid cell , is the unit direction vector of the center point of the controlled dividing line pointing to grid cell .

[0040] It should be noted that the conductivity connection table can be presented in the form of a matrix, and each element in the matrix represents the conductivity between two adjacent grid cells divided by the controlled dividing line.

[0041] Figure 4 This is a schematic diagram of the conductivity between two adjacent grid cells divided by the controlled dividing line provided by the embodiment of the present application. Figure 4 The interpretation of the relevant parameters in

[0042] S106: Use the production index to represent the production capacity of the oil production well under the first pressure difference, and use the injection index to represent the injection capacity of the injection well under the second pressure difference.

[0043] Use the production index to represent the production capacity of the oil production well under the first pressure difference, and use the injection index to represent the injection capacity of the injection well under the second pressure difference, including: The first pressure difference is: . Among them, is the first pressure difference, is the reservoir pressure, is the bottom-hole pressure. The second pressure difference is: . Wherein, is the second pressure difference.

[0044] The calculation formula of the productivity index or injection index is: . Wherein, is the productivity index or the injection index , is the unit conversion factor, is the effective angle around the wellbore, is the reservoir capacity, is the pressure equivalent radius of the grid cell, is the actual radius of the wellbore, is the skin effect factor, is the logarithmic term, representing the flow resistance from the grid cell to the wellbore.

[0045] Specifically, in reservoir development and fluid management, accurately quantifying the production capacity of production wells and the injection capacity of injection wells is crucial. For this purpose, this application introduces the productivity index and the injection index these two key indicators. These two indices not only reflect the fluid handling capacity of the wellbore under a specific pressure difference, but also contain the influence of reservoir geological characteristics and wellbore conditions on fluid flow. The first pressure difference is the main driving force for driving reservoir fluids to flow into the wellbore. The second pressure difference is the key driving force for pushing fluids from the wellbore into the reservoir. It should be noted that the calculation formulas of the productivity index and the injection index are the same.

[0046] S107: For all adjacent two grid cells divided by the controlled dividing line, construct a mass conservation equation according to the finite volume method, and iteratively solve the pressure vector and the total fluid saturation until the convergence condition is met.

[0047] The mass conservation equation is: . Wherein, is the fluid residual, is the grid cell porosity, is the grid cell density of the fluid in, is the grid cell fluid in saturation, is the time step change amount within, is the grid cell volume, For the mass flow rate flowing from grid cell into grid cell , is the set of grid cells adjacent to grid cell , and

[0048] is the source / sink term. . Among them, is the relative permeability of fluid in grid cell , is the viscosity of fluid in grid cell , is the productivity index .

[0049] When the grid cell is a production well, the source / sink term is: , is the injection index .

[0050] The expression for the mass flow rate flowing from grid cell into the adjacent grid cell is: . Among them, is the conductivity between two adjacent grid cells divided by the control boundary, is the average value of the mobility parameter between grid cell and grid cell , is the density of fluid , is the relative permeability of fluid , is the viscosity of fluid , represents the gravity term, is the gravitational acceleration, is the height difference between adjacent grid cells, is the pressure difference between adjacent grid cells.

[0051] . Among them, is the fluid compressibility, is the fluid density at the corresponding reference pressure , is the pressure of the current grid cell, is the reference pressure, is the base of the natural logarithm.

[0052] . Among them, is the fluid mobility parameter in the grid cell , is the fluid mobility parameter in the grid cell .

[0053] The pressure vector includes the pressures of all grid cells, and the total fluid saturation includes the saturations of the fluids in all grid cells.

[0054] Furthermore, the pressure vector and the total fluid saturation are iteratively solved until the convergence condition is satisfied, including the following steps.

[0055] Before the iterative solution, a series of initial conditions need to be set, including the initial pressure field , the saturation of the initial fluid , and the time step and the convergence threshold . These conditions provide a starting point for the iterative process and a criterion for judging convergence.

[0056] According to the current pressure field (k represents the number of iterations), calculate the residual of the mass conservation equation. The residual includes the oil-phase residual and the water-phase residual , corresponding to each grid cell respectively. Within the time step , assuming that variables such as the density, viscosity, and grid porosity of the fluid remain unchanged, take the partial derivatives of the residual with respect to the variables to be solved, namely the pressure of the grid cell and the saturation of the fluid in the grid cell. These partial derivatives will be used to construct the Jacobian matrix .

[0057] Specifically, the calculation process of taking the partial derivatives with respect to the variables to be solved, namely the pressure of the grid cell and the saturation of the fluid in the grid cell, is as follows.

[0058] The partial derivative of the oil-phase residual with respect to the oil-phase saturation is: . It should be noted that the subscripts all represent the relevant parameters of the oil phase in the grid cell . is the density of the oil phase in the grid cell , is the viscosity of the oil phase in the grid cell , is the relative permeability of the oil phase . is the symbol of partial derivative, is the relative permeability of the oil phase in the grid cell and the grid cell ; is the density of the oil phase, is the viscosity of the oil phase.

[0059] The oil phase residual with respect to the pressure of the grid cell has the partial derivative as follows: .

[0060] The oil phase residual with respect to the pressure of the adjacent grid cell has the partial derivative as follows: .

[0061] The water phase residual with respect to the water phase saturation has the partial derivative as follows: . It should be noted that the subscript all represents the relevant parameters of the water phase in the grid cell ; is the density of the water phase in the grid cell ; is the viscosity of the water phase in the grid cell ; is the relative permeability of the water phase ; is the relative permeability of the water phase in the grid cell and the grid cell ; is the density of the water phase, is the viscosity of the water phase.

[0062] The water phase residual with respect to the pressure of the grid cell has the partial derivative as follows: .

[0063] The water phase residual with respect to the pressure of the adjacent grid cell has the partial derivative as follows: .

[0064] Furthermore, using the previously calculated partial derivatives, the Jacobian matrix 。The Jacobian matrix is a square matrix whose number of rows and columns is equal to the number of variables to be solved. In this application, it includes the pressure of each grid cell and the saturation of the fluid in the grid cell. Constructing the Jacobian matrix Each term of is: Among them, is the element corresponding to the th row and the th column in the Jacobian matrix, is the th element of the residual vector , is the th element of the variable to be solved . , is the pressure vector, is the total fluid saturation. , is the pressure of different grid cells, is the total number of grid cells, is the transpose operation, , is the saturation of the fluid in different grid cells. , , is the oil-phase residual vector, is the oil-phase residual of different grid cells, , is the water-phase residual vector, is the water-phase residual of different grid cells.

[0065] Solve the linear equations: Among them, . is the change in the variable to be solved, which refers to the difference between the th iteration and the th iteration, is the variable to be solved at the th iteration, is the variable to be solved at the th iteration.

[0066] Set a small positive number as the convergence threshold. When , the convergence condition is met, and the solution is completed. Record the pressure of each grid cell, the saturation of the fluid in the grid cell, and the flow rate between any adjacent grid cells. If the convergence condition is not met, continue the iteration.

[0067] S108: Traverse the conductivity connection table. For all adjacent two grid cells separated by the controlled dividing line, calculate the flow component between the injection well and the production well, accumulate all the flow components, and calculate the total flow of the fluid between the injection well and the production well.

[0068] The calculation formula for the flow component is: . Wherein, is the flow component between the adjacent two grid cells between the injection well and the production well.

[0069] The calculation formula for the total flow is: . Wherein, is the total flow of the fluid between the injection well and the production well, is the control area of the injection well , is the control area of the production well , is the inter-well connection coefficient, is the average pressure within the control area of the injection well , is the average pressure within the control area of the production well .

[0070] S109: Derive the expression of the inter-well connection coefficient from the calculation formula of the total flow, and calculate the inter-well connection coefficient according to the expression of the inter-well connection coefficient.

[0071] The expression of the inter-well connection coefficient is: ; wherein, , , is the pressure of the grid cell , is the pressure of the grid cell , is the volume of the grid cell , is the volume of the grid cell , is the total volume of the control area of the injection well , is the total volume of the control area of the production well .

[0072] Figure 5 This is the calculation data graph of the inter-well connection coefficient provided by the embodiment of the present application. As Figure 5 shown, the abscissa represents the well name, and the ordinate is the inter-well connection coefficient.

[0073] Figure 6 This is the comparison result graph of the reservoir production capacity prediction between the method based on numerical simulation and the method based on the material balance principle provided by the embodiment of the present application. AsFigure 6 As shown, the abscissa represents the prediction time, and the ordinate represents the daily crude oil production of the production well. The red broken line in the figure represents the prediction result of numerical simulation, and the green broken line represents the prediction result of material balance. The method used in this application is based on the principle of material balance. From Figure 6 it can be seen that for Well W15-17, the broken line graphs of the method based on numerical simulation and the method based on the principle of material balance in this application basically coincide, and the accuracy rate of predicted production is maintained at 80%, which proves the reliability of the well connectivity quantification method in this application.

[0074] The embodiment of this application also provides a well connectivity quantification and analysis device 700, as Figure 7 shown. The device includes: an acquisition distribution information module 701, a grid cell division module 702, a control area division module 703, a determination module 704, a conductivity acquisition module 705, a representation module 706, an iterative solution module 707, a total flow calculation module 708, and a well connectivity coefficient calculation module 709.

[0075] The acquisition distribution information module 701 is used to perform interlayer division on the original reservoir model to obtain the distribution information of injection wells and production wells in the horizontal layer system.

[0076] The grid cell division module 702 is used to process the distribution information of injection wells and production wells in the horizontal layer system by using the geometric segmentation method, divide the reservoir into multiple interconnected grid cells, and make the injection wells and production wells located within the grid cells.

[0077] The control area division module 703 is used to use the division rules to divide the control area for each injection well and production well. Among them, the control area of a well is a polygonal area surrounded by the intersection points of the control dividing lines of six pairs of injection wells or production wells around it.

[0078] The determination module 704 is used to traverse all grid cells, and determine whether the grid cell belongs to the control area of the injection well or the production well by comparing the Euclidean distance between the center point of each grid cell and all injection wells and production wells.

[0079] The conductivity acquisition module 705 is used to acquire the conductivity between two adjacent grid cells divided by the control dividing line, and record the conductivity between all two adjacent grid cells divided by the control dividing line using a conductivity connection table.

[0080] The representation module 706 is used to represent the production capacity of the production well under the first pressure difference using the production index, and represent the injection capacity of the injection well under the second pressure difference using the injection index.

[0081] The iterative solution module 707 is used to construct a mass conservation equation for two adjacent grid cells divided by all controlled dividing lines according to the finite volume method, and iteratively solve the pressure vector and the total fluid saturation until the convergence condition is met.

[0082] The total flow rate calculation module 708 is used to traverse the conductivity connection table, calculate the flow rate components between the injection well and the production well for all two adjacent grid cells divided by the controlled dividing lines, accumulate all the flow rate components, and calculate the total flow rate of the fluid between the injection well and the production well.

[0083] The inter-well connectivity coefficient calculation module 709 is used to obtain the expression of the inter-well connectivity coefficient from the calculation formula of the total flow rate, and calculate the inter-well connectivity coefficient according to the expression of the inter-well connectivity coefficient.

[0084] Some modules in the device described in this application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform specific tasks or implement specific abstract data types. This application can also be practiced in a distributed computing environment, where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media including storage devices.

[0085] The device or module illustrated in the above application embodiments can be specifically implemented by a computer chip or entity, or by a product with a certain function. For the convenience of description, the above device is described by dividing it into various modules according to functions. When implementing the embodiments of this application, the functions of each module can be implemented in the same or multiple software and / or hardware. Of course, the module that implements a certain function can also be implemented by combining multiple sub-modules or sub-units.

[0086] The methods, devices or modules described in this application can be implemented in the form of computer-readable program codes. The controller can be implemented in any appropriate manner. For example, the controller can take the form of, for example, a microprocessor or a processor, a computer-readable medium storing computer-readable program codes (such as software or firmware) executable by the (micro)processor, logic gates, switches, application specific integrated circuits (ASICs), programmable logic controllers, and embedded microcontrollers. Examples of the controller include, but are not limited to, the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320. The memory controller can also be implemented as part of the control logic of the memory. Those skilled in the art also know that, in addition to implementing the controller in the form of pure computer-readable program codes, the method steps can be logically programmed to enable the controller to be implemented in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, and embedded microcontrollers to achieve the same function. Therefore, such a controller can be considered a hardware component, and the devices included therein for implementing various functions can also be regarded as the structures within the hardware component. Or even, the devices for implementing various functions can be regarded as either software modules for implementing the method or structures within the hardware component.

[0087] As Figure 8 shown, an inter-well connectivity quantification and analysis server is further provided in an embodiment of this application, including a memory 801 and a processor 802; the memory 801 is used to store computer-executable instructions; the processor 802 is used to execute the computer-executable instructions to implement an inter-well connectivity quantification and analysis method described above in the embodiment of this application.

[0088] An embodiment of this application further provides a computer-readable storage medium, which stores executable instructions. When a computer executes the executable instructions, it can implement an inter-well connectivity quantification and analysis method described above in the embodiment of this application.

[0089] As can be seen from the description of the above embodiments, those skilled in the art can clearly understand that the present application can be implemented by means of software plus necessary hardware. Based on such an understanding, the technical solution of the present application, in essence or the part that contributes to the prior art, can be embodied in the form of a software product, or can also be embodied in the implementation process of data migration. The computer software product can be stored in a storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions for causing a computer device (which can be a personal computer, a mobile terminal, a server, or a network device, etc.) to execute the method described in the embodiments of the present application.

[0090] The various embodiments in this specification are described in a progressive manner. For the same or similar parts between the various embodiments, reference can be made to each other. Each embodiment focuses on the differences from other embodiments. All or part of the present application can be used in many general or special computer system environments or configurations.

[0091] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit the present application. Although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the present application.

Claims

1. A method for quantifying and analyzing well connectivity, characterized in that: include: Divide the original reservoir model into layers to obtain the distribution information of water injection wells and oil production wells in the horizontal layer system; The geometric segmentation method is used to process the distribution information of water injection wells and oil production wells in the horizontal layer system, and the oil reservoir is divided into multiple interconnected grid units, so that the water injection wells and oil production wells are all located in the grid units; Use the division rule to divide each water injection well and oil production well into its control area; the control area of ​​a well is a polygonal area surrounded by the intersection points of the control boundary lines of the six pairs of water injection wells or oil production wells around it; Traverse all grid cells, and determine whether the grid cell belongs to the control area of ​​the water injection well or the control area of ​​the oil production well by comparing the Euclidean distance between the center point of each grid cell and all water injection wells and oil production wells; Obtain the conductivity between two adjacent grid cells separated by the control boundary line, and use the conductivity connectivity table to record the conductivity between all two adjacent grid cells separated by the control boundary line; The production index is used to represent the production capacity of the oil production well under the first pressure difference, and the water injection index is used to represent the water injection capacity of the water injection well under the second pressure difference; For all two adjacent grid cells separated by the control boundary, the mass conservation equation is constructed according to the finite volume method, and the pressure vector and total fluid saturation are iteratively solved until the convergence condition is met; Traverse the conductivity connectivity table, calculate the flow components between the water injection well and the oil production well for all two adjacent grid cells separated by the control boundary line, accumulate all flow components, and calculate the total flow of the fluid between the water injection well and the oil production well; The expression of the inter-well connectivity coefficient is obtained from the calculation formula of the total flow, and the inter-well connectivity coefficient is calculated based on the expression of the inter-well connectivity coefficient.

2. The well connectivity quantification and analysis method according to claim 1, characterized in that: The division rules include: For any pair of injection wells and production wells, calculate the coordinates of their midpoints and the slope of the connecting line, and calculate the slope of the perpendicular bisector; Obtain the equation of the perpendicular bisector according to its slope, and use the perpendicular bisector as the control dividing line; The coordinates of the midpoint are calculated as: ;in, Midpoint The coordinates of For oil wells The coordinates of For water injection well The coordinates of The slope of the connecting line is calculated as: ;in, is the slope of the connecting line; The formula for calculating the slope of the perpendicular bisector is: ;in, is the slope of the perpendicular bisector; The equation of the perpendicular bisector is: ;in, and are the variables in the equation for the perpendicular bisector, is the horizontal coordinate of any point in the horizontal layer system, is the vertical coordinate of any point in the horizontal layer system.

3. The method for quantifying and analyzing well connectivity according to claim 2, characterized in that: The traversing of all grid cells and determining whether the grid cell belongs to the control area of ​​the water injection well or the control area of ​​the oil production well by comparing the Euclidean distance between the center point of each grid cell and all water injection wells and oil production wells includes: If the Euclidean distance from the center point of a grid cell to a water injection well is less than the Euclidean distance from the center point to all oil production wells, the grid cell is determined to belong to the control area of ​​the water injection well; If the Euclidean distance from the center point of a grid cell to an oil well is less than the Euclidean distance from the center point to all water injection wells, the grid cell is determined to belong to the control area of ​​the oil well.

4. The method for quantifying and analyzing well connectivity according to claim 1, characterized in that: The step of obtaining the conductivity between two adjacent grid cells separated by the control boundary line comprises: The calculation formula of the conductivity between two adjacent grid cells separated by the control boundary is: , , ;in, is the conductivity between two adjacent grid cells separated by the control boundary, and is the intermediate variable, For grid cells The permeability, For grid cells The permeability, For grid cells The area of For grid cells The area of For grid cells The distance from the center point to the control boundary line, For grid cells The distance from the center point to the control boundary line, To control the direction of the dividing line to the grid cell The unit normal vector of To control the direction of the dividing line to the grid cell The unit normal vector of To control the center point of the dividing line to point to the grid unit The unit direction vector of To control the center point of the dividing line to point to the grid unit The unit direction vector of .

5. The method for quantifying and analyzing well connectivity according to claim 4, characterized in that: The method of using the production index to represent the production capacity of the oil production well under the first pressure difference and using the water injection index to represent the water injection capacity of the water injection well under the second pressure difference includes: The first pressure difference is: ;in, is the first pressure difference, is the reservoir pressure, is the bottom hole pressure; The second pressure difference is: ;in, is the second pressure difference; The calculation formula of production index or water injection index is: ;in, Production Index or water injection index , is the unit conversion factor, is the effective angle around the wellbore, is the reservoir capacity, is the pressure equivalent radius of the grid cell, is the actual radius of the wellbore, is the epidermal effect factor, It is a logarithmic term, which represents the flow resistance from the grid cell to the wellbore.

6. The method for quantifying and analyzing well connectivity according to claim 5, characterized in that: The method of constructing a mass conservation equation for all two adjacent grid cells separated by the control boundary line according to the finite volume method, and iteratively solving the pressure vector and the total fluid saturation until the convergence condition is met includes: The mass conservation equation is: ;in, is the fluid residual, For grid cells The porosity, For grid cells The density of the fluid, For grid cells Medium fluid The saturation of is the time step The amount of change within For grid cells The volume of From the grid unit Flow into grid cell The mass flow rate, For the grid unit The set of adjacent grid cells, is the source-sink term; When the grid cell contains an oil well, the source and sink items are: ;in, For fluid In the grid cell The relative permeability, For fluid In the grid cell The viscosity of Production Index ; When there is an injection well in the grid cell, the source and sink terms are: , Water injection index ; From the grid cell Flowing into adjacent grid cells Mass flow rate The expression is: ;in, is the conductivity between two adjacent grid cells separated by the control boundary, For the grid cells and grid cells The average value of the liquidity parameter between For fluid The density of For fluid The relative permeability, For fluid The viscosity of represents the gravity term, is the acceleration due to gravity, is the height difference between adjacent grid cells, is the pressure difference between adjacent grid cells; ;in, is the fluid compressibility coefficient, Corresponding to the reference pressure The fluid density under is the pressure of the current grid cell, is the reference pressure, is the base of natural logarithm; ;in, For the grid cells The fluidity parameters of the fluid, For the grid cells The fluidity parameters of the medium; The pressure vector includes the pressure of all grid cells, and the total fluid saturation includes the saturation of the fluid in all grid cells.

7. The method for quantifying and analyzing well connectivity according to claim 6, characterized in that: The traversal conductivity connectivity table calculates the flow components between the water injection well and the oil production well for all two adjacent grid cells separated by the control boundary line, accumulates all flow components, and calculates the total flow of the fluid between the water injection well and the oil production well, including: The flow rate component is calculated as: ;in, is the flow component between the water injection well and the oil production well in two adjacent grid cells; The total flow rate is calculated as: ;in, is the total flow rate of fluid between the injection well and the production well, For water injection well control area, For oil wells control area, is the well connectivity coefficient, For water injection well The average pressure in the control area, For oil wells The average pressure in the control area.

8. The method for quantifying and analyzing well connectivity according to claim 7, characterized in that: The expression of the inter-well connectivity coefficient is obtained from the calculation formula of the total flow, and the inter-well connectivity coefficient is calculated according to the expression of the inter-well connectivity coefficient, including: The expression of the well connectivity coefficient is: ;in, , , For grid cells pressure, For grid cells pressure, For grid cells The volume of For grid cells The volume of For water injection well The total volume of the control area, For oil wells The total volume of the control area.

9. A device for quantifying and analyzing well connectivity, characterized in that: include: The module for obtaining distribution information is used to divide the original reservoir model into layers and obtain the distribution information of water injection wells and oil production wells in the horizontal layer system; The grid unit division module is used to process the distribution information of water injection wells and oil production wells in the horizontal layer system by using the geometric segmentation method, and divide the oil reservoir into multiple interconnected grid units so that the water injection wells and oil production wells are all located in the grid units; The control area division module is used to divide the control area for each water injection well and oil production well using the division rule; wherein the control area of ​​a well is a polygonal area surrounded by the intersection points of the control boundary lines of six pairs of water injection wells or oil production wells around it; A determination module is used to traverse all grid cells and determine whether the grid cell belongs to the control area of ​​the water injection well or the control area of ​​the oil production well by comparing the Euclidean distance between the center point of each grid cell and all water injection wells and oil production wells; A conductivity acquisition module is used to acquire the conductivity between two adjacent grid cells divided by the control boundary line, and use a conductivity connectivity table to record the conductivity between all two adjacent grid cells divided by the control boundary line; A representation module, used for representing the production capacity of the oil production well under the first pressure difference by using the production index, and representing the water injection capacity of the water injection well under the second pressure difference by using the water injection index; An iterative solution module is used to construct the mass conservation equation for all two adjacent grid cells separated by the control boundary line according to the finite volume method, and iteratively solve the pressure vector and the total fluid saturation until the convergence condition is met; The total flow calculation module is used to traverse the conductivity connectivity table, calculate the flow components between the water injection well and the oil production well for all two adjacent grid cells separated by the control boundary line, accumulate all flow components, and calculate the total flow of the fluid between the water injection well and the oil production well; The module for calculating the inter-well connectivity coefficient is used to obtain the expression of the inter-well connectivity coefficient from the calculation formula of the total flow rate, and calculate the inter-well connectivity coefficient according to the expression of the inter-well connectivity coefficient.

Images