Polymer flooding numerical simulation method based on high-order WENO method

Abstract

The invention discloses a polymer flooding numerical simulation method based on a high-order WENO method, and the method comprises the steps: fixing the polymer concentration at the current moment, and carrying out the integration through a finite volume method. A polynomial function meeting the r-order precision is found in the WENO method, and the numerical flux on the boundary is reconstructed. The convex combination of a plurality of templates is adopted to solve the problem of discontinuous units, and the smoothness factor is utilized to measure the smoothness of each template. After the water saturation at the next moment is obtained through the WENO method, the water saturation is fixed, and for a seepage differential equation of water components, the polymer solution concentration at the next moment is obtained through the WENO method. And alternately solving to obtain a distribution rule of parameters such as polymer concentration and water saturation at all moments. And finally, compared with the traditional first-order windward difference method, the result proves that the method provided by the invention has higher calculation precision on the saturation and the concentration front edge. The method is used for high-precision calculation of the polymer flooding concentration front and provides guidance for polymer flooding numerical simulation.

Description

Technical Field

[0001] The invention belongs to the technical field of numerical simulation of reservoir polymer flooding, and in particular relates to a polymer flooding numerical simulation method based on a high-order WENO method. Background Art

[0002] Polymer flooding technology is one of the important methods to improve oil recovery. In this process, by injecting high-viscosity polymer solutions, the flow characteristics of the fluid can be improved and the displacement efficiency of oil and water can be enhanced. However, its effect mainly depends on the accurate simulation and prediction of the flow behavior of the polymer solution in the formation.

[0003] Traditional low-order numerical methods such as first-order finite difference or finite volume method may introduce numerical dispersion, resulting in unclear calculated concentration fronts. Numerical oscillations can cause instability in the calculation results, especially in high-gradient regions such as concentration fronts.

[0004] Therefore, establishing a polymer numerical simulation method that can capture detailed features with high accuracy and good stability is of great significance for accurately describing the polymer driving process and improving the efficiency of oilfield development. Summary of the invention

[0005] The problem to be solved by the present invention is to provide a polymer flooding numerical simulation method based on the high-order WENO method, which aims to use the WENO method to stably and accurately calculate the polymer flooding front concentration and saturation, and provide theoretical guidance for accurately describing the polymer flooding process.

[0006] In order to solve the above technical problems, the technical solution adopted by the present invention is: a polymer flooding numerical simulation method based on a high-order WENO method, comprising the following steps:

[0007] S1: Viscosity correction for polymer solution;

[0008] S2: describes the ratio of oil-water relative permeability, which is substituted into the fractional flow equation to calculate the water cut;

[0009] S3: Solve the seepage differential equation of polymer solution using WENO method;

[0010] S4: Compare with the upwind difference format to verify its computational performance.

[0011] Further, in S1, the polymer solution changes the oil-water mobility ratio by viscosity increase, and the viscosity correction is performed using the following cubic polynomial:

[0012] μ w =μ wi (1+rC p +sC p2 +tC p 3 )

[0013] Among them, μ wi and μ w are the viscosity of the polymer solution before and after adding the polymer; C p is the polymer concentration; r, s, t are the equation parameters measured experimentally.

[0014] Further, the S2 comprises the following steps:

[0015] S21: Assuming that the ratio of oil-water relative permeability varies linearly with water saturation, the simplified expression based on Archie's formula is:

[0016]

[0017] Among them, a and b are constant parameters to be determined, which are obtained by fitting experimental data; S w is the water saturation;

[0018] S22: Considering the influence of residual oil saturation and irreducible water saturation, the modified model expression is:

[0019]

[0020] Among them, S wc is the bound water saturation, S or is the residual oil saturation;

[0021] S23: Ignoring the effects of capillary force and gravity, the polymer solution is regarded as the water phase, and the modified viscosity equation and relative permeability are brought in. The flow equation expression is:

[0022]

[0023] Further, the S3 comprises the following steps:

[0024] S31: The polymer flooding process is regarded as a one-dimensional oil-water two-phase flow. The differential equation for the seepage of the polymer solution is expressed as:

[0025]

[0026] S32: After fixing the polymer concentration, the water content is considered as a function of water saturation as follows,

[0027] f w (S w ,C p (n) );

[0028] S33: For the water phase seepage differential equation, the finite volume method is used in [x i-Δx / 2 ,x i+Δx / 2 ], we can get the following formula:

[0029]

[0030] S34: Let f(x, t) = f w (x, t), dividing both sides by Δx, we get,

[0031]

[0032] Among them, S wi It is expressed as the water saturation on the i unit; S wi is the integration domain [x i-Δx / 2 ,x i+Δx / 2 ], that is, f(x,t) represents the numerical flux of the polymer solution;

[0033] S35: Solve the water saturation equation using the WENO method;

[0034] S36: In the WENO method, the smoothing factor β k To measure the smoothness of each template;

[0035] S37: After using the WENO method to calculate the water saturation at time n+1, the seepage differential equation for the water component in the polymer solution is as follows:

[0036]

[0037] S38: The polymer concentration at time n+1 is calculated using the WENO method. Based on the polymer concentration and water saturation at time n+1, the calculation is repeated to obtain the water saturation and polymer concentration distribution at time n+2 and thereafter, and the parameter distribution at the outlet is obtained.

[0038] Further, the S35 includes the following steps:

[0039] S351: Use the combination of three templates S0, S1 and S2 to establish the WENO discrete format.

[0040]

[0041] S352: To ensure stability and consistency, the weight ω r The following conditions are met,

[0042]

[0043] S353: When all templates S r When all are in the smooth region, that is, any unit S in the template r (i) are not discontinuous units, then let ω r =y r ,get,

[0044]

[0045] S354: When template S r When located in the discontinuous region, ω r The calculation form of is as follows,

[0046]

[0047] The ε coefficient is introduced to prevent the denominator from being zero.

[0048] Furthermore, in S36, the smoothness calculation formula of each template is as follows:

[0049]

[0050] Further, the S4 comprises the following steps:

[0051] S41: In the traditional upwind difference method, the one-dimensional convection equation is as follows:

[0052]

[0053] Among them, S w To substitute variables, Q, A, and φ are constants;

[0054] S42: The concentration of the polymer solution at a fixed time n, the water content is a function of the water saturation f w (S w );

[0055] S43: Discretize the one-dimensional convection equation in space and time, using a grid of points x i and time step t n , the grid spacing is Δx, the time step is Δt, then the discretization form of the first-order upwind format is,

[0056]

[0057] S44: Use the formula of S43 to perform iterative calculation to obtain the water saturation at time n+1;

[0058] S45: using the water saturation at time n+1 to alternately solve and perform iterative calculations to obtain the polymer solution concentration and its water content distribution at time n+1;

[0059] S46: After calculating the polymer solution concentration at each fixed time, repeat the above steps to obtain the water saturation, polymer solution concentration and water content distribution at all times.

[0060] Furthermore, the present invention also provides a device for running any one of the data processing methods described above.

[0061] Furthermore, the present invention also provides a device, comprising a memory, a processor, and an algorithm stored in the memory and executable on the processor, wherein the processor implements the above-mentioned data processing method when executing the computer program.

[0062] Furthermore, the present invention also provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer algorithm, and the computer algorithm implements the above-mentioned data processing when executed by a processor.

[0063] The advantages and positive effects of the present invention are:

[0064] The present invention first fixes the polymer concentration at the current moment, and integrates the seepage differential equation of the polymer solution using the finite volume method. In the WENO method, a polynomial function that satisfies the r-order accuracy is found, and the numerical flux on the boundary is reconstructed. In order to eliminate the numerical oscillation caused by the discontinuous units contained in the fixed template, a convex combination of multiple templates is used to solve the problem of discontinuous units, and the smoothness of each template is measured using a smooth factor. After the water saturation at the next moment is obtained by the WENO method, it is fixed, and the polymer solution is regarded as consisting of two components, polymer and water. For the seepage differential equation of the water component, the WENO method is used to solve the polymer solution concentration at the next moment. In this way, the distribution law of parameters such as polymer concentration and water saturation at all times is obtained by alternating solutions. Finally, it is compared with the traditional first-order upwind difference method, and the results verify that the method proposed by the present invention has higher calculation accuracy for saturation and concentration front. The present invention is used to calculate the concentration front of polymer flooding with high precision, and provides guidance for polymer flooding numerical simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 It is a schematic diagram of the overall process of an embodiment of the present invention.

[0066] Figure 2 Graph 4 is a graph of a single template S and a polynomial function p(x) according to an embodiment of the present invention.

[0067] Figure 3 is a multi-template S of an embodiment of the present invention i The values ​​and their combination diagrams.

[0068] Figure 4 It is a water saturation distribution diagram calculated by different methods in Example 1 of the present invention.

[0069] Figure 5 1 is a polymer concentration distribution diagram calculated by different methods in Example 1 of the present invention.

[0070] Figure 6 It is a water saturation distribution diagram calculated by different methods in Example 2 of the present invention.

[0071] Figure 7 1 is a polymer concentration distribution diagram calculated by different methods in Example 2 of the present invention. DETAILED DESCRIPTION

[0072] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0073] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicating the orientation or positional relationship, are based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance.

[0074] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0075] The embodiments of the present invention are further described below in conjunction with the accompanying drawings:

[0076] like Figure 1 As shown, a polymer flooding numerical simulation method based on a high-order WENO method includes the following steps.

[0077] S1: Viscosity correction for polymer solution.

[0078] Assuming that the changes in oil phase viscosity and relative permeability are small during polymer solution flooding and are considered constant, the water content is mainly affected by the concentration of the polymer solution. In addition, the water content is also related to the water saturation. The polymer solution changes the oil-water mobility ratio by viscosity increase, and a cubic polynomial is used to correct the viscosity. The specific formula is as follows:

[0079] μ w =μ wi (1+rC p +sC p 2 +tC p 3 ) (1)

[0080] Among them, μ wi and μ w are the viscosity of the polymer solution before and after adding the polymer; C p is the polymer concentration; r, s, t are the equation parameters measured experimentally.

[0081] S2: Based on the mechanism, an expression for the viscosity of the polymer solution is obtained to describe the ratio of the relative permeabilities of oil and water, and it is substituted into the fractional flow equation to calculate the water cut.

[0082] Specifically, the oil-water relative permeability ratio is usually expressed as a function of water saturation, which is used to describe the flow behavior of multiphase fluids in porous media. The general form of the oil-water relative permeability ratio changing with water saturation can be obtained by fitting experimental data. Commonly used models include empirical models based on simplified assumptions and Corey models.

[0083] Assuming that the ratio of oil-water relative permeability varies linearly with water saturation, the simplified form based on Archie's formula can be expressed as:

[0084]

[0085] Among them, a and b are unknown constant parameters that need to be obtained by fitting experimental data; Sw is water saturation.

[0086] Considering the influence of residual oil saturation and irreducible water saturation, the modified exponential function model can more accurately describe the influence of water saturation on the ratio of oil-water relative permeability. In this modified model, it is assumed that the relative permeabilities of the oil phase and the water phase change exponentially with the water saturation, but only work within the effective saturation range (i.e., between irreducible water saturation and residual oil saturation). The modified model can be expressed as:

[0087]

[0088] Among them, S wc is the bound water saturation, Sor is the residual oil saturation.

[0089] Ignoring the effects of capillary force and gravity, the polymer solution is regarded as the water phase, which is equivalent to oil-water two-phase flow. Substituting the corrected viscosity equation and relative permeability, the flow equation can be expressed as:

[0090]

[0091] S3: The WENO method is used to solve the percolation differential equation of polymer solution.

[0092] Considering polymer flooding as oil-water two-phase flow in one-dimensional direction, the differential equation for the seepage of polymer solution (water phase) can be written as:

[0093]

[0094] At this time, if the polymer concentration at time n is fixed, the water content can actually be regarded as a function of water saturation: f w (S w ,C p (n) ).

[0095] For the water phase seepage differential equation, the finite volume method is used in [x i-Δx / 2 ,x i+Δx / 2 ], we can get:

[0096]

[0097] Let f(x,t)=f w (x, t), dividing both sides by Δx, we get:

[0098]

[0099] Among them, S wi It is expressed as the water saturation on unit i. It can be understood as S wi is the integration domain [x i-Δx / 2 ,x i+Δx / 2 ], f(x,t) is denoted as the numerical flux of the polymer solution (aqueous phase).

[0100] The exact solution of the water phase seepage differential equation is shown in equation (7). However, x i The numerical flux f(x i-1 / 2 ,t) and f(x i+1 / 2 ,t) is difficult to obtain, and a high-order reconstruction of the numerical flux f on the boundary is required.

[0101] In the WENO method, we need to find a polynomial function p(x) such that the polynomial function satisfies f(x i ,t)=p(x i ), and has r-order accuracy. It can be expressed as:

[0102]

[0103] As long as we know the specific form of the polynomial function p(x), we can obtain the numerical flux approximation on the boundary with sufficient accuracy. Right now:

[0104]

[0105] Therefore, combined with formula (17), it can be converted into:

[0106]

[0107] In order to obtain the polynomial function p(x), it is necessary to establish a template S(i) with r-1 grids near the i node. Assuming that there are m grids on the left side of the i node and n grids on the right side, S(i) is expressed as:

[0108] S(i)={I i-m ,…,I i ,…,I i+n}, m+n=r-1 (11)

[0109] Under the condition that the function value on the template S(i) is known (given by the initial condition at t = 0), the Lagrange interpolation method is used to represent the polynomial function p(x), then the numerical flux at the boundary is It is expressed as:

[0110]

[0111] Among them, the polynomial coefficient l i+k It can be expressed as:

[0112]

[0113] This embodiment uses the fifth-order WENO format to solve the hyperbolic partial differential equation. Therefore, for the template S(x), three numerical flux units need to be selected to solve the polynomial function p(x), such as Figure 2 In addition, the present invention uses the WENO method to solve the water saturation equation in formula (7). In order to eliminate the numerical oscillation caused by the discontinuous units contained in the fixed template, the WENO format uses multiple templates S iThe convex combination of S0, S1 and S2 is used to solve the problem of discontinuous units. In the present invention, three template combinations S0, S1 and S2 are used to establish the WENO discrete format, such as Figure 3 shown.

[0114] In the 5th order WENO format, a single template S1 = {I i-1 , I i , I i+1}, combining formula (21) and formula (22), the expression of the polynomial function p(x) can be solved as follows:

[0115]

[0116] in:

[0117]

[0118] After the polynomial p(x) is known, the reconstructed The form is expressed as:

[0119]

[0120] Similarly, the reconstructed numerical flux in templates S2 and S3 can be solved. The WENO method requires all templates S i Convex combinations are created in the following form:

[0121]

[0122] To ensure stability and consistency, the weight ω r The following conditions must be met:

[0123]

[0124] When all templates S r When all are in the smooth region, that is, any unit S in the template r (i) None of them are discontinuous units. In this case, let ω r =y r ,available:

[0125]

[0126] When the template S r In the discontinuous area, in order to reduce the numerical oscillation, the template weight ω r It should be small enough to reduce the impact of intermittent values ​​on the calculation.

[0127] Under discontinuous conditions, Jiang and Shu (1996) gave ω through a large number of experiments. r The calculation form is:

[0128]

[0129] In formula (22), the ε coefficient is introduced to prevent the denominator from being 0. In this embodiment, ε=10 -6 .

[0130] In the WENO method, the smoothing factor β k To measure the smoothness of each template. This embodiment directly adopts the calculation form given by Jiang and Shu (1996):

[0131]

[0132] After using WENO to solve the water saturation at time n+1, the water saturation at time n+1 is fixed, and the polymer solution is regarded as two components, polymer and water. For the water component in the polymer solution, its seepage differential equation can be written as:

[0133]

[0134] Then use the WENO method to solve the polymer concentration at time n+1. In this way, the polymer concentration and water saturation at time n+1 are obtained, and the water saturation and polymer concentration distribution at time n+2... are obtained repeatedly, and the parameter distribution at the outlet end can also be obtained.

[0135] S4: Compare with the upwind difference format to verify the computational performance of the new method.

[0136] In the traditional first-order upwind difference method, the one-dimensional convection equation is considered:

[0137]

[0138] Where S w Substitute variables for variables, Q, A, and φ are constants.

[0139] At this time, the concentration of the polymer solution at time n is fixed, and the water content can actually be regarded as a function of water saturation: f w (S w ).

[0140] Discretize the above equation in space and time, using the grid points x i and time step t n , the grid spacing is Δx, and the time step is Δt.

[0141] The discretized form of the first-order upwind format is:

[0142]

[0143] After sorting, we get:

[0144]

[0145] By iteratively calculating using the above formula, the water saturation at time n+1 can be obtained.

[0146] At this time, the water saturation at time n+1 is used to alternately solve the problem, and the polymer solution is regarded as two components, polymer and water. For the water component in the polymer solution, its seepage differential equation can be written as:

[0147]

[0148] Q w =Q·f w

[0149] The discretized form of the first-order upwind format is:

[0150]

[0151] After sorting, we get:

[0152]

[0153] By iteratively calculating using the above formula, the concentration and water content distribution of the polymer solution at time n+1 can be obtained.

[0154] At this time, the concentration of the polymer solution at time n+1 is fixed, and the above steps are repeated to calculate the water saturation, polymer solution concentration, and water content distribution at all times.

[0155] In order to verify the computational performance of the new method proposed in the present invention, parameters were selected and two groups of embodiments were designed to compare the high-order WENO method with the first-order upwind scheme.

[0156] The present invention will be specifically described below in conjunction with specific embodiments:

[0157] Example 1

[0158] S1: Viscosity correction for polymer solution.

[0159] S2: Describe the ratio of oil-water relative permeability, substitute it into the fractional flow equation to calculate the water content. Specifically, the water content function is selected as a simplified expression based on Archie's formula, where a = -0.21, b = -1.43, and other physical parameters are shown in Table 1:

[0160] Table 1 Physical properties

[0161]

[0162] S3: Use the WENO method to solve the seepage differential equation of the polymer solution. Specifically, the initial and boundary conditions are set as:

[0163]

[0164] S4: The water saturation distribution and polymer concentration distribution calculated by the high-order WENO method and the first-order upwind finite difference method are shown in Figure 4. Figure 4 , Figure 5 shown.

[0165] Example 2

[0166] The moisture content function is selected as the modified model, and the same physical property parameters as in Table 1 are used to calculate the water saturation and component concentration using the high-order WENO method and the first-order upwind difference method. Figure 6 , Figure 7 shown.

[0167] The advantages and positive effects of the present invention are:

[0168] The present invention first fixes the polymer concentration at the current moment, and integrates the seepage differential equation of the polymer solution using the finite volume method. In the WENO method, a polynomial function that satisfies the r-order accuracy is found, and the numerical flux on the boundary is reconstructed. In order to eliminate the numerical oscillation caused by the discontinuous units contained in the fixed template, a convex combination of multiple templates is used to solve the problem of discontinuous units, and the smoothness of each template is measured using a smooth factor. After the water saturation at the next moment is obtained by the WENO method, it is fixed, and the polymer solution is regarded as consisting of two components, polymer and water. For the seepage differential equation of the water component, the WENO method is used to solve the polymer solution concentration at the next moment. In this way, the distribution law of parameters such as polymer concentration and water saturation at all times is obtained by alternating solutions. Finally, it is compared with the traditional first-order upwind difference method, and the results verify that the method proposed by the present invention has higher calculation accuracy for saturation and concentration front. The present invention is used to calculate the concentration front of polymer flooding with high precision, and provides guidance for polymer flooding numerical simulation.

[0169] The above is a detailed description of an embodiment of the present invention, but the content is only a preferred embodiment of the present invention and cannot be considered to limit the scope of implementation of the present invention. All equivalent changes and improvements made within the scope of the present invention should still fall within the scope of the patent coverage of the present invention.

Claims

1. A polymer flooding numerical simulation method based on a high-order WENO method, characterized in that: The following steps are included: S1: Viscosity correction for polymer solution; S2: describes the ratio of oil-water relative permeability, which is substituted into the fractional flow equation to calculate the water cut; S3: Solve the seepage differential equation of polymer solution using WENO method; S4: Compare with the upwind difference format to verify its computational performance.

2. The polymer flooding numerical simulation method based on the high-order WENO method according to claim 1, characterized in that: In S1, the polymer solution changes the oil-water mobility ratio by viscosity increase, and the viscosity correction is performed using the following cubic polynomial, m w =μ wi (1+rC p +sC p 2 +tC p 3 ) Among them, μ wi and μ w are the viscosity of the polymer solution before and after adding the polymer; C p is the polymer concentration; r, s, t are the equation parameters measured experimentally.

3. A polymer flooding numerical simulation method based on a high-order WENO method according to claim 1 or 2, characterized in that: The S2 comprises the following steps, S21: Assuming that the ratio of oil-water relative permeability varies linearly with water saturation, the simplified expression based on Archie's formula is: Among them, a and b are constant parameters to be determined, which are obtained by fitting experimental data; S w is the water saturation; S22: Considering the influence of residual oil saturation and irreducible water saturation, the modified model expression is: Among them, S wc is the bound water saturation, S or is the residual oil saturation; S23: Ignoring the effects of capillary force and gravity, the polymer solution is regarded as the water phase, and the modified viscosity equation and relative permeability are brought in. The flow equation expression is:

4. A polymer flooding numerical simulation method based on a high-order WENO method according to claim 1 or 2, characterized in that: The S3 comprises the following steps, S31: The polymer flooding process is regarded as a one-dimensional oil-water two-phase flow. The differential equation for the seepage of the polymer solution is expressed as: S32: After fixing the polymer concentration, the water content is considered as a function of water saturation as follows, f w (S w ,C p (n) ); S33: For the water phase seepage differential equation, the finite volume method is used in [x i-Δx / 2 ,x i+Δx / 2 ], we can get the following formula: S34: Let f(x, t) = f w (x, t), dividing both sides by Δx, we get, Among them, S wi It is expressed as the water saturation on the i unit; S wi is the domain of integration [x i-Δx / 2 ,x i+Δx / 2 ], that is, f(x,t) represents the numerical flux of the polymer solution; S35: Solve the water saturation equation using the WENO method; S36: In the WENO method, the smoothing factor β k To measure the smoothness of each template; S37: After using the WENO method to calculate the water saturation at time n+1, the seepage differential equation for the water component in the polymer solution is as follows: S38: The polymer concentration at time n+1 is calculated using the WENO method. Based on the polymer concentration and water saturation at time n+1, the calculation is repeated to obtain the water saturation and polymer concentration distribution at time n+2 and thereafter, and the parameter distribution at the outlet is obtained.

5. A polymer flooding numerical simulation method based on a high-order WENO method according to claim 4, characterized in that: The S35 comprises the following steps: S351: Use the combination of three templates S0, S1 and S2 to establish the WENO discrete format. S352: To ensure stability and consistency, the weight ω r The following conditions are met, S353: When all templates S r When all are in the smooth region, that is, any unit S in the template r (i) are not discontinuous units, then let ω r =y r ,get, S354: When template S r When located in the discontinuous region, ω r The calculation form of is as follows, The ε coefficient is introduced to prevent the denominator from being zero.

6. A polymer flooding numerical simulation method based on a high-order WENO method according to claim 4, characterized in that: In S36, the smoothness calculation formula of each template is as follows:

7. A polymer flooding numerical simulation method based on a high-order WENO method according to claim 1 or 2, characterized in that: The S4 comprises the following steps, S41: In the traditional upwind difference method, the one-dimensional convection equation is as follows: Among them, S w To substitute variables, Q, A, and φ are constants; S42: The concentration of the polymer solution at a fixed time n, the water content is a function of the water saturation f w (S w ); S43: Discretize the one-dimensional convection equation in space and time, using a grid of points x i and time step t n , the grid spacing is Δx, the time step is Δt, then the discretization form of the first-order upwind format is, S44: Use the formula of S43 to perform iterative calculation to obtain the water saturation at time n+1; S45: using the water saturation at time n+1 to alternately solve and perform iterative calculations to obtain the polymer solution concentration and its water content distribution at time n+1; S46: After calculating the polymer solution concentration at each fixed time, repeat the above steps to obtain the water saturation, polymer solution concentration and water content distribution at all times.

8. A device, characterized in that: Run the data processing method according to any one of claims 1 to 7.

9. A device comprising a memory, a processor, and an algorithm stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the data processing method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium storing a computer algorithm, characterized in that: When the computer algorithm is executed by a processor, the data processing according to any one of claims 1 to 7 is implemented.

Images