CN110593863A - Identification method and identification system for water-consuming zone of high-water-cut oil reservoir - Google Patents
Identification method and identification system for water-consuming zone of high-water-cut oil reservoir Download PDFInfo
- Publication number
- CN110593863A CN110593863A CN201910870315.4A CN201910870315A CN110593863A CN 110593863 A CN110593863 A CN 110593863A CN 201910870315 A CN201910870315 A CN 201910870315A CN 110593863 A CN110593863 A CN 110593863A
- Authority
- CN
- China
- Prior art keywords
- water
- well
- injection
- zone
- production
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 585
- 238000004519 manufacturing process Methods 0.000 claims abstract description 298
- 239000003129 oil well Substances 0.000 claims abstract description 77
- 238000011161 development Methods 0.000 claims abstract description 68
- 238000004088 simulation Methods 0.000 claims abstract description 45
- 239000003921 oil Substances 0.000 claims description 124
- 238000010521 absorption reaction Methods 0.000 claims description 108
- 238000002347 injection Methods 0.000 claims description 62
- 239000007924 injection Substances 0.000 claims description 62
- 239000010779 crude oil Substances 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000035699 permeability Effects 0.000 claims description 22
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 238000013461 design Methods 0.000 abstract description 7
- 239000010410 layer Substances 0.000 description 97
- 230000008569 process Effects 0.000 description 13
- 230000001186 cumulative effect Effects 0.000 description 11
- 230000000694 effects Effects 0.000 description 11
- 238000009826 distribution Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000004005 microsphere Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000006187 pill Substances 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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-consuming zone of an oil reservoir in a high water-containing period. The identification method comprises the following steps: fitting the production dynamics of the oil well and the water well by using a numerical reservoir simulator based on the geological data of a target area in the high water cut oil reservoir and the production dynamics data of the oil well and the water well in the target area to obtain a numerical reservoir simulation model; calculating the identification coefficient of a water consumption zone between each water well and an oil well around the water well based on the numerical reservoir simulation model; and identifying a development level of the water-consuming layer zone based on the identification coefficient of the water-consuming layer zone. The method can quickly judge and identify the development level of the water-consuming zone and quantitatively characterize the development level, thereby effectively identifying the development direction of the high-water-consuming zone and playing an effective guiding role in the design of a regulation and control scheme in the subsequent oil field development stage.
Description
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-consuming zone of an oil reservoir in a high water-containing period.
Background
After many years of water flooding development, most oil fields in China enter a high/ultrahigh water content development stage, and part of the oil fields even enter an ultrahigh water content later stage (the water content is higher than 95%), so that a high water consumption zone generally develops. The method has the advantages that a large amount of injected water is circulated inefficiently and inefficiently, the dynamic heterogeneity of a reservoir is obviously enhanced, the contradiction between an inner layer and an outer layer is enlarged, the water flooding efficiency of the oil field is reduced, production increasing measures are difficult to implement, the development cost is increased, and the economic recovery ratio of the oil field is finally reduced. Therefore, the method effectively identifies the high water-consumption zone and controls the high water-consumption zone by adopting a reasonable regulation and control technology, and has very important significance for further improving the development effect of the oil deposit in the later period of ultrahigh water content.
In view of the above, some technical solutions have been disclosed, which identify the high water-consumption zone as a water flow dominant channel, for example, by measuring the production dynamics of an oil field production well by using a water-containing index characteristic curve to detect a cross flow channel of the oil field; or identifying the water injection dominant channel by establishing an explanation model of the radius and permeability of the throat of the reservoir by using the core data and the logging data, and preferably blocking the particle size of the microsphere.
According to the technical scheme, whether a certain well group has a high water-consumption zone can be only roughly judged, but the characterization effect of the high water-consumption zone is poor, so that 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-consumption zone of an oil reservoir in a high water-cut period, which can quickly judge and quantitatively characterize the development level of the water-consumption zone, thereby effectively identifying the development direction of the high water-consumption zone and playing an effective guiding role in the design of a regulation and control scheme in the subsequent oil field development stage.
In order to achieve the above object, an aspect of the present invention provides a method for identifying a water-consuming zone of a high water-cut reservoir, the method comprising: fitting the production dynamics of the oil well and the water well by using a numerical reservoir simulator based on the geological data of a target area in the high water cut oil reservoir and the production dynamics data of the oil well and the water well in the target area to obtain a numerical reservoir simulation model; calculating the identification coefficient of a water consumption zone between each water well and an oil well around the water well based on the numerical reservoir simulation model; and identifying a development level of the water-consuming layer zone based on the identification coefficient of the water-consuming layer zone.
Preferably, said fitting production dynamics of said oil well and said water well with a reservoir numerical simulator comprises: based on geological data of a target region in the high water-cut oil deposit, establishing a fine geological model of the target region by adopting a geological modeling algorithm; and importing the refined geological model into the reservoir numerical simulator; and adjusting characteristic parameters in the fine geological model in the oil reservoir numerical simulator by combining the production dynamic data of the oil well and the water well in the target area so as to fit the production dynamic of the oil well and the water well.
Preferably, the calculating of the identification coefficient of the water-consuming zone between each water well and the oil well around the water well comprises: based on the numerical reservoir simulation model, obtaining the effective thickness, water saturation, daily water absorption and daily oil production of each layer of each water well on each injection-production well pair; calculating the ratio of dimensionless water consumption and economic water consumption of each layer of each water well on each injection-production well pair based on the daily water absorption and daily oil production of each layer of each water well on each injection-production well pair and the economic water consumption, 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 on each injection-production well pair to the water absorption intensity of each water well on the basis of the effective thickness of each layer of each water well on each injection-production well pair, the daily water absorption capacity and the effective thickness of each water well; calculating the ratio of the water saturation of each layer of each water well on each injection-production well pair to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection-production well pair; and calculating a water consumption zone identification coefficient of each layer of each water well on each injection-production well pair based on the ratio of the dimensionless water consumption and the economic water consumption of each layer on each injection-production well pair, the ratio of the water absorption intensity of each layer on each injection-production well pair to the water absorption intensity of each water well and the ratio of the water saturation of each layer on each injection-production well pair to the average water saturation, wherein the injection-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 calculating the ratio of the dimensionless water consumption to the economic water consumption of each horizon of each water well on each injection-production well pair comprises: by the formulaCalculating the dimensionless water consumption E of the horizon i on the injection-production well pair kikWherein, IikwThe daily water absorption, Q, of the layer i on the injection-production well pair kikoThe daily oil production of the horizon i on the injection-production well pair k; by the formulaCalculating economic water consumption EwWherein P isoPrice per unit volume of crude oil, CwCost per unit water injection; and by formulaCalculating the ratio R of the dimensionless water consumption and the economic water consumption of the horizon i on the injection-production well pair kike。
Preferably, the calculating the ratio of the water absorption intensity of each horizon of each water well on each injection-production well pair to the water absorption intensity of each water well comprises: by the formulaCalculating the water absorption intensity omega of the horizon i on the injection-production well pair kikWherein, IikwThe daily water absorption, h, of the layer i on the injection-production well pair kikThe effective thickness of the horizon i on the injection-production well pair k; by the formulaCalculating the water absorption intensity omega of each water well, wherein I is the water absorption capacity of each water well,h is the effective thickness of each well; and by formulaCalculating the ratio R of the water absorption intensity of the horizon i on the injection-production well pair k to the water absorption intensity of each wellikd。
Preferably, the calculating the ratio of the water saturation of each horizon on each injection-production well pair to the average water saturation in the target area comprises: by the formulaCalculating the water saturation S of the horizon i on the injection-production well pair kwikAnd the average water saturationRatio R ofiks。
Preferably, the calculating the water-consuming zone identification coefficient of the horizon on each injection-production well pair comprises: calculating the water consumption zone identification coefficient C of the horizon i on the injection-production well pair k by the following formulaik,Cik=ln(Rike×Riks×Rikd) Wherein R isikeThe ratio of dimensionless water consumption and economic water consumption of the horizon i on the injection-production well pair k is shown; rikdThe ratio of the water absorption intensity of the horizon i on the injection-production well pair k to the water absorption intensity of each well is shown; and RiksIs at horizon iA ratio of water saturation over the injection-production well pair k to the average water saturation.
Preferably, said identifying a developmental level of said water-consuming layer band comprises: identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is greater than or equal to a first preset coefficient, the identification horizon i is a first high-water-consumption horizon on the injection-production well pair k; identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is smaller than the first preset coefficient and is larger than or equal to a second preset coefficient, the identification horizon i is a second high water consumption horizon on the injection-production well pair k; and the identification coefficient C in the water-consuming zoneikAnd under the condition that the identification horizon i is smaller than the second preset coefficient, identifying that the horizon i is a common water-consuming horizon on the injection-production well pair k.
Preferably, the identification method further comprises: and identifying the direction of an injection and production well pair k corresponding to the first high water-consumption zone and/or the second high water-consumption zone as the development direction of the high water-consumption zone, wherein the injection and production well pair refers to a reservoir between each water well and a specific oil well positioned around the water well in the control area of the water well.
Preferably, the identification method further comprises: after performing the step of identifying a developmental level of the water-consuming layer band, performing the following operations: acquiring the permeability grade difference, the water content, the crude oil viscosity and the water absorption strength of each injection-production well pair between each water well and each oil well around each water well on the basis of the numerical reservoir simulation model; and calculating the volumes of the water consumption zone 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 and production well pair, wherein the injection and production well pair refers to a reservoir between each water well and a specific oil well positioned around the water well in the control area of each water well.
Preferably, the developmental levels of the water-consuming layer zones include: ordinary water-consuming layer area reaches high water-consuming layer area, high water-consuming layer area includes: a first high water-consuming zone and a second high water-consuming zone, respectively, said calculating the volume of said water-consuming zones of different developmental levels comprising: calculating the common by the following three formulasThe volume percentages of the water passing layer zone, the second high water consumption layer zone and the first high water consumption layer zone in the injection and production well pairAnd
wherein x is1Is the permeability level difference of each injection-production well pair, mu is the viscosity of the crude oil, fwThe water content is the above; omega is the water absorption intensity of each water well; and calculating the volume V of the j-th water-consumption zone of the horizon i on the injection-production well pair k by the following formulaikj,Wherein j is 1, 2 or 3, and the class 1 water-consuming layer strip corresponds to the common water-consuming layer strip; the class 2 water-consuming layer zone corresponds to the second high water-consuming layer zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, AikThe area of a horizontal transverse plane h of the horizon i in the direction of the injection-production well pair kikIs the effective thickness of horizon i over the pair k of injection and production wells.
Compared with the prior art, the identification method of the water-consuming zone has the following advantages:
(1) fitting the production dynamics of an oil well and a water well in a target area in the high water cut oil reservoir to obtain an oil reservoir numerical simulation model; then, based on the obtained numerical reservoir simulation model, calculating the identification coefficient of a water consumption zone between each water well and the surrounding oil wells; finally, according to the identification coefficient of the water-consuming zone, different development levels are identified, so that the development levels of the water-consuming zone can be rapidly judged and quantitatively characterized, the development direction of the high-water-consuming zone can be effectively identified, and an effective guiding effect is played for the design of a regulation and control scheme in the subsequent oil field development stage.
(2) The direction of each injection well pair of each layer is taken as a research object, and the water-consumption zone of the oil reservoir in the high water-cut period is identified, so that the development layer position, the direction, the level and the volume of the water-consumption zone can be quickly and accurately judged, and the characterization effect of the water-consumption 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 zone identification coefficient, namely, the influence of economic factors such as crude oil price on the development of the oil reservoir in the high water-cut 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 zone according to the market demand condition, thereby adjusting the development strategy in time and ensuring to obtain the best economic benefit.
In a second aspect, the present invention provides a system for identifying a water-consuming zone of a high water-cut reservoir, the system comprising: the fitting device is used for fitting the production dynamics of the oil well and the water well by using a numerical reservoir simulator based on the geological data of a target area in the oil reservoir in the high water cut period and the production dynamics data of the oil well and the water well in the target area so as to obtain a numerical reservoir simulation model; the identification coefficient calculation device is used for calculating the identification coefficient of a water consumption zone between each water well and the oil well around the water well on the basis of the numerical reservoir simulation model; and an identification device for identifying the development level of the water-consuming zone based on the identification coefficient of the water-consuming zone.
The specific implementation details and effects of the identification system of the water-consumption zone of the oil reservoir in the high water cut stage are the same as those of the identification method of the water-consumption zone of the oil reservoir in the high water cut stage, and are not repeated herein.
A third aspect of the invention provides a machine-readable storage medium having stored thereon instructions for causing a machine to perform the above-described method for identifying a water-depletion zone of a high-water-cut reservoir.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 is a flow chart of a method for identifying a water-consuming zone of a high water-cut reservoir according to an embodiment of the present invention;
FIG. 2 is a flow chart of a process of 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 coefficients of a water-consuming zone between each water well and the wells surrounding the water well, according to an embodiment of the present invention;
FIG. 4 is a flow chart of a process for identifying a water-consuming zone provided by an embodiment of the present invention;
FIG. 5 is a water cut fit within a target region provided by an embodiment of the invention;
FIG. 6 is a cumulative oil production fit within a target area provided by an embodiment of the present invention;
FIG. 7A is an oil saturation distribution of the L3 layer in the target area provided by an embodiment of the present invention;
FIG. 7B is an oil saturation distribution of the L35 layer in the target area provided by an embodiment of the present invention;
FIG. 8A is a schematic illustration of the control areas between a water well I1 in a target area and four oil wells P1, P2, P3, P4 located around the water well I1 according to an embodiment of the present invention;
FIG. 8B is a schematic diagram of a pair of injection and production wells 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-cut reservoir according to an 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 detailed description of embodiments of the invention refers to the accompanying drawings. It should be understood that the detailed description and specific examples, while indicating the present invention, are given by way of illustration and explanation only, not limitation.
Because the high water-cut period (the high water-cut period in the oil reservoir engineering can be generally divided into three stages with different water contents, namely, the high water-cut period with the water content of 80-90%, the ultrahigh water-cut period with the water content of more than 90%, and the ultrahigh water-cut later stage with the water content of more than 95%, especially the ultrahigh water-cut later stage) of the oil reservoir, the high water-cut zone of the oil reservoir generally develops, so that a large amount of injected water is inefficient and ineffective in circulation, the dynamic heterogeneity of the reservoir is remarkably increased, the intra-layer and inter-layer contradiction is enlarged, the water drive efficiency of the oil field is reduced, the development cost is increased, and the economic recovery ratio of the oil field is finally reduced. In order to solve the above problems, it is necessary to efficiently identify the high water-consuming layer band. However, the prior art can only roughly judge whether a certain well group has a high water-consumption zone, but cannot confirm the development position, direction, level and volume of the high water-consumption zone, and has poor characterization effect on the high water-consumption zone. Of course, the three stages of the high water cut period in the present invention are not limited to the three specific water cut ranges, and other reasonable dividing manners are possible.
Based on this, the invention identifies the water-consuming zone by: on the basis of collecting geological data and production dynamic data of a high water cut oil reservoir, firstly, carrying out fine oil reservoir numerical simulation research based on the geological data of the high water cut oil reservoir; then, production history fitting is carried out by combining the production dynamic data of the water well (namely, a water injection well) and the oil well (namely, an oil production well), so that the data of water saturation, daily water absorption, daily oil production and the like of each small layer (namely, horizon) on each injection and production well pair can be obtained; finally, characterizing the water-consuming zone based on the obtained data, namely calculating an identification coefficient of the water-consuming zone, and quantitatively evaluating the development degree of the water-consuming zone through the identification coefficient; judging the development level and direction of a high water-consumption zone of each layer in the direction of each injection-production well pair by combining a preset identification criterion; and the volume of the water consumption layer belt with different levels can be calculated according to a simulation formula. Therefore, the method for identifying the water-consumption zone of the oil reservoir in the high water-cut period can accurately determine the position, direction, level and volume of the high water-consumption zone, so that the representation effect of the high water-consumption zone is improved.
Before describing particular embodiments of the present invention, the following references to pairs of injection and production wells are explained and illustrated. The injection-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 the water well. The control area of the water well is the sum of the planar areas of the water well and the oil wells located around the water well, which are horizontally crossed (the planar area between the water well I1 and the oil well P2 shown in fig. 8B is the area a). Wherein, the direction of each injection-production well pair indicates different water injection directions of the water well. For example, the area shown in fig. 8A is the control area of well I1.
Fig. 1 is a flow chart of a method for identifying a water-consuming zone of a high water-cut reservoir provided by the invention. As shown in fig. 1, the identification method may include the steps of: step S101, based on geological data of a target area in the high water-cut oil deposit and production dynamic data of an oil well and a water well in the target area, fitting the production dynamics of the oil well and the water well by using an oil deposit numerical simulator to obtain an oil deposit numerical simulation model; step S102, calculating the identification coefficient of a water consumption zone between each water well and an oil well around the water well based on the numerical reservoir simulation model; and a step S103 of identifying the development level of the water-consuming zone based on the identification coefficient of the water-consuming zone.
Before step S101 is executed, geological data, production dynamic data, etc. of the high water-cut reservoir need to be collected. Specifically, the collecting geological data of the target area in the high water-cut reservoir may include: and collecting a structural contour map, fault track data, a sand body thickness distribution contour map, an effective thickness distribution contour map, a porosity distribution contour map, a permeability distribution contour map, a partition distribution map and the like of different layers in the target region. The collecting production dynamics for oil and water wells within the target zone may comprise: the daily water injection amount (equivalent to daily water absorption) of each water well, the daily oil production amount and water content of each oil well around each water well, and the water absorption/liquid production profile data of the oil-water well 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 oil deposit, a geological modeling algorithm is adopted to establish a fine geological model of the target area.
A fine geological model of the target region may be established using a geological modeling algorithm in Petrel software based on the collected geological data.
And S202, importing the fine geological model into the numerical reservoir simulator.
The refined geological model of the target region may be imported into a reservoir numerical simulator Eclipse for simulation to fit production dynamics of oil and water wells within the target region in the reservoir numerical simulator Eclipse.
And S203, 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 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 and the like) by means of an oil reservoir numerical simulation method, and compare the calculated development indexes (such as pressure, yield, comprehensive water content and the like) 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, the model parameters are required to be properly adjusted, and then the calculation is carried out again by using the adjusted model parameters until the calculation result is consistent with the actual dynamic state or within an allowable error range.
Specifically, in the production history fitting process, the involved characteristic parameters include: the original geological reserves of the reservoir, the dynamic parameters of the reservoir, and the dynamic parameters of the individual wells (referred to as individual oil and water wells). The dynamic parameters of the reservoir may include: the cumulative oil production, the cumulative water production, the cumulative gas-oil ratio, the comprehensive water content, the average pressure of the stratum, the cumulative water injection and the like change along with time. The dynamic parameters of the well may include: instantaneous oil production, water content, water production, instantaneous gas-oil ratio, cumulative oil production, cumulative water production, bottom hole pressure and the like change with time. The dynamic parameters of the water well may include: instantaneous water injection, cumulative water injection, bottom hole pressure, etc. over time. In the actual production history fitting process, the original geological reserves of the oil reservoir are fitted, then the comprehensive parameters (such as cumulative oil production, cumulative water production, cumulative gas-oil ratio and the like) of the oil reservoir 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 a relative permeability curve, permeability data of an interwell oil layer, an effective reservoir thickness, and a compressibility of rock and fluid in the fine geological model.
Production history fitting requires that the fitting error of the original geological reserves of the oil reservoir is within 2 percent, and the fitting error of the comprehensive parameters of the oil reservoir is within 5 percent. For each parameter involved in the production dynamics of a single well, the accuracy of the history fit is often evaluated using the following 4 indicators: (1) maximum forward relative error Epmax: when the calculated value is higher than the actual value, the ratio of the difference obtained by subtracting the corresponding value on the actual curve from the parameter value on the calculated curve to the actual value; (2) negative maximum relative error Enmax: meterWhen the calculated value is lower than the actual value, calculating the absolute value of the ratio of the difference obtained by subtracting the corresponding value on the actual curve from the parameter value on the curve and the actual value; (3) mean value of relative error Erave: calculating an average of absolute values of errors between the values and the actual values; (4) standard deviation of absolute error Eas: the square of the average of the sum of squares of the difference between the value and the actual value is calculated. For a certain parameter involved in the production dynamics of a single well, the following condition is satisfied: epmax<10%,Enmax<10%,Erave<5%, and Eas<0.1, the parameter is considered to have been fitted. If 90% of all the individual wells meet the above criteria, the resulting refined geological model (i.e. the reservoir numerical simulation model) can be considered accurate, i.e. the fitting accuracy meets the requirements. Of course, each fitting error threshold in the present invention is not limited to the above specific data, and may be reasonably adjusted according to actual situations.
Fig. 5 and 6 are fitting results of the water content and the accumulated oil production in the target region, respectively, and it can be found that the prediction result of the numerical reservoir simulation model is well matched with the actual production dynamics of the reservoir, so that the numerical reservoir 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 around the water well may include the following steps:
and S301, acquiring the effective thickness, water saturation, daily water absorption and daily oil production of each layer of each water well on each injection-production well pair based on the numerical reservoir simulation model.
Step S302, calculating the ratio of the dimensionless water consumption and the economic water consumption of each layer of each water well on each injection-production well pair based on the daily water absorption and the daily oil production of each layer of each water well on each injection-production well pair and the economic water consumption.
Wherein the economic water consumption is related to the price of the crude oil and the cost of the water injection. In particular, it can be represented by a formulaCalculating the dimensionless water consumption E of the horizon i on the injection-production well pair kikWherein, IikwThe daily water absorption of the layer i on the injection-production well pair k, m3/d,QikoFor the daily oil production, m, of horizon i on injection-production well pair k3D; by the formulaCalculating economic water consumption EwWherein P isoIs the price per unit volume of crude oil, yuan/m3,CwCost per unit water injection, yuan/m3(ii) a And by formulaCalculating the ratio R of the dimensionless water consumption and the economic water consumption of the horizon i on the injection-production well pair kike。
The embodiment provides that the economic water consumption is taken as an important index for identifying the water consumption zone, not only can reflect the scientific fact that the relative flow capacity difference of oil and water in a high water-cut stage (including an ultra-high water-cut later stage) is increased rapidly, but also considers the influence of economic factors such as the price of crude oil on the development of the oil deposit in the high water-cut stage, and the like, thereby meeting the requirement of high-efficiency development of the oil deposit in the high water-cut stage under the condition of low oil price. The index can change along with the oil price, so that the development condition of the high-water-consumption zone can be dynamically judged by an oil field company according to the market demand condition, the development strategy is adjusted in time, and the optimal economic benefit is ensured to be obtained.
Step S303, calculating the ratio of the water absorption intensity of each layer of each well on each injection-production well pair to the water absorption intensity of each well based on the effective thickness of each layer of each well on each injection-production well pair, the daily water absorption capacity and the effective thickness of each well.
In particular, it can be represented by a formulaCalculating the water absorption intensity omega of the horizon i on the injection-production well pair kikWhich isIn (I)ikwThe daily water absorption of the layer i on the injection-production well pair k, m3/(d·m),hikThe effective thickness m of the horizon i on the injection-production well pair k; by the formulaCalculating the water absorption intensity omega, m of each well3V (d.m), wherein I is the water absorption capacity of each well, m3/d,h is the effective thickness of each well (as shown in fig. 8B), m; and by formulaCalculating the ratio R of the water absorption intensity of the horizon i on the injection-production well pair k to the water absorption intensity of each wellikd。
And S304, calculating the ratio of the water saturation of each layer of each water well on each injection-production well pair to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection-production well pair.
In particular, it can be represented by a formulaCalculating the water saturation S of the horizon i on the injection-production well pair kwikAnd the average water saturationRatio R ofiks。
And S305, calculating a water consumption zone identification coefficient of each layer of each water well on each injection-production well pair based on the ratio of the dimensionless water consumption and the economic water consumption of each layer on each injection-production well pair, the ratio of the water absorption intensity of each layer on each injection-production well pair to the water absorption intensity of each water well, and the ratio of the water saturation of each layer on each injection-production well pair to the average water saturation.
In particular, can be prepared byCalculating the water consumption zone identification coefficient C of the horizon i on the injection-production well pair k by the following formulaik,Cik=ln(Rike×Riks×Rikd) Wherein R isikeThe ratio of dimensionless water consumption and economic water consumption of the horizon i on the injection-production well pair k is shown; rikdThe ratio of the water absorption intensity of the horizon i on the injection-production well pair k to the water absorption intensity of each well is shown; and RiksAnd the ratio of the water saturation of the horizon i on the injection-production well pair k to the average water saturation is shown.
For step S103, the identifying the development level of the water-consuming layer zone based on the identification coefficient of the water-consuming layer zone may include: identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is greater than or equal to a first preset coefficient, the identification horizon i is a first high-water-consumption horizon on the injection-production well pair k; identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is smaller than the first preset coefficient and is larger than or equal to a second preset coefficient, the identification horizon i is a second high water consumption horizon on the injection-production well pair k; and the identification coefficient C in the water-consuming zoneikAnd under the condition that the identification horizon i is smaller than the second preset coefficient, the identification horizon i is a common water-consuming zone on the injection-production well pair k, wherein the water consumption capacities of a first high water-consuming zone (such as an extreme high water-consuming zone), a second high water-consuming zone (a common high water-consuming zone) and the common water-consuming zone are weakened in sequence. For example, the level of development of the water-consuming zone may be determined using the following identification criteria: if Cik<0, the layer position i is a common water-consuming layer zone on the injection-production well pair k; if 0 is less than or equal to Cik<0.8, the layer position i is a common high water consumption zone on the injection-production well pair k; if CikAnd if the average value is more than or equal to 0.8, the horizon i is an extremely high water-consumption horizon on the injection-production well pair k.
The identification method may further include: and identifying the direction of an injection-production well pair k corresponding to the first high water-consumption zone and/or the second high water-consumption zone as the development direction of the high water-consumption zone. For example, the direction of the pair k of injection and production wells of a general high water-consuming zone and an extreme high water-consuming zone can be identified as the development direction of the high water-consuming zone.
After performing step 103, embodiments of the present invention may also perform a quantitative calculation of the volume of the water-consuming zone at different developmental levels. The specific calculation process may include: acquiring the permeability grade difference, the water content, the crude oil viscosity and the water absorption strength of each injection-production well pair between each water well and each oil well around each water well on the basis of the numerical reservoir simulation model; and calculating the volumes of the water consumption 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.
Specifically, first, the volume percentages of the ordinary water-consumption zone, the second high water-consumption zone and the first high water-consumption zone in the injection-production well pair can be calculated respectively through the following three formulasAnd
wherein x is1Mu is the viscosity of the crude oil, mPa.s, f is the permeability level difference of each injection-production well pairwThe water content is the above; omega is the water absorption intensity of each water well.
Secondly, the volume V of the j-th water-consumption zone of the horizon i on the injection-production well pair k can be calculated by the following formulaikj,Wherein j is 1, 2 or 3, and the class 1 water-consuming layer strip corresponds to the common water-consuming layer strip; class 2 water-consuming layer bands correspond toThe second high water-consumption zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, AikThe area of the horizontal cross-section plane of the horizon i in the direction of the injection-production well pair k (area A shown in FIG. 8B), m2,hikIs the effective thickness, m, of horizon i over the pair of injection and production wells k.
Specifically, the process of identifying (or characterizing) water-consuming zones is explained and illustrated in detail below with respect to reservoirs as an example.
After many years of water injection development, the XX unit of a certain oil field enters the later stage of exploitation with ultra-high water content, and a high water consumption zone generally develops, so that a large amount of injected water is in low-efficiency and ineffective circulation, the water content of the oil well rises quickly, and the water drive utilization degree is low. The target area is a river facies positive rhythm oil reservoir, 2 wells (1-2-95 and 1-2-11) and 6 wells (1-2-9, 1-2-92, 1-2-94, 1-2N91, 1-2-111 and 1-2X112) are commonly used for normal production at present, the comprehensive water content is 95.6 percent, and the method is a typical representative of the oil reservoir which is completely filled at the later stage with ultrahigh water content.
As shown in fig. 4, the identification process of the water-consuming zone may include the steps of:
step S401, collecting geological data in the target area, and production dynamic data and monitoring data of oil wells and water wells in the target area.
And S402, fitting the production dynamics of the oil well and the water well.
Step S403, calculating the identification coefficient of the water-consuming layer zone.
Firstly, for each water well, reading the effective thickness, water saturation, daily water absorption, daily oil production and other data of each layer in each injection and production well pair direction from an oil reservoir numerical simulation model after history fitting to obtain a data table shown in table 1. The data are easy to obtain, the calculation process is simple, the workload of field technicians is greatly reduced, and the operability is strong.
TABLE 1 basic data table of each horizon on each injection-production well pair
As can be seen from Table 1, wells 1-2-95 correspond to 4 horizons, and for wells 1-2-95, only 1 well 1-2-9 is located at horizon L3; at horizon L21, corresponding to 3 wells, 1-2-9, 1-2-92, and 1-2-94, respectively. The 3 injection-production well pairs for wells 1-2-95 at horizon L21 are reservoirs at horizon L21 between them and the 3 wells (1-2-9, 1-2-92, and 1-2-94), respectively. The data of the daily oil production, daily water absorption and the like of the injection and production well pair corresponding to the horizon L3 water feeding well 1-2-95 and the oil well 1-2-9 are the data of the first row in the table 1.
Secondly, calculating the ratio of the dimensionless water consumption and the economic water consumption of each layer on each injection-production well pair aiming at each well.
Firstly, according to the daily water absorption I of the horizon I on the injection-production well pair k in the table 1ikwAnd daily oil production QikoCalculating the dimensionless water consumption E of the horizon i on the injection-production well pair k by using the following formulaik:And the calculated data are listed in table 2. For example, dimensionless water consumption was calculated to be 326.46 for horizon L3 on 1-2-95 and 1-2-9 injection-production well pairs, which are filled in Table 2.
TABLE 2 dimensionless water consumption, water absorption strength and water saturation values in each injection-production well pair
Second, assume that the current crude oil price is $ 50/barrel, the exchange rate of $ 6.7 yuan for Renminbi, and the cost per unit of water injection CwIs 8 yuan/m3Then the economic water consumption can be calculated using the following equation:wherein, PoIs the price per unit volume of crude oil, yuan/m3. The economic water consumption is 263.4m3。
Finally, the horizon can be calculated using the following formulai ratio of dimensionless water consumption to economic water consumption (abbreviated as water consumption ratio) R on injection-production well pair kike:And the calculated data are listed in table 3. For example, the ratio of dimensionless water consumption to economic water consumption for horizon L3 on 1-2-95 and 1-2-9 injection-production well pairs was calculated to be 1.239 and filled in Table 3.
TABLE 3 result of calculation of water-consumption zone identification coefficient of each injection-production well pair
Thirdly, calculating the ratio of the water absorption intensity of each layer on each injection-production well pair to the water absorption intensity of the whole well (namely each well) aiming at each well.
First, from the data in Table 1, the water absorption strength ω of the horizon i on the injection-production well pair k can be calculated using the following formulaik:And the calculated data are listed in table 2. For example, horizon L3 has a water absorption strength of 1.123 on the 1-2-95 and 1-2-9 pair of injection and production wells, and is filled in Table 2.
Secondly, the water absorption strength ω of the whole well can be calculated using the following formula:wherein I is the water absorption capacity of the whole well,m3d; h is the total effective thickness, m, corresponding to the whole well.
For example, from the data in Table 1, the daily water absorption of wells 1-2-95 was calculated to be 63.536m3D (note that only daily water absorption of a certain injection-production well pair is calculated on the horizon L21), and the total effective thickness corresponding to the whole well is 8.388m (which can be directly obtained from an oil reservoir numerical simulation model), so thatThe water absorption intensity of the whole well of the water wells 1-2-95 is 7.574m3/(d.m). The water absorption strength of the whole well for wells 1-2-11, which is 6.672m, can also be obtained in a similar manner3/(d.m). The water absorption intensity values for the whole well are listed in table 4.
Finally, the water absorption strength omega of the horizon i on the injection-production well pair k can be calculated by using the following formulaikRatio R to total well Water absorption Strength omega (Water absorption Strength ratio for short)ikd:And the calculated data are listed in table 3.
Fourthly, calculating the ratio of the water saturation of each layer on each injection-production well pair to the average water saturation for each well.
Calculating the water saturation of the horizon i on the injection-production well pair k by using the following formula; swikAnd average water saturationRatio of (simply referred to as water saturation ratio) Riks:And the calculated data are listed in table 3.
Specifically, the SwikAndall 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 pairs of injection and production wells on the L3 pill is 0.5465, and the average water saturation of the cell is 0.5211, so that R isiksThe value of (d) is 1.048, which is filled in table 3.
And fifthly, calculating the water consumption zone identification coefficient of each layer on each injection-production well pair aiming at each water well.
Calculating the water consumption zone identification coefficient C of the horizon i on the injection-production well pair k by using the following formulaik:Cik=ln(Rike×Riks×Rikd). Specifically, canSubstituting the data in the table 3 into the formula to calculate the water-consuming zone identification coefficient C of the horizon i on the injection-production well pair kikAnd the calculated data are listed in table 3.
And step S404, judging the development level and direction of the water-consuming zone by combining the identification criterion of the water-consuming zone.
If Cik<0, the water consumption zone of the layer i on the injection-production well pair k is a common water consumption zone;
if 0 is less than or equal to Cik<0.8, the water consumption zone of the layer i on the injection-production well pair k is a common high water consumption zone;
if CikAnd if the water consumption zone of the horizon i on the injection-production well pair k is more than or equal to 0.8, the water consumption zone is an extremely high water consumption zone.
And identifying the direction of the injection-production well pair corresponding to the general high water-consumption zone and the extreme high water-consumption zone as the development direction of the high water-consumption zone.
The identification criterion of the water-consuming zone is determined by the following method: according to the actual geological data of the oil deposit field and the production dynamic characteristics of each water well and each oil well, a large number of conceptual models are established to carry out numerical simulation research. And respectively calculating the water consumption zone identification coefficient of each layer on each injection-production well pair according to the simulation result, and performing cluster analysis on the calculation result. Research results show that the simulated 2516 group data can be divided into three categories through a clustering algorithm, and then the identification criterion of each category of water consumption zone can be determined according to clustering analysis results.
The identification results of the water-consuming layer bands in this example are shown in table 3. From Table 3, it can be seen that the horizon L21 layer within the target zone develops extremely high water-consumption zones between wells 1-2-95 and wells 1-2-9, 1-2-94; an extremely high water consumption zone develops between the water feeding well 1-2-95 at the L23 layer and the oil well 1-2-92; horizon L35 layers develop extremely high water-consumption 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 zone between the water wells 1-2-95 and the oil wells 1-2-9, as shown in FIG. 7A; horizon L35 layers developed extremely high water zones between wells 1-2-11 and wells 1-2-94, 1-2-111, and 1-2X112, as shown in figure 7B (direction of arrows), which is consistent with the results listed in table 3.
Step S405, the volumes of the water-consuming zone of different levels are calculated.
Firstly, parameters such as permeability grade difference, water content, crude oil viscosity, water absorption strength and the like in each injection and production well pair are obtained from the numerical reservoir simulation model. Table 4 lists only the permeability grade difference, water content, crude oil viscosity, water absorption strength and other parameters in each injection and production well pair of the extremely high water-consumption zone.
TABLE 4 basic data of each injection-production well pair
Secondly, according to the data in the table 4, for each water well, the volume percentage of the common water-consuming zone in each injection-production well pair is calculated by adopting the following formula
Aiming at each water well, the volume percentage of the general high water consumption zone in each injection-production well pair is calculated by adopting the following formula
Aiming at each water well, the volume percentage of the extremely high water consumption zone in each injection-production well pair is calculated by adopting the following formula
The volume percent of the different levels of water-depleted zone for each injection-production well pair direction in this example is shown in table 5.
TABLE 5 percent by volume of water-consuming layer tapes of different grades
The calculation formula of the volume percentage of the water consumption zone with different levels in each injection-production well pair is determined by adopting the following method: according to the actual geological data of the oil deposit field and the production dynamic characteristics of each water well and each oil well, a large number of conceptual models are established to carry out numerical simulation research. Respectively calculating the water consumption zone identification coefficients of each layer on each injection-production well pair according to the simulation result; determining the volumes of the water consumption zone in different levels in each injection-production well pair by using the identification criteria of the water consumption zone; and then sensitivity analysis finds that the volume of the water-consuming zone is mainly influenced by parameters such as permeability grade difference, water content, crude oil viscosity, water absorption strength and the like, so that statistical regression is carried out on the volume of the water-consuming zone obtained by simulation to obtain a volume percentage calculation formula of the water-consuming zones of different grades in the injection and production well pair.
Finally, the volume V of the j-th water-consuming zone of the horizon i on the injection-production well pair k is calculated by adopting the following formulaikj:Wherein j is 1, 2 or 3, and the class 1 water-consuming layer strip corresponds to the common water-consuming layer strip; the class 2 water-consuming layer zone corresponds to the second high water-consuming layer zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, AikThe area of the horizontal cross-section plane of the horizon i in the direction of the injection-production well pair k (area A shown in FIG. 8B), m2,hikIs the effective thickness, m, of horizon i over the pair of injection and production wells k.
The volumes of the water consumption zone zones with different levels can provide guidance for the actual profile control and water plugging design of the oil field site. According to the identification knotAccording to the method, deep profile control is carried out on a water well in a target area in 2018 and 10 months, and oil field site water absorption profile test data and an interwell tracer monitoring result show that an extremely high water consumption zone is effectively blocked, the water absorption capacity of a general high water consumption zone is effectively inhibited, residual oil in a general water consumption zone is effectively used, and the average daily oil yield of a well group is increased by 13.2m3And 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 oil field shows that the identification method can effectively identify the high water-consumption zone in the high water-containing period (including the ultra-high water-containing later period), the actual coincidence rate with the oil field reaches more than 90%, and the volume of the obtained water-consumption zone can effectively guide the design of the regulation and control measures.
In the embodiment of the invention, by collecting geological data of the oil deposit at the later stage of ultrahigh water content, and production dynamic characteristics and monitoring data of each water well and each oil well, fine oil deposit numerical simulation research is developed and production history fitting is carried out, so that data such as water saturation, daily water absorption, daily oil yield and the like of each layer on each injection-production well pair can be obtained, and basic data is provided for identification of a water-consuming zone; the development degree of the water-consuming zone can be quantitatively evaluated through the calculated identification coefficient of the water-consuming zone, and the development level and direction of each layer on each injection-production well pair can be judged by combining the identification criterion of the water-consuming layer. Furthermore, the volume of the water-consuming zone of different levels can also be calculated. Therefore, the method for identifying the water-consumption zone of the oil reservoir in the high water-cut period can accurately determine the position, direction, level and volume of the high water-consumption zone, and improves the characterization effect of the water-consumption zone.
In conclusion, the invention creatively fits the production dynamics of the oil well and the water well in the target area in the high water-cut oil deposit so as to obtain an oil deposit numerical simulation model; then, based on the obtained numerical reservoir simulation model, calculating the identification coefficient of a water consumption zone between each water well and the surrounding oil wells; finally, according to the identification coefficient of the water-consuming zone, different development levels are identified, so that the development levels of the high-water-consuming zone can be quickly judged and quantitatively characterized, the development direction of the high-water-consuming zone can be effectively identified, and an effective guiding effect is played for the design of a regulation and control scheme in the subsequent oil field development stage.
Correspondingly, the invention also provides an identification system of the water-consuming zone of the oil reservoir in the high water-cut period. As shown in fig. 9, the recognition system may include: the fitting device 10 is used for fitting the production dynamics of the oil well and the water well by using a numerical reservoir simulator based on the geological data of a target area in the oil reservoir in the high water cut period and the 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; an 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 numerical reservoir simulation model; and an identification means 30 for identifying a development level of the water-consuming zone based on the identification coefficient of the water-consuming zone.
Preferably, the fitting device 10 may include: the fine geological model building module is used for building a fine geological model of a target region by adopting a geological modeling algorithm based on geological data of the target region in the high water cut oil reservoir; the leading-in module is used for leading the fine geological model into the oil 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 of the oil well and the water well.
Preferably, the identification coefficient calculation means 20 may include: the parameter acquisition module is used for acquiring the effective thickness, water saturation, daily water absorption and daily oil production of each layer of each water well on each injection-production well pair based on the numerical reservoir simulation model; and the water consumption ratio calculation module is used for calculating the ratio of the dimensionless water consumption and the economic water consumption of each layer of each water well on each injection-production well pair based on the daily water absorption and the daily oil production of each layer of each water well on each injection-production well pair and the economic water consumption, wherein the economic water consumption isWater consumption is related to the price of crude oil and the cost of water injection; the water absorption intensity ratio calculation module is used for calculating the ratio of the water absorption intensity of each layer of each well on each injection and production well pair to the water absorption intensity of each well based on the effective thickness of each layer of each well on each injection and production well pair, the daily water absorption amount and the effective thickness of each well; the water saturation ratio calculation module is used for calculating the ratio of the water saturation of each layer of each water well on each injection-production well pair to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection-production well pair; and the identification coefficient calculation module is used for calculating the water consumption zone identification coefficient of each layer of each water well on each injection and production well pair based on the ratio of the dimensionless water consumption and the economic water consumption of each layer on each injection and production well pair, the ratio of the water absorption intensity of each layer on each injection and production well pair to the water absorption intensity of each water well and the ratio of the water saturation of each layer on each injection and production well pair to the average water saturation, 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 device 30 may include: a development level identification module to perform the following operations: identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is greater than or equal to a first preset coefficient, the identification horizon i is a first high-water-consumption horizon on the injection-production well pair k; identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is smaller than the first preset coefficient and is larger than or equal to a second preset coefficient, the identification horizon i is a second high water consumption horizon on the injection-production well pair k; and the identification coefficient C in the water-consuming zoneikAnd under the condition that the identification horizon i is smaller than the second preset coefficient, identifying that the horizon i is a common water-consuming horizon on the injection-production well pair k.
Preferably, the identification system further comprises: and the development direction identification module is used for identifying the direction of an injection and production well pair k corresponding to the first high water consumption zone and/or the second high water consumption zone as the development direction of the high water consumption zone, wherein the injection and production well pair refers to a reservoir between each water well and an oil well around the water well in the control area of each water well.
Preferably, the identification system may further include: the parameter acquisition device is used for acquiring the permeability grade difference, the water content, the crude oil viscosity and the water absorption strength of each water well and each injection and production well pair between each water well and each oil well around each water well based on the numerical reservoir simulation model; and the volume calculation device is used for calculating the volumes of the water consumption 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 each water well and a specific oil well positioned 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 oil reservoir in the high water-cut period of the embodiment of the invention are the same as those of the embodiment of the identification method, and are not repeated herein.
Accordingly, the present invention also provides a machine-readable storage medium having stored thereon instructions for causing a machine to perform the above-described method for identifying a water-depletion zone of a high-hydration reservoir.
The machine-readable storage medium includes, but is not limited to, 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, magnetic disk storage or other magnetic storage devices, and various media capable of storing program code.
The preferred embodiments of the present invention have been described in detail with reference to the accompanying drawings, however, 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 technical idea of the present invention, and these simple modifications are within the protective scope of the present invention.
It should be noted that the various features described in the above embodiments may be combined in any suitable manner without departing from the scope of the invention. The invention is not described in detail in order to avoid unnecessary repetition.
In addition, any combination of the various embodiments of the present invention is also possible, and the same should be considered as the disclosure of the present invention as long as it does not depart from the spirit of the present invention.
Claims (18)
1. A method for identifying a water-consuming zone of a high water-cut reservoir, the method comprising:
fitting the production dynamics of the oil well and the water well by using a numerical reservoir simulator based on the geological data of a target area in the high water cut oil reservoir and the production dynamics data of the oil well and the water well in the target area to obtain a numerical reservoir simulation model;
calculating the identification coefficient of a water consumption zone between each water well and an oil well around the water well based on the numerical reservoir simulation model; and
identifying a development level of the water-consuming zone based on the identification coefficient of the water-consuming zone.
2. The method of identifying water-consuming zone of a high water cut reservoir of claim 1, wherein said fitting production dynamics of said oil well and said water well with a reservoir numerical simulator comprises:
based on geological data of a target region in the high water-cut oil deposit, establishing a fine geological model of the target region by adopting a geological modeling algorithm; and
importing the refined geological model into the reservoir numerical simulator;
and adjusting characteristic parameters in the fine geological model in the oil reservoir numerical simulator by combining the production dynamic data of the oil well and the water well in the target area so as to fit the production dynamic of the oil well and the water well.
3. The method of claim 1, wherein calculating the identification coefficients of the water-consuming zone between each water well and the oil wells surrounding the water well comprises:
based on the numerical reservoir simulation model, obtaining the effective thickness, water saturation, daily water absorption and daily oil production of each layer of each water well on each injection-production well pair;
calculating the ratio of dimensionless water consumption and economic water consumption of each layer of each water well on each injection-production well pair based on the daily water absorption and daily oil production of each layer of each water well on each injection-production well pair and the economic water consumption, 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 on each injection-production well pair to the water absorption intensity of each water well on the basis of the effective thickness of each layer of each water well on each injection-production well pair, the daily water absorption capacity and the effective thickness of each water well;
calculating the ratio of the water saturation of each layer of each water well on each injection-production well pair to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection-production well pair; and
calculating the water consumption zone identification coefficient of each layer of each water well on each injection-production well pair based on the ratio of the dimensionless water consumption and the economic water consumption of each layer on each injection-production well pair, the ratio of the water absorption intensity of each layer on each injection-production well pair to the water absorption intensity of each water well, and the ratio of the water saturation of each layer on each injection-production well pair to the average water saturation,
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.
4. The method of claim 3, wherein the calculating the ratio of dimensionless water consumption to economic water consumption for each horizon for each water well on each injection-production well pair comprises:
by the formulaCalculating the dimensionless water consumption E of the horizon i on the injection-production well pair kikWherein, IikwThe daily water absorption, Q, of the layer i on the injection-production well pair kikoThe daily oil production of the horizon i on the injection-production well pair k;
by the formulaCalculating economic water consumption EwWherein P isoPrice per unit volume of crude oil, CwCost per unit water injection; and
by the formulaCalculating the ratio R of the dimensionless water consumption and the economic water consumption of the horizon i on the injection-production well pair kike。
5. The method of claim 3, wherein calculating the ratio of the water absorption intensity of each layer of each water well on each injection-production well pair to the water absorption intensity of each water well based on the effective thickness of each layer of each water well on each injection-production well pair, the daily water absorption capacity, and the effective thickness of each water well comprises:
by the formulaCalculating the water absorption intensity omega of the horizon i on the injection-production well pair kikWherein, IikwFor daily suction of horizon i on injection-production well pair kAmount of water, hikThe effective thickness of the horizon i on the injection-production well pair k;
by the formulaCalculating the water absorption intensity omega of each water well, wherein I is the water absorption capacity of each water well,h is the effective thickness of each well; and
by the formulaCalculating the ratio R of the water absorption intensity of the horizon i on the injection-production well pair k to the water absorption intensity of each wellikd。
6. The method of identifying water-consuming zone of a high water-cut reservoir of claim 3, wherein the calculating the ratio of water saturation of each horizon over each injection-production well pair to the average water saturation in the target region comprises:
by the formulaCalculating the water saturation S of the horizon i on the injection-production well pair kwikAnd the average water saturationRatio R ofiks。
7. The method for identifying water-consuming zone of high-water-cut reservoir according to claim 3, wherein the calculating the water-consuming zone identification coefficient of the horizon on each injection-production well pair comprises:
calculating the water consumption zone identification coefficient C of the horizon i on the injection-production well pair k by the following formulaik,
Cik=ln(Rike×Riks×Rikd)
Wherein R isikeThe ratio of dimensionless water consumption and economic water consumption of the horizon i on the injection-production well pair k is shown; rikdThe ratio of the water absorption intensity of the horizon i on the injection-production well pair k to the water absorption intensity of each well is shown; and RiksAnd the ratio of the water saturation of the horizon i on the injection-production well pair k to the average water saturation is shown.
8. The method of identifying a water-consuming zone of a high water-cut reservoir as claimed in claim 1, wherein the identifying the level of development of the water-consuming zone comprises:
identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is greater than or equal to a first preset coefficient, the identification horizon i is a first high-water-consumption horizon on the injection-production well pair k;
identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is smaller than the first preset coefficient and is larger than or equal to a second preset coefficient, the identification horizon i is a second high water consumption horizon on the injection-production well pair k; and
identification coefficient C in the water-consuming zoneikAnd under the condition that the identification horizon i is smaller than the second preset coefficient, identifying that the horizon i is a common water-consuming horizon on the injection-production well pair k.
9. The method of identifying water-consuming zones of a high water-cut reservoir as claimed in claim 8, further comprising:
identifying the direction of an injection-production well pair k corresponding to the first high water-consumption zone and/or the second high water-consumption zone as the development direction of the high water-consumption zone,
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.
10. The method of identifying water-consuming zone of a high water-cut reservoir as claimed in claim 1, further comprising:
after performing the step of identifying a developmental level of the water-consuming layer band, performing the following operations:
acquiring the permeability grade difference, the water content, the crude oil viscosity and the water absorption strength of each injection-production well pair between each water well and each oil well around each water well on the basis of the numerical reservoir simulation model; and
calculating the volumes of the water consumption 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,
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.
11. The method of identifying water-consuming zone of a high water-cut reservoir as claimed in claim 10, wherein the development level of the water-consuming zone comprises: ordinary water-consuming layer area reaches high water-consuming layer area, high water-consuming layer area includes: a first high water-consumption layer belt and a second high water-consumption layer belt,
accordingly, the calculating the volumes of the water-consuming zone of different developmental levels comprises:
calculating the volume percentages of the common water-consumption zone, the second high water-consumption zone and the first high water-consumption zone in the injection-production well pair respectively through the following three formulasAnd
wherein x is1Is the permeability level difference of each injection-production well pair, mu is the viscosity of the crude oil, fwThe water content is the above; omega is the water absorption intensity of each water well; and
calculating the volume V of the j-th water-consumption zone of the horizon i on the injection-production well pair k by the following formulaikj,
Wherein j is 1, 2 or 3, and the class 1 water-consuming layer strip corresponds to the common water-consuming layer strip; the class 2 water-consuming layer zone corresponds to the second high water-consuming layer zone; class 3 water-consuming zone corresponds to the first high water-consuming zone, AikThe area of a horizontal transverse plane h of the horizon i in the direction of the injection-production well pair kikIs the effective thickness of horizon i over the pair k of injection and production wells.
12. An identification system for a water-consuming zone of a high water-cut reservoir, the identification system comprising:
the fitting device is used for fitting the production dynamics of the oil well and the water well by using a numerical reservoir simulator based on the geological data of a target area in the oil reservoir in the high water cut period and the production dynamics data of the oil well and the water well in the target area so as to obtain a numerical reservoir simulation model;
the identification coefficient calculation device is used for calculating the identification coefficient of a water consumption zone between each water well and the oil well around the water well on the basis of the numerical reservoir simulation model; and
and the identification device is used for identifying the development level of the water-consuming layer band based on the identification coefficient of the water-consuming layer band.
13. The system for identifying water-consuming zone of a high water-cut reservoir as claimed in claim 12, wherein said fitting means comprises:
the fine geological model building module is used for building a fine geological model of a target region by adopting a geological modeling algorithm based on geological data of the target region in the high water cut oil reservoir; and
the leading-in module is used for leading the fine geological model into the oil 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 of the oil well and the water well.
14. The system for identifying a water-consuming zone of a high water-cut reservoir as claimed in claim 12, wherein the identification coefficient calculating means comprises:
the parameter acquisition module is used for acquiring the effective thickness, water saturation, daily water absorption and daily oil production of each layer of each water well on each injection-production well pair based on the numerical reservoir simulation model;
the water consumption ratio calculation module is used for calculating the ratio of the dimensionless water consumption and the economic water consumption of each layer of each water well on each injection and production well pair based on the daily water absorption and the daily oil yield of each layer of each water well and the economic water consumption, wherein the economic water consumption is related to the price of crude oil and the cost of water injection;
the water absorption intensity ratio calculation module is used for calculating the ratio of the water absorption intensity of each layer of each well on each injection and production well pair to the water absorption intensity of each well based on the effective thickness of each layer of each well on each injection and production well pair, the daily water absorption amount and the effective thickness of each well;
the water saturation ratio calculation module is used for calculating the ratio of the water saturation of each layer of each water well on each injection-production well pair to the average water saturation in the target area based on the water saturation of each layer of each water well on each injection-production well pair; and
an identification coefficient calculation module for calculating the water consumption zone identification coefficient of each layer of each water well on each injection-production well pair based on the ratio of the dimensionless water consumption and the economic water consumption of each layer on each injection-production well pair, the ratio of the water absorption intensity of each layer on each injection-production well pair to the water absorption intensity of each water well, and the ratio of the water saturation of each layer on each injection-production well pair to the average water saturation,
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.
15. The system for identifying a water-consuming zone of a high water-cut reservoir as claimed in claim 12, wherein the identifying means comprises:
a development level identification module to perform the following operations:
identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is greater than or equal to a first preset coefficient, the identification horizon i is a first high-water-consumption horizon on the injection-production well pair k;
identification coefficient C in the water-consuming zoneikUnder the condition that the identification horizon i is smaller than the first preset coefficient and is larger than or equal to a second preset coefficient, the identification horizon i is a second high water consumption horizon on the injection-production well pair k; and
identification coefficient C in the water-consuming zoneikAnd under the condition that the identification horizon i is smaller than the second preset coefficient, identifying that the horizon i is a common water-consuming horizon on the injection-production well pair k.
16. The system for identifying water-consuming zones of a high water-cut reservoir as claimed in claim 12, further comprising:
a development direction identification module used for identifying the direction of the injection-production well pair k corresponding to the first high water-consumption zone and/or the second high water-consumption zone as the development direction of the high water-consumption zone,
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.
17. The system for identifying water-consuming zones of a high water-cut reservoir as claimed in claim 12, further comprising:
the parameter acquisition device is used for acquiring the permeability grade difference, the water content, the crude oil viscosity and the water absorption strength of each water well and each injection and production well pair between each water well and each oil well around each water well based on the numerical reservoir simulation model;
volume calculating means for calculating the volumes of the water-consuming zone of different development levels based on the difference in permeability level, the water content, the viscosity of crude oil, and the water absorption intensity of each water well within each injection-production well pair,
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.
18. A machine-readable storage medium having stored thereon instructions for causing a machine to perform the method for identifying a water-depletion zone of a high-water-cut reservoir according to any one of claims 1-11.
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 true CN110593863A (en) | 2019-12-20 |
CN110593863B 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) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111911135A (en) * | 2020-07-21 | 2020-11-10 | 中国石油大学(华东) | Dynamic description method for high water consumption strip of water-drive reservoir |
CN112632864A (en) * | 2020-12-14 | 2021-04-09 | 中国石油化工股份有限公司胜利油田分公司孤岛采油厂 | Method for identifying high-water-consumption zone in mine by dynamic and static combination |
CN113297740A (en) * | 2021-05-26 | 2021-08-24 | 中国石油大学(华东) | Horizontal well gel profile control parameter optimization method |
CN113505507A (en) * | 2021-06-29 | 2021-10-15 | 中国石油化工股份有限公司 | Method for identifying high water consumption zone between injection and production wells |
CN113537585A (en) * | 2021-07-09 | 2021-10-22 | 中海石油(中国)有限公司天津分公司 | Oil field production increasing measure recommendation method based on random forest and gradient lifting decision tree |
CN114763736A (en) * | 2021-01-13 | 2022-07-19 | 中国石油化工股份有限公司 | Method for treating high water consumption of fault block oil reservoir |
CN114810006A (en) * | 2021-01-27 | 2022-07-29 | 中国石油化工股份有限公司 | Potential evaluation method for regulating and controlling high-water-consumption zone by separate-layer water injection after heterogeneous flooding |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017212056A1 (en) * | 2016-06-10 | 2017-12-14 | 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 |
CN108301813A (en) * | 2017-12-20 | 2018-07-20 | 中国石油化工股份有限公司 | The multilayer sandstone oil reservoir modification scenario method of scattered sand body development |
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 |
-
2019
- 2019-09-16 CN CN201910870315.4A patent/CN110593863B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017212056A1 (en) * | 2016-06-10 | 2017-12-14 | 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 |
CN108301813A (en) * | 2017-12-20 | 2018-07-20 | 中国石油化工股份有限公司 | The multilayer sandstone oil reservoir modification scenario method of scattered sand body development |
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 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111911135A (en) * | 2020-07-21 | 2020-11-10 | 中国石油大学(华东) | Dynamic description method for high water consumption strip of water-drive reservoir |
CN111911135B (en) * | 2020-07-21 | 2021-06-25 | 中国石油大学(华东) | Dynamic description method for high water consumption strip of water-drive reservoir |
CN112632864A (en) * | 2020-12-14 | 2021-04-09 | 中国石油化工股份有限公司胜利油田分公司孤岛采油厂 | Method for identifying high-water-consumption zone in mine by dynamic and static combination |
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 |
CN114810006A (en) * | 2021-01-27 | 2022-07-29 | 中国石油化工股份有限公司 | Potential evaluation method for regulating and controlling high-water-consumption zone by separate-layer water injection after heterogeneous flooding |
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 |
CN113505507A (en) * | 2021-06-29 | 2021-10-15 | 中国石油化工股份有限公司 | Method for identifying high water consumption zone between injection and production wells |
CN113505507B (en) * | 2021-06-29 | 2022-09-09 | 中国石油化工股份有限公司 | Method for identifying high water consumption zone between injection and production wells |
CN113537585A (en) * | 2021-07-09 | 2021-10-22 | 中海石油(中国)有限公司天津分公司 | Oil field production increasing measure recommendation method based on random forest and gradient lifting decision tree |
Also Published As
Publication number | Publication date |
---|---|
CN110593863B (en) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110593863A (en) | Identification method and identification system for water-consuming zone of high-water-cut oil reservoir | |
CN104790926B (en) | A kind of fracture-pore reservoir waterflooding development effect evaluation method | |
CN102041995B (en) | System for monitoring complicated oil deposit flooding conditions | |
CN105134191B (en) | Method for evaluating reserves of tight oil well | |
CN111695303A (en) | Method for evaluating water filling strength of coal seam roof sandstone aquifer | |
CN103075142A (en) | Method for selecting water plugging oil well during development of oilfield through water injection | |
CN107832574B (en) | Horizontal well water flooded layer interpretation method based on logging while drilling | |
CN106372810A (en) | TBM (tunnel boring machine) quantitative type selection method by use of efficacy coefficient method | |
CN105404972A (en) | Method for reservoir development uncertainty research and risk control | |
CN107130955B (en) | Method for determining bottom hole flowing pressure and method for determining natural energy of reservoir | |
CN108194053B (en) | method and device for determining matrix acidizing acid liquid dosage of open-hole horizontal well | |
CN113803029B (en) | Method and device for adjusting coalbed methane well pattern and computer readable storage medium | |
CN108843300B (en) | Method and device for determining type of main flow channel in complex porous medium | |
CN114897436B (en) | Stratum pressure balancing method and controller for gas storage | |
CN113901546B (en) | Tunnel water inflow calculation method and device | |
CN113032996B (en) | Water channeling channel identification method for hypotonic fractured reservoir horizontal well | |
CN109973068B (en) | Method and device for identifying oil reservoir water injection induced cracks | |
US11574083B2 (en) | Methods and systems for selecting inflow control device design simulations based on case selection factor determinations | |
CN114060015A (en) | Method and device for evaluating gas content of tight sandstone | |
Zakirov et al. | In-situ determination of displacement efficiency and oil and water relative permeability curves through integrated well test study at exploration-to-pilot stage of the oilfield development project | |
CN113704950B (en) | Method and device for standardized evaluation of reservoir gas content through gas logging of all-hydrocarbon | |
CN105528732B (en) | Method for predicting productivity of gas testing well | |
CN116950637B (en) | Single well group profile control method and system based on micro-pressure drive and readable storage medium | |
CN117432399B (en) | Method for calculating recoverable resource of geothermal fluid with pulse heat storage | |
CN110442950B (en) | Multi-factor fusion residual oil distribution prediction graph compiling method |
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 |