Calculation method and system for effective communication of injection-production system under intelligent identification constraint

Abstract

The invention relates to the field of oil reservoir engineering communication technology research, and discloses a calculation method and system for effective communication of an injection-production system under intelligent identification constraint, and the method comprises the steps: obtaining oil field injection-production comprehensive data, carrying out the standardized preprocessing of the obtained oil field injection-production comprehensive data, achieving the index statistics of the oil field injection-production comprehensive data, and achieving the effective communication of the oil field injection-production system. The injection-production relation between the injection-production well and the injection-production well pattern on each sand body plane is calculated through a ground effective injection-production well pattern and an underground effective injection-production well pattern in sequence according to the oilfield injection-production comprehensive data subjected to standardization pretreatment, and a one-element-one-convergence-element injection-production system mathematical model is established; according to the injection-production relationship between an injection-production well and an injection-production well pattern on each sand body plane, the attributes of an element injection-production system are determined under the constraint of an identification structure and lithology, the effective communication relationship and effective action time of the injection-production system are intelligently identified, and then parameters such as a water drive theoretical sweep area and a water drive reserve utilization area are calculated. The per capita working efficiency can be greatly improved, and the labor time-consuming cost is reduced.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil reservoir engineering communication, in particular to a calculation method and system for intelligently identifying effective connectivity of injection-production systems under constraints. BACKGROUND

[0002] The water drive control degree refers to the ratio of the reserves in the oil-bearing area that can be affected by the injected water under the existing well pattern to the total producing geological reserves of the reservoir, and the water drive control degree is a reflection of the water injection volume sweep efficiency, which is affected by factors such as structure, sedimentation, sand body distribution, and artificial well pattern. At present, the most commonly used and most popular simplified "thickness method" defined by the industry standard is used for calculation, that is, the ratio of the effective thickness of the production well that is connected with the injection well to the total effective thickness of the production well in the well group. This method actually only considers the vertical injection-production correspondence and does not reflect the concept of reserves.

[0003] The water drive producing degree is defined as the ratio of the water drive producing reserves to the geological reserves of the reservoir, and the simplified standard algorithm in the industry is to calculate the ratio of the total water absorption thickness to the total connected thickness of the injection well, or the ratio of the total liquid production thickness to the total connected thickness of the oil well, according to the water absorption profile of all the tested water wells and the liquid production profile of all the tested oil wells in a year. Similarly, the concept of reserves is not reflected.

[0004] As described above, the statistical methods of water drive control degree and producing degree are both simplified "thickness methods", but even so, for large work areas or thin interbedded reservoirs with long well sections, the development unit and the number of oil and water wells are large, resulting in a large amount of manual statistical work and a long time-consuming period. In addition, this algorithm ignores the constraints of the actual geological background (such as fault cutting, sand body pinch-out, and heterogeneity influence, etc.), causing misjudgment of injection-production connectivity to occur from time to time, so the calculation results often differ greatly from the true situation, and even the producing degree value is larger than the control degree value. SUMMARY

[0005] The present application aims to provide a calculation method and system for intelligently identifying effective connectivity of injection-production systems under constraints, to solve the technical problems of large manual repetitive calculation workload and large calculation structure error in the prior art.

[0006] The present application is achieved by the following technical solutions:

[0007] In a first aspect, the present application provides a calculation method for intelligently identifying effective connectivity of injection-production systems under constraints, comprising:

[0008] acquiring oilfield injection-production comprehensive data, and performing standardized preprocessing on the acquired oilfield injection-production comprehensive data;

[0009] The standardized pretreated oilfield injection-production comprehensive data is sequentially subjected to ground effective injection-production well pattern and subsurface effective injection-production well pattern to calculate injection-production relationship of injection-production well and injection-production well pattern on each sand body plane;

[0010] A unary one-hub injection-production system mathematical model is established, and injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is used to determine injection-production system attribute under the identification of structure and lithology constraint;

[0011] According to the injection-production system attribute, effective swept superimposed area and effective producing superimposed area in single sand body are obtained, and water drive control degree and water drive producing degree of single sand body are calculated according to area method, and single sand body swept area boundary coordinates are extracted according to water drive control degree and water drive producing degree of single sand body.

[0012] The water drive control degree and water drive producing degree of all single sand bodies contained in the reservoir and development unit are summed up, and the calculation of injection-production system effective connectivity under intelligent identification constraint is completed.

[0013] Preferably, the oilfield injection-production comprehensive data includes drilling and completion data, injection-production relationship data, sand body boundary data, structure boundary data, well trajectory data, stratigraphic correlation data, interpretation conclusion data, layer opening and closing well history data, and production and suction profile data.

[0014] Preferably, the standardization pretreatment process is as follows: unification of well name in English, unification of static and dynamic layer numbers, unification of injection-production effective connectivity determination standard, unification of sand body boundary and oil-bearing area boundary envelope, and unification of well trajectory geological engineering intersection point coordinate calculation method.

[0015] Preferably, the specific process of sequentially performing ground effective injection-production well pattern and subsurface effective injection-production well pattern on the standardized pretreated oilfield injection-production comprehensive data to calculate injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is as follows:

[0016] S1, the standardized pretreated oilfield injection-production comprehensive data is subjected to effective injection-production well pattern calculation to obtain ground injection-production relationship, and the ground injection-production relationship is stored in the db_surface data model;

[0017] S2, according to the ground injection-production relationship, the subsurface effective injection-production well pattern is calculated, the intermediate calculation results are expanded in three dimensions according to the wellbore, the envelope relationship with the outer boundary and the inner boundary of single sand body is determined through well trajectory coordinates layer by layer, the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is generated, the corresponding conclusion is output according to the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane, and the conclusion is stored in the db_under data model.

[0018] Preferably, a unit injection-production system mathematical model is established, and the process of calculating the unit injection-production system attribute under the constraints of structure and lithology according to the injection-production relationship of injection-production wells and injection-production well networks on each sand body plane is as follows

[0019] According to the seepage field theory, an orthogonal coordinate system is established, the coordinates of the unit injection-production well points and the coordinates of the sand body inner and outer boundaries are called, and the breakpoint coordinates and the sand body boundary coordinates in the same sand body plane are called, and the actual two-dimensional unit injection-production system plane graph affected by structure and lithology is projected on the orthogonal coordinate system;

[0020] The coordinates of the actual two-dimensional unit injection-production system plane graph are translated and rotated as a whole around the water well, so as to become an ellipse with the water well at the origin of the orthogonal coordinate system and located at the left focus point and the oil well at the right vertex, and the effective swept area under the unit injection-production system is calculated by calculating the area of the ellipse;

[0021] The determination condition threshold range is set to determine the unit injection-production system effective swept area identification mechanism model, the degree of cutting of the unit injection-production system swept shape by the fault line and the sand body pinchout line is judged to obtain the unit injection-production system attribute, if the plane graph combination relationship meets the condition threshold range, it is determined that the injection-production connectivity of the unit injection-production system attribute is objective and effective, and is valued as 1; otherwise, it is determined that the injection-production connectivity of the unit injection-production system attribute is invalid, and is valued as 0.

[0022] Further, the unit injection-production system effective swept area identification mechanism model includes an Ic type unit injection-production system model, an If type unit injection-production system model, an nc type unit injection-production system model, an mc type unit injection-production system model, an nf type unit injection-production system model and an mf type unit injection-production system model;

[0023] When the condition threshold range is that the container arrSand count and the container arrFault count are both zero, it is an Ic type unit injection-production system model;

[0024] When the condition threshold range is that the container arrSand count is zero and the container arrFault count is greater than zero, it is an If type unit injection-production system model;

[0025] When the container arrSand count is greater than zero and the container arrFault count is zero, it is judged whether there is an intersection between the sand body edge line and the X positive half axis, when there is an intersection, it is an nc type unit injection-production system model, otherwise it is an mc type unit injection-production system model;

[0026] When the container arrSand count and the container arrFault count are both greater than zero, it is judged whether there is an intersection between the sand body edge line and the X positive half axis, when there is an intersection, it is an nf type unit injection-production system model, otherwise it is an mf type unit injection-production system model.

[0027] Further, the unit injection-production system attribute includes address oiliness, fault blocking, effective connectivity and effective dynamic property;

[0028] Among them, if the oil well interpretation conclusion in the oil-bearing layer is oil layer, oil-water layer, oil-gas layer, poor oil layer, then it is geologically oil-bearing, and the value is 1, otherwise the value is 0;

[0029] In the fault blocking, if the intersection of the fault and the X-axis in the orthogonal coordinate system satisfies greater than 0 and less than the distance between the oil well and the water well, it is considered as fault blocking, and the value is 1, otherwise the value is 0;

[0030] In the effective connectivity, the perforation state of the oil well and the water well is counted per year, if the oil well and the water well are in the perforation production state within the counted year, the value is 1, otherwise the value is 0;

[0031] In the effective production, for medium-high permeability or complete liquid production data of the reservoir, if there is corresponding production or absorption, or corresponding effect due to micro test, the value is 1, otherwise the value is 0; for low permeability or missing liquid production data of the reservoir, the water absorption of the water well is used for judgment, the test shows that the water absorption or the water absorption is not tested but the KH ratio is greater than 0.25, the value is 1, otherwise the value is 0.

[0032] Preferably, the calculation process of the water drive control degree of a single sand body is as follows:

[0033] The oil saturation parameters, the crude oil density parameters and the volume coefficient parameters of each injection-production system in the same single sand body are considered to be consistent, and according to the volumetric reserve calculation formula, the single sand body water drive reserve control degree expression of the reserve definition method can be simplified as single sand body water drive control area / single sand body oil-bearing area;

[0034] The calculation process of the water drive production degree of a single sand body is as follows:

[0035] For quantitative calculation of single sand body water drive production reserves, based on quantitative calculation of single sand body water drive control reserves, the effective injection-production production area of the injection-production system marked in the previous module is called, and then the single sand body water drive reserve production degree expression can be simplified as single sand body water drive sweep area / single sand body oil-bearing area.

[0036] Preferably, the water drive control degree calculation formula of all single sand bodies contained in the reservoir and the development unit is as follows:

[0037]

[0038] The water drive production degree calculation formula of all single sand bodies contained in the reservoir and the development unit is as follows:

[0039]

[0040] Wherein, R is the number of oil-bearing single sand bodies of the development unit; N oi is the geological reserves of each sand body, t; N diis the water drive development degree of each single sand body, t; E d is the water drive development degree of the development unit, %; N Wi is the water drive control degree of the development unit; N i is the water drive swept reserve of each sand body.

[0041] In the second aspect, the application provides a computing system for intelligently identifying the effective connectivity of a constrained injection-production system, which is used to realize the computing method for intelligently identifying the effective connectivity of the constrained injection-production system, and comprises the following steps:

[0042] A data preprocessing module is configured to acquire oilfield injection-production comprehensive data and perform standardization preprocessing on the acquired oilfield injection-production comprehensive data.

[0043] A first data calculation module is configured to sequentially calculate the injection-production relationship of injection-production wells and injection-production well networks on the plane of each sand body by using the oilfield injection-production comprehensive data after standardization preprocessing and by using the ground effective injection-production well network and the underground effective injection-production well network.

[0044] An attribute determination module is configured to establish a one-element-one-sink-element injection-production system mathematical model, and determine the attribute of the element injection-production system under the recognition constraints of structure and lithology according to the injection-production relationship of injection-production wells and injection-production well networks on the plane of each sand body.

[0045] A second data calculation module is configured to obtain the effective swept superimposed area and the effective development superimposed area in a single sand body according to the attribute of the element injection-production system, calculate the water drive control degree and the water drive development degree of the single sand body by using the area method, and extract the swept area boundary coordinates of the single sand body according to the water drive control degree and the water drive development degree of the single sand body.

[0046] A third data calculation module is configured to sum up and aggregate the water drive control degree and the water drive development degree of all single sand bodies contained in the reservoir and the development unit, and complete the calculation of the effective connectivity of the constrained injection-production system.

[0047] Compared with the prior art, the application has the following beneficial technical effects:

[0048] The application provides a calculation method for intelligently identifying effective connectivity of injection-production systems, and the method comprises the following steps: obtaining oilfield injection-production comprehensive data, performing standardization preprocessing on the obtained oilfield injection-production comprehensive data, realizing index statistics of the oilfield injection-production comprehensive data, and sequentially calculating injection-production relationships of injection-production wells and injection-production well networks on each sand body plane by using the oilfield injection-production comprehensive data after standardization preprocessing, establishing a one-in-one-out injection-production system mathematical model, determining injection-production system properties under the constraints of structure and lithology according to the injection-production relationships of injection-production wells and injection-production well networks on each sand body plane, intelligently identifying effective connectivity and effective action time of the injection-production system, and further calculating parameters such as water drive theoretical swept area and water drive reserve producing area, and finally obtaining water drive control degree and producing degree of single sand body, different development units and the whole reservoir, which can greatly improve the work efficiency per capita, reduce the labor time cost, increase the cost benefit of enterprises, provide more reliable and real water drive reserve control degree and producing degree for reservoir development decision makers, and improve the calculation accuracy of water drive reserve control degree and producing degree.

[0049] Further, the effective swept superimposed area and the effective producing superimposed area in the single sand body are obtained according to the injection-production system properties, the water drive control degree and the water drive producing degree of the single sand body are obtained by using the area method, the water drive control degree and the water drive producing degree of the single sand body are extracted to obtain the swept area boundary coordinates, the automation of the water drive reserve control and producing reserve calculation is realized, the constraint conditions for determining the calculation accuracy of the reserves are limited to the actual underground injection-production well spacing, the pre-adjusted injection-production well spacing, the single sand body drawing precision and the structure interpretation precision, the method is more inclined to the guidance of geological results compared with the traditional method, and the scientific research and analysis efficiency is improved.

[0050] Further, the boundary filtering, coordinate translation and rotation, injection-production system effective area parameter extraction, injection-production system type identification, injection-production system effective boundary splicing, reverse rotation and translation to restore coordinate points and injection-production system representation property calculation are involved in the process of determining the injection-production system properties, and the accuracy of the water drive control degree and the water drive producing degree of the injection-production relationships of the injection-production well networks on each sand body plane in the injection-production system under the constraints of structure and lithology is improved.

[0051] The application further provides a computing system for intelligently identifying effective connectivity of injection-production systems under constraints, which realizes automation of calculation of water-drive controlled reserves and produced reserves by means of a data preprocessing module, a first data calculation module, an attribute determination module, a second data calculation module and a third data calculation module, intelligently identifies effective connectivity and effective action time of the injection-production systems through the computing system, and further calculates parameters such as water-drive theoretical swept area and water-drive produced area, and finally obtains water-drive control degree and production degree of single sand body, different development units and the whole reservoir, thereby ensuring calculation accuracy of the water-drive control degree and the production degree. BRIEF DESCRIPTION OF DRAWINGS

[0052] Figure 1 A flow chart of the computing method for intelligently identifying effective connectivity of injection-production systems under constraints in the application;

[0053] Figure 2 A principle diagram of the computing system for intelligently identifying effective connectivity of injection-production systems under constraints in the application;

[0054] Figure 3 A flow chart of attribute determination of the injection-production system in the application;

[0055] Figure 4 A schematic diagram of the mathematical model of the injection-production system in the application;

[0056] Figure 5 A schematic diagram of coordinate offset rotation calculation and restoration model in the application;

[0057] Figure 6 A schematic diagram of the mechanism model structure of effective swept area identification of the injection-production system in the application;

[0058] In the figure: 1-data preprocessing module; 2-first data calculation module; 3-attribute determination module; 4-second data calculation module; 5-third data calculation module. DETAILED DESCRIPTION

[0059] In order to make the personnel in the technical field better understand the application scheme, the technical solutions in the embodiments of the application will be described clearly and completely in combination with the drawings in the embodiments of the application. Obviously, the described embodiments are only a part of the embodiments of the application, rather than all the embodiments. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative labor should belong to the protection scope of the application.

[0060] It is to be understood that the terms "first", "second", and the like, used in the description and the claims of the application and the above description of the drawings, are used to distinguish between similar objects and are not necessarily used to describe a particular sequential or chronological order. It is to be understood that the use of data so designated is not to be construed as limiting of the embodiments of the application described herein to only those embodiments absolutely recited in the description and the claims of the application. Further, the terms "include," "includes," and "including" and any variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, system, product or apparatus that comprises a list of steps or units are not necessarily limited to those steps or units that are expressly listed, but can include other steps or units not expressly listed or inherent to such process, method, product or apparatus.

[0061] The application will be further described in detail below with reference to the accompanying drawings:

[0062] The purpose of the present application is to provide a kind of intelligent identification constraint under the betting system effective communication calculation method and system, to solve the technical problems of large amount of repetitive calculation work and large calculation structure error in prior art by artificial.

[0063] Embodiment 1

[0064] Referring to Figure 1 In an embodiment of the present application, an intelligent identification constraint under the betting system effective communication calculation method is provided, which belongs to the algorithm packaging part of the key technology in large software program structure, has a communication interface of lower database architecture, and the calculation process of key part is the specific design and implementation of the embodiment. The calculation method can be compiled using a general-purpose computer program design language, integrated in the server backend of B / S architecture or the client in C / S architecture, and the specific process is as follows:

[0065] Step 1, obtaining oilfield injection and production comprehensive data, standardizing and preprocessing the obtained oilfield injection and production comprehensive data;

[0066] Specifically, the pre-computation condition preparation is performed, the comprehensive data in the database, including drilling and completion data, injection and production relationship data, sand body boundary data, structure boundary data, well trajectory data, stratigraphic correlation data, interpretation conclusion data, layer opening and closing well history data, production and suction profile data, etc. are standardized and preprocessed, and the data source is standardized and preprocessed, which specifically includes five data quality standard unified processing: well name English unification, layer number dynamic and static unification, injection and production effective communication judgment standard unification, sand body boundary and oil-bearing area boundary envelope unification, and well trajectory geological engineering intersection coordinate calculation method unification. The processed comprehensive data is stored in the corresponding database table model for subsequent calculation step calling.

[0067] Step 2, the normalized pretreated oilfield injection-production comprehensive data is sequentially subjected to ground effective injection-production well pattern and subsurface effective injection-production well pattern to calculate injection-production relationship of injection-production well and injection-production well pattern on each sand body plane;

[0068] Specifically, the specific process of sequentially subjecting the normalized pretreated oilfield injection-production comprehensive data to ground effective injection-production well pattern and subsurface effective injection-production well pattern to calculate injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is as follows:

[0069] S1, the normalized pretreated oilfield injection-production comprehensive data is subjected to effective injection-production well pattern calculation to obtain ground injection-production relationship, and the ground injection-production relationship is stored in the db_surface data model;

[0070] S2, according to the ground injection-production relationship, the subsurface effective injection-production well pattern is calculated, the intermediate calculation results are expanded in three dimensions according to the wellbore, the envelope relationship with the outer boundary and the inner boundary of the single sand body is determined layer by layer through the well trajectory coordinates, the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is generated, the corresponding conclusion is output according to the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane, and the conclusion is stored in the db_under data model.

[0071] Step 3, a unary one-hub injection-production system mathematical model is established, the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is used to determine the injection-production system attribute under the identification of structure and lithology constraints, and the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane is used to determine the injection-production relationship of injection-production well and injection-production well pattern on each sand body plane. Figure 3 As shown in the specific process as follows:

[0072] Step 31, taking the midpoint of the oil well and the water well as the center of the circle, defining the filter radius as the radius of the circle, taking the sand body boundary and the fault boundary corresponding to the sand body number in the db_under data model;

[0073] Traverse the sand body boundary points, calculate the distance from the point to the center of the circle, if the distance is less than the radius of the circle, it is the record value, and store this sand body boundary point in the container arrSand,

[0074] Traverse the fault boundary points, calculate the distance from the point to the center of the circle, if the distance is less than the radius of the circle, it is the record value, and store this fault boundary point in the container arrFault;

[0075] Specifically, in this embodiment, the unary one-hub injection-production system on the single sand body is defined as the injection-production system, that is, the minimum injection-production system analysis unit, according to the percolation field theory, the unary one-hub percolation flow line formed by injection-production is a cluster of circles with the center moving on the y-axis in the homogeneous ideal model, considering that in actual production, it will be disturbed by two or more percolations, therefore, the effective swept range corresponding to the injection-production system is approximately elliptical, and the elliptical area is approximately equivalent to the effective swept area under the injection-production system;

[0076] Step 32, all point data are translated along the vector of the water well, and all points are rotated around the origin with the angle β between the oil well and the positive half of the X axis.

[0077] Step 33, given a reasonable φ value, the meta injection-production system is constructed to determine the parameters such as the major and minor axes of the ellipse.

[0078] Specifically, by calling the coordinates of the meta injection-production well points and the coordinates of the inner and outer boundaries of the sand body in the db_under data model, and simultaneously calling the coordinates of the breakpoints and the sand body boundaries in the same sand body plane, the actual two-dimensional meta injection-production system plane map affected by structure and lithology can be projected. For easy calculation, the coordinates of the projected meta injection-production system need to be translated and rotated around the water well, wherein the meta injection-production system includes the well site, the sand body boundary, the oil sand boundary, and the breakpoint, so that it becomes an ellipse with the water well at the origin of the orthogonal coordinate system and located at the left focus and the oil well at the right vertex. The relevant equations are constructed for solving.

[0079] In the mono-hydraulic injection-production system with a well spacing D, according to the percolation field theory, the water well is set at the left focus and the oil well is set at the right vertex, the well spacing D = a + c, c = φa, and then the standard equation of the ellipse can be established, as shown in the following formula. Figure 4 The complete ellipse is not affected by factors such as faults and pinch-outs; and the water drive effective swept area can be expressed as:

[0080]

[0081] In the formula, φ is the conversion relationship between the permeability, the starting pressure difference, and the major and minor axes of the ellipse in the meta injection-production system (the value range of φ is [0.5, 1]); D is the well spacing, km; Au is the effective swept area of the meta injection-production system, km 2 .

[0082] The calculation of the effective swept area of the meta injection-production system can be further optimized, and the expression is as follows:

[0083]

[0084] In the formula, ω is the three relationships between the sand body pinch-out line and the meta effective swept area in the corresponding figure (ω = φ in the figures lc and lf; ω = 0.9φ in the figures mc and mf; ω = 0.8φ in the figures nc and nf); Δx i is the grid step in the x direction; Δy i is the grid step in the y direction; M is the number of grid cells of the ineffective swept range in the shale; and N is the number of grid cells of the ineffective swept range in the sand body cut by the fault.

[0085] Step 34, the type of the meta injection-production system effective swept area identification mechanism model is determined by setting a determination condition threshold range.

[0086] Specifically, in the embodiment, the sand body boundary is depicted based on the actual data of the drilled well, and thus the injection-production system located at the sand body boundary will have a relative increase in the swept area on the other side due to the compression of the single-side flow space, such as the convex bank deposition of meandering river, the convex bank dam deposition of braided river, and the two-wing tip of the beach bar deposition. In more complex geological conditions, the formed reservoir will develop different levels of faults in the tectonic movement, and thus the injection-production system will be cut by the faults, resulting in the further compression and distortion of the effective swept area form. It is found that there are six states of the injection-production system.

[0087] According to Figure 6 As shown in the figure, the effective swept area identification mechanism model of the injection-production system includes an Ic-type injection-production system model, an If-type injection-production system model, an nc-type injection-production system model, an mc-type injection-production system model, an nf-type injection-production system model, and an mf-type injection-production system model.

[0088] When the condition threshold range is that the container arrSand count and the container arrFault count are both zero, it is the Ic-type injection-production system model.

[0089] When the condition threshold range is that the container arrSand count is zero and the container arrFault count is greater than zero, it is the If-type injection-production system model.

[0090] When the container arrSand count is greater than zero and the container arrFault count is zero, it is determined whether the sand body edge line intersects with the X positive half-axis. When there is an intersection, it is the nc-type injection-production system model, otherwise it is the mc-type injection-production system model.

[0091] When the container arrSand count and the container arrFault count are both greater than zero, it is determined whether the sand body edge line intersects with the X positive half-axis. When there is an intersection, it is the nf-type injection-production system model, otherwise it is the mf-type injection-production system model.

[0092] Step 35, after determining the injection-production system, taking the water well as the coordinate origin and the oil well as the right vertex, 44 boundary streamline points are created by using the elliptical difference value, and the envelope relationship is judged with the sand body boundary and the fault boundary to generate the effective swept area peripheral boundary points of the injection-production system.

[0093] Step 36, taking the water well as the vector, translating all point data according to the vector, rotating all points around the origin at a translation-β angle, recording the restored effective boundary of the injection-production system, and recording the data in the container.

[0094] In order to solve the unified meta injection model conveniently, the oil well is rotated by a certain angle around the water well in the meta injection system. Assuming that the oil well is point A (x, y), and after rotating around the water well B (a, b), it becomes point C (c, d). Assuming that the angle of point A before rotation is θ, the angle after rotating counterclockwise to point C is θ+β. The distance between the oil well and the water well r is unchanged, as shown in Figure 5 Therefore, we have

[0095] x / cos(θ)=y / sin(θ)=c / cos(β+θ)=d / sin(β+θ)

[0096] The coordinates of point C (c, d) after rotation can be calculated as

[0097] c=r*cos(β+θ)=r*cos(β)cos(θ)-r*sin(β)sin(θ)=xcos(β)-ysin(β)

[0098] d=r*sin(β+θ)=r*sin(β)cos(θ)+r*cos(β)sin(θ)=ycos(β)+xsin(β)

[0099] From the formula derivation, it can be seen that the coordinates c and d after rotation are only related to the coordinates x and y before rotation and the angle β of rotation. In addition, the coordinates can be restored by rotating a negative angle in the opposite direction. According to the odd-even nature of trigonometric functions, the coordinates of point A after rotating clockwise by an angle β become

[0100] x1=xcos(β)+ysin(β)

[0101] y1=ycos(β)-xsin(β)。

[0102] Step 37, calculate the meta injection system attribute, the meta injection system attribute includes address oiliness, fault blocking property, effective connectivity and effective productivity;

[0103] Among them, if the oil well interpretation conclusion is oil layer, oil-water layer, oil-gas layer, poor oil layer, then it is geologically oil-bearing, and the value is 1, otherwise the value is 0;

[0104] In the fault blocking property, if the intersection of the fault and the X-axis in the orthogonal coordinate system satisfies greater than 0 and less than the distance between the oil well and the water well, it is considered as fault blocking, and the value is 1, otherwise the value is 0;

[0105] In the effective connectivity, the perforation state of the oil well and the water well is counted per year. If the oil well and the water well are in perforation production state within the statistical year, the value is 1, otherwise the value is 0;

[0106] The effective utilization is 1 if the corresponding production and absorption or the corresponding test result is effective, and 0 otherwise for the middle-high permeability reservoir or the reservoir with complete production data. For the low permeability reservoir or the reservoir with missing production data, the water absorption of the water well is used for judgment. If the test result shows water absorption or the test result does not show water absorption but the KH proportion is greater than 0.25, the value is 1, otherwise the value is 0.

[0107] The oil-bearing property is used to determine whether the oil well is in the oil-bearing area. The value is 1 for subsequent calculation, otherwise it is marked as 0 and does not participate in the calculation of swept area.

[0108] The control effectiveness is calculated only when the injection-production connectivity objective effectiveness value is 1 and the oil-bearing property value is 1. Otherwise, the value is directly assigned as 0. The effective range of the corresponding number axis is composed of the perforation and plugging time of the oil and water wells. The control effectiveness of the injection-production system is calculated year by year with the water well betting time as the starting point and the step length of 1 year. If the year meets the effective perforation range of the oil well and the water well at the same time, it is determined to be controlled effectively, and the value is 1, otherwise the value is 0. The annual calculation results of the injection-production system composed of the same pair of oil and water wells are stored in db_el.

[0109] The injection-production utilization state is calculated only when the control effectiveness value is 1. Otherwise, the value is directly assigned as 0. The production and absorption profile test results of the oil and water wells are extracted by computer algorithm. The injection-production system with the same sand body number and effective production and absorption test corresponding elements is determined to be injection-production connected or injection-production utilization effective, and the state result is stored in db_el.

[0110] Step 38, record all process parameters and attributes in DB-eil data model.

[0111] Step 4, calculate the effective swept and superimposed area and the effective utilization superimposed area in the single sand body according to the injection-production system attribute. The water drive control degree and the water drive utilization degree of the single sand body are calculated according to the area method. The swept area boundary coordinates of the single sand body are extracted according to the water drive control degree and the water drive utilization degree of the single sand body.

[0112] Specifically, by calling each element of the effective swept area stored in db_el contained in a single sand body, further development of each element of the swept area superposition calculation of injection-production system in single sand body is carried out. This method is based on the disclosed computer two-dimensional geometric calculation method, and the maximum water drive swept area in multi-direction water drive is error corrected by eliminating the grid number of the overlapping swept area of multiple injection-production systems. For each element of the injection-production system in the same single sand body, the oil saturation, oil density, volume coefficient and other parameters are considered to be consistent. According to the volumetric method of reserve calculation formula, the single sand body water drive reserve control degree expression of reserve definition method can be simplified as single sand body water drive control area / single sand body oil-bearing area. For the quantitative calculation of single sand body water drive producing reserves, the same is based on the quantitative calculation of single sand body water drive control reserves, and the effective injection-production area of the element injection-production system marked in the previous module is called, so that the expression of the degree of single sand body water drive reserve production can be simplified as single sand body water drive swept area / single sand body oil-bearing area, and the two-degree results of single sand body are stored in db_ewed.

[0113] Step 5, the water drive control degree and water drive production degree of all single sand bodies contained in the reservoir and development unit are summed up, and the calculation work of effective connectivity of injection-production system under intelligent identification constraint is completed.

[0114] Specifically, the water drive control degree calculation formula of all single sand bodies contained in the reservoir and development unit is as follows:

[0115]

[0116] The water drive production degree calculation formula of all single sand bodies contained in the reservoir and development unit is as follows:

[0117]

[0118] Wherein, R is the number of oil-bearing single sand bodies of the development unit; N oi is the geological reserves of each sand body, t; N di is the water drive production degree of each single sand body, t; E d is the water drive reserve production degree of the development unit, %; N Wi is the water drive reserve control degree of the development unit; N i is the water drive swept reserves of each sand body.

[0119] The embodiment aims to provide a convenient and accurate water drive "two-degree" index statistics and calculation tool for reservoir development management and researchers. The tool can intelligently identify the effective connection relationship and effective action time of the injection-production system based on computer program algorithm, and then calculate parameters such as water drive theoretical swept area and water drive reserve producing area, and finally obtain the water drive control degree and producing degree of single sand body, different development units and the whole reservoir. The tool can replace the traditional manual calculation mode, not only can greatly improve the work efficiency per capita, reduce the cost of manual time, increase the cost benefit for the enterprise, but also can provide more reliable and real water drive reserve control degree and producing degree for the reservoir development decision makers, and automatically provide the corresponding oil and water well adjustment and potential tapping workload suggestion.

[0120] Embodiment 2

[0121] According to Figure 2 The embodiment provides a calculation system for intelligently identifying the effective connection of injection-production system under constraints, which is used for realizing the calculation method for intelligently identifying the effective connection of injection-production system under constraints, and comprises the following steps:

[0122] A data preprocessing module 1 is configured to acquire oilfield injection-production comprehensive data, and perform standardization preprocessing on the acquired oilfield injection-production comprehensive data.

[0123] A first data calculation module 2 is configured to sequentially perform ground effective injection-production well pattern and underground effective injection-production well pattern calculation on the oilfield injection-production comprehensive data after standardization preprocessing, so as to calculate the injection-production relationship of injection-production wells and injection-production well pattern on each sand body plane.

[0124] An attribute determination module 3 is configured to establish a one-element-one-hub-element injection-production system mathematical model, and determine the element injection-production system attribute under the recognition constraints of structure and lithology according to the injection-production relationship of injection-production wells and injection-production well pattern on each sand body plane.

[0125] A second data calculation module 4 is configured to obtain the effective swept superposition area and effective producing superposition area in a single sand body according to the element injection-production system attribute, calculate the water drive control degree and water drive producing degree of the single sand body according to the area method, and extract the swept area boundary coordinates of the single sand body according to the water drive control degree and water drive producing degree of the single sand body.

[0126] A third data calculation module 5 is configured to sum and aggregate the water drive control degree and water drive producing degree of all single sand bodies contained in the reservoir and development unit, and complete the calculation of the effective connection of injection-production system under constraints.

[0127] The automation of the water drive controlled reserve and the producing reserve calculation is realized by the data preprocessing module, the first data calculation module, the attribute determination module, the second data calculation module and the third data calculation module in the embodiment, the effective connection relationship and the effective action time of the injection-production system are intelligently identified by the calculation system, and then the water drive theoretical swept area, the water drive producing area and other parameters are calculated, and finally the water drive control degree and the producing degree of the single sand body, different development units and the whole reservoir are obtained, so that the calculation accuracy of the water drive controlled reserve and the producing reserve is ensured.

[0128] To sum up, the application provides a calculation method and system for intelligently identifying the effective connection of the injection-production system under constraints, realizes the automation of the water drive controlled reserve and the producing reserve calculation, and is tested in Wunan Oilfield in Qaidam Basin. The constraint conditions that determine the calculation accuracy of the reserve definition are limited to the actual underground injection-production well spacing, the pre-adjusted injection-production well spacing, the single sand body drawing precision and the structure interpretation precision, which is more inclined to the guidance of geological results compared with the traditional method. From the aspects of reservoir potential tapping and development contradictions, the calculation results of the reserve definition method are compared, and it is found that the actual Ew value calculated by the reserve definition method is reduced by 3.8 percentage points, but Ed is reduced by 11.5 percentage points, which is more in line with the actual situation of serious casing damage and uneven vertical water absorption in the reservoir development. The whole calculation greatly optimizes the repetitive work flow of the previous manual "thickness method" calculation, improves the research and analysis work efficiency, and the unified calculation based on the comprehensive basic data of geology, engineering, production and dynamic monitoring takes about 15 hours, the new well maintenance takes 1-5 minutes, and the "water drive two-degree calculation process" takes 7 hours, and the local update takes 5-10 minutes. For the in-depth statistical analysis of the intermediate results of step four, the remaining work amount of the reservoir injection-production system improvement and adjustment can be obtained, which helps the technical personnel to arrange the scheme work amount faster. In summary, the application technology can complete the work amount of 3-4 months of traditional manual work through about 30 hours of automatic calculation, and has significant benefits for the business support of reservoir development research.

[0129] The database ORM mapping technology of EFCore and the C# language programming technology are used in the application, and the back-end server program under B / S is realized. The method of the application runs through the main program of the server client, wherein the key step 3 is encapsulated as a.dll dynamic support library, and the step 1 of configuration automation is encapsulated as a.dll dynamic support library.

[0130] It should be pointed out finally that the above embodiments are only used to illustrate the technical solutions of the present application but not to limit it. Although the present application has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that the specific embodiments of the present application can be modified or replaced equivalently without departing from the spirit and scope of the present application, and any modification or equivalent replacement without departing from the spirit and scope of the present application should be covered in the protection scope of the claims of the present application.

Claims

1. A computing method for intelligently identifying valid connectivity of a bet placement system under constraints, characterized in that, The method comprises the following steps: obtaining oilfield injection-production comprehensive data, and performing standardization preprocessing on the obtained oilfield injection-production comprehensive data; performing ground effective injection-production well pattern and underground effective injection-production well pattern calculation on the oilfield injection-production comprehensive data after standardization preprocessing to obtain injection-production relationship of injection-production wells and injection-production well patterns on each sand body plane; establishing a one-to-one injection-production system mathematical model, and determining injection-production system attributes under the constraints of structure and lithology according to injection-production relationship of injection-production wells and injection-production well patterns on each sand body plane; obtaining effective swept superimposed area and effective produced superimposed area in a single sand body according to injection-production system attributes, calculating water drive control degree and water drive production degree of the single sand body according to area method, and extracting single sand body swept area boundary coordinates according to water drive control degree and water drive production degree of the single sand body; summing up water drive control degree and water drive production degree of all single sand bodies contained in a reservoir and a development unit to complete calculation of effective connectivity of injection-production system under intelligent identification constraints.

2. The method of claim 1, wherein the method is effective for communication of the system under the constraint of intelligent identification and betting, characterized in that, The oilfield injection-production comprehensive data comprises drilling and completion data, injection-production relationship data, sand body boundary data, structure boundary data, well trajectory data, stratum correlation data, interpretation conclusion data, layer opening and closing well history data, and production and suction profile data.

3. The method of claim 1, wherein the method further comprises: The standardization preprocessing process comprises the following steps in sequence: unifying well names in English, unifying static and dynamic layer numbers, unifying injection-production effective connectivity determination standards, unifying sand body boundary and oil-bearing area boundary envelope, and unifying well trajectory geological engineering intersection point coordinate calculation methods.

4. The method of claim 1, wherein the method is effective for communication of the system under the smart identification constraint betting system. The specific process of performing ground effective injection-production well pattern and underground effective injection-production well pattern calculation on the oilfield injection-production comprehensive data after standardization preprocessing to obtain injection-production relationship of injection-production wells and injection-production well patterns on each sand body plane is as follows: S1, performing effective injection-production well pattern calculation on the oilfield injection-production comprehensive data after standardization preprocessing to obtain ground injection-production relationship, and storing the ground injection-production relationship into a db_surface data model; S2, performing underground effective injection-production well pattern calculation according to the ground injection-production relationship, calculating intermediate results in three-dimensional expansion from a wellbore, determining envelope relationship with single sand body outer boundary and inner boundary through well trajectory coordinates layer by layer, generating injection-production relationship of injection-production wells and injection-production well patterns on each sand body plane, outputting corresponding conclusions according to the injection-production relationship of injection-production wells and injection-production well patterns on each sand body plane, and storing the conclusions into a db_under data model.

5. The method of claim 1, wherein the method further comprises: The process of establishing a one-to-one injection-production system mathematical model and calculating injection-production system attributes under the constraints of structure and lithology according to injection-production relationship of injection-production wells and injection-production well patterns on each sand body plane is as follows an orthogonal coordinate system is established according to percolation field theory, injection-production well point coordinates and sand body inner and outer boundary coordinates are called, and breakpoint coordinates and sand body boundary coordinates in the same sand body plane are called, and an actual two-dimensional injection-production system plane affected by structure and lithology is projected in the orthogonal coordinate system; the coordinates of the actual two-dimensional injection-production system plane are translated and rotated as a whole around a water well, so that the water well is at the origin of the orthogonal coordinate system and located at the left focus point, and the oil well is at the right vertex of an ellipse, and the effective swept area under the injection-production system is calculated by calculating the area of the ellipse; The set determination condition threshold range determines the effective swept area of the injection-production system to identify the mechanism model, judges the degree of the swept form of the injection-production system being cut by the fault line and the sand body pinch-out line to obtain the injection-production system attribute, and if the planar graph combination relationship meets the condition threshold range, it is determined that the injection-production connectivity of the injection-production system attribute is objectively effective, and is assigned a value of 1; Otherwise, it is determined that the injection-production connectivity of the injection-production system attribute is invalid, and is assigned a value of 0.

6. The method of claim 5, wherein the method further comprises: The mechanism model for identifying the effective swept area of the injection-production system includes an Ic-type injection-production system model, an If-type injection-production system model, an nc-type injection-production system model, an mc-type injection-production system model, an nf-type injection-production system model, and an mf-type injection-production system model; When the condition threshold range is that the container arrSand count and the container arrFault count are both zero, it is the Ic-type injection-production system model; When the condition threshold range is that the container arrSand count is zero and the container arrFault count is greater than zero, it is the If-type injection-production system model; When the container arrSand count is greater than zero and the container arrFault count is zero, it is determined whether there is an intersection between the sand body edge line and the X positive half axis, and when there is an intersection, it is the nc-type injection-production system model, otherwise it is the mc-type injection-production system model; When the container arrSand count and the container arrFault count are both greater than zero, it is determined whether there is an intersection between the sand body edge line and the X positive half axis, and when there is an intersection, it is the nf-type injection-production system model, otherwise it is the mf-type injection-production system model.

7. The method of claim 5, wherein the method further comprises: The injection-production system attribute includes oil-bearing property, fault blocking property, effective connectivity, and effective dynamic property; In the oil-bearing property, if the oil well interpretation conclusion is oil layer, oil-water layer, oil-gas layer, or poor oil layer, it is considered to be geologically oil-bearing, and is assigned a value of 1, otherwise it is assigned a value of 0; In the fault blocking property, if the intersection between the fault and the X axis in the orthogonal coordinate system satisfies greater than 0 and less than the distance between the oil well and the water well, it is considered to be fault blocking, and is assigned a value of 1, otherwise it is assigned a value of 0; In the effective connectivity, the perforation state of the oil well and the water well is counted per year, and if both the oil well and the water well are in the perforation production state within the counted year, it is assigned a value of 1, otherwise it is assigned a value of 0; In the effective dynamic property, for a medium-high permeability reservoir or a reservoir with complete liquid production data, if there is corresponding production absorption or corresponding micro-test, it is assigned a value of 1, otherwise it is assigned a value of 0; for a low-permeability reservoir or a reservoir with missing liquid production data, the water absorption of the water well is used for judgment, and if the test shows water absorption or the test does not show water absorption but the KH ratio is greater than 0.25, it is assigned a value of 1, otherwise it is assigned a value of 0.

8. The method of claim 1, wherein the system is effective to communicate under the constraint of intelligent identification. The calculation process of the water drive control degree of a single sand body is as follows: The oil saturation parameters, the crude oil density parameters, and the volume coefficient parameters of each injection-production system in the same single sand body are considered to be consistent, and according to the volumetric method reserve calculation formula, the single sand body water drive reserve control degree expression of the reserve definition method can be simplified as single sand body water drive control area / single sand body oil-bearing area; The calculation process of the water drive dynamic degree of a single sand body is as follows: For single sand body water drive dynamic reserves quantification calculation, also based on single sand body water drive control reserves quantification calculation, through calling the marked effective injection-production system in the former module, the single sand body water drive reserves producing degree expression can be simplified as single sand body water drive swept area / single sand body oil-bearing area.

9. The calculation method for effective connectivity of an injection-progression system under intelligent identification constraints according to claim 1, characterized in that, The water drive control degree calculation formula of all single sand bodies contained in the reservoir and development unit is as follows: The water drive dynamic degree calculation formula of all single sand bodies contained in the reservoir and development unit is as follows: Wherein, R is the number of oil-bearing single sand body of development unit; N oi is the geological reserves of each sand body, t; N di is the water drive displacement degree of each single sand body, t; E d is the water drive reserves displacement degree of development unit, %; N Wi is the water drive reserves control degree of development unit; N i is the water drive swept reserves of each sand body.

10. A computing system for intelligently identifying valid connectivity of a restricted betting system, for implementing the method of intelligently identifying valid connectivity of a restricted betting system according to any one of claims 1-9, characterized in that, It includes: The data preprocessing module (1) is used for obtaining oilfield injection-production comprehensive data, and standardizing and preprocessing the obtained oilfield injection-production comprehensive data; The first data calculation module (2) is used for calculating the injection-production relationship of injection wells and injection-production well pattern on each sand body plane by sequentially performing ground effective injection-production well pattern and underground effective injection-production well pattern on the standardized and preprocessed oilfield injection-production comprehensive data; The attribute determination module (3) is used for establishing a one-to-one injection-production system mathematical model, and determining the injection-production system attribute under the recognition constraint of structure and lithology according to the injection-production relationship of injection wells and injection-production well pattern on each sand body plane; The second data calculation module (4) is used for obtaining the effective swept superimposed area and effective producing superimposed area in the single sand body, calculating the water drive control degree and water drive dynamic degree of the single sand body according to the area method, and extracting the swept area boundary coordinates of the single sand body according to the water drive control degree and water drive dynamic degree of the single sand body; The third data calculation module (5) is used for summing and collecting the water drive control degree and water drive dynamic degree of all single sand bodies contained in the reservoir and development unit, and completing the calculation work of the effective connectivity of the injection-production system under the intelligent recognition constraint.

Images