Method, system and equipment for calculating optimal dosage of profile control agent and readable storage medium

Abstract

The invention relates to the technical field of oil exploitation, and discloses a method, a system and equipment for calculating the optimal dosage of a profile control agent and a readable storage medium, the method comprises the following steps: calculating the water absorption index of each permeable layer at the t-th time step according to the profile control agent plugging radius of each permeable layer at the t-th time step, the water absorption index Lorentz coefficient corresponding to the tth time step is calculated, a time-varying curve of the water absorption index Lorentz coefficient is drawn, and the water absorption index Lorentz coefficient is used for reflecting the non-uniform degree of the layered water absorption index; based on the curve that the water absorption index Lorentz coefficient changes along with time, the profile control agent dosage corresponding to the time step where the minimum water absorption index Lorentz coefficient is located is selected as the optimal profile control agent dosage. The method aims at minimizing the water absorption index Lorentz coefficient, the optimal dosage of the chemical profile control agent can be accurately predicted, and technical support is provided for profile control process decision making.

Description

Technical Field

[0001] This invention relates to the field of oil extraction technology, and in particular to a method, system, equipment, and readable storage medium for calculating the optimal dosage of profile control agent. Background Technology

[0002] In the later stages of waterflooding development in oilfields, the pore structure, physical parameters, and seepage characteristics of the reservoir will change due to the long-term scouring of injected water. The vertical and planar heterogeneity of the oil layer will intensify, and the waterflooding development effect will continuously deteriorate.

[0003] Water injection well profile control technology improves reservoir water absorption profile by blocking high-permeability zones, thereby increasing the sweep efficiency and water drive efficiency of injected water. The effectiveness of profile control and production enhancement measures largely depends on the scientific and reasonable optimization design of injection parameters. Existing technologies have established selective injection volume and selective injection pressure models for profile control agents using reservoir engineering methods. However, these methods do not consider the distribution radius of the profile control agent in the reservoir at different time steps, thus limiting their application. Summary of the Invention

[0004] This invention takes into account the different plugging radii of the profile control agent in each permeable layer at different time steps, and provides a method that can accurately predict the optimal dosage of chemical profile control agent with the goal of minimizing the water absorption index and Lorentz coefficient.

[0005] To achieve the above objectives, the first aspect of the present invention provides a method for calculating the optimal dosage of a profile control agent, comprising:

[0006] Calculate the water absorption index of each permeable layer at time step t based on the profile control agent sealing radius of each permeable layer at time step t.

[0007] Based on the water absorption index of each permeable layer at time step t, calculate the Lorentz coefficient of the water absorption index corresponding to time step t, and plot the curve of the Lorentz coefficient of the water absorption index changing with time. The Lorentz coefficient of the water absorption index is used to reflect the degree of non-uniformity of the water absorption index of the layers.

[0008] Based on the curve of the water absorption index and Lorentz coefficient changing with time, the amount of profile control agent corresponding to the time step with the minimum water absorption index and Lorentz coefficient is selected as the optimal amount of profile control agent.

[0009] Furthermore, based on the profile control agent sealing radius of each permeable layer at time step t, the water absorption index of each permeable layer at time step t is calculated, including:

[0010] Step S10: Let t = 0, no profile control agent flows into the permeable layer, and the profile control agent sealing radius r f (i,0)=0, calculate the water absorption index J(i,0) of the i-th permeable layer;

[0011]

[0012] Step S11: When t>0, determine the sealing radius r of the profile control agent in the i-th permeable layer at time step t. f (i,t), calculate the water absorption index J(i,t) of the i-th permeable layer at time step t;

[0013]

[0014] Where k(i) is the average permeability of the i-th permeable layer, in μm 2 ;

[0015] h(i) is the effective thickness of the i-th permeation layer, in meters;

[0016] 'a' is the unit conversion factor, and its value is a constant.

[0017] μ w The viscosity of formation water is expressed in mPa·s.

[0018] r w The radius of the injection well shaft, in meters (m).

[0019] r e The control radius of the injection well, in meters;

[0020] RRF is the residual drag coefficient of the profile control agent, and its value is a constant;

[0021] Step S12: When t = t + Δt, determine the sealing radius r of the profile control agent in the i-th permeable layer at time step t. f Calculate the profile control agent's cleavage component in the i-th permeable layer within time intervals (i,t) and Δt, and calculate the profile control agent's plugging radius r in each permeable layer at time step t+1. f (i,t+1), where Δt is the interval between two time steps;

[0022] Step S13: Repeat steps S11-S12 above to obtain the water absorption index of each permeable layer corresponding to the time step.

[0023] Furthermore, when t = t + Δt, calculate the profile control agent sealing radius r of each permeable layer at time step t+1. f (i,t+1):

[0024]

[0025] ΔQ(i,t) represents the splitting component of the profile control agent in the i-th permeable layer during the time interval Δt, in meters. 3 ;

[0026]

[0027] Where q is the injection rate of the profile control agent;

[0028] Δt is the interval between two time steps;

[0029] φ(i) represents the porosity of the i-th permeable layer;

[0030] s wc To bind water saturation;

[0031] n represents the total number of permeable layers.

[0032] Furthermore, based on the water absorption index of each permeable layer at time step t, the Lorentz coefficient of the water absorption index corresponding to time step t is calculated, including:

[0033] Arrange the water absorption index of each permeable layer in descending order at time step t, and calculate the cumulative proportion of the water absorption index contribution value and the cumulative proportion of the ordinal contribution value of each permeable layer at time step t.

[0034] Among them, the cumulative proportion Y of the water absorption index contribution value of the i-th permeable layer i Represented as:

[0035]

[0036] The cumulative proportion of the ordinal contribution value of the i-th permeation layer X i Represented as:

[0037]

[0038] Where m(j) is the ordinal number of the permeable layer, and n is the total number of permeable layers;

[0039] Plot the Lorentz curve of water absorption index at time step t with the cumulative proportion of the ordinal contribution value of the permeable layer as the horizontal axis and the cumulative proportion of the water absorption index contribution value as the vertical axis.

[0040] The ratio of the area enclosed by the Lorentz curve of the water absorption index at time step t to the area enclosed by the Lorentz curve of the water absorption index in extremely heterogeneous reservoirs is taken as the Lorentz coefficient of the water absorption index at time step t.

[0041] Furthermore, the Lorentz coefficient of the water absorption index at time step t ranges from [0,1].

[0042] Furthermore, based on the curve of the water absorption index (Lorentz coefficient) changing over time, the amount of profile control agent corresponding to the time step with the minimum water absorption index (Lorentz coefficient) is selected as the optimal amount of profile control agent, including:

[0043] On the curve of the water absorption index Lorentz coefficient changing with time, the time step t corresponding to the minimum Lorentz coefficient is selected. m ;

[0044] Based on time step t m Calculate the optimal dosage Q of the profile control agent. m Unit m 3 ;

[0045] Q m =qt m Formula (7) for Δt

[0046] Where q is the injection rate of the profile control agent;

[0047] Δt is the interval between two time steps.

[0048] Furthermore, the method also includes:

[0049] Calculate the profile control agent dosage corresponding to the time step;

[0050] Using the Lorentz coefficient of water absorption index as the ordinate and the amount of profile control agent as the abscissa, a curve showing the relationship between the Lorentz coefficient of water absorption index and the amount of profile control agent at different time steps was established.

[0051] The optimal amount of profile control agent is determined by the Lorentz coefficient corresponding to the minimum water absorption index.

[0052] A second aspect of the present invention also provides a system for calculating the optimal dosage of a profile control agent, comprising:

[0053] The calculation unit is used to calculate the water absorption index of each permeable layer at time step t based on the profile control agent sealing radius of each permeable layer at time step t.

[0054] The image plotting unit is used to calculate the Lorentz coefficient of the water absorption index corresponding to the time step t based on the water absorption index of each permeable layer at the time step t, and plot the curve of the Lorentz coefficient of the water absorption index changing with time. The Lorentz coefficient of the water absorption index is used to reflect the degree of non-uniformity of the water absorption index of the layers.

[0055] The calculation unit is used to select the amount of profile control agent corresponding to the time step with the minimum Lorentz coefficient of water absorption index as the optimal amount of profile control agent, based on the curve of the Lorentz coefficient of water absorption index changing over time.

[0056] Based on the same inventive concept, the present invention also provides an electronic device, including: a memory and a processor; the processor is used to read and execute a computer program stored in the memory to implement the aforementioned method for calculating the optimal dosage of the profile control agent.

[0057] Based on the same inventive concept, the present invention also provides a computer storage medium storing computer-executable instructions, which, when executed, implement the aforementioned method for calculating the optimal dosage of profile control agent.

[0058] The technical effects and advantages of this invention are as follows: Based on the general mechanism of profile control in multi-layer reservoirs, and considering the different sealing radii of profile control agents in each permeable layer at different time steps, this invention calculates the water absorption index of each permeable layer at different time steps and the Lorentz coefficient of the water absorption index corresponding to the time step. It uses a time step iterative calculation method to simulate the injection process of the profile control agent. This invention aims to minimize the water absorption index and the Lorentz coefficient, and can accurately predict the optimal dosage of chemical profile control agents, providing technical support for profile control process decisions.

[0059] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures pointed out in the description, claims and drawings. Attached Figure Description

[0060] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0061] Figure 1 This is a flowchart illustrating the method for calculating the optimal dosage of the profile control agent in an embodiment of the present invention.

[0062] Figure 2 This is the Lorentz curve of the water absorption index at time step t in this embodiment of the invention;

[0063] Figure 3 This is a curve showing the relationship between the water absorption index (Lorentz coefficient) and the amount of profile control agent in an embodiment of the present invention.

[0064] Figure 4 This is a schematic diagram of the structure of the calculation system for the optimal dosage of profile control agent in an embodiment of the present invention;

[0065] Calculation unit 100; image drawing unit 200; calculation unit 300. Detailed Implementation

[0066] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0067] In multi-layered syngas reservoirs, due to the vertical heterogeneity of the reservoir, the injection of profile control agents exhibits a natural selection effect, meaning that the amount of profile control agent absorbed varies across different permeable layers. In layers with higher water absorption indices, profile control agents are more widely distributed, while in layers with lower water absorption indices, they are less distributed or not distributed at all. The water absorption index refers to the daily water injection rate of an injection well under a unit injection pressure differential, reflecting the water injection capacity of the well and the water absorption capacity of the reservoir. It can be used to analyze the operating status of the injection well and changes in the water absorption capacity of the reservoir. The distribution of profile control agents is closely related to the layered water absorption indices.

[0068] like Figure 1 As shown, this embodiment of the invention provides a method for calculating the optimal dosage of profile control agent, including the following steps:

[0069] Step S1: Calculate the water absorption index of each permeable layer at time step t based on the profile control agent sealing radius of each permeable layer at time step t.

[0070] Step S2: Based on the water absorption index of each permeable layer at time step t, calculate the Lorentz coefficient of the water absorption index corresponding to time step t, and plot the curve of the Lorentz coefficient of the water absorption index changing with time. The Lorentz coefficient of the water absorption index is used to reflect the degree of non-uniformity of the water absorption index of the layers.

[0071] Step S3: Based on the curve of the water absorption index Lorentz coefficient changing with time, select the profile control agent dosage corresponding to the time step with the minimum water absorption index Lorentz coefficient as the optimal dosage of profile control agent.

[0072] In this embodiment of the invention, based on the general mechanism of profile control in multi-layer reservoirs, and considering that the sealing radius of the profile control agent in each permeable layer is different at different time steps, the water absorption index of each permeable layer at different time steps is calculated, and the Lorentz coefficient of the water absorption index corresponding to the time step is calculated. The injection process of the profile control agent is simulated by using a time step iterative calculation method. The present invention aims to minimize the water absorption index and the Lorentz coefficient, and can accurately predict the optimal dosage of chemical profile control agent, providing technical support for profile control process decision-making.

[0073] According to a specific implementation method, the present invention uses a time-step iterative calculation method to simulate the entire injection process of the profile control agent. Based on the profile control agent plugging radius and the performance of the profile control agent in each permeable layer at different time steps, the stratification water absorption index of each permeable layer is calculated.

[0074] Step S1, which involves calculating the water absorption index of each permeable layer at time step t based on the sealing radius of the profile control agent at time step t, includes the following steps:

[0075] Step S10: Let t = 0, no profile control agent flows into the permeable layer, and the profile control agent sealing radius r f(i,0)=0, calculate the water absorption index J(i,0) of the i-th permeable layer;

[0076]

[0077] Step S11: When t>0, determine the sealing radius r of the profile control agent in the i-th permeable layer at time step t. f (i,t), calculate the water absorption index J(i,t) of the i-th permeable layer at time step t;

[0078]

[0079] Where k(i) is the average permeability of the i-th permeable layer, in μm 2 ;

[0080] h(i) is the effective thickness of the i-th permeation layer, in meters;

[0081] 'a' is the unit conversion factor, and its value is a constant.

[0082] μ w The viscosity of formation water is expressed in mPa·s.

[0083] r w The radius of the injection well shaft, in meters (m).

[0084] r e The control radius of the injection well, in meters;

[0085] RRF is the residual drag coefficient of the profile control agent, and its value is a constant;

[0086] Step S12: When t = t + Δt, determine the sealing radius r of the profile control agent in the i-th permeable layer at time step t. f Calculate the profile control agent's cleavage component in the i-th permeable layer within time intervals (i,t) and Δt, and calculate the profile control agent's plugging radius r in each permeable layer at time step t+1. f (i,t+1), where Δt is the interval between two time steps;

[0087] Step S13: Repeat steps S11-S12 above to obtain the water absorption index of each permeable layer corresponding to the time step.

[0088] Furthermore, in step S12 above, when t = t + Δt, the profile control agent sealing radius r of each permeable layer at time step t+1 is calculated. f (i,t+1):

[0089]

[0090] ΔQ(i,t) represents the splitting component of the profile control agent in the i-th permeable layer during the time interval Δt, in meters. 3;

[0091]

[0092] Where q is the injection rate of the profile control agent;

[0093] Δt is the interval between two time steps;

[0094] φ(i) represents the porosity of the i-th permeable layer;

[0095] s wc To bind water saturation;

[0096] n represents the total number of permeable layers.

[0097] It should be noted that the aforementioned time step t, time step t+1, or other time steps arranged in sequence can all be considered as time step t. The present invention does not impose any restrictions on the value of t. According to the above method, the permeable layer water absorption index corresponding to each time step is obtained.

[0098] In this embodiment of the invention, step S2 above, which involves calculating the Lorentz coefficient of the water absorption index corresponding to time step t based on the water absorption index of each permeable layer at time step t, includes the following steps:

[0099] Step S20: Arrange the water absorption index of each permeable layer in descending order at time step t, and calculate the cumulative proportion of the water absorption index contribution value and the cumulative proportion of the ordinal contribution value of each permeable layer at time step t.

[0100] Among them, the cumulative proportion Y of the water absorption index contribution value of the i-th permeable layer i Represented as:

[0101]

[0102] The cumulative proportion of the ordinal contribution value of the i-th permeation layer X i Represented as:

[0103]

[0104] m(j) is the ordinal number of the permeable layer, for example, the j-th permeable layer arranged in sequence, and n is the total number of permeable layers; X i and Y i The value range of is [0,1];

[0105] Step S21: Plot the Lorentz curve of water absorption index at time step t, with the cumulative proportion of the ordinal contribution value of the permeable layer as the horizontal axis and the cumulative proportion of the water absorption index contribution value as the vertical axis.

[0106] In practical applications, if the reservoir's fluid production (water absorption) profile is uniform, it means that the reservoir's water absorption index is uniform, such as... Figure 2 As shown, the Lorenz curve for the water absorption index of a homogeneous reservoir is a straight line with a slope of 45°. Figure 2 The line segment AC in the diagram is called the "completely homogeneous line"; if all the production (water absorption) of the reservoir is concentrated at one point, it means that the water absorption index of the reservoir is extremely heterogeneous, and the Lorenz curve of the water absorption index of an extremely heterogeneous reservoir is a broken line ( Figure 2 The broken line ABC in the diagram is called a "completely heterogeneous line". Generally, the Lorenz curve of the water absorption index in conventional heterogeneous reservoirs falls between these two cases, and is a "convex" curve. Figure 2 The curve ADC is a curve in which the distance from any point other than the two endpoints A and C to the two coordinate axes is not equal.

[0107] In this embodiment, the Lorentz curve of the water absorption index at time step t represents the Lorentz curve of the water absorption index of a conventional heterogeneous reservoir, which is used to reflect the degree of variation of the reservoir water absorption index. The closer it is to the "completely homogeneous line", the more uniform the production (water absorption) profile is.

[0108] Step S22: The ratio of the area enclosed by the Lorentz curve of the water absorption index at time step t to the area enclosed by the Lorentz curve of the water absorption index of the extremely heterogeneous reservoir is taken as the Lorentz coefficient of the water absorption index at time step t.

[0109] Specifically, the water absorption index and Lorentz coefficient D at time step t. J The calculation method for (t) is as follows:

[0110]

[0111] Among them, S ADC S represents the area of ​​the curvilinear ADC, and S is the dimensionless degree. ABC Let be the area of ​​triangle ABC, and be a dimensionless number.

[0112] In this embodiment of the invention, following the definition of the Gini coefficient in economics, the following is used: Figure 2 The area enclosed by the curve ADC in the middle (i.e. Figure 2 The ratio of the shaded area in the diagram to the area enclosed by the broken line ABC is called the Lorentz coefficient, which reflects the degree of non-uniformity of the water absorption index in stratification. The value of the Lorentz coefficient is used to quantitatively describe its non-uniformity.

[0113] According to a specific embodiment of the present invention, the Lorentz coefficient of the water absorption index at time step t ranges from [0,1]. Specifically, the Lorentz coefficient of the water absorption index is 0 for a completely homogeneous reservoir, 1 for an extremely heterogeneous reservoir, and between 0 and 1 for a conventional reservoir. A smaller Lorentz coefficient indicates a more uniform water production profile in the reservoir, and vice versa.

[0114] Furthermore, based on the Lorentz coefficient of water absorption index at different time steps, a curve showing the change of the Lorentz coefficient of water absorption index over time is plotted. On this curve, the amount of profile control agent corresponding to the time step with the lowest Lorentz coefficient of water absorption index is selected as the optimal amount of profile control agent.

[0115] According to a specific implementation method, on the curve of the water absorption index Lorentz coefficient changing with time, the time step t corresponding to the minimum Lorentz coefficient is selected. m According to time step t m Calculate the optimal dosage Q of the profile control agent. m Unit m 3 ;

[0116] Q m =qt m Formula (7) for Δt

[0117] Where q is the injection rate of the profile control agent; Δt is the interval between two time steps.

[0118] In this embodiment of the invention, the Lorentz curve is used to analyze the non-uniformity of the reservoir water absorption index, which can intuitively and clearly show the non-uniformity of the production and absorption profile. The Lorentz coefficient quantitatively describes the non-uniformity of the water absorption index in each stratification. In this invention, by calculating the Lorentz coefficient of the water absorption index at different time steps, the amount of profile control agent corresponding to the time node with the lowest Lorentz coefficient is selected as the optimal amount, thereby accurately predicting the optimal amount of chemical profile control agent and providing technical support for profile control process decisions.

[0119] Based on the above method, according to a preferred embodiment, the method for calculating the optimal dosage of the profile control agent further includes:

[0120] Calculate the profile control agent dosage corresponding to the time step; establish the relationship curve between the water absorption index Lorentz coefficient and the profile control agent dosage at different time steps with the water absorption index Lorentz coefficient as the ordinate and the profile control agent dosage as the abscissa; take the profile control agent dosage corresponding to the minimum water absorption index Lorentz coefficient as the optimal dosage of profile control agent.

[0121] Example

[0122] A certain water injection well has a depth of 2500m and an oil drainage radius of 250m. The current average daily water injection volume is 100m³. 3 / d. The formation where the injection well is located is highly heterogeneous, and the injected water rushes along a single layer, resulting in severe ineffective water circulation. After comprehensive analysis, it was decided to conduct a general profile adjustment.

[0123] First, reservoir geological parameters and development dynamic data are collected, for example, the density of the injected chemical profiler is 1.02 g / cm³. 3The drag coefficient is 29, the residual drag coefficient is 16, and the small-layer physical property data of the injection well are as follows:

[0124]

[0125] Iterative calculations were performed at different time steps to obtain the total amount of profile control agent used and the water absorption index (Lorentz coefficient) at different time steps, such as... Figure 3 As shown, a curve relating the water absorption index (Lorentz coefficient) to the amount of profile control agent was plotted.

[0126] from Figure 3 It can be seen that when the injection volume of the chemical profile control agent is 1205m... 3 At that time, the water absorption index and Lorentz coefficient of the reservoir reached their minimum after profile control measures. Figure 3 The formation's water absorption profile was improved to the maximum extent (shown as 0.368). Therefore, the optimal dosage of profile control agent for this injection well is 1205 m³. 3 .

[0127] Based on the same inventive concept, embodiments of the present invention also provide a system for calculating the optimal dosage of profile control agent, such as... Figure 4 As shown, the system includes:

[0128] The calculation unit 100 is used to calculate the water absorption index of each permeable layer at time step t based on the profile control agent sealing radius of each permeable layer at time step t.

[0129] The image rendering unit 200 is used to calculate the Lorentz coefficient of the water absorption index corresponding to the time step t based on the water absorption index of each permeable layer at the time step t, and to plot the curve of the Lorentz coefficient of the water absorption index changing with time.

[0130] Among them, the Lorentz coefficient of water absorption index is used to reflect the degree of non-uniformity of the water absorption index of stratification;

[0131] The calculation unit 300 is used to select the amount of profile control agent corresponding to the time step with the minimum Lorentz coefficient of water absorption index as the optimal amount of profile control agent based on the curve of the Lorentz coefficient of water absorption index changing with time.

[0132] Regarding the system in the above embodiments, the specific manner in which each unit module performs operations has been described in detail in the embodiments related to the method, and will not be elaborated here.

[0133] Based on the same inventive concept, embodiments of the present invention also provide an electronic device, including: a memory and a processor, wherein the processor is used to read and execute a computer program stored in the memory to implement the aforementioned method for calculating the optimal dosage of profile control agent.

[0134] Based on the same inventive concept, embodiments of the present invention also provide a computer storage medium storing computer-executable instructions, which, when executed, implement the aforementioned method for calculating the optimal dosage of profile control agent.

[0135] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of modules is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple modules or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between apparatuses or modules may be electrical, mechanical, or other forms.

[0136] The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical modules; that is, they may be located in one place or distributed across multiple network modules. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Furthermore, the functional modules in the various embodiments of this invention can be integrated into one processing module, or each module can exist physically separately, or two or more modules can be integrated into one module. The integrated modules described above can be implemented in hardware or as software functional modules.

[0137] If the integrated module is implemented as a software functional module and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0138] It should be noted that, for the sake of simplicity, the foregoing method embodiments are all described as a series of actions. However, those skilled in the art should understand that the present invention is not limited to the described order of actions, because according to the present invention, some steps can be performed in other orders or simultaneously. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are preferred embodiments, and the actions and modules involved are not necessarily essential to the present invention.

[0139] In the above embodiments, the descriptions of each embodiment have different focuses. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0140] Finally, it should be noted that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for calculating the optimum amount of profile control agent, characterized in that, The method comprises: According to the profile control agent plugging radius of each permeable layer at the t time step, the water absorption index of each permeable layer at the t time step is calculated; According to the water absorption index of each permeable layer at the t time step, the water absorption index Lorenz coefficient corresponding to the t time step is calculated, and a curve of the water absorption index Lorenz coefficient changing with time is drawn, wherein the water absorption index Lorenz coefficient is used to reflect the uneven degree of the layered water absorption index; Based on the curve of the water absorption index Lorenz coefficient changing with time, the profile control agent dosage corresponding to the time step where the minimum water absorption index Lorenz coefficient is located is selected as the optimal profile control agent dosage.

2. The method of claim 1, wherein, The water absorption index of each permeable layer at the t time step is calculated according to the profile control agent plugging radius of each permeable layer at the t time step, comprising: Step S10, let t = 0, the permeable layer has no profile control agent flowing in, the profile control agent plugging radius r f (i,0) = 0, calculate the water absorption index J(i,0) of the i permeable layer; Step S11, when t>0, calculating the profile control agent plugging radius r of the i th permeable layer at the t th time step according to the formula f (i,t), calculating the water absorption index J(i,t) of the i th permeable layer at the t th time step; wherein k(i) is the average permeability of the i-th permeation layer, in pm 2 ; h(i) is the effective thickness of the i permeable layer, with a unit of m; a is a unit conversion coefficient, which is a constant; μ w Viscosity of the formation water, in mPa-s; r w R is the wellbore radius for the injection well in meters; r e R is the control radius of the injection well, in meters; RRF is the residual resistance factor of the profile control agent, which is a constant; Step S12, when t=t+Δt, the profile control agent plugging radius r of the i-th permeable layer at the t time step is calculated according to the profile control agent plugging radius r f (i,t) and the split quantity of the profile control agent in the i-th permeable layer in the Δt time period, the profile control agent plugging radius r of each permeable layer at the t+1 time step is calculated f (i,t+1), wherein Δt is the interval of two time steps. Step S13, repeat the above steps S11-S12 to obtain the water absorption index of each permeable layer corresponding to the time step.

3. The method of claim 2, wherein, When t = t + Δt, the profile control agent plugging radius r of each permeable layer at the t + 1 time step is calculated f (i,t+1): ΔQ(i, t) is the split quantity of the profile control agent in the ith permeable layer in the time period of Δt, unit m 3 ; Wherein q is the injection displacement of the profile control agent; Δt is the interval of two time steps; φ(i) is the porosity of the i permeable layer; s wc Sw is the irreducible water saturation; n is the total number of permeable layers.

4. The method of claim 3, wherein, The water absorption index Lorenz coefficient corresponding to the t time step is calculated according to the water absorption index of each permeable layer at the t time step, comprising: Arranging the water absorption index of each permeable layer at the t time step in descending order, and calculating the water absorption index contribution value cumulative proportion and the ordinal contribution value cumulative proportion of each permeable layer at the t time step respectively; wherein the cumulative proportion Y of the water absorption index contribution value of the i-th permeation layer i is represented as: Cumulative proportion of ordinal contribution value of the i-th permeation layer X i is represented as: Wherein m(j) is the ordinal of the permeable layer, and n is the total number of permeable layers; Taking the ordinal contribution value cumulative proportion of the permeable layer as the abscissa and the water absorption index contribution value cumulative proportion as the ordinate, the water absorption index Lorenz curve at the t time step is drawn; The ratio of the area surrounded by the water absorption index Lorenz curve at the t time step to the area surrounded by the water absorption index Lorenz curve of the extremely heterogeneous reservoir is taken as the water absorption index Lorenz coefficient at the t time step.

5. The method of claim 4, wherein, The value range of the water absorption index Lorenz coefficient at the t time step is [0, 1].

6. The method of claim 1, wherein, Based on the curve of the water absorption index Lorenz coefficient changing with time, the profile control agent dosage corresponding to the time step where the minimum water absorption index Lorenz coefficient is located is selected as the optimal profile control agent dosage, comprising: on the curve of the water absorption index Lorentz coefficient changing with time, the time step t corresponding to the minimum Lorentz coefficient is selected m ; According to time step t m Calculating the best amount of profile control agent Q m , unit m 3 ; Q m = qt m Δt Formula (7) Wherein q is the injection displacement of the profile control agent; Δt is the interval of two time steps.

7. The method of claim 1, wherein, The method further comprises: Calculating the profile control agent dosage corresponding to the time step; Taking the water absorption index Lorenz coefficient as the ordinate and the profile control agent dosage as the abscissa, a relationship curve of the water absorption index Lorenz coefficient and the profile control agent dosage at different time steps is established; The profile control agent dosage corresponding to the minimum water absorption index Lorenz coefficient is taken as the optimal profile control agent dosage.

8. A system for calculating the optimum amount of profile control agent, characterized in that, The system comprises: A calculation unit is configured to calculate the water absorption index of each permeable layer at the t time step according to the profile control agent plugging radius of each permeable layer at the t time step. The image drawing unit is used for calculating the water absorption index Lorenz coefficient corresponding to the t time step according to the water absorption index of each permeable layer at the t time step, and drawing the curve of the water absorption index Lorenz coefficient changing with time, wherein the water absorption index Lorenz coefficient is used for reflecting the uneven degree of layered water absorption index. The estimation unit is used for selecting the profile control agent dosage corresponding to the time step of the minimum water absorption index Lorenz coefficient as the optimal dosage of the profile control agent based on the curve of the water absorption index Lorenz coefficient changing with time.

9. An electronic device, comprising: Comprise: Memory, processor; The processor is used for reading and executing the computer program stored in the memory, so as to realize the method in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer executable instructions, and the computer executable instructions are executed to realize the method in any one of claims 1-7.

Images