CN110593863B - Identification method and identification system for water consumption zone of high water-cut oil reservoir - Google Patents

Identification method and identification system for water consumption zone of high water-cut oil reservoir Download PDF

Info

Publication number
CN110593863B
CN110593863B CN201910870315.4A CN201910870315A CN110593863B CN 110593863 B CN110593863 B CN 110593863B CN 201910870315 A CN201910870315 A CN 201910870315A CN 110593863 B CN110593863 B CN 110593863B
Authority
CN
China
Prior art keywords
water
well
layer
injection
zone
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
CN201910870315.4A
Other languages
Chinese (zh)
Other versions
CN110593863A (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 East China
Original Assignee
China University of Petroleum East China
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 East China filed Critical China University of Petroleum East China
Priority to CN201910870315.4A priority Critical patent/CN110593863B/en
Publication of CN110593863A publication Critical patent/CN110593863A/en
Application granted granted Critical
Publication of CN110593863B publication Critical patent/CN110593863B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to the technical field of oil and gas field development, and discloses a method and a system for identifying a water consumption layer zone of a high water-containing period oil reservoir. The identification method comprises the following steps: fitting the production dynamics of the oil well and the water well by using an oil reservoir numerical simulator based on geological data of a target area in the oil reservoir in the high water-cut period and production dynamics data of the oil well and the water well in the target area so as to obtain an oil reservoir numerical simulation model; based on the oil reservoir numerical simulation model, calculating an identification coefficient of a water consumption layer zone between each water well and an oil well around the water well; and identifying a development level of the water-consuming zone based on the identification coefficient of the water-consuming zone. The invention can rapidly judge and identify the development level of the water-consuming layer and quantitatively characterize the water-consuming layer, thereby effectively identifying the development direction of the high-water-consuming layer and playing an effective guiding role in the design of the regulation and control scheme in the subsequent oil field development stage.

Description

Identification method and identification system for water consumption zone of high water-cut oil reservoir
Technical Field
The invention relates to the technical field of oil and gas field development, in particular to a method and a system for identifying a water consumption layer zone of a high water-containing period oil reservoir.
Background
Through years of water injection development, most of oil fields in China enter a high/ultra-high water content development stage, and part of oil fields even enter an ultra-high water content later stage (the water content is higher than 95%), so that a high water consumption layer is generally developed. The method has the advantages that a large amount of injected water is recycled in an inefficient and ineffective way, the dynamic heterogeneity of the reservoir is remarkably aggravated, the contradiction between layers is enlarged, the water flooding efficiency of the oil field is further reduced, the implementation of the yield increasing measures is difficult, the development cost is increased, and finally the economic recovery ratio of the oil field is reduced. Therefore, the method effectively identifies the high water consumption zone and adopts a reasonable regulation and control technology to treat the high water consumption zone, and has very important significance for further improving the development effect of the oil reservoir in the ultra-high water content later period.
In view of the foregoing, there have been disclosed some solutions for identifying a high water-consumption zone that approximates a water-flow dominant channel, for example, by measuring the production dynamics of an oilfield production well using a water-cut index profile to detect a channeling channel of an oilfield; or, by utilizing core data and logging data, identifying a water injection dominant channel by establishing an interpretation model of the radius and permeability of the throat of the reservoir, and preferably plugging the particle size of the microspheres.
The technical scheme only can roughly judge whether a certain well group has a high water consumption zone or not, but has poor characterization effect on the high water consumption zone, so the technical scheme cannot effectively guide the design of a regulation and control scheme in the subsequent oil field development stage.
Disclosure of Invention
The invention aims to provide a method and a system for identifying a water-consuming layer of a high-water-content oil reservoir, which can rapidly judge and identify the development level of the water-consuming layer and quantitatively characterize the development level of the water-consuming layer, so that the development direction of the high-water-consuming layer can be effectively identified, and the design of a regulation and control scheme in the subsequent oil field development stage can be effectively guided.
In order to achieve the above object, according to one aspect of the present invention, there is provided a method for identifying a water-consuming zone of a high water-containing reservoir, the method comprising: fitting the production dynamics of the oil well and the water well by using an oil reservoir numerical simulator based on geological data of a target area in the oil reservoir in the high water-cut period and production dynamics data of the oil well and the water well in the target area so as to obtain an oil reservoir numerical simulation model; based on the oil reservoir numerical simulation model, calculating an identification coefficient of a water consumption layer zone between each water well and an oil well around the water well; and identifying a development level of the water-consuming zone based on the identification coefficient of the water-consuming zone.
Preferably, said fitting the production dynamics of the oil well and the water well using a reservoir numerical simulator comprises: establishing a fine geological model of the target area by adopting a geological modeling algorithm based on geological data of the target area in the oil reservoir in the high water-cut period; importing the fine geologic model into the reservoir numerical simulator; and adjusting characteristic parameters in the fine geological model in the oil reservoir numerical simulator by combining production dynamic data of the oil well and the water well in the target area so as to fit the production dynamic data of the oil well and the water well.
Preferably, the calculating the identification coefficient of the water-consuming zone between each well and the wells located around the well comprises: based on the oil reservoir numerical simulation model, obtaining the effective thickness, the water saturation, the daily water absorption and the daily oil production of each horizon of each water well on each injection well pair; calculating the ratio of dimensionless water consumption to economic water consumption of each layer of each water injection well pair based on the daily water absorption capacity and daily oil production capacity and economic water consumption of each layer of each water injection well pair, wherein the economic water consumption is related to the price of crude oil and the cost of water injection; calculating the ratio of the water absorption intensity of each layer of each water well to the water absorption intensity of each water injection well pair based on the effective thickness of each layer of each water injection well pair, the daily water absorption amount and the effective thickness of each water well; calculating the ratio of the water saturation of each layer of each water well to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection well pair; and calculating a water consumption zone identification coefficient for each layer of each well on each injection well pair based on a ratio of dimensionless water consumption to economic water consumption for each layer on each injection well pair, a ratio of water absorption intensity for each layer on each injection well pair to water absorption intensity for each well, and a ratio of water saturation for each layer on each injection well pair to average water saturation, wherein the injection well pair refers to a reservoir between the well and a particular well located around the well within a control area for each well.
Preferably, the calculating the ratio of dimensionless water consumption to economic water consumption on each injection well pair for each layer of each well comprises: by the formula
Figure BDA0002202609930000031
Calculating dimensionless water consumption E of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, Q iko Daily oil production for horizon i on injection and production well pair k; by the formula
Figure BDA0002202609930000032
Calculating economic water consumption E w Wherein P is o Price per unit volume of crude oil, C w Cost per unit water injection amount; by the formula->
Figure BDA0002202609930000033
Calculating the ratio R of dimensionless water consumption to economic water consumption of horizon i on injection well pair k ike
Preferably, the calculating the ratio of the water absorption intensity of each layer of each water well to the water absorption intensity of each water well on each injection well pair comprises: by the formula
Figure BDA0002202609930000034
Calculating the water absorption intensity omega of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, h ik The effective thickness of the horizon i on the injection well pair k is the effective thickness of the horizon i on the injection well pair k; by the formula->
Figure BDA0002202609930000035
Calculating the water absorption intensity omega of each water well, wherein I is the water absorption capacity of each water well, < > and->
Figure BDA0002202609930000036
h is the effective thickness of each water well; by the formula->
Figure BDA0002202609930000037
Calculating the ratio R of the water absorption intensity of the horizon i on the injection and production well pair k to the water absorption intensity of each water well ikd
Preferably, the calculating of the water saturation on each injection well pair for each horizon and the target zoneThe ratio of the average water saturation within includes: by the formula
Figure BDA0002202609930000041
Calculating the water saturation S of horizon i on injection and production well pair k wik Is +.about.the average water saturation>
Figure BDA0002202609930000042
Ratio R of (2) iks
Preferably, the calculating the water loss zone identification coefficients for the horizons on each injection and production well pair comprises: calculating the water consumption layer zone identification coefficient C on the injection well pair k by the following formula ik ,C ik =ln(R ike ×R iks ×R ikd ) Wherein R is ike The ratio of dimensionless water consumption to economic water consumption of the horizon i on the injection well pair k is the ratio of dimensionless water consumption to economic water consumption of the horizon i; r is R ikd The ratio of the water absorption intensity of the horizon i on the injection well pair k to the water absorption intensity of each water well is given; r is as follows iks Is the ratio of the water saturation of horizon i over the injection well pair k to the average water saturation.
Preferably, said identifying the development level of the water-consuming zone comprises: identification coefficient C of the water-consuming belt ik If the identification layer i is larger than or equal to a first preset coefficient, the identification layer i is a first high water consumption layer zone on the injection well pair k; identification coefficient C of the water-consuming belt ik Under the condition that the first preset coefficient is smaller than the first preset coefficient and larger than or equal to the second preset coefficient, the identification horizon i is a second high water consumption zone on the injection well pair k; identification coefficient C of the water-consuming belt ik And under the condition that the identification layer position i is smaller than the second preset coefficient, the identification layer position i is a common water-consuming layer zone on the injection well pair k.
Preferably, the identification method further comprises: and identifying the development direction of the high water consumption zone in the direction of the injection and production well pair k corresponding to the first high water consumption zone and/or the second high water consumption zone, wherein the injection and production well pair refers to a reservoir between each water well and a specific oil well around the water well in the control area of each water well.
Preferably, the identification method further comprises: after performing the step of identifying the development level of the water-consuming zone, performing the following operations: based on the oil reservoir numerical simulation model, obtaining the permeability level difference, the water content, the crude oil viscosity and the water absorption strength in each injection well pair between each water well and each oil well around each water well; and calculating the volume of the water-consuming layer zone with different development levels based on the permeability level difference, the water content, the crude oil viscosity and the water absorption intensity of each water well in each injection and production well pair, wherein the injection and production well pair refers to a reservoir between the water well and a specific oil well around the water well in the control area of each water well.
Preferably, the development level of the water-consuming zone comprises: ordinary water consumption layer area reaches high water consumption layer area, high water consumption layer area includes: the first high water-consuming layer zone and the second high water-consuming layer zone, respectively, the calculating the volumes of the water-consuming layer zones of different development levels comprises: calculating the volume percentage of the common water consumption layer zone, the second high water consumption layer zone and the first high water consumption layer zone in the injection and production well pair by the following three formulas respectively
Figure BDA0002202609930000051
Is->
Figure BDA0002202609930000052
Figure BDA0002202609930000053
Figure BDA0002202609930000054
Figure BDA0002202609930000055
Wherein x is 1 Permeability in pairs for each of the injection and production wellsGrade difference, mu is the viscosity of the crude oil, f w The water content is the water content; omega is the water absorption intensity of each water well; and calculating the volume V of the j-th class water-consuming zone of the horizon i on the injection and production well pair k by the following formula ikj
Figure BDA0002202609930000056
Wherein j=1, 2 or 3, the class 1 water-consuming zone corresponds to the normal water-consuming zone; a class 2 water-consuming zone corresponds to the second high water-consuming zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, A ik Is the plane area of the horizon i which is transverse to the horizontal direction in the direction of the injection well pair k and h ik Is the effective thickness of horizon i over the injection well pair k.
Compared with the prior art, the water consumption layer identification method has the following advantages:
(1) Fitting production dynamics of an oil well and a water well of a target area in the oil reservoir in a high water-cut period so as to obtain an oil reservoir numerical simulation model; then, based on the obtained oil reservoir numerical simulation model, calculating an identification coefficient of a water consumption layer zone between each water well and surrounding oil wells; finally, different development levels are identified according to the identification coefficient of the water-consuming layer, so that the development level of the water-consuming layer can be rapidly identified and quantitatively characterized, the development direction of the high-water-consumption layer can be effectively identified, and the design of a regulation and control scheme in the subsequent oil field development stage can be effectively guided.
(2) The water consumption layer zone of the oil reservoir in the high water-containing period is identified by taking the direction of each injection well of each horizon as a research object, so that the development horizons, directions, levels and volumes of the water consumption layer zone can be rapidly and accurately judged, and the characterization effect of the water consumption layer zone is effectively improved. All parameters required by the identification process are easy to obtain, the calculation process is very simple, the workload of field technicians is greatly reduced, and the operability is high.
(3) The economic water consumption is used as an important index of the water consumption layer with the identification coefficient, that is, the influence of economic factors such as crude oil price and the like on the development of the oil deposit in the high water-containing period is considered. With the change of the oil price, the invention can help oil field companies to dynamically judge the development condition of the high water consumption layer according to market demand conditions, thereby timely adjusting the development strategy and ensuring the best economic benefit.
In a second aspect, the present invention provides a system for identifying a water-consuming zone of a high water-bearing reservoir, the system comprising: the fitting device is used for fitting the production dynamics of the oil well and the water well by using an oil reservoir numerical simulator based on geological data of a target area in the oil reservoir in the high water content period and production dynamics data of the oil well and the water well in the target area so as to obtain an oil reservoir numerical simulation model; the identification coefficient calculation device is used for calculating the identification coefficient of the water consumption layer zone between each water well and the oil well around the water well based on the oil reservoir numerical simulation model; and the identification device is used for identifying the development level of the water-consuming layer belt based on the identification coefficient of the water-consuming layer belt.
The specific implementation details and effects of the identification system of the water consumption zone of the high water-containing period oil reservoir are the same as those of the identification method of the water consumption zone of the high water-containing period oil reservoir, and are not repeated here.
A third aspect of the present invention provides a machine-readable storage medium having instructions stored thereon for causing a machine to perform the above-described method of identifying a water-consuming zone of a high water-bearing reservoir.
Drawings
The accompanying drawings are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate the invention and together with the description serve to explain, without limitation, the invention. In the drawings:
FIG. 1 is a flow chart of a method for identifying water-consuming zones of a high water-bearing reservoir according to one embodiment of the present invention;
FIG. 2 is a flow chart of a process for fitting production dynamics of the oil well and the water well provided by an embodiment of the present invention;
FIG. 3 is a flow chart of a process for calculating the identification coefficient of the water loss zone between each well and the wells located around the well according to one embodiment of the present invention;
FIG. 4 is a flow chart of a process for identifying a water-loss zone according to one embodiment of the present invention;
FIG. 5 is a graph showing the fit of the water content in a target area according to an embodiment of the present invention;
FIG. 6 is a graph showing a fitting result of cumulative oil production within a target zone according to an embodiment of the present invention;
FIG. 7A is an oil saturation distribution of an L3 layer within a target area provided by an embodiment of the present invention;
FIG. 7B is an oil saturation distribution of the L35 layer within the target area provided by an embodiment of the present invention;
FIG. 8A is a schematic illustration of the control area between a well I1 in a target area and four wells P1, P2, P3, P4 located around the well I1 provided by an embodiment of the present invention;
FIG. 8B is a schematic diagram of an injection well pair formed by well I1 and well P2 according to an embodiment of the present invention; and
FIG. 9 is a block diagram of a system for identifying water-consuming zones of a high water-bearing reservoir, according to one embodiment of the present invention.
Description of the reference numerals
10 fitting device 20 identification coefficient calculating device
30 identification device
Detailed Description
The following describes specific embodiments of the present invention in detail with reference to the drawings. It should be understood that the detailed description and specific examples, while indicating and illustrating the invention, are not intended to limit the invention.
The high water-consumption layer of the oil reservoir generally develops in a high water-containing period (the high water-containing period in the oil reservoir engineering can be divided into three stages with different water contents, namely a high water-containing period with the water content of 80% -90%, an extra-high water-containing period with the water content of more than 90%, and an extra-high water-containing later stage with the water content of more than 95%, especially the extra-high water-containing later stage), so that a large amount of injected water is in low-efficiency and invalid circulation, dynamic heterogeneity of the reservoir is obviously aggravated, in-layer and interlayer contradiction is enlarged, further, the water drive efficiency of an oil field is reduced, the development cost is increased, and finally the economic recovery ratio of the oil field is reduced. In order to solve the above problems, it is necessary to efficiently identify a high water-consuming zone. However, the prior art only can roughly judge whether a certain well group has a high water-consuming zone, but cannot confirm the developed horizon, direction, level and volume of the high water-consuming zone, and has poor characterization effect on the high water-consuming zone. Of course, the three stages of the high water content period in the present invention are not limited to the stages corresponding to the three specific water content ranges, and other reasonable dividing methods are possible.
Based on this, the present invention identifies the water-consuming zones by: on the basis of collecting geological data and production dynamic data of a high-water-content oil reservoir, firstly, carrying out fine oil reservoir numerical simulation research based on the geological data of the high-water-content oil reservoir; then, carrying out production history fitting by combining production dynamic data of a water well (namely a water injection well) and an oil well (namely a production well), so that data such as water saturation, daily water absorption, daily oil production and the like of each small layer (namely a horizon) on each injection and production well pair can be obtained; finally, characterizing the water-consuming layer belt based on the acquired data, namely, calculating the identification coefficient of the water-consuming layer belt, and quantitatively evaluating the development degree of the water-consuming layer belt through the identification coefficient; judging the development level and the development direction of the high water consumption zone in the direction of each injection well by combining with a preset recognition criterion; and the volumes of the water-consuming zones of different levels can be calculated according to the simulation formula. Therefore, the identification method of the water-consuming layer zone of the oil reservoir in the high water-containing period can accurately determine the layer position, the direction, the level and the volume of the water-consuming layer zone, so that the characterization effect of the water-consuming layer zone is improved.
Before describing particular embodiments of the present invention, injection and production well pairs referred to below are explained and illustrated. Injection and production well pairs refer to reservoirs between each well and a particular well located around that well within the control area of that well. The control area of a well refers to the sum of the horizontal cross-sectional planar areas between the well and the wells located around the well (the planar area between well I1 and well P2 is shown in fig. 8B as area a). The direction of each injection well pair indicates different water injection directions of the well. For example, the area shown in fig. 8A is the control area of the water well I1.
FIG. 1 is a flow chart of a method for identifying a water-consuming zone of a high water-bearing reservoir provided by the invention. As shown in fig. 1, the identification method may include the steps of: step S101, fitting production dynamics of an oil well and a water well in the target area by using an oil reservoir numerical simulator based on geological data of the target area in the oil reservoir in the high water content period and production dynamics data of the oil well and the water well in the target area so as to obtain an oil reservoir numerical simulation model; step S102, calculating the identification coefficient of a water consumption layer zone between each water well and the oil wells around the water well based on the oil reservoir numerical simulation model; and step S103, identifying the development level of the water-consuming layer belt based on the identification coefficient of the water-consuming layer belt.
Before executing step S101, it is also necessary to collect geological data, production dynamic data, etc. of the high water-cut reservoir. Specifically, the collecting geological data of the target region within the high water reservoir may include: and collecting a structural contour map, fault track data, sand thickness distribution contour map, effective thickness distribution contour map, porosity distribution contour map, permeability distribution contour map, interlayer distribution map and the like of different layers in the target area. The collecting production dynamics of oil and water wells within the target area may include: daily water injection rate (corresponding to daily water absorption rate) of each water well, daily oil production rate and water content of each oil well around each water well, water absorption/liquid production profile data of the oil-water well and the like are collected.
For step S101, as shown in fig. 2, the process of fitting the production dynamics of the oil well and the water well may include the steps of:
step S201, based on geological data of a target area in the high water-cut reservoir, a fine geological model of the target area is built by adopting a geological modeling algorithm.
A fine geologic model of the target region may be established using a geologic modeling algorithm in Petrel software based on the collected geologic data.
Step S202, the fine geologic model is imported into the reservoir numerical simulator.
The fine geologic model of the target zone may be imported into a reservoir numerical simulator Eclipse for simulation to fit production dynamics of the oil and water wells within the target zone in the reservoir numerical simulator Eclipse.
And step S203, adjusting characteristic parameters in the fine geological model in the oil reservoir numerical simulator in combination with the production dynamic data of the oil well and the water well in the target area so as to fit the production dynamic data of the oil well and the water well.
Fitting the production dynamics of the oil well and the water well refers to fitting the production history of the oil well and the water well. The production history fitting is to simulate the development history of the oil field by using the existing oil reservoir parameters (such as permeability, porosity, relative permeability curve, etc.) by means of an oil reservoir numerical simulation method, and compare the calculated development indexes (such as pressure, yield, comprehensive water content, etc.) with the actual development dynamics of the oil reservoir. If the calculation result is inconsistent with the actual dynamic state, the geological knowledge of the oil reservoir is not clear, the model parameters are inconsistent with the actual parameters, proper adjustment is needed, and the calculation is performed again by using the adjusted model parameters until the calculation result is consistent with the actual dynamic state or within an allowable error range.
In particular, in the production history fitting process, the feature parameters involved include: original geological reserves of the oil reservoir, dynamic parameters of the oil reservoir, and dynamic parameters of a single well (referring to each oil well and water well). Dynamic parameters of the reservoir may include: the accumulated oil production, the accumulated water production, the accumulated gas-oil ratio, the comprehensive water content, the average stratum pressure, the accumulated water injection amount and the like are changed along with time. Dynamic parameters of the well may include: the instantaneous oil production, water content, water production, instantaneous gas-oil ratio, cumulative oil production, cumulative water production, bottom hole pressure, etc. over time. The dynamic parameters of the well may include: the instantaneous water injection rate, the accumulated water injection rate, the bottom hole pressure, etc. change with time. In the actual production history fitting process, the original geological reserves of the oil reservoirs are fitted first, then the comprehensive parameters (such as accumulated oil production, accumulated water production, accumulated gas-oil ratio and the like) of the oil reservoirs are fitted, and finally the dynamic parameters of the single well are fitted.
In this embodiment, the production dynamics of the oil well and the water well are fitted mainly by adjusting characteristic parameters such as the relative permeability curve in the fine geologic model, the permeability data of the oil layer between the wells, the effective thickness of the reservoir, and the compression coefficients of the rock and the fluid.
Production history fitting requires that the fitting error of the original geological reserves of the oil reservoir be within 2% and the fitting error of the comprehensive parameters of the oil reservoir be within 5%. For each parameter related to the production dynamics of a single well, the following 4 indexes are often adopted to evaluate the accuracy of history fit: (1) Forward maximum relative error E pmax : when the calculated value is higher than the actual value, the parameter value on the calculated curve is subtracted from the corresponding value on the actual curve, and the ratio of the obtained difference to the actual value is calculated; (2) Negative maximum relative error E nmax : when the calculated value is lower than the actual value, calculating the absolute value of the ratio of the obtained difference to the actual value after subtracting the corresponding value on the actual curve from the parameter value on the curve; (3) Average value of relative error E rave : calculating an average of absolute values of errors between the numerical values and the actual numerical values; (4) Standard deviation of absolute error E as : the square of the sum of squares of the differences between the values and the actual values is calculated. For a parameter involved in the production dynamics of a single well, if the following conditions are met: e (E) pmax <10%,E nmax <10%,E rave <5, and E as <0.1, then the parameter is considered to have been fitted. If 90% of all individual wells meet the above criteria, the resulting fine geologic model (i.e., reservoir numerical simulation model) may be considered accurate, i.e., the fitting accuracy meets the requirements. Of course, the fitting error thresholds in the present invention are not limited to the specific data, and can be reasonably adjusted according to practical situations.
Fig. 5 and fig. 6 are fitting results of the water content and the accumulated oil yield in the target area, respectively, and it can be found that the prediction result of the oil reservoir numerical simulation model is better matched with the actual production dynamics of the oil reservoir, so that the oil reservoir numerical simulation model can better reflect the underground condition.
For step S102, as shown in fig. 3, the process of calculating the identification coefficient of the water-consuming zone between each water well and the oil well located around the water well may include the steps of:
and step S301, based on the oil reservoir numerical simulation model, acquiring the effective thickness, the water saturation, the daily water absorption and the daily oil production of each layer of each water well on each injection and production well pair.
Step S302, calculating the ratio of dimensionless water consumption to economic water consumption of each layer of each water injection well pair based on the daily water absorption capacity and daily oil production capacity and economic water consumption of each layer of each water injection well pair.
Wherein the economic water consumption is related to the price of crude oil and the cost of water injection. Specifically, the formula can be passed
Figure BDA0002202609930000121
Calculating dimensionless water consumption E of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, m 3 /d,Q iko Daily oil production for horizon i on injection and production well pair k, m 3 /d; by the formula->
Figure BDA0002202609930000122
Calculating economic water consumption E w Wherein P is o Price per unit volume of crude oil, yuan/m 3 ,C w Cost per unit water injection quantity, yuan/m 3 The method comprises the steps of carrying out a first treatment on the surface of the By the formula->
Figure BDA0002202609930000123
Calculating the ratio R of dimensionless water consumption to economic water consumption of horizon i on injection well pair k ike
The embodiment provides an important index for identifying a water consumption zone by taking economic water consumption as an important index, not only can reflect the scientific fact that the difference of oil-water relative flow capacity in a high water content period (including an ultra-high water content later period) is greatly increased, but also considers the influence of economic factors such as crude oil price and the like on the development of the oil deposit in the high water content period, and meets the requirement of high-efficiency development of the oil deposit in the high water content period under the condition of low oil price. The index can be changed along with the oil price, so that an oil field company can be helped to dynamically judge the development condition of the high-water-consumption layer according to the market demand condition, the development strategy can be timely adjusted, and the best economic benefit can be ensured.
Step S303, calculating the ratio of the water absorption intensity of each layer of each water injection well pair to the water absorption intensity of each water well based on the effective thickness of each layer of each water injection well pair, the daily water absorption amount and the effective thickness of each water well.
Specifically, the formula can be passed
Figure BDA0002202609930000124
Calculating the water absorption intensity omega of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, m 3 /(d·m),h ik The effective thickness of the horizon i on the injection well pair k is m; by the formula->
Figure BDA0002202609930000125
Calculating the water absorption intensity omega, m of each water well 3 (d.m), wherein I is the water absorption capacity of each well, m 3 /d,/>
Figure BDA0002202609930000126
h is the effective thickness of each well (as shown in fig. 8B), m; by the formula->
Figure BDA0002202609930000127
Calculating the ratio R of the water absorption intensity of the horizon i on the injection and production well pair k to the water absorption intensity of each water well ikd
Step S304, calculating the ratio of the water saturation of each layer of each water well to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection well pair.
Specifically, the formula can be passed
Figure BDA0002202609930000131
Calculating the water saturation S of horizon i on injection and production well pair k wik Is +.about.the average water saturation>
Figure BDA0002202609930000132
Ratio R of (2) iks
Step S305, calculating the water consumption layer zone identification coefficient of each layer of each injection well pair based on the ratio of the dimensionless water consumption to the economic water consumption of each layer on each injection well pair, the ratio of the water absorption intensity of each layer on each injection well pair to the water absorption intensity of each well, and the ratio of the water saturation of each layer on each injection well pair to the average water saturation.
Specifically, the water loss zone identification coefficient C of horizon i on injection and production well pair k can be calculated by the following formula ik ,C ik =ln(R ike ×R iks ×R ikd ) Wherein R is ike The ratio of dimensionless water consumption to economic water consumption of the horizon i on the injection well pair k is the ratio of dimensionless water consumption to economic water consumption of the horizon i; r is R ikd The ratio of the water absorption intensity of the horizon i on the injection well pair k to the water absorption intensity of each water well is given; r is as follows iks Is the ratio of the water saturation of horizon i over the injection well pair k to the average water saturation.
For step S103, the identifying the development level of the water-consuming zone based on the identification coefficient of the water-consuming zone may include: identification coefficient C of the water-consuming belt ik If the identification layer i is larger than or equal to a first preset coefficient, the identification layer i is a first high water consumption layer zone on the injection well pair k; identification coefficient C of the water-consuming belt ik Under the condition that the first preset coefficient is smaller than the first preset coefficient and larger than or equal to the second preset coefficient, identifying the horizon i on the injection well pair kA second high water-consuming layer; identification coefficient C of the water-consuming belt ik And under the condition that the identification layer position i is smaller than the second preset coefficient, the identification layer position i is a common water-consuming layer zone on the injection well pair k, wherein the water-consuming capacities of the first high water-consuming layer zone (such as an extremely high water-consuming layer zone), the second high water-consuming layer zone (common high water-consuming layer zone) and the common water-consuming layer zone are sequentially weakened. For example, the development level of the water-depleted zone may be determined using the following identification criteria: if C ik <0, the horizon i is a common water-consuming zone on the injection well pair k; if 0 is less than or equal to C ik <0.8, the horizon i is a general high water consumption zone on the injection well pair k; if C ik And if the water consumption of the horizon i is more than or equal to 0.8, the horizon i is an extremely high water consumption zone on the injection well pair k.
The identification method may further include: and identifying the development direction of the high water consumption layer zone in the direction of the injection well pair k corresponding to the first high water consumption layer zone and/or the second high water consumption layer zone. For example, the direction in which the pairs k of injection and production wells of a general high water-loss zone and an extremely high water-loss zone are located can be identified as the development direction of the high water-loss zone.
After step 103 is performed, the embodiment of the present invention may also quantitatively calculate the volumes of the water-consuming zones of different development levels. The specific calculation process may include: based on the oil reservoir numerical simulation model, obtaining the permeability level difference, the water content, the crude oil viscosity and the water absorption strength in each injection well pair between each water well and each oil well around each water well; and calculating the volumes of the water consumption zones with different development levels based on the permeability level difference, the water content, the crude oil viscosity and the water absorption strength of each water well in each injection well pair.
Specifically, firstly, the volume percentage of the common water-consuming layer, the second high water-consuming layer and the first high water-consuming layer in the injection well pair can be calculated through the following three formulas respectively
Figure BDA0002202609930000141
Is->
Figure BDA0002202609930000142
Figure BDA0002202609930000143
Figure BDA0002202609930000144
/>
Figure BDA0002202609930000145
Wherein x is 1 For the permeability level difference in each injection well pair, mu is the viscosity of the crude oil, mPa.s, f w The water content is the water content; omega is the water absorption intensity of each water well.
Next, the volume V of the j-th class water-consuming zone of horizon i on the injection and production well pair k can be calculated by the following formula ikj
Figure BDA0002202609930000146
Wherein j=1, 2 or 3, the class 1 water-consuming zone corresponds to the normal water-consuming zone; a class 2 water-consuming zone corresponds to the second high water-consuming zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, A ik For the horizontal cross-sectional plane area (area A shown in FIG. 8B) of horizon i in the direction of injection well pair k, m 2 ,h ik Is the effective thickness of horizon i over the injection well pair k, m.
Specifically, the process of identifying (or characterizing) the water-loss zone is explained and illustrated in detail using the following reservoir example.
Through years of water injection development, a certain oilfield XX unit enters an extra-high water content later mining stage, and a high water consumption layer is generally developed, so that a large amount of injected water is subjected to inefficient and invalid circulation, the water content of an oil well rises rapidly, and the water driving utilization degree is low. The target area is a river positive rhythm oil reservoir, currently, 2 wells (1-2-95 and 1-2-11) and 6 oil wells (1-2-9, 1-2-92, 1-2-94, 1-2N91, 1-2-111 and 1-2X 112) are commonly produced, the comprehensive water content is 95.6%, and the method is a typical representation of the ultra-high water content later-stage integral oil reservoir.
As shown in fig. 4, the process of identifying the water-consuming zone may include the steps of:
and S401, collecting geological data in the target area, production dynamic data and monitoring data of oil wells and water wells in the target area.
Step S402, fitting production dynamics of the oil well and the water well.
Step S403, calculating the identification coefficient of the water consumption layer.
First, for each well, data such as effective thickness, water saturation, daily water absorption, daily oil production and the like of each horizon in the direction of each injection well are read from the oil reservoir numerical simulation model after history fitting, and a data table shown in table 1 is obtained. The data are easy to obtain, the calculation process is simple, the workload of field technicians is greatly reduced, and the operability is high.
TABLE 1 basic data sheet for each layer on each injection well pair
Figure BDA0002202609930000161
As can be seen from Table 1, water wells 1-2-95 correspond to 4 horizons, for water wells 1-2-95, only 1 well 1-2-9 is on horizon L3; corresponding to 3 wells at horizon L21 are 1-2-9, 1-2-92 and 1-2-94, respectively. The 3 injection and production well pairs of wells 1-2-95 on horizon L21 are reservoirs on horizon L21 between them and the 3 wells (1-2-9, 1-2-92 and 1-2-94), respectively. The data of daily oil production, daily water absorption and the like of the injection and production well pair corresponding to the water feeding well 1-2-95 and the oil well 1-2-9 on the horizon L3 are the first row of data in the table 1.
Second, for each well, the ratio of dimensionless water consumption to economic water consumption on each injection and production well pair is calculated for each layer.
First, according to the daily water absorption I of horizon I on injection well pair k in Table 1 ikw And daily oil production Q iko Calculating horizon i in injection and production by using the following methodDimensionless water consumption E on well pair k ik
Figure BDA0002202609930000162
And the calculated data are listed in table 2. For example, the dimensionless water consumption of horizon L3 over 1-2-95 and 1-2-9 injection and production well pairs is calculated to be 326.46, which is filled into Table 2.
Table 2 dimensionless water consumption, water absorption strength and water saturation values in each injection and production well pair
Figure BDA0002202609930000171
Second, assuming that the current crude oil price is 50 dollars/barrel, the dollar RMB exchange rate is 6.7 yuan, cost per unit water injection rate C w Is 8 yuan/m 3 The economic water consumption can be calculated using the following formula:
Figure BDA0002202609930000172
wherein P is o Price per unit volume of crude oil, yuan/m 3 . The economic water consumption is 263.4m 3
Finally, the ratio (simply called water consumption ratio) R of the dimensionless water consumption to the economic water consumption of the horizon i on the injection well pair k can be calculated by using the following formula ike
Figure BDA0002202609930000173
And the calculated data are listed in table 3. For example, the ratio of dimensionless water consumption to economic water consumption of horizon L3 over 1-2-95 and 1-2-9 injection well pairs was calculated to be 1.239 and filled into Table 3.
Table 3 calculation results of Water-loss zone identification coefficients for each injection well pair
Figure BDA0002202609930000174
Thirdly, for each well, the ratio of the water absorption intensity of each layer on each injection well pair to the water absorption intensity of the whole well (i.e. each well) is calculated.
First, from the data in Table 1, the water absorption intensity ω of horizon i over the injection and production well pair k can be calculated using the following formula ik
Figure BDA0002202609930000181
And the calculated data are listed in table 2. For example, horizon L3 has a water absorption strength of 1.123 over 1-2-95 and 1-2-9 injection well pairs and is filled in Table 2.
Second, the water absorption intensity ω of the whole well can be calculated using the following formula:
Figure BDA0002202609930000182
wherein I is the water absorption capacity of the whole well,
Figure BDA0002202609930000183
m 3 /d; h is the total effective thickness corresponding to the whole well, and m.
For example, as can be calculated from the data in Table 1, the daily water absorption of wells 1-2-95 is 63.536m 3 D (note that only the daily water absorption capacity of a certain injection well pair is calculated on the horizon L21), the total effective thickness of the whole well is 8.388m (can be directly obtained from the oil reservoir numerical simulation model), so the water absorption intensity of the whole well of the water wells 1-2-95 is 7.574m 3 /(d.m). The water absorption strength of the whole well for wells 1-2-11 can also be obtained in a similar manner, which is 6.672m 3 /(d.m). The water absorption intensity values for the whole well are listed in table 4.
Finally, the water absorption intensity omega of horizon i on injection and production well pair k can be calculated by using the following formula ik The ratio of the water absorption intensity omega to the whole well (simply called as the water absorption intensity ratio) R ikd
Figure BDA0002202609930000184
And the calculated data are listed in table 3.
Fourth, for each well, the ratio of water saturation to average water saturation for each zone over the respective injection well pair is calculated.
The water saturation of horizon i on injection well pair k can be calculated using the following equation; s is S wik And average water saturation
Figure BDA0002202609930000185
Is referred to as the water saturation ratio R iks :/>
Figure BDA0002202609930000186
And the calculated data are listed in table 3.
Specifically, the S wik And (3) with
Figure BDA0002202609930000187
All can be obtained from table 1. For example, as can be seen from Table 1, the water saturation in the 1-2-95 and 1-2-9 injection well pairs on the L3 small layer is 0.5465 and the average water saturation of the unit is 0.5211, thus R iks Is 1.048, which is filled in table 3.
Fifth, for each water well, calculate the water-loss zone identification coefficient on each injection well pair for each horizon.
The water consumption layer zone identification coefficient C on the injection well pair k can be calculated by using the following formula ik :C ik =ln(R ike ×R iks ×R ikd ). Specifically, the data in table 3 can be substituted into the above formula to calculate the water-consuming layer identification coefficient C of horizon i on injection well pair k ik And the calculated data are listed in table 3.
Step S404, the development level and direction of the water-consuming layer are judged by combining the identification criteria of the water-consuming layer.
If C ik <0, the water-consuming layer zone on the injection well pair k is a common water-consuming layer zone;
if 0 is less than or equal to C ik <0.8, the water-consuming layer zone of the horizon i on the injection well pair k is a general high water-consuming layer zone;
if C ik And if the water consumption layer area of the horizon i on the injection well pair k is more than or equal to 0.8, the water consumption layer area of the horizon i on the injection well pair k is an extremely high water consumption layer area.
And identifying the development direction of the high water-consuming zone in the direction of the injection well pair corresponding to the general high water-consuming zone and the extremely high water-consuming zone.
The identification criteria of the water-consuming belt are determined by the following method: and (3) establishing a large number of conceptual models according to actual geological data of an oil reservoir site and production dynamic characteristics of each water well and each oil well, and carrying out numerical simulation research. And respectively calculating the water consumption layer identification coefficient of each layer on each injection well pair according to the simulation result, and carrying out cluster analysis on the calculation result. The research result shows that the simulated 2516 data can be divided into three major categories through a clustering algorithm, and then the recognition criterion of each category of water consumption layer can be determined according to the clustering analysis result.
The results of the identification of the water-consuming zones in this example are shown in table 3. It can be seen from Table 3 that horizon L21 within the target zone develops an extremely high water-loss zone between wells 1-2-95 and wells 1-2-9, 1-2-94; an extremely high water consumption layer is developed between the L23 layer water feeding well 1-2-95 and the oil well 1-2-92; layer L35 develops extremely high water-loss zones between wells 1-2-11 and wells 1-2-94, 1-2-111 and 1-2X 112. The oil reservoir numerical simulation result shows that the horizon L3 layer is a common water consumption layer between the water wells 1-2-95 and the oil wells 1-2-9, as shown in FIG. 7A; horizon L35 develops extremely high water-loss zones between wells 1-2-11 and wells 1-2-94, 1-2-111 and 1-2X112, as shown in FIG. 7B (direction of arrow), which is consistent with the results set forth in Table 3.
In step S405, the volumes of the water-consuming zones of different levels are calculated.
Firstly, parameters such as permeability level difference, water content, crude oil viscosity, water absorption intensity and the like in each injection well pair are obtained from the oil reservoir numerical simulation model. Table 4 only lists the permeability level differences, water contents, crude oil viscosity, water absorption strength, etc. within each injection and production well pair for extremely high water loss zones.
Table 4 base data for each injection well pair
Figure BDA0002202609930000201
Second, according to the data in Table 4, for Each water well adopts the following formula to calculate the volume percentage of the common water consumption layer in each injection well pair
Figure BDA0002202609930000202
Figure BDA0002202609930000203
For each water well, the following method is adopted to calculate the volume percentage of the general high water consumption layer in each injection well pair
Figure BDA0002202609930000204
Figure BDA0002202609930000205
For each water well, the volume percentage of the extremely high water consumption layer in each injection and production well pair is calculated by adopting the following formula
Figure BDA0002202609930000206
Figure BDA0002202609930000207
The volume percentages of the different levels of water-loss zones in the direction of each injection well in this example are shown in table 5.
TABLE 5 volume percent of different grades of water-consuming bands
Figure BDA0002202609930000211
The calculation formulas of the volume percentages of the water consumption zones with different levels in each injection well pair are determined by adopting the following methods: and (3) establishing a large number of conceptual models according to actual geological data of an oil reservoir site and production dynamic characteristics of each water well and each oil well, and carrying out numerical simulation research. According to the simulation result, respectively calculating the water consumption layer zone identification coefficient of each injection well pair; determining the volumes of the water consumption zones of different levels in each injection and production well pair by using the recognition criteria of the water consumption zones; and then, the sensitivity analysis finds that the volume of the water-consuming layer is mainly influenced by parameters such as the permeability level difference, the water content, the crude oil viscosity, the water absorption strength and the like, so that the volume of the water-consuming layer obtained by simulation is subjected to statistical regression, and a calculation formula of the volume percentage of the water-consuming layer with different levels in the injection and production well pair is obtained.
Finally, calculating the volume V of the j-th water-consuming zone of the horizon i on the injection and production well pair k by adopting the following method ikj
Figure BDA0002202609930000212
Wherein j=1, 2 or 3, the class 1 water-consuming zone corresponds to the normal water-consuming zone; a class 2 water-consuming zone corresponds to the second high water-consuming zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, A ik For the horizontal cross-sectional plane area (area A shown in FIG. 8B) of horizon i in the direction of injection well pair k, m 2 ,h ik Is the effective thickness of horizon i over the injection well pair k, m.
The volumes of the water consumption zones with different levels can provide guidance for the actual profile control and water shutoff design of the oilfield site. According to the identification result, deep profile control is carried out on the water well in the target area in 10 months 2018, and the oil field on-site water absorption profile test data and the monitoring result of the tracer between the wells show that the extremely high water consumption layer is effectively plugged, the water absorption capacity of the high water consumption layer is effectively inhibited, the residual oil in the common water consumption layer is effectively utilized, and the average daily oil yield of the well group is increased by 13.2m 3 And/d, the water content is reduced by 10.4 percent, and the predicted recovery ratio is improved by 2.87 percent.
The practical application effect of the embodiment of the invention in the oilfield field shows that the identification method can effectively identify the high water consumption zone in the high water content period (including the ultra-high water content later period), the actual coincidence rate with the oilfield field is more than 90%, and the volume of the obtained water consumption zone can effectively guide the design of regulation measures.
In the embodiment of the invention, the data such as water saturation, daily water absorption, daily oil production and the like of each layer on each injection and production well pair can be obtained by collecting geological data of the extra-high water-containing later-stage oil reservoir, production dynamic characteristics and monitoring data of each well and each oil well, carrying out fine oil reservoir numerical simulation research and carrying out production history fitting, and basic data is provided for identifying water consumption zones; the development degree of the water-consuming layer zone can be quantitatively evaluated through the calculated identification coefficient of the water-consuming layer zone, and the development level and the development direction of the water-consuming layer zone on each injection and production well pair can be judged by combining the identification criterion of the water-consuming layer. And the volumes of the water-consuming zones of different levels can also be calculated. Therefore, the identification method of the water-consuming layer zone of the oil reservoir in the high water-containing period can accurately determine the layer position, the direction, the level and the volume of the water-consuming layer zone, and improves the characterization effect of the water-consuming layer zone.
In summary, the invention creatively fits the production dynamics of the oil well and the water well of the target area in the oil reservoir in the high water-cut period, thereby obtaining the oil reservoir numerical simulation model; then, based on the obtained oil reservoir numerical simulation model, calculating an identification coefficient of a water consumption layer zone between each water well and the surrounding oil wells; finally, different development levels are identified according to the identification coefficient of the water-consuming layer, so that the development level of the high-water-consuming layer can be rapidly identified and quantitatively characterized, the development direction of the high-water-consuming layer can be effectively identified, and the design of a regulation and control scheme in the subsequent oil field development stage can be effectively guided.
Correspondingly, the invention also provides a system for identifying the water consumption zone of the oil reservoir in the high water-containing period. As shown in fig. 9, the identification system may include: fitting means 10 for fitting production dynamics of the oil well and the water well by using an oil reservoir numerical simulator based on geological data of a target area in the high water-cut period oil reservoir and production dynamics data of the oil well and the water well in the target area to obtain an oil reservoir numerical simulation model; identification coefficient calculation means 20 for calculating an identification coefficient of a water-consuming zone between each well and an oil well located around the well based on the reservoir numerical simulation model; and identification means 30 for identifying the development level of the water-consuming zone based on the identification coefficient of the water-consuming zone.
Preferably, the fitting device 10 may comprise: the fine geologic model building module is used for building a fine geologic model of the target area by adopting a geologic modeling algorithm based on geologic data of the target area in the high-water-cut oil reservoir; and an importing module for importing the fine geologic model into the reservoir numerical simulator; and the fitting module is used for combining the production dynamic data of the oil well and the water well in the target area, and adjusting the characteristic parameters in the fine geological model in the oil reservoir numerical simulator so as to fit the production dynamic data of the oil well and the water well.
Preferably, the identification coefficient calculating means 20 may include: the parameter acquisition module is used for acquiring the effective thickness, the water saturation, the daily water absorption and the daily oil production of each layer of each injection well pair based on the oil reservoir numerical simulation model; the water consumption ratio calculation module is used for calculating the ratio of the dimensionless water consumption to the economic water consumption of each layer of each water injection and production well pair based on the daily water absorption and daily oil production and the economic water consumption of each layer of each water injection and production well pair, wherein the economic water consumption is related to the price of crude oil and the cost of water injection; the water absorption intensity ratio calculating module is used for calculating the ratio of the water absorption intensity of each layer of each water injection well pair to the water absorption intensity of each water injection well based on the effective thickness and daily water absorption amount of each layer of each water injection well pair and the effective thickness of each water injection well; the water saturation ratio calculating module is used for calculating the ratio of the water saturation of each layer of each water well to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection well pair; and an identification coefficient calculation module for each layer based on each layer The ratio of dimensionless water consumption to economic water consumption on each injection well pair, the ratio of the water absorption intensity of each layer on each injection well pair to the water absorption intensity of each well, and the ratio of the water saturation of each layer on each injection well pair to the average water saturation are calculated, and the water consumption layer on each injection well pair is provided with an identification coefficient, wherein the injection well pair refers to a reservoir between a specific well around the well in the control area of each well. Preferably, the recognition device 30 may include: the development level identification module is used for executing the following operations: identification coefficient C of the water-consuming belt ik If the identification layer i is larger than or equal to a first preset coefficient, the identification layer i is a first high water consumption layer zone on the injection well pair k; identification coefficient C of the water-consuming belt ik Under the condition that the first preset coefficient is smaller than the first preset coefficient and larger than or equal to the second preset coefficient, the identification horizon i is a second high water consumption zone on the injection well pair k; identification coefficient C of the water-consuming belt ik And under the condition that the identification layer position i is smaller than the second preset coefficient, the identification layer position i is a common water-consuming layer zone on the injection well pair k.
Preferably, the identification system further comprises: and the development direction identification module is used for identifying the development direction of the high water consumption zone in the direction of the injection well pair k corresponding to the first high water consumption zone and/or the second high water consumption zone, wherein the injection well pair is in the control area of each water well, and a reservoir between the water well and the oil wells positioned around the water well.
Preferably, the identification system may further comprise: the parameter obtaining device is used for obtaining the permeability level difference, the water content, the crude oil viscosity and the water absorption intensity in each injection well pair between each water well and each oil well around each water well based on the oil reservoir numerical simulation model; and the volume calculating device is used for calculating the volumes of the water consumption zones with different development levels based on the permeability level difference, the water content, the crude oil viscosity and the water absorption intensity of each water well in each injection and production well pair, wherein the injection and production well pair refers to a reservoir between each water well and a specific oil well around the water well in the control area of each water well.
The specific implementation details and effects of the identification system of the water-consuming zone of the high water-cut oil reservoir in the embodiment of the present invention are the same as those of the embodiment of the identification method, and are not repeated here.
Correspondingly, the invention also provides a machine-readable storage medium, and the machine-readable storage medium is stored with instructions, and the instructions are used for enabling a machine to execute the identification method of the water-consuming zone of the high water-bearing reservoir.
The machine-readable storage medium includes, but is not limited to, phase-change Memory (abbreviation for phase-change random access Memory, phase Change Random Access Memory, PRAM, also known as RCM/PCRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash Memory (Flash Memory) or other Memory technology, compact disc read only Memory (CD-ROM), digital Versatile Disc (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices, and the like.
The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings, but the present invention is not limited to the specific details of the above embodiments, and various simple modifications can be made to the technical solution of the present invention within the scope of the technical concept of the present invention, and all the simple modifications belong to the protection scope of the present invention.
In addition, the specific features described in the above embodiments may be combined in any suitable manner without contradiction. The various possible combinations of the invention are not described in detail in order to avoid unnecessary repetition.
Moreover, any combination of the various embodiments of the invention can be made without departing from the spirit of the invention, which should also be considered as disclosed herein.

Claims (13)

1. The method for identifying the water consumption zone of the oil reservoir in the high water-containing period is characterized by comprising the following steps of:
fitting the production dynamics of the oil well and the water well by using an oil reservoir numerical simulator based on geological data of a target area in the oil reservoir in the high water-cut period and production dynamics data of the oil well and the water well in the target area so as to obtain an oil reservoir numerical simulation model;
based on the oil reservoir numerical simulation model, calculating an identification coefficient of a water consumption layer zone between each water well and an oil well around the water well; and
identifying the development level of the water-consuming zone based on the identification coefficient of the water-consuming zone,
wherein, calculate the discernment coefficient of the water consumption layer area between every well and the oil well that is located around this well includes:
Based on the oil reservoir numerical simulation model, obtaining the effective thickness, the water saturation, the daily water absorption and the daily oil production of each horizon of each water well on each injection well pair;
calculating the ratio of dimensionless water consumption to economic water consumption of each layer of each water injection well pair based on the daily water absorption capacity and daily oil production capacity and economic water consumption of each layer of each water injection well pair, wherein the economic water consumption is related to the price of crude oil and the cost of water injection;
calculating the ratio of the water absorption intensity of each layer of each water well to the water absorption intensity of each water injection well pair based on the effective thickness of each layer of each water injection well pair, the daily water absorption amount and the effective thickness of each water well;
calculating the ratio of the water saturation of each layer of each water well to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection well pair; and
calculating the water consumption zone identification coefficient of each layer of each injection well pair based on the ratio of dimensionless water consumption to economic water consumption of each layer on each injection well pair, the ratio of water absorption intensity of each layer on each injection well pair to water absorption intensity of each well, and the ratio of water saturation of each layer on each injection well pair to average water saturation,
Wherein the injection well pair refers to a reservoir between the well and a specific well located around the well within a control area of each well,
the calculating of the ratio of dimensionless water consumption to economic water consumption on each injection well pair of each layer of each water well comprises the following steps:
by the formula
Figure FDA0004038651300000021
Calculating dimensionless water consumption E of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, Q iko Daily oil production for horizon i on injection and production well pair k;
by the formula
Figure FDA0004038651300000022
Calculating economic water consumption E w Wherein P is o Price per unit volume of crude oil, C w Cost per unit water injection amount; and
by the formula
Figure FDA0004038651300000023
Calculating the ratio R of dimensionless water consumption to economic water consumption of horizon i on injection well pair k ike
Wherein, based on the effective thickness, daily water absorption of each layer on each injection well pair of each well, and the effective thickness of each well, the ratio of the water absorption intensity of each layer on each injection well pair to the water absorption intensity of each well is calculated to include:
by the formula
Figure FDA0004038651300000024
Calculating the water absorption intensity omega of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, h ik The effective thickness of the horizon i on the injection well pair k is the effective thickness of the horizon i on the injection well pair k;
by the formula
Figure FDA0004038651300000025
Calculating the water absorption intensity omega of each water well, wherein I is the water absorption capacity of each water well, < > and->
Figure FDA0004038651300000026
h is the effective thickness of each water well; and
by the formula
Figure FDA0004038651300000027
Calculating the ratio R of the water absorption intensity of the horizon i on the injection and production well pair k to the water absorption intensity of each water well ikd
The calculating of the water consumption layer zone identification coefficient of each layer of each water well on each injection well pair comprises the following steps:
calculating the water consumption layer zone identification coefficient C on the injection well pair k by the following formula ik ,C ik =ln(R ike ×R iks ×R ikd ),
Wherein R is ike The ratio of dimensionless water consumption to economic water consumption of the horizon i on the injection well pair k is the ratio of dimensionless water consumption to economic water consumption of the horizon i; r is R ikd The ratio of the water absorption intensity of the horizon i on the injection well pair k to the water absorption intensity of each water well is given; r is as follows iks Is the ratio of the water saturation of horizon i over the injection well pair k to the average water saturation.
2. The method of identifying water-consuming zones of a high water-bearing reservoir of claim 1, wherein said fitting production dynamics of the oil well and the water well using a reservoir numerical simulator comprises:
establishing a fine geological model of the target area by adopting a geological modeling algorithm based on geological data of the target area in the oil reservoir in the high water-cut period; and
Importing the fine geologic model into the reservoir numerical simulator;
and adjusting characteristic parameters in the fine geological model in the oil reservoir numerical simulator by combining production dynamic data of the oil well and the water well in the target area so as to fit the production dynamic data of the oil well and the water well.
3. The method of identifying water-consuming zones of a high water-bearing reservoir of claim 1, wherein the calculating of the ratio of the water saturation on each injection well pair to the average water saturation in the target zone for each horizon of each well comprises:
by the formula
Figure FDA0004038651300000031
Calculating the water saturation S of horizon i on injection and production well pair k wik And the average water saturation
Figure FDA0004038651300000041
Ratio R of (2) iks
4. The method of identifying a water-consuming zone of a high water-bearing reservoir of claim 1, wherein the identifying a level of development of the water-consuming zone comprises:
identification coefficient C of the water-consuming belt ik If the identification layer i is larger than or equal to a first preset coefficient, the identification layer i is a first high water consumption layer zone on the injection well pair k;
identification coefficient C of the water-consuming belt ik Under the condition that the first preset coefficient is smaller than the first preset coefficient and larger than or equal to the second preset coefficient, the identification horizon i is a second high water consumption zone on the injection well pair k; and
Identification coefficient C of the water-consuming belt ik Under the condition that the identification level i is smaller than the second preset coefficient, identifying the horizon i in the injection well pair kThe upper part is a common water consumption layer belt.
5. The method of identifying a water-consuming zone of a high water-bearing reservoir of claim 4, further comprising:
identifying the development direction of the high water-consuming layer zone in the direction of the injection well pair k corresponding to the first high water-consuming layer zone and/or the second high water-consuming layer zone,
the injection and production well pair refers to a reservoir between each well and a specific oil well around the well in the control area of the well.
6. The method of identifying a water-consuming zone of a high water-bearing reservoir of claim 1, further comprising:
after performing the step of identifying the development level of the water-consuming zone, performing the following operations:
based on the oil reservoir numerical simulation model, obtaining the permeability level difference, the water content, the crude oil viscosity and the water absorption strength in each injection well pair between each water well and each oil well around each water well; and
calculating the volumes of the water-consuming zones with different development levels based on the permeability level difference, the water content, the crude oil viscosity and the water absorption intensity of each water well in each injection and production well pair,
The injection and production well pair refers to a reservoir between each well and a specific oil well around the well in the control area of the well.
7. The method of claim 6, wherein the development level of the water-consuming zone comprises: ordinary water consumption area and high water consumption area, high water consumption area includes: a first high water-consuming layer belt and a second high water-consuming layer belt,
correspondingly, the calculating the volume of the water-consuming zone at different developmental levels comprises:
respectively by the following threeCalculating the volume percentage of the common water consumption layer zone, the second high water consumption layer zone and the first high water consumption layer zone in the injection and production well pair according to a formula
Figure FDA0004038651300000051
Is->
Figure FDA0004038651300000052
Figure FDA0004038651300000053
Figure FDA0004038651300000054
Figure FDA0004038651300000055
Wherein x is 1 For the permeability level difference in each injection well pair, mu is the viscosity of the crude oil, f w The water content is the water content; omega is the water absorption intensity of each water well; and
calculating the volume V of the j-th water-consuming zone on the injection well pair k by the following formula ikj
Figure FDA0004038651300000056
Wherein j=1, 2 or 3, the class 1 water-consuming zone corresponds to the normal water-consuming zone; a class 2 water-consuming zone corresponds to the second high water-consuming zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, A ik Is the plane area of the horizon i which is transverse to the horizontal direction in the direction of the injection well pair k and h ik Is the effective thickness of horizon i over the injection well pair k.
8. A system for identifying a water-consuming zone of a high water-bearing reservoir, the system comprising:
the fitting device is used for fitting the production dynamics of the oil well and the water well by using an oil reservoir numerical simulator based on geological data of a target area in the oil reservoir in the high water content period and production dynamics data of the oil well and the water well in the target area so as to obtain an oil reservoir numerical simulation model;
the identification coefficient calculation device is used for calculating the identification coefficient of the water consumption layer zone between each water well and the oil well around the water well based on the oil reservoir numerical simulation model; and
identification means for identifying a development level of the water-consuming zone based on the identification coefficient of the water-consuming zone,
wherein the identification coefficient calculating means includes:
the parameter acquisition module is used for acquiring the effective thickness, the water saturation, the daily water absorption and the daily oil production of each layer of each injection well pair based on the oil reservoir numerical simulation model;
the water consumption ratio calculation module is used for calculating the ratio of the dimensionless water consumption to the economic water consumption of each layer of each water injection and production well pair based on the daily water absorption and daily oil production and the economic water consumption of each layer of each water injection and production well pair, wherein the economic water consumption is related to the price of crude oil and the cost of water injection;
The water absorption intensity ratio calculating module is used for calculating the ratio of the water absorption intensity of each layer of each water injection well pair to the water absorption intensity of each water injection well based on the effective thickness and daily water absorption amount of each layer of each water injection well pair and the effective thickness of each water injection well;
the water saturation ratio calculating module is used for calculating the ratio of the water saturation of each layer of each water well to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection well pair; and
an identification coefficient calculation module for calculating the identification coefficient of the water consumption layer of each horizon on each injection well pair based on the ratio of the dimensionless water consumption to the economic water consumption of each horizon on each injection well pair, the ratio of the water absorption intensity of each injection well pair to the water absorption intensity of each water well, and the ratio of the water saturation of each horizon on each injection well pair to the average water saturation,
wherein the injection well pair refers to a reservoir between the well and a specific well located around the well within a control area of each well,
The calculating of the ratio of dimensionless water consumption to economic water consumption on each injection well pair of each layer of each water well comprises the following steps:
by the formula
Figure FDA0004038651300000071
Calculating dimensionless water consumption E of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, Q iko Daily oil production for horizon i on injection and production well pair k;
by the formula
Figure FDA0004038651300000072
Calculating economic water consumption E w Wherein P is o Price per unit volume of crude oil, C w Cost per unit water injection amount; and
by the formula
Figure FDA0004038651300000073
Calculating the ratio R of dimensionless water consumption to economic water consumption of horizon i on injection well pair k ike
Wherein, based on the effective thickness, daily water absorption of each layer on each injection well pair of each well, and the effective thickness of each well, the ratio of the water absorption intensity of each layer on each injection well pair to the water absorption intensity of each well is calculated to include:
by the formula
Figure FDA0004038651300000074
Calculating the water absorption intensity omega of horizon i on injection well pair k ik Wherein I ikw For the daily water absorption of horizon i on injection well pair k, h ik The effective thickness of the horizon i on the injection well pair k is the effective thickness of the horizon i on the injection well pair k;
by the formula
Figure FDA0004038651300000075
Calculating the water absorption intensity omega of each water well, wherein I is the water absorption capacity of each water well, < > and- >
Figure FDA0004038651300000076
h is the effective thickness of each water well; and
by the formula
Figure FDA0004038651300000081
Calculating the ratio R of the water absorption intensity of the horizon i on the injection and production well pair k to the water absorption intensity of each water well ikd
The calculating of the water consumption layer zone identification coefficient of each layer of each water well on each injection well pair comprises the following steps:
calculating the water consumption layer zone identification coefficient C on the injection well pair k by the following formula ik ,C ik =ln(R ike ×R iks ×R ikd ),
Wherein R is ike The ratio of dimensionless water consumption to economic water consumption of the horizon i on the injection well pair k is the ratio of dimensionless water consumption to economic water consumption of the horizon i; r is R ikd The ratio of the water absorption intensity of the horizon i on the injection well pair k to the water absorption intensity of each water well is given; r is as follows iks Is the ratio of the water saturation of horizon i over the injection well pair k to the average water saturation.
9. The system for identifying a water-consuming zone of a high water reservoir of claim 8, wherein the fitting means comprises:
the fine geologic model building module is used for building a fine geologic model of the target area by adopting a geologic modeling algorithm based on geologic data of the target area in the high-water-cut reservoir; and
an importing module for importing the fine geologic model into the reservoir numerical simulator;
and the fitting module is used for combining the production dynamic data of the oil well and the water well in the target area, and adjusting the characteristic parameters in the fine geological model in the oil reservoir numerical simulator so as to fit the production dynamic data of the oil well and the water well.
10. The system for identifying a water-consuming zone of a high water reservoir of claim 8, wherein the identifying means comprises:
the development level identification module is used for executing the following operations:
identification coefficient C of the water-consuming belt ik If the identification layer i is larger than or equal to a first preset coefficient, the identification layer i is a first high water consumption layer zone on the injection well pair k;
identification coefficient C of the water-consuming belt ik Under the condition that the first preset coefficient is smaller than the first preset coefficient and larger than or equal to the second preset coefficient, the identification horizon i is a second high water consumption zone on the injection well pair k; and
identification coefficient C of the water-consuming belt ik And under the condition that the identification layer position i is smaller than the second preset coefficient, the identification layer position i is a common water-consuming layer zone on the injection well pair k.
11. The system for identifying a water-consuming zone of a high water reservoir of claim 8, further comprising:
a development direction identification module for identifying the development direction of the high water-consuming layer zone in the direction of the injection well pair k corresponding to the first high water-consuming layer zone and/or the second high water-consuming layer zone,
the injection and production well pair refers to a reservoir between each well and a specific oil well around the well in the control area of the well.
12. The system for identifying a water-consuming zone of a high water reservoir of claim 8, further comprising:
the parameter obtaining device is used for obtaining the permeability level difference, the water content, the crude oil viscosity and the water absorption intensity in each injection well pair between each water well and each oil well around each water well based on the oil reservoir numerical simulation model;
a volume calculating device for calculating the volumes of the water-consuming zones with different development levels based on the permeability level difference, the water content, the crude oil viscosity and the water absorption strength of each water well in each injection well pair,
the injection and production well pair refers to a reservoir between each well and a specific oil well around the well in the control area of the well.
13. A machine-readable storage medium having instructions stored thereon for causing a machine to perform the method of identifying water-depleted zones of a high water reservoir according to any one of claims 1 to 7.
CN201910870315.4A 2019-09-16 2019-09-16 Identification method and identification system for water consumption zone of high water-cut oil reservoir Active CN110593863B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910870315.4A CN110593863B (en) 2019-09-16 2019-09-16 Identification method and identification system for water consumption zone of high water-cut oil reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910870315.4A CN110593863B (en) 2019-09-16 2019-09-16 Identification method and identification system for water consumption zone of high water-cut oil reservoir

Publications (2)

Publication Number Publication Date
CN110593863A CN110593863A (en) 2019-12-20
CN110593863B true CN110593863B (en) 2023-05-09

Family

ID=68859767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910870315.4A Active CN110593863B (en) 2019-09-16 2019-09-16 Identification method and identification system for water consumption zone of high water-cut oil reservoir

Country Status (1)

Country Link
CN (1) CN110593863B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111911135B (en) * 2020-07-21 2021-06-25 中国石油大学(华东) Dynamic description method for high water consumption strip of water-drive reservoir
CN112632864B (en) * 2020-12-14 2023-08-04 中国石油化工股份有限公司胜利油田分公司孤岛采油厂 Ore field method for identifying high-water-consumption zone by dynamic and static combination
CN114763736A (en) * 2021-01-13 2022-07-19 中国石油化工股份有限公司 Method for treating high water consumption of fault block oil reservoir
CN114810006B (en) * 2021-01-27 2024-05-31 中国石油化工股份有限公司 Heterogeneous flooding post-treatment separate layer water injection regulation and control high water consumption layer potential evaluation method
CN113297740A (en) * 2021-05-26 2021-08-24 中国石油大学(华东) Horizontal well gel profile control parameter optimization method
CN113505507B (en) * 2021-06-29 2022-09-09 中国石油化工股份有限公司 Method for identifying high water consumption zone between injection and production wells
CN113537585B (en) * 2021-07-09 2023-04-07 中海石油(中国)有限公司天津分公司 Oil field production increasing measure recommendation method based on random forest and gradient lifting decision tree

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190145226A1 (en) * 2016-06-10 2019-05-16 Total Sa A method and device for optimizing a modelling of flows within a reservoir for optimization of oil/gas production
CN107674661A (en) * 2017-11-09 2018-02-09 中国石油化工股份有限公司 It is a kind of to utilize blocking agent of Oil Field Waste Solid thing production and preparation method thereof, application
CN108301813B (en) * 2017-12-20 2020-01-03 中国石油化工股份有限公司 Method for developing and adjusting multilayer sandstone oil reservoir for development of scattered sand bodies
CN109441437A (en) * 2018-09-07 2019-03-08 中国石油化工股份有限公司 A kind of recognition methods in the highly water intensive channel of ultra-high water cut reservoir

Also Published As

Publication number Publication date
CN110593863A (en) 2019-12-20

Similar Documents

Publication Publication Date Title
CN110593863B (en) Identification method and identification system for water consumption zone of high water-cut oil reservoir
CN109577929B (en) Quantitative evaluation method for establishing effective displacement of ultra-low permeability tight reservoir horizontal well
CN103075142B (en) A kind of waterflooding oil field water blockoff oil well well choosing method
CN103902827A (en) Flow unit division method of carbonate-rock horizontal wells
WO2020063603A1 (en) Dynamic data processing method for oilfield development and production
CN103912248A (en) Method for predicting water contents of water-drive oilfields
CN107130955B (en) Method for determining bottom hole flowing pressure and method for determining natural energy of reservoir
CN108194053A (en) A kind of determining method and device of open-hole horizontal well matrix acidifying acid solution dosage
Boah et al. Decline curve analysis and production forecast studies for oil well performance prediction: a case study of reservoir X
Sheely Description of field tests to determine residual oil saturation by single-well tracer method
CN109670729A (en) A kind of top plate aquifer water well evaluation method
CN207829870U (en) A kind of oil-field flooding fouling experimental provision
Zulqarnain et al. Field‐scale well leakage risk assessment using reduced‐order models
CN112709547A (en) Method for judging water plugging time of water producing gas well
CN116950637B (en) Single well group profile control method and system based on micro-pressure drive and readable storage medium
Bánki et al. Aiding Reservoir Simulation and Maintenance Using analytical Formulae
CN109033575A (en) A kind of water saturation calculation method of Complicated Pore Structures reservoir
CN109033508A (en) A kind of determination method and apparatus of oil reservoir moisture percentage and water_bearing escalating rate
CN109063228A (en) A kind of determination method and device of water-drive pool water_bearing escalating rate variation
Zhou et al. Recovery method and parameter optimization of a pilot test for conformance control flooding and thermal recovery in the offshore heavy oilfield
CN115807657A (en) Method and device for determining residual oil distribution
CN110442950B (en) Multi-factor fusion residual oil distribution prediction graph compiling method
Zhang et al. Effects of structural fracture and in situ stress combination characteristics on the production of coalbed methane wells: A case in the Mabidong block, China
CN118134041A (en) Prediction method for liquid yield of horizontal well sub-tank exploitation
Zhang et al. Evaluation of Well Pattern Adaptability in Low-Permeability Carbonate Gas Reservoir

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
CB02 Change of applicant information

Address after: 266580 No. 66 Changjiang West Road, Huangdao District, Qingdao, Shandong.

Applicant after: China University of Petroleum (East China)

Address before: 266555 No. 66 Changjiang West Road, Qingdao economic and Technological Development Zone, Shandong

Applicant before: China University of Petroleum (East China)

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant