CN111080789B - Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area - Google Patents

Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area Download PDF

Info

Publication number
CN111080789B
CN111080789B CN201911367075.2A CN201911367075A CN111080789B CN 111080789 B CN111080789 B CN 111080789B CN 201911367075 A CN201911367075 A CN 201911367075A CN 111080789 B CN111080789 B CN 111080789B
Authority
CN
China
Prior art keywords
well
oil
control area
area
grid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911367075.2A
Other languages
Chinese (zh)
Other versions
CN111080789A (en
Inventor
张波
刘慧卿
叶继根
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201911367075.2A priority Critical patent/CN111080789B/en
Publication of CN111080789A publication Critical patent/CN111080789A/en
Application granted granted Critical
Publication of CN111080789B publication Critical patent/CN111080789B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/30—Specific pattern of wells, e.g. optimizing the spacing of wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/30—Specific pattern of wells, e.g. optimizing the spacing of wells
    • E21B43/305—Specific pattern of wells, e.g. optimizing the spacing of wells comprising at least one inclined or horizontal well
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining

Abstract

The application discloses a method and a device for determining a well position of a ciphered well in a complex fault block oil reservoir exploitation area, wherein the method comprises the following steps: acquiring a three-dimensional geological model of a complex fault block oil reservoir exploitation area; aiming at each non-well point grid in the three-dimensional geological model, determining the number of oil layers encountered by single well drilling and the actual control area of the single well in each oil layer when the well is drilled in the non-well point grid according to geological conditions; determining the remaining movable oil reserves of the single well in each oil layer according to the number of grids contained in the actual control area of the single well in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserve of the single well; and determining the well-point-free grid corresponding to the maximum total remaining movable oil reserves as the deployment well position of the encrypted well. The well position of the encryption well can be determined through a reasonable method, and the reserve utilization degree of the complex fault block oil reservoir exploitation area is improved.

Description

Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area
Technical Field
The application relates to the technical field of petroleum exploration and development, in particular to a method and a device for determining a well position of an encrypted well in a complex fault block oil reservoir exploitation area.
Background
This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
The complex fault block oil reservoir has large reserves when being stored in the world oil reserves, the number of faults in the complex fault block oil reservoir is large and complex, and the oil reservoir is cut into a plurality of fault blocks with different sand body forms. The complex fault block oil reservoir has small and complex oil-containing area and is provided with a plurality of sets of oil-water systems, each set of oil-water system has different reservoir characteristics and stratum structures, and generally, the natural energy is insufficient. In the complex fault block oil reservoir, the multilayer complex fault block oil reservoir occupies a very high specific gravity, the complex fault block oil reservoir develops a plurality of oil layers in the longitudinal direction, the interlayer heterogeneity is large, and the reservoir thickness and the geological reserves of different sand bodies are different.
Domestic and foreign researches show that the oil reservoir with sufficient natural energy is suitable to be exploited in a well pattern encryption mode. The well pattern encryption is to reduce the well distance and increase the density of the well pattern on the basis of fully analyzing the reserve consumption degree, to recombine the well pattern for the reservoir which is not used and the oil layer which has poor consumption degree, and to adjust the well ratio of the unit oil-containing area, thereby realizing the increase of the recoverable reserve. For oil reservoirs with insufficient natural energy, especially complex fault-block oil reservoirs in the later development stage and with seriously decreased yield, injection wells are additionally added to form an injection-production well pattern, so that the formation energy can be increased, and the recovery efficiency is improved. Therefore, the well position of the encryption well is optimized, the technology recoverable reserves can be effectively increased, the reserve utilization degree is improved, the oil deposit recovery degree is improved, and the development effect of the oil field is improved.
At present, the conventional method for selecting the well position of the encrypted well is to select the well position of the encrypted well according to the type of the well pattern to be deployed, for example, the well position of the encrypted well to be deployed is determined according to the well pattern of the nine-point method, the five-point method, and the like. However, the area of the exploitation region of the complex fault block oil reservoir is usually small, and generally no well pattern is clear, so that the well position of the infill well is difficult to determine by adopting a conventional method. Even if an obvious well pattern exists in a complex fault block oil reservoir exploitation area, a plurality of oil layers and a plurality of sets of oil-water systems are distributed in the underground of the area, and the purpose of improving the reserve utilization degree cannot be well achieved by selecting the deployed encryption well through a conventional method.
Disclosure of Invention
The embodiment of the application provides a method for determining a well position of a ciphered well in a complex fault block oil reservoir exploitation area, which is used for determining the well position of the ciphered well through a reasonable method and improving the reserve utilization degree of the complex fault block oil reservoir exploitation area, and comprises the following steps:
acquiring a three-dimensional geological model established according to the geological condition of the complex fault block oil reservoir exploitation area, wherein the three-dimensional geological model is composed of grids, the grids are associated with well arrangement conditions, and the residual oil saturation and pore volume of the corresponding complex fault block oil reservoir exploitation area are obtained; aiming at each non-well point grid in the three-dimensional geological model, determining the number of oil layers encountered by single well drilling and the actual control area of the single well in each oil layer when the well is drilled in the non-well point grid according to geological conditions; determining the remaining movable oil reserves of the single well in each oil layer according to the number of grids contained in the actual control area of the single well in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserve of the single well; and comparing the total residual movable oil reserves of the single wells drilled in all the non-well point grids, and determining the non-well point grid corresponding to the maximum total residual movable oil reserve as the deployment well position of the encrypted well.
The embodiment of the present application still provides a complicated fault block oil reservoir exploitation region encrypted well position determining device for confirm the well position of encrypted well through reasonable method, improve complicated fault block oil reservoir exploitation region's reserves and use the degree, the device includes:
the acquisition module is used for acquiring a three-dimensional geological model established according to the geological condition of the complex fault block oil reservoir exploitation area, wherein the three-dimensional geological model is composed of grids which are associated with well arrangement conditions, and the residual oil saturation and pore volume of the corresponding complex fault block oil reservoir exploitation area; the determining module is used for determining the number of oil layers encountered by single well drilling and the actual control area of the single well in each oil layer when the well is drilled in the well-free point grid according to the geological condition aiming at each well-free point grid in the three-dimensional geological model acquired by the acquiring module; determining the remaining movable oil reserves of the single well in each oil layer according to the number of grids contained in the actual control area of the single well in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserve of the single well; the determining module is further used for comparing the total remaining movable oil reserves of the single wells drilled in all the non-well point grids, and determining the non-well point grid corresponding to the maximum total remaining movable oil reserves as the deployment well position of the encrypted well.
In the embodiment of the application, when a single well is deployed in all grids without well points, the residual movable oil reserves of the single well are calculated, and the grid with the maximum residual movable oil reserves is selected as the deployed well position of the single well. By the quantitative calculation method, the purpose of improving the reserve utilization degree can be achieved by definitely determining the mesh in which the single well is deployed and obtaining the largest oil quantity by utilizing the single well; in addition, the method has important significance for recognizing the particularity of the complex fault block oil reservoir and researching the complex fault block oil reservoir, particularly for researching how to economically and reasonably develop the complex fault block oil reservoir which is in the later development stage and has seriously reduced yield.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a flow chart of a method for determining a well position of a infill well in a complex fault block reservoir production area according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a grid in a three-dimensional geological model according to an embodiment of the present application;
FIG. 3(a) is a schematic representation of a well trajectory of a deviated well in an embodiment of the present application;
FIG. 3(b) is a schematic representation of a well trajectory of another deviated well according to an embodiment of the present application;
FIG. 4(a) is a schematic view of a horizontal well section of a horizontal well in an embodiment of the present disclosure;
FIG. 4(b) is a schematic view of another horizontal well section in an embodiment of the present disclosure;
FIG. 4(c) is a schematic view of another horizontal well section in an embodiment of the present disclosure;
FIG. 4(d) is a schematic view of another horizontal well section in an embodiment of the present disclosure;
FIG. 4(e) is a schematic view of another horizontal well section in an embodiment of the present disclosure;
FIG. 5(a) is a schematic illustration of an injection-production well pattern in an embodiment of the present application;
FIG. 5(b) is a schematic diagram of an irregular spot-like water injection well pattern according to an embodiment of the present application;
FIG. 5(c) is a schematic view of a drainage pattern in an embodiment of the present application;
FIG. 5(d) is a schematic illustration of a cut injection and production well pattern in an embodiment of the present application;
FIG. 5(e) is a schematic illustration of depleted single well site production in an embodiment of the present application;
fig. 6 is a schematic structural diagram of a device for determining a position of a infill well in a complex fault block reservoir production area in an embodiment of the present application.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present application more apparent, the embodiments of the present application are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present application are provided herein to explain the present application and not to limit the present application.
At present, the well position of the encryption well of the complex fault block oil reservoir is selected without the research of a general quantitative detailed calculation method, so that the residual oil saturation of part of the oil reservoir of the complex fault block oil reservoir is often larger, the residual oil saturation of other oil reservoirs is smaller or the complex fault block oil reservoir is already flooded, the development is uneven, the difficulty of subsequent conversion of a development mode is large, and the economic benefit is low. Therefore, the conventional well location selection method is not completely suitable for multilayer complex fault block oil reservoirs, and aims at the scale characteristics of different fault blocks, and adopts a non-uniform well arrangement mode to reasonably arrange encryption wells and perform well pattern optimization adjustment according to the shapes of sand bodies and residual oil enrichment areas of single or combined oil layers.
The main factors influencing the well position determination of the encrypted well comprise the following factors:
(1) the remaining mobile oil reserve. The residual movable oil reserves are the material basis of the oil well yield, and generally speaking, under the similar reservoir physical property conditions, the larger the residual movable oil reserves of the oil reservoir are, the more the oil production of the oil well is, and the better the economic benefit is. The infill well should control as large a reserve of mobile oil remaining as possible.
(2) Physical properties of the oil layer. Generally, the better the reservoir properties, the more oil production from the well, but the higher the water content of the well may be quickly and adversely affecting the well production.
(3) Degree of flooding of oil layer. The flooding degree of an oil layer is one of important indexes reflecting the water production rate of the oil layer, and generally speaking, the higher the flooding degree of the oil layer is, the higher the water production rate of an oil well is. The reservoir flooding controlled by the infill wells should be as low as possible.
(4) Number of oil layers. Under the condition that the controlled residual movable oil reserves are the same, the more the layers of oil encountered by the encryption well are, the less the residual movable oil reserves of each layer are. Therefore, the more the number of layers of oil to be encountered is, the better the well position of the well is.
(5) Maximum reservoir depth. The depth of the oil layer directly determines the depth of the shaft, and the larger the depth of the oil layer is, the larger the well depth of the encrypted well is, and the higher the drilling cost is.
(6) Faults, and the like. Infill wells, whether vertical, deviated or horizontal, the well track should not cross the fault as the reservoirs on either side of the fault generally do not belong to the same pressure system and oil and water system.
The residual movable oil reserves are the most main factors influencing the well position selection of the encryption well, and generally, the larger the residual movable oil reserves controlled by the encryption well, the more the accumulated oil production. Therefore, the basic principle of the well location selection of the encrypted well is as follows: and taking a certain well-point-free grid in the geological model as a center, calculating the residual movable oil reserves of the encryption well in the control area of each penetrated oil layer, simultaneously recording parameters such as permeability, oil saturation and the like of the penetrated oil layer of the oil well, overlapping the residual movable oil reserves of all the penetrated oil layers of the oil well in a certain set of strata, determining the residual movable oil reserves of all the penetrated oil layers of the encryption well, and selecting the well-point-free grid with the largest residual movable oil reserve to arrange the encryption well.
According to the basic principle of the infill well position selection, the embodiment of the application provides a infill well position determination method for a complex fault block oil reservoir production area, as shown in fig. 1, the method includes steps 101 to 103:
step 101, obtaining a three-dimensional geological model established according to the geological condition of the complex fault block oil reservoir exploitation region.
The three-dimensional geological model is formed by grids, the grids are related to well arrangement conditions, and the residual oil saturation and pore volume of the corresponding complex fault block oil reservoir exploitation area are obtained.
The three-dimensional geological model is established according to the geological condition of the complex fault block oil reservoir exploitation area, and the well arrangement condition, the residual saturation and the pore volume of the related grid in each grid need to be obtained from oil reservoir numerical simulation data. Specifically, a three-dimensional geological model is established according to the geological condition of the complex fault block oil reservoir exploitation area and the well arrangement condition; acquiring numerical reservoir simulation data, wherein the numerical reservoir simulation data comprises well arrangement conditions, residual oil saturation and pore volume at different relative coordinates; establishing a corresponding relation between geodetic coordinates of a three-dimensional geological model grid and relative coordinates of numerical reservoir simulation data; and associating each grid with the well placement, remaining oil saturation and pore volume at the corresponding relative coordinates.
The method for establishing the corresponding relation between the geodetic coordinates of the three-dimensional geological model grid and the relative coordinates of the numerical reservoir simulation data comprises the following steps: acquiring a proportional relation between a set three-dimensional geological model and a complex fault block oil reservoir exploitation region; determining the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate according to the proportional relation; and determining the relative coordinates corresponding to the abscissa and the ordinate of the three-dimensional geological model grid through interpolation calculation according to the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate.
The process of correspondence between the relative coordinates of the three-dimensional geological model grid and the geodetic coordinates of the numerical reservoir simulation data, and the process of associating the grid with the numerical reservoir simulation data will be described in detail below with reference to the accompanying drawings.
(1) Establishment of corresponding relation between geodetic coordinates and relative coordinates of three-dimensional geological model grid
Petrel is widely applied as three-dimensional modeling software of a complex fault block oil reservoir and Eclipse is widely applied as numerical reservoir simulation software, and the Petrel and Eclipse respectively become mainstream tool software for geological modeling and numerical reservoir simulation at present. The grids in the three-dimensional geological model data file generated by Petrel are all represented by geodetic coordinates, and the coordinates used in numerical reservoir simulation are relative coordinates. And (3) carrying out well position optimization on the encrypted well, and firstly establishing a corresponding relation between the three-dimensional geological model grid geodetic coordinates and the relative coordinates in the numerical reservoir simulation data.
The corner point grid model in the three-dimensional geological model comprises two keywords of COORD and ZCORN. Wherein the data stored in COORD is one line of coordinates, totaling (NDIVIX +1) × (NDIVIY +1) lines. Each coordinate line is defined by 2 different points, each of which is a three-dimensional data point corresponding to geodetic coordinates, i.e. one coordinate line is composed of 6 data points. The first three data of each line are vertex coordinates, and the last three data are bottom coordinates. And the three-dimensional model grid defines geodetic coordinates (x, y, z) corresponding to all grid corner points (i, j, k) through COORD keywords.
Stored in the ZCORN key are depth values (k values) of the corner models. Each grid block has 8 corners, each corner corresponding to a depth value, so there are 8NDIVIX NDIVIY x NDIVIZ data points in the key. And the coordinate data of the i, j direction of each mesh can be obtained by interpolation with the depth value of each mesh.
Each grid point corresponds to four lines and 8 corner points, and as shown in fig. 2, grid (i, j, k) corresponds to four lines, each controlled by two points whose specific coordinate values are already given in the COORD key.
All the relative coordinate values (24 data points in total) of the 8 angular points of the grid can be calculated through interpolation according to the depth values of the 8 angular points and the known relative coordinate values of the top points and the bottom points of the 4 lines, so that the one-to-one corresponding relation between the relative coordinates (i, j, k) of the grid and the ground coordinates of the grid is established.
(2) Reading of relevant attributes in a three-dimensional geological model grid
The change in saturation field in the three-dimensional model is directly related to the remaining oil distribution and the remaining mobile oil reserves. After the oil reservoir numerical simulation history fitting is completed, the oil reservoir numerical simulation data file generated by Eclipse software contains saturation field information before and after the history fitting. By utilizing the functions of Eclipse software, an oil saturation field is derived from a generated data file to form a new data file, wherein oil saturation data are sequentially stored in the x direction, the y direction and the z direction in sequence and are in one-to-one correspondence with grids in a three-dimensional geological model, and thus, the one-to-one correspondence relation between oil saturation (Soil) and grid coordinates (i, j, k) is established.
The correlation between the permeability (PERMX, PERMY, PERMZ), pore volume (PORV), and other properties of the grid and the grid is established similarly to the reading of the oil saturation.
For river facies complex fault block oil fields, faults and abandoned river channels develop, and a single sand layer is divided into a plurality of single sand bodies (independent oil and gas reservoirs). Different single sand bodies are distinguished by adopting a method of dividing a balance area, and one single sand body corresponds to one balance area. In a similar way, a correspondence between the partition data and the grid coordinates (i, j, k) can be established.
(3) Establishment of incidence relation between layer position and grid of existing oil-water well
The position coordinate of the existing oil-water well is an indispensable parameter for calculating the distance between the oil-water well and the encrypted well, and the three-dimensional geological model does not contain the position coordinate information of the existing oil-water well, so the position coordinate information needs to be loaded into the three-dimensional geological model.
The oil-water well horizon coordinate provided by geology personnel generally comprises a well name, a horizon name, a vertical depth, an X coordinate and a Y coordinate. If the position coordinates of the oil-water well are directly loaded into the model, when the distance between the oil-water well and the encrypted well is calculated, because one oil well or water well often passes through a plurality of single sand layers, each single sand layer needs to be subjected to one comparison operation, and the calculation amount is large. The position coordinates of the oil-water wells are carefully observed, the plane distance between the well mouths and the well bottoms of most wells is usually less than the length of half a grid, so that the coordinates of all the passing positions of the oil-water wells are averaged, namely the position coordinates of any oil layer penetrated by one oil-water well are the same, and the position coordinates are corresponding to the grid of the three-dimensional geological model, so that the calculation is greatly simplified.
Step 102, aiming at each non-well point grid in the three-dimensional geological model, determining the number of oil layers encountered by single well drilling and the actual control area of a single well in each oil layer when the well is drilled in the non-well point grid according to geological conditions; determining the remaining movable oil reserves of the single well in each oil layer according to the number of grids contained in the actual control area of the single well in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; and overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserve of the single well.
It should be noted that the well point-free grid is a grid in which no single well is deployed in the three-dimensional geological model, and the grid is a grid corresponding to the ground in the three-dimensional geological model.
The single well has three well types of a vertical well, an inclined well and a horizontal well, and if the well type of the encrypted well is determined before the well position of the encrypted well is selected, the residual movable yield of the single well can be directly calculated according to the step 102. If the well type of the encryption well is not determined, the residual movable production rate of the encryption well with different well types is calculated for each non-well point grid respectively, and the well type with the highest residual movable oil production rate and the non-well point grid are selected as the well type and the well position of the deployed encryption well.
Specifically, for each non-well point grid in the three-dimensional geological model, the number of oil layers encountered by the different types of single wells when the different types of single wells are drilled in the non-well point grid and the actual control areas of the different types of single wells in each oil layer are determined according to geological conditions; determining the remaining movable oil reserves of the single wells of different types in each oil layer according to the number of grids contained in the actual control area of the single wells of different types in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserves of the single wells of different types; and comparing the total residual movable oil production of the single wells with different types drilled in all the non-well point grids, and respectively determining the type of the single well corresponding to the maximum total residual movable oil production and the non-well point grids as the deployment type and the deployment position of the encrypted well.
The actual control areas of the wells are different according to the well types of the encrypted wells, and a method for calculating the actual control area of each well type of the individual well will be described in detail below.
(1) Vertical well
And when the residual movable oil reserves of the vertical well are calculated, the control radius of the encryption well is determined to be r according to the actual condition of oil field development. Assuming the number of grids in the three-dimensional geological model is i x j x k, Gij(with planar coordinates (i, j)) is one of the well-free grids, denoted GijCenter point of (1)ij(with geodetic coordinates x)i，yj) As the well drilling position of the encrypted well, r is taken as the radius, a circle is determined, and the area S covered by the circle is taken asijAs the theoretical control area size of the infill well.
Because an oil field generally consists of a plurality of oil layers, the encryption well usually passes through a plurality of oil layers when being drilled. The theoretical control area of each oil layer is the area of a circle which is formed by taking the grid of the oil layer penetrated by the encryption well as the center and taking the control radius r as the radius.
However, the different oil layer injection and production structures and sand bodies of the actual oil fieldThe boundary and the like have great difference, and the actual control area of each oil layer encryption well is often different from the theoretical control area. In the actual three-dimensional geological model, the sizes of sand bodies are greatly changed due to natural conditions such as faults, bottom water, lithologic boundaries and the like, and the theoretical control area region S of the infill wellijThe area may contain only one sand body or may contain a plurality of sand bodies. When S isijWhen the area spans a plurality of sand bodies, the fault and the impermeable strip form the area SijDivided into two parts, SijTwo sand bodies are covered. In this case, the judgment can be made by using the equilibrium region. Suppose that the sand body of the oil well belongs to the balance area 1, and the other partial area belongs to the balance area 2. The actual control area of the infill well is not the whole SijThe area of the zone is the area of the sand body where the vertical well is located in two or more sand bodies divided by the fault and the like, and the area of the sand body where the vertical well is located is used as the actual control area of the current oil layer vertical well, namely the area of the balance zone 1 is the actual control area of the encryption well. When calculating the remaining mobile oil reserves of the infill well at this layer, it is only necessary to calculate the remaining mobile oil reserves in the equilibrium zone 1.
When S isijWhen the area is contained in a sand body, the theoretical control area can be used as the actual control area of the current oil layer vertical well. However, due to the geological conditions of oil wells, water wells, faults, abandoned river channels or strong edge water on sand bodies, the actual control area of the vertical well may not be equal to the theoretical control area of the vertical well. When the following two conditions exist, the actual control area of the vertical well is determined as follows.
a. When the infill is in a regular pattern
The actual control area of the regular well pattern encryption well such as the inverse nine-point method, the five-point method, the seven-point method and the like can be determined through the range of the shunt line control. Each oil well equally divides the area of the whole water injection unit, and if the area of the water injection unit is 1, the actual control area of the encryption well is 1/n.
b. When infill wells are in an irregular pattern
The irregular well pattern of the complex fault block oil reservoir is encrypted, and the well position, the well type and the boundary attribute of the existing oil well need to be considered. When a water injection well exists, determining the actual control area of the encryption well according to the relation between the encryption well and the oil production well, the water injection well and the boundary, and calculating the residual recoverable reserve in the actual control area; when no water injection well exists, the injection-production relation needs to be analyzed and determined according to the sand body area and the reserve, otherwise, the ultimate recovery factor of the production in failure is taken. The calculation shows that when the distance L between the encryption well and the oil well is 2r, the actual control area of the encryption well is larger than that when the distance is less than 2r, and when L is 2r, the actual control area of the encryption well is the largest. The area of a circle which is formed by taking the grid penetrated by the encryption well as the center of the circle and taking 2r as the radius can be selected as the actual control area of the encryption well.
(2) Inclined shaft
The inclined shaft is one of well types common in oil field sites, and is widely applied in complex fault block oil reservoir exploitation. The reason is that the stacking relation of the single sand bodies of the complex fault block oil reservoir on the section is irregular, the residual movable oil reserves controlled by the inclined shaft are possibly larger, and the inclined shaft can be drilled with better economic benefit.
The method for optimizing the well position of the inclined well is similar to that of a vertical well in the aspects of establishing the corresponding relation between the three-dimensional geological model grid relative coordinates and the geodetic coordinates, calculating the residual movable oil reserves and the like, and the biggest difference from the vertical well is that a well inclination angle and a well inclination azimuth angle exist. In general, the well trajectory of a deviated well generally has different well inclination angles and well inclination azimuths at different depths in the oil reservoir, so that the actual well trajectory of the deviated well is a three-dimensional space curve and not a straight line, as shown in fig. 3 (a). If well position optimization is carried out according to an actual inclined well track model on site, the calculation amount is very large, and the optimization is difficult to realize.
The well type is changed into the inclined well on the basis of the optimized vertical well, the well positions of the inclined well penetrating through different oil layers are determined according to the inclination of the inclined well, and the actual control area range of the inclined well for the different oil layers and the residual recoverable reserves controlled by the inclined well for the single well are calculated. And optimizing the well position of the inclined well by comparing the residual recoverable reserves controlled by the inclined wells with different inclination angles.
Therefore, the following assumptions are made for the inclined well trajectory in the application: firstly, from the ground to the top surface of an oil layer, a well track is straight, namely straight line, and the inclination is started from the top surface of the first oil layer in contact; the angle alpha of the well inclination angle of the inclined well is a constant, and as shown in figure 3(b), the value of the well inclination angle alpha in the text does not exceed 30 degrees; and thirdly, fixing the values of the well deviation azimuth angles, taking 8 inclination directions for the well track of one well position, performing trial calculation on the corresponding well deviation azimuth angles according to the 8 inclination directions, calculating the actual control areas of the inclined wells in the 8 inclination directions, and selecting the direction with the largest residual movable oil storage amount in the 8 directions as the well track direction of the well position. The deviated well trajectory obtained according to the above assumptions is a straight line. Wherein the 8 tilt directions include east, south, west, north, southeast, northeast, southwest, and northwest.
During specific implementation, the inclination direction of the inclined well is changed, the grid of each oil layer which is penetrated by the inclined well in each inclination direction is determined, the grid penetrated by each oil layer inclined well is taken as the center of a circle, the control radius of a single well is taken as the radius, and the area of the circle formed by the center of the circle and the radius is taken as the current inclination direction and the theoretical control area of the current oil layer inclined well; . Similar to the determination method of the actual control area of the vertical well, if a sand body exists in the theoretical control area, the theoretical control area is used as the actual control area of the current oil layer inclined well in the current inclined direction; and if a plurality of sand bodies exist in the theoretical control area, taking the area of the sand body where the inclined shaft is located as the actual control area of the current oil layer inclined shaft in the current inclined direction.
(3) Horizontal well
Because the single well has large controlled reserve, high yield and better economic benefit, the horizontal well is more and more widely applied in recent years, and the drilling of the horizontal well also becomes one of the important measures for improving the recovery ratio in the secondary development of the old oil field, and is widely applied in China.
The horizontal well position optimization method is similar to a vertical well in the aspects of establishing the corresponding relation between the three-dimensional geological model grid relative coordinates and the geodetic coordinates, calculating the residual movable oil reserves and the like, and is mainly different from the vertical well in the aspect that the actual control area is different from that of the vertical well.
The actual control area of the horizontal well is calculated on the basis of the following premise: firstly, a horizontal section of a horizontal well only passes through one oil layer; and secondly, the azimuth angle value of the horizontal section is fixed, the horizontal section of one well position takes 8 extending directions, the horizontal section of each well position is calculated by trial according to the 8 extending directions, and the extending direction of the horizontal section with the largest residual movable oil storage amount is selected as the well track direction of the well point. Wherein, the extending direction includes east, south, west, north, southeast, northeast, southwest and northwest.
When the remaining movable oil reserve of the horizontal well is calculated, an oil layer with a larger remaining movable oil reserve is preferably selected as an oil layer for deploying the horizontal section of the horizontal well.
When the residual movable oil reserves of single-well control are calculated in a single sand layer, and when no interference exists, namely if a lithologic boundary and a strong edge water boundary do not exist in the theoretical control area of the horizontal well, and the vertical distance between the theoretical control area and the horizontal section of the horizontal well is greater than the distance between the control radius and is less than or equal to 2 times of the control radius, no oil well or water well exists, the theoretical control area of the horizontal well is determined as the actual control area. The theoretical control area of the horizontal well is the area of two semicircles formed by taking the preset control radius of a single well as the radius and taking the starting point and the ending point of the horizontal section of the horizontal well as the circle center, and is added with the area of a rectangle formed by taking the length of the horizontal section of the horizontal well as the width of 2 times of the control radius, such as the area of two semicircles and the area of a rectangle in fig. 4(a), wherein r is the control radius.
Due to the existence of faults, edge water, abandoned river channels, lithologic boundaries and other interferences, the actual control surface of the horizontal well is discussed in several cases:
if a lithologic boundary, a fault or a waste river channel exists in the theoretical control area of the horizontal well, and the theoretical control area of the horizontal well is divided into two parts by the lithologic boundary, the fault or the waste river channel, the area of the part where the horizontal well is located in the two parts is determined as the actual control area of the horizontal well, and as shown in fig. 4(b), the actual control area of the horizontal well is the area of the polygonal CDFE.
And secondly, if the vertical distance between the oil well and the horizontal well is larger than the control radius and is less than or equal to 2 times of the control radius, namely when the distance L between the horizontal well and the oil well meets the condition that r < L < > is 2r, the oil well can influence the actual control area of the horizontal well, taking half of a connecting line of midpoints of horizontal sections of the oil well and the horizontal well as the midpoint of the connecting line as the radius, making a circle by taking the midpoint of the connecting line as the center of the circle, overlapping the circle with the theoretical control area, and determining the difference value between the theoretical control area and the area of the overlapped part as the actual control area of the horizontal well. As shown in fig. 4(C), point C is taken as the midpoint of the connecting line between the horizontal section midpoint O of the horizontal well and the well point F of the oil well, CF is taken as the radius to make a circle, and the intersection rectangle is at point D, E, so that the actual control area of the horizontal well is the difference between the areas of the two semicircles and the rectangle and the bow-shaped area CDE.
And thirdly, if the vertical distance between the water well and the horizontal well is larger than the control radius and is less than or equal to 2 times of the control radius, namely when the distance L between the horizontal well and the existing water well meets r < L < -2 r, taking the water well as an end point to two semicircles of the theoretical control area as tangent lines, forming a polygon by the two tangent lines and the boundary of the theoretical control area of the horizontal well, determining the sum of the theoretical control area and the polygon area as the actual control area of the horizontal well, and taking the sum of the areas of the two semicircles and the rectangle and the area of the polygon CDF as the actual control area as shown in figure 4(d), wherein the CD and the CF are tangent lines from the well point F of the water well to the semicircles.
If a strong edge water boundary parallel to the horizontal section of the horizontal well exists in the theoretical control area of the horizontal well, the theoretical control area is divided into two parts by the strong edge water boundary, according to the seepage mechanics mirror image reflection principle, the sum of the length of the horizontal section of the horizontal well and the length of 2 times of the control radius is taken as the length, the vertical distance between the strong edge water boundary and the horizontal section of the horizontal well is a wide rectangular shape, the half of the control area which the strong edge water boundary does not penetrate and the rectangular shape form a polygon, and the area of the formed polygon is determined as the actual control area of the horizontal well. As shown in fig. 4(e), the area of the polygon CEGHFD is determined as the actual control area of the horizontal well.
And calculating the control area of the horizontal well on a certain single sand layer, and then calculating the residual exploitation reserve controlled by the horizontal well. When the well position of the horizontal well is optimized, the direction with the largest residual movable oil reserve is taken as the extension direction of the horizontal section of the horizontal well, a series of parameters corresponding to the better horizontal well are obtained along with the movement of the well point, the analogy is repeated, the well position of the horizontal well is optimized on other single sand layers, all the obtained well positions of all the single sand layers suitable for drilling the horizontal well are compared, and the well position with the largest residual movable oil reserve is preferably selected as the well position of the horizontal well.
After determining the actual control area of the individual well of each well type in each oil layer, the remaining movable oil reserves N of the individual well in the individual oil layer are calculated according to the following formular(i,j)：
Figure BDA0002338706480000111
Wherein s iso(xi,yj,zk) Represents the remaining oil saturation at grid (i, j, k); sorRepresents residual oil saturation, as a constant; rhooRepresents the crude oil density; vφRepresents the pore volume at grid (i, j, k); n represents the grid number contained in the actual control area of the current oil layer; l denotes the l-th mesh among the meshes contained in the current reservoir actual control area.
Then, the remaining mobile oil reserves N of the individual wells in all the oil layers are calculated according to the following formulaGeneral assembly：
Figure BDA0002338706480000112
Wherein p represents the p-th reservoir; layers represent the total number of layers a single well passes through.
It should be noted, however, that there are constraints on the calculations using the above equations. Let (x)a，yb，zc) The coordinates of the well position of a certain oil production well or water injection well existing in the same single sand layer, the existing well position and the encryption well OijThe distance between the encryption well and the old well should be larger than the limit well distance, that is:
(xa-xi)2+(yb-yi)2≥r2
and 103, comparing the total residual movable oil reserves of the single wells drilled in all the non-well point grids, and determining the non-well point grid corresponding to the maximum total residual movable oil reserve as the deployment well position of the encrypted well.
The Sudan 124 area is in a high water content stage at the later development stage, and when a gas injection and water injection pilot test is carried out on part of complex fault block oil reservoirs, the method disclosed by the application is adopted, so that a good application effect is obtained.
In the embodiment of the application, when a single well is deployed in all grids without well points, the residual movable oil reserves of the single well are calculated, and the grid with the maximum residual movable oil reserves is selected as the deployed well position of the single well. By the quantitative calculation method, the purpose of improving the reserve utilization degree can be achieved by definitely determining the mesh in which the single well is deployed and obtaining the largest oil quantity by utilizing the single well; in addition, the method has important significance for recognizing the particularity of the complex fault block oil reservoir and researching the complex fault block oil reservoir, particularly for researching how to economically and reasonably develop the complex fault block oil reservoir which is in the later development stage and has seriously reduced yield.
In addition, on the basis of analyzing the residual movable oil reserves and the extraction degree, the layer with larger residual movable oil reserves and lower utilization degree is selected, the encryption well is arranged on the oil reservoir with the highest residual movable oil reserves, and the oil layer with concentrated residual oil is perforated. The well location of the new infill well is further optimized based on several aspects, mainly:
(1) a movable effective thickness profile of the oil reservoir;
(2) seismic sections;
(3) well types (vertical, horizontal and deviated);
(4) a secondary perforation location;
(5) and analyzing the encrypted well by combining the measure effect.
On the basis of the optimization of the well position of the encrypted oil well, the well position of the water injection well preferably follows the following principle by taking the main force single sand body as a basic unit:
(1) generally at a relatively low location in the formation or adjacent to the well;
(2) generally located in an area with relatively high flooding degree;
(3) the water injection well does not need to be bundled;
(4) a water injection well is not generally arranged near strong edge water, and a water injection well is not generally arranged in an oil reservoir with sufficient bottom water energy;
(5) the injection well ratio is not lower than that before the injection well pattern is optimized.
The main methods adopted during the optimal configuration of the injection-production well pattern comprise the following steps:
(1) arranging the wells which are selected according with the limit well spacing condition and the single-well control limit residual movable oil reserves according to the size of the single-well control movable oil reserves, and selecting the most favorable encryption well position;
(2) and (4) preferably selecting the well position of the water injection well on the basis of the main force single sand body and the main force encryption well position. The injection-production well ratio of the single sand bodies with different scales reaches a certain numerical value, and when the numerical value cannot be met, a transfer injection well is preferably selected from an encryption well and the existing oil production well, so that the maximum storage capacity of the residual movable oil after transfer injection is ensured;
(3) and determining the well pattern form according to the type, the structural form and the oil layer distribution characteristics of the oil reservoir. Aiming at the geological characteristics of the complex fault block oil reservoir, the well pattern is preferably a triangular well pattern. This is because: a. the triangular well pattern well rows are distributed in a staggered mode, so that the triangular well pattern well rows are suitable for irregular complex fault oil reservoirs, and are also beneficial to implementing small faults and mastering the distribution of lens sand bodies; b. the complex fault block oil reservoir is broken, irregular dotted area water injection is adopted, a more complete injection and production system is more easily formed by a triangular well pattern, and the water injection sweep coefficient is improved.
In addition, according to the research findings at home and abroad, the design principle of a plurality of complex fault block oil reservoir well patterns is based on the description result of a single sand body and the distribution characteristics of residual oil, on the basis of the existing injection and production wells, the well patterns are reasonably arranged according to the form of the single sand body and a residual oil enrichment area, the uneven well distribution mode is adopted, the water flooding wave and the volume are furthest expanded, the foundation is laid for the deep profile control and flooding of a main force block, and the following well pattern modes are designed by taking a certain foreign oil field as an example:
(1) one-injection one-production well net
As shown in fig. 5(a), for strip-shaped or triangular sand bodies with small oil-containing area and low reserves, a new well can not be drilled and a well pattern can not be perfected for the sand bodies, and the new well of the main sand body can be used for realizing the purpose of realizing directional well, so that a well pattern form of one injection and one production is formed, the sand bodies are used to the maximum extent, and the residual oil of the sand bodies is fully excavated.
(2) Irregular punctiform water injection well pattern
For irregular sand bodies with slightly large oil-containing area and slightly high reserve, an irregular punctiform water injection well pattern is formed by combining with other new sand body wells and adopting measures of multiple target points and the like, one injection and two extraction or one injection and multiple extraction are carried out, the injected water wave and volume are expanded as much as possible, and residual oil at the edge of the sand body is excavated, as shown in fig. 5 (b).
(3) Row-shaped injection-production well pattern
For some fault-controlled forward-constructed sand bodies, production wells are concentrated at the high part of the fault edge, water injection wells are arranged at the low part to form a drainage water injection well pattern, and residual oil at the root of the fault is deeply developed by combining with natural edge water propulsion, as shown in fig. 5 (c).
(4) Broken block internal cutting injection-production well pattern
For a large block of integrally-installed sand body with serious middle flooding, low oil saturation and high oil saturation at the edge, water is injected into the sand body and oil is extracted from the edge by combining the current situation of the existing oil-water well, an internal cutting injection-production well network is formed, the original water injection benefiting direction is changed, and residual oil in a digging well is excavated to the maximum extent, as shown in fig. 5 (d).
(5) Exhausted single well point mining
For some sand bodies with small oil-containing areas and only controlled by a single well point, an injection and production well pattern cannot be formed, and only conventional hole-filling measures and other measures can be adopted to carry out exhaustive mining by utilizing natural energy, as shown in fig. 5 (e).
The embodiment of the application also provides a device for determining the well position of the infill well in the exploitation region of the complex fault block oil reservoir, as shown in fig. 6, the device 600 comprises an obtaining module 601 and a determining module 602.
The obtaining module 601 is configured to obtain a three-dimensional geological model established according to the geological condition of the complex fault block oil reservoir exploitation area, where the three-dimensional geological model is formed by grids, and the grids are associated with the well arrangement condition, and the residual oil saturation and pore volume of the complex fault block oil reservoir exploitation area are corresponding to the grid.
A determining module 602, configured to determine, for each grid without a well point in the three-dimensional geological model acquired by the acquiring module 601, according to a geological condition, the number of oil layers encountered by a single well during drilling in the grid without a well point, and an actual control area of the single well in each oil layer; determining the remaining movable oil reserves of the single well in each oil layer according to the number of grids contained in the actual control area of the single well in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserve of the single well;
the determining module 602 is further configured to compare the total remaining movable oil reserves of the single wells drilled in all the non-well point grids, and determine the non-well point grid corresponding to the maximum total remaining movable oil reserve as the deployed well position of the infill well.
In an implementation manner of the embodiment of the present application, the apparatus 600 further includes:
and the model establishing module 603 is used for establishing a three-dimensional geological model according to the geological condition of the complex fault block oil reservoir exploitation area and the well arrangement condition.
The obtaining module 601 is further configured to obtain reservoir numerical simulation data, where the reservoir numerical simulation data includes well arrangement conditions, remaining oil saturation, and pore volume at different relative coordinates;
the model establishing module 603 is further configured to establish a correspondence between geodetic coordinates of the three-dimensional geological model grid and relative coordinates of the numerical reservoir simulation data acquired by the acquiring module.
The model building module 603 is further configured to associate each grid with the well placement, remaining oil saturation, and pore volume at the corresponding relative coordinates.
In an implementation manner of the embodiment of the present application, the model establishing module 603 is configured to:
acquiring a proportional relation between a set three-dimensional geological model and a complex fault block oil reservoir exploitation region;
determining the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate according to the proportional relation;
and determining the relative coordinates corresponding to the abscissa and the ordinate of the three-dimensional geological model grid through interpolation calculation according to the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate.
In an implementation manner of the embodiment of the present application, the determining module 602 is configured to:
aiming at each non-well point grid in the three-dimensional geological model, respectively determining the number of oil layers encountered by the single wells of different types when the single wells of different types are drilled in the non-well point grid and the actual control areas of the single wells of different types in each oil layer according to geological conditions; determining the remaining movable oil reserves of the single wells of different types in each oil layer according to the number of grids contained in the actual control area of the single wells of different types in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserves of the single wells of different types; the single well type comprises a vertical well, an inclined well and a horizontal well;
and comparing the total residual movable oil production of the single wells with different types drilled in all the non-well point grids, and respectively determining the type of the single well corresponding to the maximum total residual movable oil production and the non-well point grids as the deployment type and the deployment position of the encrypted well.
In an implementation manner of the embodiment of the present application, the determining module 602 is configured to:
aiming at each oil layer penetrated by the vertical well, calculating the area of a circle which takes the control radius as the radius and takes the grid where the vertical well is located as the center of the circle according to the preset control radius of the single well, and taking the area as the theoretical control area of the vertical well;
judging whether a plurality of sand bodies exist in the theoretical control area of the vertical well;
if a sand body exists in the theoretical control area, taking the theoretical control area as the actual control area of the current oil layer vertical well;
and if a plurality of sand bodies exist in the theoretical control area, taking the area of the sand body where the vertical well is located as the actual control area of the current oil layer vertical well.
In an implementation manner of the embodiment of the present application, the determining module 602 is configured to:
changing the well deviation azimuth angle of the inclined well, and determining the grid of each oil layer penetrated by the inclined well in each well deviation azimuth angle direction, wherein the well deviation azimuth angle directions comprise east, south, west, north, south-east, north-east, south-west and north-west;
aiming at each oil layer penetrated by the inclined shaft, calculating the area of a circle which takes the control radius as the radius and takes a grid penetrated by the inclined shaft as the center of the circle in the current inclined direction as the theoretical control area of the inclined shaft according to the preset control radius of the single shaft;
if a sand body exists in the theoretical control area, taking the theoretical control area as the actual control area of the current oil layer inclined well in the current inclined direction;
and if a plurality of sand bodies exist in the theoretical control area, taking the area of the sand body where the inclined shaft is located as the actual control area of the current oil layer inclined shaft in the current inclined direction.
In an implementation manner of the embodiment of the present application, the determining module 602 is configured to:
selecting an oil layer with the most residual movable oil reserves as a deployed oil layer of the horizontal well, changing the extension direction of the horizontal well in the deployed oil layer, and determining grids passed by the horizontal well in each extension direction of each oil layer; wherein, the extending direction comprises east, south, west, north, southeast, northeast, southwest and northwest;
aiming at each oil layer where a horizontal well is located, the area of two semicircles formed by taking the preset control radius of the single well as the radius and taking the starting point and the ending point of the horizontal section of the horizontal well as the circle center is added with the area of a rectangle formed by taking the length of the horizontal section of the horizontal well as the width of 2 times of the control radius to obtain the theoretical control area of the horizontal well;
if a lithologic boundary exists in the theoretical control area of the horizontal well, the lithologic boundary divides the theoretical control area of the horizontal well into two parts, and the area of the part where the horizontal well is located in the two parts is determined as the actual control area of the horizontal well;
if the vertical distance between the oil well and the horizontal well is larger than the control radius and is smaller than or equal to 2 times of the control radius, taking a half of a connecting line of midpoints of horizontal sections of the oil well and the horizontal well as the center of the connecting line as a radius, making a circle by taking the midpoint of the connecting line as the center of the circle, overlapping the made circle with the theoretical control area, and determining the difference value between the theoretical control area and the area of the overlapped part as the actual control area of the horizontal well;
if the vertical distance between the water well and the horizontal well is larger than the control radius and smaller than or equal to 2 times of the control radius, taking the water well as an end point to make tangent lines to two semicircles of the theoretical control area, forming a polygon on the boundary of the two tangent lines and the theoretical control area of the horizontal well, and determining the sum of the theoretical control area and the polygon area as the actual control area of the horizontal well;
if a strong edge water boundary parallel to the horizontal section of the horizontal well exists in the theoretical control area of the horizontal well, the theoretical control area is divided into two parts by the strong edge water boundary, the sum of the length of the horizontal section of the horizontal well and the length of 2 times of the control radius is taken as the length, the vertical distance between the strong edge water boundary and the horizontal section of the horizontal well is taken as the width to form a rectangle, the half of the control area which is not penetrated by the strong edge water boundary and the rectangle form a polygon, and the area of the formed polygon is determined as the actual control area of the horizontal well;
and if no lithologic boundary and strong edge water boundary exist in the theoretical control area of the horizontal well, no oil well or water well exists in the distance that the vertical distance between the horizontal section of the horizontal well and the control radius is less than or equal to 2 times of the control radius, determining the theoretical control area of the horizontal well as the actual control area.
In an implementation manner of the embodiment of the present application, the determining module 602 is configured to:
according to
Figure BDA0002338706480000161
Calculating the remaining mobile oil reserves N of the individual wells in the individual oil reservoirr(i,j)；
Wherein s iso(xi,yj,zk) Represents the remaining oil saturation at grid (i, j, k); sorRepresents residual oil saturation, as a constant; rhooRepresents the crude oil density; vφRepresents the pore volume at grid (i, j, k); n represents the grid number contained in the actual control area of the current oil layer; l denotes the l-th mesh among the meshes contained in the current reservoir actual control area.
In the embodiment of the application, when a single well is deployed in all grids without well points, the residual movable oil reserves of the single well are calculated, and the grid with the maximum residual movable oil reserves is selected as the deployed well position of the single well. By the quantitative calculation method, the purpose of improving the reserve utilization degree can be achieved by definitely determining the mesh in which the single well is deployed and obtaining the largest oil quantity by utilizing the single well; in addition, the method has important significance for recognizing the particularity of the complex fault block oil reservoir and researching the complex fault block oil reservoir, particularly for researching how to economically and reasonably develop the complex fault block oil reservoir which is in the later development stage and has seriously reduced yield.
The embodiment of the present application further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, any one of the methods in step 101 to step 103 is implemented.
An embodiment of the present application further provides a computer-readable storage medium, where a computer program for executing any one of the methods in step 101 to step 103 is stored in the computer-readable storage medium.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are further described in detail for the purpose of illustrating the invention, and it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (12)

1. A well position determination method for a ciphered well in a complex fault block oil reservoir exploitation region is characterized by comprising the following steps:
acquiring a three-dimensional geological model established according to the geological condition of the complex fault block oil reservoir exploitation area, wherein the three-dimensional geological model is composed of grids, the grids are associated with well arrangement conditions, and the residual oil saturation and pore volume of the corresponding complex fault block oil reservoir exploitation area are obtained;
aiming at each non-well point grid in the three-dimensional geological model, respectively determining the number of oil layers encountered by the single wells of different types when the single wells of different types are drilled in the non-well point grid and the actual control areas of the single wells of different types in each oil layer according to geological conditions; determining the remaining movable oil reserves of the single wells of different types in each oil layer according to the number of grids contained in the actual control area of the single wells of different types in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserves of the single wells of different types; the single well type comprises a vertical well, an inclined well and a horizontal well; the angle of the well inclination angle of the inclined well does not exceed 30 degrees;
comparing the total residual movable oil production of the single wells with different types drilled in all the non-well point grids, and respectively determining the type of the single well corresponding to the maximum total residual movable oil production and the non-well point grids as the deployment type and the deployment well position of the encrypted well;
wherein, when determining to drill a horizontal well in the grid without well points according to the geological conditions, the actual control area of the horizontal well in each oil layer comprises:
selecting an oil layer with the most residual movable oil reserves as a deployed oil layer of the horizontal well, changing the extension direction of the horizontal well in the deployed oil layer, and determining grids passed by the horizontal well in each extension direction of each oil layer; wherein, the extending direction comprises east, south, west, north, southeast, northeast, southwest and northwest;
aiming at each oil layer where a horizontal well is located, the area of two semicircles formed by taking the preset control radius of the single well as the radius and taking the starting point and the ending point of the horizontal section of the horizontal well as the circle center is added with the area of a rectangle formed by taking the length of the horizontal section of the horizontal well as the width of 2 times of the control radius to obtain the theoretical control area of the horizontal well;
if a lithologic boundary exists in the theoretical control area of the horizontal well, the lithologic boundary divides the theoretical control area of the horizontal well into two parts, and the area of the part where the horizontal well is located in the two parts is determined as the actual control area of the horizontal well;
if the vertical distance between the oil well and the horizontal well is larger than the control radius and is smaller than or equal to 2 times of the control radius, taking a half of a connecting line of midpoints of horizontal sections of the oil well and the horizontal well as the center of the connecting line as a radius, making a circle by taking the midpoint of the connecting line as the center of the circle, overlapping the made circle with the theoretical control area, and determining the difference value between the theoretical control area and the area of the overlapped part as the actual control area of the horizontal well;
if the vertical distance between the water well and the horizontal well is larger than the control radius and smaller than or equal to 2 times of the control radius, taking the water well as an end point to make tangent lines to two semicircles of the theoretical control area, forming a polygon on the boundary of the two tangent lines and the theoretical control area of the horizontal well, and determining the sum of the theoretical control area and the polygon area as the actual control area of the horizontal well;
if a strong edge water boundary parallel to the horizontal section of the horizontal well exists in the theoretical control area of the horizontal well, the theoretical control area is divided into two parts by the strong edge water boundary, the sum of the length of the horizontal section of the horizontal well and the length of 2 times of the control radius is taken as the length, the vertical distance between the strong edge water boundary and the horizontal section of the horizontal well is taken as the width to form a rectangle, the half of the control area which is not penetrated by the strong edge water boundary and the rectangle form a polygon, and the area of the formed polygon is determined as the actual control area of the horizontal well;
and if no lithologic boundary and strong edge water boundary exist in the theoretical control area of the horizontal well, no oil well or water well exists in the distance that the vertical distance between the horizontal section of the horizontal well and the control radius is less than or equal to 2 times of the control radius, determining the theoretical control area of the horizontal well as the actual control area.
2. The method of claim 1, wherein prior to obtaining the three-dimensional geological model based on complex fault reservoir production zone geology and well placement, the method further comprises:
establishing a three-dimensional geological model according to the geological condition and well arrangement condition of the complex fault block oil reservoir exploitation area;
acquiring numerical reservoir simulation data, wherein the numerical reservoir simulation data comprises well arrangement conditions, residual oil saturation and pore volume at different relative coordinates;
establishing a corresponding relation between geodetic coordinates of a three-dimensional geological model grid and relative coordinates of numerical reservoir simulation data;
and associating each grid with the well placement, remaining oil saturation and pore volume at the corresponding relative coordinates.
3. The method of claim 2, wherein establishing a correspondence between geodetic coordinates of the three-dimensional geological model grid and relative coordinates of the reservoir numerical simulation data comprises:
acquiring a proportional relation between a set three-dimensional geological model and a complex fault block oil reservoir exploitation region;
determining the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate according to the proportional relation;
and determining the relative coordinates corresponding to the abscissa and the ordinate of the three-dimensional geological model grid through interpolation calculation according to the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate.
4. The method of claim 1, wherein determining an actual control area of vertical wells in each reservoir when drilling vertical wells in the wellless grid based on geology comprises:
aiming at each oil layer penetrated by the vertical well, calculating the area of a circle which takes the control radius as the radius and takes the grid where the vertical well is located as the center of the circle according to the preset control radius of the single well, and taking the area as the theoretical control area of the vertical well;
judging whether a plurality of sand bodies exist in the theoretical control area of the vertical well;
if a sand body exists in the theoretical control area, taking the theoretical control area as the actual control area of the current oil layer vertical well;
and if a plurality of sand bodies exist in the theoretical control area, taking the area of the sand body where the vertical well is located as the actual control area of the current oil layer vertical well.
5. The method of claim 1, wherein determining an actual control area of the slant wells in each reservoir when drilling the slant wells in the wellless grid based on the geology comprises:
changing the well deviation azimuth angle of the inclined well, and determining the grid of each oil layer penetrated by the inclined well in each well deviation azimuth angle direction, wherein the well deviation azimuth angle directions comprise east, south, west, north, south-east, north-east, south-west and north-west;
aiming at each oil layer penetrated by the inclined shaft, calculating the area of a circle which takes the control radius as the radius and takes a grid penetrated by the inclined shaft as the center of the circle in the current inclined direction as the theoretical control area of the inclined shaft according to the preset control radius of the single shaft;
if a sand body exists in the theoretical control area, taking the theoretical control area as the actual control area of the current oil layer inclined well in the current inclined direction;
and if a plurality of sand bodies exist in the theoretical control area, taking the area of the sand body where the inclined shaft is located as the actual control area of the current oil layer inclined shaft in the current inclined direction.
6. The method of any one of claims 1 to 5, wherein determining the remaining mobile oil reserves for the individual well in each reservoir based on the number of grids contained by the individual well within the actual control area in each reservoir, the remaining oil saturation for each grid, the residual oil saturation, the crude oil density, and the pore volume comprises:
according to
Figure FDA0003333434160000031
Calculating the remaining mobile oil reserves N of the individual wells in the individual oil reservoirr(i,j)；
Wherein s iso(xi,yj,zk) Indicating remaining oil saturation at grid (i, j, k)Degree; sorRepresents residual oil saturation, as a constant; rhooRepresents the crude oil density; vφRepresents the pore volume at grid (i, j, k); n represents the grid number contained in the actual control area of the current oil layer; l denotes the l-th mesh among the meshes contained in the current reservoir actual control area.
7. A well position determination device for a ciphered well in a complex fault block oil deposit exploitation region is characterized by comprising:
the acquisition module is used for acquiring a three-dimensional geological model established according to the geological condition of the complex fault block oil reservoir exploitation area, wherein the three-dimensional geological model is composed of grids which are associated with well arrangement conditions, and the residual oil saturation and pore volume of the corresponding complex fault block oil reservoir exploitation area;
the determining module is used for determining the number of oil layers encountered by the single wells of different types when the single wells of different types are drilled in the well-free grids and the actual control areas of the single wells of different types in each oil layer according to geological conditions aiming at each well-free grid in the three-dimensional geological model; determining the remaining movable oil reserves of the single wells of different types in each oil layer according to the number of grids contained in the actual control area of the single wells of different types in each oil layer, the remaining oil saturation of each grid, the remaining oil saturation, the crude oil density and the pore volume; overlapping the residual movable oil reserves of all oil layers to obtain the total residual movable oil reserves of the single wells of different types; the single well type comprises a vertical well, an inclined well and a horizontal well; the angle of the well inclination angle of the inclined well does not exceed 30 degrees;
the determining module is further used for comparing the total residual movable oil production of the single wells with different types drilled in all the non-well point grids, and determining the type of the single well corresponding to the maximum total residual movable oil production and the non-well point grids as the deployment type and the deployment well position of the encrypted well respectively;
the determining module is used for selecting an oil layer with the most residual movable oil reserves as a deployed oil layer of the horizontal well, changing the extending direction of the horizontal well in the deployed oil layer and determining grids which the horizontal well passes through in each extending direction of each oil layer; wherein, the extending direction comprises east, south, west, north, southeast, northeast, southwest and northwest; aiming at each oil layer where a horizontal well is located, the area of two semicircles formed by taking the preset control radius of the single well as the radius and taking the starting point and the ending point of the horizontal section of the horizontal well as the circle center is added with the area of a rectangle formed by taking the length of the horizontal section of the horizontal well as the width of 2 times of the control radius to obtain the theoretical control area of the horizontal well; if a lithologic boundary exists in the theoretical control area of the horizontal well, the lithologic boundary divides the theoretical control area of the horizontal well into two parts, and the area of the part where the horizontal well is located in the two parts is determined as the actual control area of the horizontal well; if the vertical distance between the oil well and the horizontal well is larger than the control radius and is smaller than or equal to 2 times of the control radius, taking a half of a connecting line of midpoints of horizontal sections of the oil well and the horizontal well as the center of the connecting line as a radius, making a circle by taking the midpoint of the connecting line as the center of the circle, overlapping the made circle with the theoretical control area, and determining the difference value between the theoretical control area and the area of the overlapped part as the actual control area of the horizontal well; if the vertical distance between the water well and the horizontal well is larger than the control radius and smaller than or equal to 2 times of the control radius, taking the water well as an end point to make tangent lines to two semicircles of the theoretical control area, forming a polygon on the boundary of the two tangent lines and the theoretical control area of the horizontal well, and determining the sum of the theoretical control area and the polygon area as the actual control area of the horizontal well; if a strong edge water boundary parallel to the horizontal section of the horizontal well exists in the theoretical control area of the horizontal well, the theoretical control area is divided into two parts by the strong edge water boundary, the sum of the length of the horizontal section of the horizontal well and the length of 2 times of the control radius is taken as the length, the vertical distance between the strong edge water boundary and the horizontal section of the horizontal well is taken as the width to form a rectangle, the half of the control area which is not penetrated by the strong edge water boundary and the rectangle form a polygon, and the area of the formed polygon is determined as the actual control area of the horizontal well; and if no lithologic boundary and strong edge water boundary exist in the theoretical control area of the horizontal well, no oil well or water well exists in the distance that the vertical distance between the horizontal section of the horizontal well and the control radius is less than or equal to 2 times of the control radius, determining the theoretical control area of the horizontal well as the actual control area.
8. The apparatus of claim 7, further comprising:
the model building module is used for building a three-dimensional geological model according to the geological condition of the complex fault block oil reservoir exploitation area and the well arrangement condition;
the acquisition module is further used for acquiring numerical reservoir simulation data, wherein the numerical reservoir simulation data comprises well arrangement conditions, residual oil saturation and pore volume at different relative coordinates;
the model establishing module is also used for establishing a corresponding relation between the geodetic coordinates of the three-dimensional geological model grid and the relative coordinates of the numerical reservoir simulation data acquired by the acquiring module;
and the model establishing module is also used for associating the well arrangement condition, the residual oil saturation and the pore volume of each grid and the corresponding relative coordinate position.
9. The apparatus of claim 8, wherein the model building module is configured to:
acquiring a proportional relation between a set three-dimensional geological model and a complex fault block oil reservoir exploitation region;
determining the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate according to the proportional relation;
and determining the relative coordinates corresponding to the abscissa and the ordinate of the three-dimensional geological model grid through interpolation calculation according to the corresponding relation between the depth coordinate of the three-dimensional geological model grid and the depth coordinate of the relative coordinate.
10. The apparatus of any one of claims 7 to 9, wherein the determining module is configured to:
according to
Figure FDA0003333434160000051
Calculating the remaining mobile oil reserves N of the individual wells in the individual oil reservoirr(i,j)；
Wherein s iso(xi,yj,zk) Represents the remaining oil saturation at grid (i, j, k); sorRepresents residual oil saturation, as a constant; rhooRepresents the crude oil density; vφRepresents the pore volume at grid (i, j, k); n represents the grid number contained in the actual control area of the current oil layer; l denotes the l-th mesh among the meshes contained in the current reservoir actual control area.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 6 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 6.
CN201911367075.2A 2019-12-26 2019-12-26 Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area Active CN111080789B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911367075.2A CN111080789B (en) 2019-12-26 2019-12-26 Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911367075.2A CN111080789B (en) 2019-12-26 2019-12-26 Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area

Publications (2)

Publication Number Publication Date
CN111080789A CN111080789A (en) 2020-04-28
CN111080789B true CN111080789B (en) 2022-04-19

Family

ID=70318187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911367075.2A Active CN111080789B (en) 2019-12-26 2019-12-26 Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area

Country Status (1)

Country Link
CN (1) CN111080789B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111594133B (en) * 2020-07-08 2022-03-11 西南石油大学 Woven well pattern for developing multilayer low-permeability oil and gas reservoir based on multi-horizontal-seam bow-shaped well
CN111594132B (en) * 2020-07-08 2022-03-11 西南石油大学 Woven displacement well pattern for fluid injection development of huge thick or multilayer oil and gas reservoir
CN112943230A (en) * 2021-02-05 2021-06-11 中国石油化工股份有限公司 Residual oil distribution prediction method for common heavy oil reservoir
CN113158275B (en) * 2021-05-08 2023-11-21 北京金阳普泰石油技术股份有限公司 Method and system for modeling pinch-out sand body patterns
CN115434686B (en) * 2022-09-21 2023-08-08 常州大学 Method for inhibiting side bottom water coning of heterogeneous broken block heavy oil reservoir and application
CN117371069B (en) * 2023-12-07 2024-03-08 中国石油大学(华东) Method and system for optimizing filling scheme of single-layer-drive streamline regulator of vertical and inclined well group

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2228433C2 (en) * 2001-04-05 2004-05-10 Открытое акционерное общество "Татнефть" им.В.Д.Шашина Method for oil extraction from watering wells and device realizing said method
CN106837289A (en) * 2016-12-29 2017-06-13 中国石油天然气集团公司 A kind of method and apparatus for determining well location
CN109753671A (en) * 2017-11-07 2019-05-14 中国石油大学(华东) A kind of accurate well location optimization method of oil reservoir based on fish-swarm algorithm

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2228433C2 (en) * 2001-04-05 2004-05-10 Открытое акционерное общество "Татнефть" им.В.Д.Шашина Method for oil extraction from watering wells and device realizing said method
CN106837289A (en) * 2016-12-29 2017-06-13 中国石油天然气集团公司 A kind of method and apparatus for determining well location
CN109753671A (en) * 2017-11-07 2019-05-14 中国石油大学(华东) A kind of accurate well location optimization method of oil reservoir based on fish-swarm algorithm

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
An Improved Spatial Clustering for the Optimization of Infill Well Placement in High-Water-Cut Oilfield;Haijun Fan等;《IEEE Xplore》;20120403;第776-778页 *
基于球体的双分支水平井地质储量计算;赵发寿等;《第十三届宁夏青年科学家论坛石化专题论坛论文集》;20170731;第2节 *
数值模拟剩余油分布规律在文中地区的应用;杨勇等;《西部探矿工程》;20051130(第11期);第88-90页 *
油田勘探开发地质对象三维可视化关键技术研究;申龙斌;《中国博士学位论文全文数据库 基础科学辑》;20120615;第2012年卷(第6期);第4章 *
高含水复杂断块油田加密井井位智能优选方法及其应用;鲍敬伟等;《石油学报》;20170430;第38卷(第4期);第0-3节 *
鲍敬伟等.高含水复杂断块油田加密井井位智能优选方法及其应用.《石油学报》.2017,第38卷(第4期),第444-452、484页. *

Also Published As

Publication number Publication date
CN111080789A (en) 2020-04-28

Similar Documents

Publication Publication Date Title
CN111080789B (en) Method and device for determining well position of encrypted well in complex fault block oil reservoir exploitation area
Tao et al. Current status and geological conditions for the applicability of CBM drilling technologies in China: A review
US10365387B1 (en) Method for secondary exploration of old oil area in fault subsidence basin
CN111750822B (en) Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method
US20160376885A1 (en) Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs
CN107808413A (en) A kind of three-dimensional geological modeling method based on GOCAD
CN104992468A (en) Fracture-cavern type carbonate hydrocarbon reservoir three-dimensional geological modeling method
CN103726481B (en) Repeatedly mined surface deformation predicting method
CN109870719B (en) Well position layout method, device and system for carbonate rock tight thin reservoir
CN104809277A (en) Geological modeling method for horizontal wells in ultra-low permeability tight reservoirs
CN110838175B (en) Geological model building method for gas injection development oil reservoir
CN104747180A (en) Fracture-cavity type oil deposit analyzing method for water-flooding development and application thereof
CN104453834A (en) Injection-production relation optimizing and adjusting method for well group
CN106875471A (en) Coal measures contains or water barrier Visualization Modeling method
CN110244021B (en) Stratum layering method based on anisotropic interpolation
CN109138837A (en) Reservoir geology guidance method
CN103758519A (en) Thick and large mineral deposit stage subsequent filling method mining piecewise optimization design and implementation method
Bo et al. Trend judgment of abandoned channels and fine architecture characterization in meandering river reservoirs: a case study of Neogene Minhuazhen Formation NmIII2 layer in Shijiutuo bulge, Chengning uplift, Bohai Bay Basin, East China
CN113156505B (en) Method for identifying reservoir of sunk lake basin reef beach by progressive constraint of three elements of&#39; substrate structure-ancient landform-seismic facies
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
CN110703334B (en) Method for quickly explaining complex fault block oil field horizon
Ciftci et al. Outcrop-based three-dimensional modeling of the Tensleep sandstone at Alkali Creek, Bighorn Basin, Wyoming
CN106154316A (en) A kind of land restoration method and apparatus
CN112282725B (en) Staggered seam-laying design method for horizontal well
CN109472867B (en) Method for quantifying influence range of drilling position information

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant