CN112329232B - Fracture-cavity type oil reservoir production dynamic characterization method, device, equipment and storage medium - Google Patents
Fracture-cavity type oil reservoir production dynamic characterization method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN112329232B CN112329232B CN202011214880.4A CN202011214880A CN112329232B CN 112329232 B CN112329232 B CN 112329232B CN 202011214880 A CN202011214880 A CN 202011214880A CN 112329232 B CN112329232 B CN 112329232B
- Authority
- CN
- China
- Prior art keywords
- oil
- fracture
- cavity
- water injection
- periodic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 59
- 238000012512 characterization method Methods 0.000 title claims abstract description 23
- 238000003860 storage Methods 0.000 title claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 115
- 239000007924 injection Substances 0.000 claims abstract description 62
- 238000002347 injection Methods 0.000 claims abstract description 62
- 230000000737 periodic effect Effects 0.000 claims abstract description 42
- 238000011084 recovery Methods 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 27
- 239000003129 oil well Substances 0.000 claims abstract description 23
- 230000006870 function Effects 0.000 claims description 46
- 230000015654 memory Effects 0.000 claims description 29
- 238000004590 computer program Methods 0.000 claims description 18
- 239000013598 vector Substances 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 8
- 239000003921 oil Substances 0.000 description 94
- 238000010586 diagram Methods 0.000 description 15
- 239000010779 crude oil Substances 0.000 description 10
- 238000011161 development Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000011158 quantitative evaluation Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 3
- 235000019994 cava Nutrition 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000084 colloidal system Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 230000020477 pH reduction Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Human Resources & Organizations (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the specification provides a fracture-cavity type oil reservoir production dynamic characterization method, a device, equipment and a storage medium, wherein the method comprises the following steps: determining the type of a fracture-cavity unit to which the target single well belongs; establishing a relation model of the periodic comprehensive water content and the periodic enhanced oil recovery amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized Darcy's law; constructing a minimized objective function according to the relation model; solving the minimized objective function by taking the observed value of the periodic comprehensive water content of the target single well as input to obtain a determined relation model; and taking the relation curve corresponding to the determined relation model as a production dynamic rule of the oil well corresponding to the type of the fracture-cavity unit under the condition of multiple times of water injection for oil replacement. The embodiment of the specification can improve the accuracy of dynamic representation of the fracture-cavity oil reservoir by performing multiple rounds of water injection instead of oil production.
Description
Technical Field
The specification relates to the technical field of fracture-cavity carbonate reservoir (hereinafter referred to as fracture-cavity reservoir) development, in particular to a method, a device, equipment and a storage medium for dynamic characterization of fracture-cavity reservoir production.
Background
Fracture-cavity reservoirs are generally special types of reservoirs with large karst caves and fractures as the main reservoir spaces. The fracture-cavity oil reservoir has various reservoir space types, large difference of single well control reserves, serious heterogeneity and abnormal and complex oil-water relationship.
The front-stage exploitation of the fracture-cavity oil reservoir mainly depends on elastic energy and natural water flooding development. Aiming at the problems of rapid yield decrease and low recovery rate amplitude caused by insufficient energy and rising water content of part of fracture-vug reservoirs, the prior art forms a later-stage exploitation technology mainly using multiple rounds of water injection for replacing oil so as to improve the recovery rate of crude oil. However, due to the different adaptability to different types of fracture-cavity reservoirs, the water injection and oil increase effects are gradually deteriorated with the increase of water injection amount, the number of failure wells is obviously increased, and a large amount of residual oil is still remained at the positions around the wells and is not used. Research shows that the dynamic law of oil replacement production by multiple rounds of water injection is accurately represented, and the method has important significance for making reasonable optimization adjustment strategies subsequently, restoring the productivity of oil wells and fully excavating used geological reserves. However, the lack of a quantitative representation of the dynamic law of oil replacement by water injection for multiple rounds has become an urgent technical scheme at present, so that the accurate representation of the dynamic law of oil replacement by water injection for multiple rounds is difficult.
Disclosure of Invention
The embodiment of the specification aims to provide a fracture-cavity oil reservoir production dynamic representation method, device, equipment and storage medium, so as to improve the accuracy of the fracture-cavity oil reservoir production dynamic representation by multiple rounds of water injection for oil replacement.
In order to achieve the above object, in one aspect, the present specification provides a fracture-cavity reservoir production dynamic characterization method, including:
determining the type of a fracture-cavity unit to which the target single well belongs;
establishing a relation model of the periodic comprehensive water content and the periodic enhanced oil recovery amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized Darcy's law;
constructing a minimized objective function according to the relation model;
taking the observed value of the periodic comprehensive water content of the target single well as an input, solving the minimized objective function, and obtaining a determined relation model;
and taking the relation curve corresponding to the determined relation model as a production dynamic rule of the oil well corresponding to the type of the fracture-cavity unit under the condition of multiple times of water injection for oil replacement.
In one method embodiment of the present description, the relational model is represented as
Wherein f is wp The water content is periodically integrated, delta R is the periodic enhanced recovery ratio amplitude, a, b, c and d are parameters to be solved, and e is a natural number.
In one method embodiment of the present specification, the minimization objective function is expressed as
Wherein J represents an objective function; u represents a control parameter vector to be solved;
respectively represents t in j-th water injection oil replacement k The predicted value and the observed value of the water content are synthesized in the period of time; n is a radical of t Representing the water injection oil replacement turns.
In one embodiment of the method of the present specification, wherein the ith component u of u i Is shown as
Wherein v is i Represents u i The value of the control parameter variable after the logarithmic transformation,
respectively represents u i Minimum and maximum values of (d), u i Satisfy the requirement of
i＝1,2,…,N u ，N u Representing the dimension of u.
In another aspect, an embodiment of the present specification further provides a fracture-cavity reservoir production dynamic characterization device, including:
the type determining module is used for determining the type of the fracture-cavity unit to which the target single well belongs;
the model establishing module is used for establishing a relation model of the periodic comprehensive water content and the periodic enhanced oil recovery amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized Darcy's law;
the function building module is used for building a minimized objective function according to the relation model;
the function solving module is used for solving the minimum objective function by taking the observed value of the periodic comprehensive water content of the target single well as input to obtain a determined relation model;
and the rule determining module is used for taking the relation curve corresponding to the determined relation model as a dynamic production rule of the oil well corresponding to the type of the fracture-cavity unit under the condition of replacing oil by multiple times of water injection.
In one embodiment of the apparatus of the present specification, the relational model is expressed as
Wherein f is wp The water content is periodically integrated, delta R is the period enhanced recovery ratio amplitude, a, b, c and d are parameters to be solved, and e is a natural number.
In one embodiment of the apparatus of the present specification, the minimization objective function is expressed as
Wherein J represents an objective function; u tableDisplaying a control parameter vector to be solved;
respectively represents t in j-th water injection oil replacement k The predicted value and the observed value of the water content are synthesized in the period of time; n is a radical of hydrogen t Representing the water injection oil replacement turns.
In an embodiment of the present specification, wherein the ith component u of u i Is shown as
Wherein v is i Represents u i The value of the control parameter variable after the logarithmic transformation,
respectively represents u i Minimum and maximum values of (d), u i Satisfy the requirement of
i＝1,2,…,N u ，N u Representing the dimension of u.
In another aspect, the embodiments of the present specification further provide a computer device, which includes a memory, a processor, and a computer program stored on the memory, and when the computer program is executed by the processor, the computer program executes the instructions of the above method.
In another aspect, the present specification further provides a computer storage medium, on which a computer program is stored, and the computer program is executed by a processor to execute the instructions of the method.
It can be seen from the technical solutions provided in the embodiments of the present specification that, in the embodiments of the present specification, the relationship curve reflects a quantitative relationship between the periodic comprehensive water content and the periodic enhanced oil recovery rate amplitude, and the quantitative relationship between the periodic comprehensive water content and the periodic enhanced oil recovery rate amplitude quantitatively represents a dynamic production rule of an oil well corresponding to the type of the fracture-cavity unit under oil replacement by multiple times of water injection, and subsequently, for any oil well in the later stage of exploitation in the fracture-cavity oil reservoir, when the type of the fracture-cavity unit is the same as the type of the fracture-cavity unit to which the target single well belongs, the dynamic production rule of the oil well corresponding to the target single well under oil replacement by multiple times of water injection can be represented by the relationship curve corresponding to the target single well. Therefore, for the fracture-cavity type oil reservoir, quantitative evaluation of the production dynamic rule of the oil well corresponding to each type of fracture-cavity unit under multiple times of water injection for oil replacement can be realized by adopting the embodiment of the specification, so that the accuracy of the multiple times of water injection for oil replacement production dynamic representation of the fracture-cavity type oil reservoir is improved. And then, applying the quantitative evaluation result to field practice, and guiding the multi-round water injection oil replacement differential comprehensive treatment of different types of fracture-cavity units and the formulation of measures for improving the recovery ratio, so that the potential of residual development is exploited to the maximum extent, and the crude oil recovery ratio of the fracture-cavity oil reservoir is improved.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a flow diagram of a method for dynamic characterization of fracture-cavity reservoir production in some embodiments of the present description;
FIG. 2 is a schematic diagram illustrating the dynamic law of production of cave-type slot-hole units in an exemplary embodiment of the present disclosure;
FIG. 3 is a schematic diagram illustrating the dynamic behavior of hole-cavity slot-hole unit production in an exemplary embodiment of the present disclosure;
FIG. 4 is a schematic diagram illustrating the dynamic behavior of hole-cavity slot-hole unit production in another exemplary embodiment of the present disclosure;
FIG. 5 is a schematic diagram illustrating the production dynamics of a fracture-hole type fracture-hole unit in an exemplary embodiment of the present disclosure;
FIG. 6 is a block diagram of a dynamic characterization apparatus for fracture-cavity reservoir production in some embodiments of the present description;
FIG. 7 is a block diagram of a computer device in some embodiments of the present description.
[ description of reference ]
61. A type determination module;
62. a model building module;
63. a function building module;
64. a function solving module;
65. a rule determining module;
702. a computer device;
704. a processor;
706. a memory;
708. a drive mechanism;
710. an input/output module;
712. an input device;
714. an output device;
716. a presentation device;
718. a graphical user interface;
720. a network interface;
722. a communication link;
724. a communication bus.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in the embodiments of the present specification will be clearly and completely described below with reference to the drawings in the embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments of the present specification, and not all of the embodiments. All other embodiments obtained by a person skilled in the art based on the embodiments in the present specification without any inventive step should fall within the scope of protection of the present specification.
Referring to fig. 1, a fracture-cavity reservoir production dynamic characterization method of some embodiments of the present description may include the following steps:
s101, determining the type of a fracture-cavity unit to which the target single well belongs.
For fracture-cavity oil reservoirs, the water drive production dynamic rules of different types of fracture-cavity units are obviously different. Therefore, in order to be beneficial to accurately characterizing the production dynamics of the fracture-cavity type oil reservoir, the fracture-cavity unit type of each single well needs to be determined so as to characterize the corresponding reservoir body type.
In some embodiments of the present disclosure, the fracture-cavity unit types (e.g., cavity type, hole-cavity type, fracture-hole type, etc.) may be classified and evaluated according to dynamic and static data of multiple rounds of water injection replacement wells in the fracture-cavity oil reservoir. Static data may include, for example, log response, seismic reflection, well test interpretation, production instability analysis, and depleted production data, among others. The dynamic data can comprise dynamic data such as the cycle exploitation time, daily oil production, peak water content, cumulative oil production, cumulative water production, dynamic geological reserve and water body size of different types of fracture-cave units during multiple rounds of water injection for oil replacement. For example, in one embodiment, the slot cell type division criteria shown in table 1 below may be formed according to the dynamic and static data.
TABLE 1
In the embodiment of the specification, the target single well is a multi-round water injection replacement oil well in the fracture-cavity oil reservoir, namely the target single well is an oil well in the fracture-cavity oil reservoir in the later exploitation stage. After the slot and hole units of different types are divided according to the dynamic data, the comprehensive water content f of the period can be calculated according to the dynamic data wp And a periodic enhanced recovery margin Δ R. And successively recurrently, determining the types of the fracture-cavity units to which all the oil wells in the later mining period in the working area belong, and determining the corresponding periodic comprehensive water content f wp And periodic enhanced recovery amplitude Δ R data.
S102, establishing a relation model of the periodic comprehensive water content and the periodic enhanced oil recovery amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized Darcy' S law.
In some embodiments of the description, an improved Craft empirical model can be used for describing the functional relationship between the oil-water relative permeability ratio and the water saturation of the water flooding oil reservoir (the core of the improved Craft empirical model is that the logarithm of the ratio of the oil-water relative permeability can be approximated by a cubic polynomial function of the water saturation), and the generalized Darcy's law and the physical meaning of the flow split (the flow split is equivalent to the water content, and the physical meaning is the percentage of the water yield in the water yield) are combined to respectively deduce the periodic comprehensive water content f of the fracture-cavity oil reservoir multi-round water flooding oil well wp And the calculation expression of the periodic enhanced recovery factor amplitude delta R can be specifically described as follows:
in the formula, Q wp 、Q op Respectively the periodic water yield and the periodic oil yield, K ro 、K rw Relative permeability of oil phase and water phase respectively; mu.s o 、μ w Oil phase and water phase viscosity respectively; v p Is the pore volume; rho o Is the crude oil density; n is dynamic geological reserves; b is o Is the volume factor of the crude oil under the formation condition; s w0 Water saturation before oil is replaced by single well water injection; s. the w And (5) the water saturation of the single well after water injection and oil replacement.
Research shows that the periodic increase of recovery ratio amplitude delta R and the change S of water saturation w -S w0 The positive correlation is formed, for a specific slot hole unit, the physical parameters of other reservoirs are fixed and unchanged, and the comprehensive water content f of a multi-cycle water injection oil replacement cycle is established according to the fixed and unchanged physical parameters wp And (3) a relation model of the periodic enhanced oil recovery amplitude delta R, wherein the expression is as follows:
wherein f is wp The water content is periodically integrated, delta R is the period enhanced recovery ratio amplitude, a, b, c and d are parameters to be solved, and e is a natural number.
As will be understood by those skilled in the art, the cycle integrated moisture content refers to: average moisture content in the current period; the periodic enhanced recovery ratio amplitude is as follows: the magnitude of change in recovery from the current cycle relative to the recovery from the previous cycle. In some embodiments of the present disclosure, the period may be a day, in which case the periodic enhanced recovery rate is the percentage increase of the recovery rate of the day relative to the recovery rate of the day before. The periodic comprehensive water content is the average water content of the current day. Of course, in other embodiments of the present disclosure, the period may also be any other suitable time period (e.g., week, month, etc.).
S103, constructing a minimized objective function according to the relation model.
In order to realize quantitative evaluation of the dynamic rule of oil replacement production by multiple rounds of water injection of a fracture-cavity oil reservoir, in some embodiments of the present description, a minimum objective function reflecting the sum of squares of errors between a predicted value and an observed value under the influence of hybrid constraints can be established according to a least square theory that a predicted value of dynamic data is matched with an observed value (i.e., an actual measured value), and an expression of the minimum objective function is as follows:
and satisfies the constraint conditions:
i＝1,2,…,N u
wherein J represents an objective function; u represents a control parameter vector to be solved;
respectively represents t in j-th water injection oil replacement k The predicted value and the observed value of the water content are synthesized in the period of time (the observed value is the actually observed value); n is a radical of t Representing the water injection oil replacement turn; u. of i The ith component of u;
respectively represents u i Minimum and maximum values of; n is a radical of hydrogen u Denotes the dimension of u, when u is [ a, b, c, d ]]When, N u ＝4。
It should be understood that the above-mentioned minimization objective function is only an example, and in other embodiments of the present disclosure, the minimization objective function may be used to reflect the difference between the predicted value and the observed value. The specification is not limited to this, and may be selected as needed.
In other embodiments of the present disclosure, the above-mentioned minimization objective function may further introduce a logarithmic transformation to eliminate the boundary constraint of the control parameter vector, so that the hybrid constraint optimization problem may be converted into an unconstrained optimization problem. The method of logarithmic transformation may be as follows:
in the formula, v i Represents u i And (4) carrying out logarithmic transformation on the control parameter variable values.
And S104, solving the minimized objective function by taking the observed value of the periodic comprehensive water content of the target single well as input to obtain a determined relation model.
Any suitable optimization algorithm may be used in some embodiments of the present description to solve the minimized objective function to fit the determined relational model.
For example, in an embodiment, a random Approximate Gradient (StoSAG) optimization algorithm may be used to perform automatic history fitting on dynamic data of oil replacement production by multiple rounds of water injection, to obtain a control parameter vector to be solved, that is, values of a, b, c, and d in the relationship model are determined by fitting, and when a, b, c, and d in the relationship model are determined, the relationship model also forms a determined relationship curve. Compared with gradient algorithms, the StoSAG algorithm has unique advantages in the aspects of calculation time consumption, avoidance of trapping in local optimal solution, processing constraint conditions and the like. When the StoSAG optimization solution is carried out, judging whether a convergence condition is met or not in real time, and if yes, terminating the circulation; otherwise, updating the control parameters and continuing the next iteration. The optimal control vector expression of the StoSAG iterative computation is as follows:
in the formula u k 、u k+1 Respectively representing the optimal control vectors of the kth iteration and the kth +1 th iteration; a is k Is a linear search step length; | d | k || ∞ Is the infinite norm of the vector; d k For controlling the parameter vector u k And the objective function J (u) k ) For approximating the gradient values of the objective function.
And S105, taking the relation curve corresponding to the determined relation model as a dynamic production rule of the oil well corresponding to the type of the fracture-cavity unit under the condition of replacing oil by multiple times of water injection.
It has been clarified above that when a, b, c and d are determined in the relation model, the relation model also forms a determined relation curve. For example, when a is 1, b is 1, c is 1, and d is 4, the relationship can be obtained as:
it is obvious that the relation reflects f wp And Δ R, and f wp And the quantitative relation between the delta R and the delta R quantitatively represents the production dynamic rule of the oil well corresponding to the fracture-cavity unit type under the condition of multiple times of water injection for replacing oilAnd subsequently, for any oil well in the fracture-cavity oil reservoir at the later exploitation stage, when the fracture-cavity unit type of the oil well is the same as that of the target single well, representing the production dynamic rule of the oil well under the condition of replacing oil by multiple times of water injection by using a corresponding relation curve of the target single well. Therefore, for the fracture-cavity type oil reservoir, quantitative evaluation of the production dynamic rule of the oil well corresponding to each type of fracture-cavity unit under multiple times of water injection for oil replacement can be realized by adopting the embodiment of the specification, so that the accuracy of the multiple times of water injection for oil replacement production dynamic representation of the fracture-cavity type oil reservoir is improved.
And then, applying the quantitative evaluation result to field practice, and guiding the multi-round water injection oil replacement differential comprehensive treatment of different types of fracture-cavity units and the formulation of measures for improving the recovery ratio, so that the potential of residual development is exploited to the maximum extent, and the crude oil recovery ratio of the fracture-cavity oil reservoir is improved.
For ease of understanding, an application example of the fracture-cavity reservoir production dynamic characterization method of the present specification is described below.
The main effective storage space of the M-fracture-cavity type carbonate rock oil reservoir of the Tarim basin comprises large-scale caves, secondary corrosion caves, cracks and the like, the oil reservoir belongs to an unsaturated oil reservoir, the earth saturation pressure difference is large (approximately 35-66 MPa), the average pressure is 50.5MPa, the formation crude oil mainly comprises light crude oil with low viscosity, low sulfur, medium colloid asphaltene content and high wax content, and the density is 0.82-1.10 g/cm 3 The volume coefficient is 1.04-1.66, the average formation pressure is 73.3MPa, the average formation temperature is 154.4 ℃, and the system is a normal temperature-pressure system. According to the dynamic and static data of 30 oil wells with the number of times of oil replacement by water injection of the M-fracture-cavity type oil reservoir being more than 5, on the basis of accurately dividing the type of a reservoir space (namely the type of a fracture-cavity unit), as shown in the table 1, the dynamic law of production by water injection for oil replacement for multiple times of different types of fracture-cavity units can be described by applying the production dynamic quantitative characterization method (for example, as shown in fig. 2 to 5).
As shown in fig. 2, for cavernous reservoirs (i.e., cavernous slotted cells), f wp The relation curve between delta R and delta R is mostly characterized as a concave characteristic, so a reasonable working system is formulated to prolong the anhydrous oil extraction period as much as possible; once the cover is closedWhen water is leaked, the development effect of replacing oil by water injection is rapidly deteriorated. And at the later stage, nitrogen gas injection from the top can be recommended to displace attic oil so as to further improve the crude oil recovery.
For cavern-cavern reservoirs (i.e. cavern-cavern slot cells), f wp The Δ R curve is characterized as either an "S" type (as shown in FIG. 3) or a "step up" (as shown in FIG. 4). For such reservoirs, a long anhydrous oil production period exists, and a large amount of residual oil is produced after water breakthrough; the residual development potential is fully excavated by combining measures such as drainage oil extraction, well closing and cone pressing and the like and by technologies such as multiple rounds of water injection oil replacement, nitrogen injection huff and puff and the like.
When f is wp When the relation curve of the delta R and the delta R is expressed as a hole-cave type reservoir body (namely a hole-cave type slot hole unit) with an S-shaped characteristic, for the reservoir body, the waterless oil production period is longer, and a large amount of residual oil is produced after the water content is broken through; therefore, measures such as drainage oil extraction, well closing and cone pressing are combined, and residual development potential is fully excavated through technologies such as multiple times of water injection for oil replacement, nitrogen injection for huff and puff and the like, so that the crude oil recovery rate is improved.
When f is shown in FIG. 5 wp When the relation curve between the delta R and the delta R is expressed as a crack-hole type reservoir body (namely a crack-hole type crack-hole unit) with a downward convex characteristic, for the reservoir body, the waterless oil extraction period is short, the effect of improving the recovery ratio is poor, and the development effect is improved by considering such excavation potential measures as repeated acidification, upward returning layer modification, hydraulic expansion, gas injection huff and puff and the like, so that the crude oil recovery ratio is improved.
According to the application of the method to each large fracture-cavity carbonate rock oil and gas field in the Tarim basin, the evaluation result of the dynamic characterization method for the fracture-cavity oil reservoir production in the embodiment of the specification is high in reliability, and theoretical basis can be provided for making different comprehensive treatment and recovery rate improvement countermeasures for different types of fracture-cavity units by means of multi-turn water injection replacement oil wells.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
Corresponding to the fracture-cavity oil reservoir production dynamic characterization method, the specification also provides an embodiment of a fracture-cavity oil reservoir production dynamic characterization device. Referring to fig. 6, in some embodiments of the present disclosure, a fracture-cavity reservoir production dynamic characterization apparatus may include:
the type determining module 61 may be configured to determine a type of a fracture-cavity unit to which the target single well belongs.
The model establishing module 62 may be configured to establish a relation model between the periodic comprehensive water content and the periodic enhanced oil recovery rate amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized darcy law.
A function building module 63, which may be configured to build a minimization objective function according to the relational model.
And the function solving module 64 may be configured to solve the minimized objective function by using the observed value of the periodic comprehensive water content of the target single well as an input, so as to obtain a determined relationship model.
The rule determining module 65 may be configured to use the relationship curve corresponding to the determined relationship model as a dynamic rule for production of the oil well corresponding to the type of the fracture-cavity unit under the condition of oil replacement by multiple rounds of water injection.
In the dynamic characterization device for fracture-cavity reservoir production of some embodiments of the present specification, the relationship model is expressed as
Wherein f is wp The water content is periodically integrated, delta R is the period enhanced recovery ratio amplitude, a, b, c and d are parameters to be solved, and e is a natural number.
In the dynamic characterization device for fracture-cavity reservoir production of some embodiments of the present specification, the minimization objective function is expressed as
Wherein J represents an objective function; u represents a control parameter vector to be solved;
respectively represents t in j-th water injection oil replacement k The predicted value and the observed value of the water content are synthesized in the period of time; n is a radical of t Representing the water injection oil replacement turns.
In the fracture-cavity reservoir production dynamic characterization device of some embodiments of the present specification, wherein the ith component u of u i Is shown as
Wherein v is i Represents u i The value of the control parameter variable after the logarithmic transformation,
respectively represents u i Minimum and maximum values of (d), u i Satisfy the requirement of
i＝1,2,…,N u ，N u Representing the dimension of u.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the various elements may be implemented in the same one or more pieces of software and/or hardware in the practice of this description.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
Corresponding to the fracture-cavity oil reservoir production dynamic characterization method, the specification also provides computer equipment. As shown in fig. 7, the computer device 702 may include one or more processors 704, such as one or more Central Processing Units (CPUs) or Graphics Processors (GPUs), each of which may implement one or more hardware threads. The computer device 702 may also include any memory 706 for storing any kind of information, such as code, settings, data, etc., and in a particular embodiment, a computer program on the memory 706 and executable on the processor 704, which computer program when executed by the processor 704 may perform instructions according to the above-described method. For example, and without limitation, the memory 706 can include any one or more of the following in combination: any type of RAM, any type of ROM, flash memory devices, hard disks, optical disks, etc. More generally, any memory may use any technology to store information. Further, any memory may provide volatile or non-volatile retention of information. Further, any memory may represent fixed or removable components of computer device 702. In one case, when the processor 704 executes associated instructions that are stored in any memory or combination of memories, the computer device 702 can perform any of the operations of the associated instructions. The computer device 702 also includes one or more drive mechanisms 708, such as a hard disk drive mechanism, an optical disk drive mechanism, or the like, for interacting with any memory.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the specification. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), 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 or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
As will be appreciated by one skilled in the art, embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, embodiments of the present description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, embodiments of the present description may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and so forth) having computer-usable program code embodied therein.
The embodiments of this specification may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The described embodiments may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of an embodiment of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (6)
1. A fracture-cavity oil reservoir production dynamic characterization method is characterized by comprising the following steps:
determining the type of a fracture-cave unit to which the target single well belongs;
establishing a relation model of the periodic comprehensive water content and the periodic enhanced oil recovery amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized Darcy's law;
constructing a minimized objective function according to the relation model;
solving the minimized objective function by taking the observed value of the periodic comprehensive water content of the target single well as input to obtain a determined relation model;
taking the relation curve corresponding to the determined relation model as a production dynamic rule of the oil well corresponding to the type of the fracture-cave unit under the condition of multiple rounds of water injection and oil replacement;
wherein the relational model is represented as
Wherein f is wp The periodic comprehensive water content is obtained, delta R is the periodic enhanced oil recovery rate amplitude, a, b, c and d are parameters to be solved, and e is a natural number;
wherein the minimization objective function is expressed as
Wherein J represents an objective function; u represents a control parameter vector to be solved;
respectively represents t in j-th water injection oil replacement k The predicted value and the observed value of the water content are synthesized in the period of time; n is a radical of t Representing the water injection oil replacement turns.
2. The fracture-cavity reservoir production dynamic characterization method of claim 1, wherein the ith component u of u, u i Is shown as
Wherein v is i Represents u i The value of the control parameter variable after the logarithmic transformation,
respectively represent u i Minimum and maximum values of (d), u i Satisfy the requirement of
N u ，N u Dimension representing u。
3. A fracture-cavity type oil reservoir production dynamic characterization device is characterized by comprising:
the type determining module is used for determining the type of the fracture-cavity unit to which the target single well belongs;
the model establishing module is used for establishing a relation model of the periodic comprehensive water content and the periodic enhanced oil recovery amplitude of the target single well under multiple rounds of water injection and oil replacement according to the generalized Darcy's law;
the function building module is used for building a minimized objective function according to the relation model;
the function solving module is used for solving the minimum objective function by taking the observed value of the periodic comprehensive water content of the target single well as input to obtain a determined relation model;
the rule determining module is used for taking the relation curve corresponding to the determined relation model as a dynamic production rule of the oil well corresponding to the type of the fracture-cavity unit under the condition of replacing oil by multiple times of water injection;
wherein the relational model is represented as
Wherein, f wp The method is characterized in that the method is a periodic comprehensive water content, delta R is a periodic enhanced recovery ratio amplitude, a, b, c and d are parameters to be solved, and e is a natural number;
wherein the minimization objective function is expressed as
Wherein J represents an objective function; u represents a control parameter vector to be solved;
respectively represents t in j-th water injection oil replacement k The predicted value and the observed value of the water content are synthesized in the period of time; n is a radical of t Representing the water injection oil replacement turns.
4. The fracture-cavity reservoir production dynamic characterization device of claim 3 wherein the ith component u of u, u i Is shown as
Wherein v is i Represents u i The value of the control parameter variable after the logarithmic transformation,
respectively represents u i Minimum and maximum values of (d), u i Satisfy the requirement of
N u ，N u Representing the dimension of u.
5. A computer device comprising a memory, a processor, and a computer program stored on the memory, characterized in that the computer program, when executed by the processor, executes instructions of the method according to claim 1 or 2.
6. A computer storage medium on which a computer program is stored, characterized in that the computer program, when being executed by a processor, executes instructions for a method according to claim 1 or 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011214880.4A CN112329232B (en) | 2020-11-04 | 2020-11-04 | Fracture-cavity type oil reservoir production dynamic characterization method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011214880.4A CN112329232B (en) | 2020-11-04 | 2020-11-04 | Fracture-cavity type oil reservoir production dynamic characterization method, device, equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112329232A CN112329232A (en) | 2021-02-05 |
| CN112329232B true CN112329232B (en) | 2022-09-20 |
Family
ID=74323599
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011214880.4A Active CN112329232B (en) | 2020-11-04 | 2020-11-04 | Fracture-cavity type oil reservoir production dynamic characterization method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112329232B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115577562A (en) * | 2022-11-09 | 2023-01-06 | 中国石油大学(华东) | Fractured reservoir well position optimization method |
| CN118094832A (en) * | 2024-02-02 | 2024-05-28 | 中国石油大学(北京) | Oil reservoir development evaluation method, device and storage medium based on data driving |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104747180A (en) * | 2014-02-26 | 2015-07-01 | 中国石油化工股份有限公司 | Fracture-cavity type oil deposit analyzing method for water-flooding development and application thereof |
| CN105626009A (en) * | 2015-12-17 | 2016-06-01 | 中国石油天然气股份有限公司 | Quantitative evaluation method for single-well water injection oil replacement effect of fracture-cave carbonate reservoir |
| CN105719097A (en) * | 2016-01-27 | 2016-06-29 | 中国石油化工股份有限公司 | Dynamic analysis and waterflooding management system of fractured-vuggy type reservoir |
-
2020
- 2020-11-04 CN CN202011214880.4A patent/CN112329232B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104747180A (en) * | 2014-02-26 | 2015-07-01 | 中国石油化工股份有限公司 | Fracture-cavity type oil deposit analyzing method for water-flooding development and application thereof |
| CN105626009A (en) * | 2015-12-17 | 2016-06-01 | 中国石油天然气股份有限公司 | Quantitative evaluation method for single-well water injection oil replacement effect of fracture-cave carbonate reservoir |
| CN105719097A (en) * | 2016-01-27 | 2016-06-29 | 中国石油化工股份有限公司 | Dynamic analysis and waterflooding management system of fractured-vuggy type reservoir |
Non-Patent Citations (1)
| Title |
|---|
| 缝洞型碳酸盐岩油藏衰竭开采特征综合聚类评价;王代刚等;《新疆石油地质》;20160415(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112329232A (en) | 2021-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Le et al. | An adaptive ensemble smoother with multiple data assimilation for assisted history matching | |
| CN109034447B (en) | Method and device for predicting water-drive reservoir water-cut rate-of-rise change | |
| CN112329232B (en) | Fracture-cavity type oil reservoir production dynamic characterization method, device, equipment and storage medium | |
| CN111814347B (en) | Method and system for predicting gas channeling channel in oil reservoir | |
| Hong et al. | Robust production optimization with capacitance-resistance model as precursor | |
| CN111101929A (en) | Method, device and system for calculating average formation pressure of oil and gas reservoir | |
| CN109829217A (en) | Pressure break Fractured Reservoir productivity simulation method and device | |
| CN110821453B (en) | Gas injection oil reservoir development scheme design method based on three-dimensional geological numerical model | |
| US11905817B2 (en) | Method and system for managing carbon dioxide supplies using machine learning | |
| CA3003701A1 (en) | Automated upscaling of relative permeability using fractional flow in systems comprising disparate rock types | |
| CN111677486A (en) | Compact oil carbon dioxide huff and puff simulation method and device and storage medium | |
| Zheng et al. | A Deep Learning-Based Surrogate Model for Rapid Assessment of Geomechanical Risks in Geologic CO2 Storage | |
| CN113392372B (en) | Gas well full life cycle non-blocking flow determination method and system | |
| Zhou et al. | Multiobjective Operation Optimization of Depleted Gas Storage Based on a Reference Vector–Guided Evolutionary Algorithm | |
| CN113010501B (en) | Recovery ratio prediction model acquisition method, recovery ratio prediction method and product | |
| CN114818514A (en) | Method and processor for determining the reservoir capacity of a gas reservoir | |
| CN112114357B (en) | Fracture-cavity type reservoir connectivity prediction method and device | |
| CN108920781B (en) | Method and device for determining water content and water content increase rate of oil reservoir | |
| CN108593877B (en) | Physical simulation experiment method for step-by-step depressurization exploitation mode of gas reservoir | |
| CN109063228B (en) | Method and device for determining water-drive reservoir water-containing rate-of-rise change | |
| CN109033508B (en) | Method and device for determining water content and water content increase rate of oil reservoir | |
| Capolei et al. | Offset risk minimization for open-loop optimal control of oil reservoirs | |
| WO2020229785A1 (en) | A method of optimizing production from a hydrocarbon reservoir | |
| CN105528732B (en) | Method for predicting productivity of gas testing well | |
| CN112052631B (en) | Method and device for determining initial unobstructed flow parameters of target well section of new drilled well |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |