Quantitative treatment method for relative permeability curve

Abstract

The invention relates to the technical field of oil and gas reservoir engineering, in particular to a quantitative treatment method for a relative permeability curve. The method comprises the following steps: deriving a water content formula according to a relative permeability curve and a Darcy law; obtaining an oil content formula according to the obtained water content formula; s2, fitting oil production error coefficient parameters are introduced, and the corrected oil content is calculated in cooperation with the oil content formula obtained in the step S2 under the condition that the liquid production is not changed; determining the corrected water content according to the obtained corrected oil content, and establishing a relative permeability curve water phase permeability correction equation by combining the Darcy law; and S4, compiling a VBA relative permeability adjustment program according to an excel vba technology, and performing quantitative processing programming on the data calculated by the correction equation obtained in the step S4 through the adjustment program. According to the method provided by the invention, the relative permeability block is quantitatively adjusted by applying the method, and the working difficulty of historical fitting of traditional numerical simulation is solved.

Description

Technical Field

[0001] The invention relates to the technical field of oil and gas reservoir engineering, and in particular to a quantitative processing method for a phase permeability curve. Background Art

[0002] The phase permeability curve is an important parameter reflecting the flow of oil and water in the reservoir. After more than 60 years of development, the reservoir properties of Daqing Oilfield have changed to varying degrees. The phase permeability curve established based on early cores cannot fully reflect the current actual underground seepage conditions. Therefore, in the numerical simulation process, adjusting the phase permeability curve is an important means to achieve historical fitting accuracy.

[0003] At present, the phase permeability adjustment work faces the following main problems: First, there are many phase permeability curves and the working time is long. Taking the simulation area using 10 phase permeability curves as an example, the one-time manual adjustment of the phase permeability curve, including phase permeability data import, phase permeability data processing, phase permeability data detection and phase permeability data export, requires about 4 hours of working time; second, the phase permeability curve adjustment is qualitative and cannot be quantitative. Since the numerical simulation of the reservoir was carried out, the numerical simulation personnel have carried out the adjustment of the phase permeability curve from a qualitative perspective based on the results of the previous numerical simulation to meet the consistency between the adjustment direction of the phase permeability curve and the difference in the simulation results. The qualitative adjustment is directional and the adjustment amplitude is also blind. From the perspective of literature, the predecessors have not carried out the quantitative adjustment of the phase permeability curve; third, there is no tool software for the quantitative adjustment of the phase permeability. Since the application of numerical simulation, there has been no programmed work on the quantitative adjustment of the phase permeability curve. The root cause is that there is no theoretical formula of the phase permeability curve itself, and there is no theoretical formula between the phase permeability curve and the actual numerical simulation results to support it, which leads to the inability to achieve efficient processing of the phase permeability curve adjustment through programming. Therefore, in view of the above shortcomings, a quantitative processing method of phase permeability curve is proposed. Summary of the invention

[0004] 1. Technical issues to be resolved

[0005] The invention provides a quantitative processing method for a phase permeability curve, so as to overcome the problems in the prior art of taking a long time to adjust the qualitative adjustment curve, a large number of numerical simulation calculations and low numerical simulation efficiency in the process of numerical simulation of the phase permeability curve.

[0006] (II) Technical solution

[0007] In order to solve the above problems, the present invention provides a quantitative processing method of a phase permeability curve, comprising:

[0008] Step S1: derive the water content formula according to the phase permeability curve and Darcy's law;

[0009] Step S2: Obtain an oil content formula according to the water content formula obtained in step S1;

[0010] Step S3: introducing the fitted oil production error coefficient parameter, and calculating the corrected oil content in combination with the oil content formula obtained in step S2 under the condition that the liquid production remains unchanged;

[0011] Step S4: determining the corrected water content through the corrected oil content obtained in step S3, and establishing a water phase permeability correction equation of the phase permeability curve in combination with Darcy's law;

[0012] Step S5: compile a VBA phase permeability adjustment program based on Excel VBA technology, and program the data calculated by the correction equation obtained in step S4 to be quantitatively processed through the adjustment program.

[0013] Preferably, in step S1, the calculation formula for the water content is:

[0014]

[0015] Among them, f w , containing water; K rw , water phase permeability; K row , oil phase permeability; μ o , oil phase viscosity; μ w , water phase viscosity.

[0016] Preferably, in step S1, the calculation formula of Darcy's law is:

[0017] Q=KFI

[0018] Where Q is the seepage rate per unit time, F is the water flow section, I is the hydraulic gradient, and K is the permeability coefficient.

[0019] Preferably, in step S2, the water content and the oil content add up to 1, and the calculation formula of the oil content is:

[0020] f r =1-f w

[0021] Among them, f r , oil content; f w , containing water.

[0022] Preferably, in step S3, the fitted oil production error coefficient parameter is the ratio of the difference between the calculated value of the cumulative oil production according to the numerical simulation and the actual value of the cumulative oil production in the development history to the cumulative oil production in the development history.

[0023] Preferably, in step S4, the calculation formula of the water phase permeability correction equation of the phase permeability curve is:

[0024]

[0025] Among them, K rw , corrected water phase permeability; K rw , water phase permeability; K row , oil phase permeability; α, error coefficient of fitted oil production; μ o , oil phase viscosity; μ w , water phase viscosity.

[0026] Preferably, in step S5, the VBA phase permeability adjustment program includes three functional modules, namely, an initialization data module, a quantitative processing module and an auxiliary function module. The quantitative processing program process is to pass the acquired data through the initialization data module, the quantitative processing module and the auxiliary function module in sequence.

[0027] Preferably, the initialization data module inputs the data of oil phase viscosity and water phase viscosity; the quantitative processing module fills the fitted oil production error coefficient value into the corresponding cell based on the numerical simulation results, and executes the quantitative processing function to obtain the specific values ​​of the oil phase and water phase permeability of the modified quantitative processing phase permeability curve; the auxiliary function module converts the monotonicity detection and results of the processed data into a txt file.

[0028] (III) Beneficial effects

[0029] The quantitative processing method of the phase permeability curve provided by the present invention can quantitatively adjust the phase permeability block by applying the method, shorten the processing time of the phase permeability curve, make the phase permeability curve smooth, reduce the number of numerical simulations, and achieve a high degree of historical fitting compliance, thus solving the difficulties in traditional numerical simulation historical fitting work. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 It is a flow chart of a quantitative processing method of a phase permeability curve according to an embodiment of the present invention;

[0031] Figure 2 It is a system diagram of the VBA phase penetration adjustment program of an embodiment of the present invention. DETAILED DESCRIPTION

[0032] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. 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.

[0033] Figure 1 Flow chart of the quantitative processing method of the phase permeability curve of the embodiment of the present invention, as shown in FIG. Figure 1As shown, the present invention provides a quantitative processing method of a phase permeability curve, comprising:

[0034] Step S1: derive the water content formula according to the phase permeability curve and Darcy's law;

[0035] Step S2: Obtain an oil content formula according to the water content formula obtained in step S1;

[0036] Step S3: introducing the fitted oil production error coefficient parameter, and calculating the corrected oil content in combination with the oil content formula obtained in step S2 under the condition that the liquid production remains unchanged;

[0037] Step S4: determining the corrected water content through the corrected oil content obtained in step S3, and establishing a water phase permeability correction equation of the phase permeability curve in combination with Darcy's law;

[0038] Step S5: compile a VBA phase permeability adjustment program based on Excel VBA technology, and program the data calculated by the correction equation obtained in step S4 to be quantitatively processed through the adjustment program.

[0039] In practical applications, in step S1, the water content is calculated as:

[0040]

[0041] Among them, f w , containing water; K rw , water phase permeability; K row , oil phase permeability; μ o , oil phase viscosity; μ w , water phase viscosity.

[0042] In this processing method, in step S1, the calculation formula of Darcy's law is:

[0043] Q=KFI

[0044] Where Q is the seepage rate per unit time, F is the water flow section, I is the hydraulic gradient, and K is the permeability coefficient.

[0045] Preferably, the calculation formula of hydraulic gradient is:

[0046] I=h / L

[0047] Among them, I is the hydraulic gradient, h is the total head loss, and L is the length of the seepage path.

[0048] In practical applications, in step S2, since the sum of the water content and the oil content is 1, the calculation formula for the oil content is:

[0049] f r =1-f w

[0050] Among them, f r , oil content; f w , containing water.

[0051] In the processing method, in step S3, the fitted oil production error coefficient parameter is the ratio of the difference between the calculated value of the numerical simulation cumulative oil production and the actual value of the development history cumulative oil production to the development history cumulative oil production.

[0052] In practical applications, in step S4, the calculation formula of the water phase permeability correction equation of the phase permeability curve is:

[0053]

[0054] Among them, K rw , corrected water phase permeability; K rw , water phase permeability; K row , oil phase permeability; α, error coefficient of fitted oil production; μ o , oil phase viscosity; μ w , water phase viscosity.

[0055] In this treatment method, Figure 2 4 is a system diagram of the VBA phase permeation adjustment program of the embodiment of the present invention, as shown in FIG. Figure 2 As shown, in step S5, the VBA phase permeability adjustment program includes three functional modules, namely, an initialization data module, a quantitative processing module and an auxiliary function module. The quantitative processing program process is to pass the acquired data through the initialization data module, the quantitative processing module and the auxiliary function module in sequence.

[0056] It should be noted that in actual applications, the initialization data module inputs the data of oil phase viscosity and water phase viscosity; the quantitative processing module fills the fitted oil production error coefficient value into the corresponding cell based on the numerical simulation results, and executes the quantitative processing function to obtain the specific values ​​of the oil phase and water phase permeability of the modified quantitative processing phase permeability curve; the auxiliary function module converts the monotonicity test and results of the processed data into a txt file.

[0057] The quantitative processing method of the phase permeability curve provided by the present invention can quantitatively adjust the phase permeability block by applying the method, shorten the processing time of the phase permeability curve, make the phase permeability curve smooth, reduce the number of numerical simulations, and achieve a high degree of historical fitting compliance, thus solving the difficulties in traditional numerical simulation historical fitting work.

[0058] Block A is located in the Changyuan structure of Daqing Oilfield, with a total of 406 wells. It was put into development in 1960. The following is a detailed description of the implementation process of this method in combination with the quantitative processing method of the phase permeability curve:

[0059] Step 1: Based on the phase permeability curve and Darcy's law, derive the water content formula;

[0060] Step 2: Obtain the oil content formula according to the water content formula obtained in step S1;

[0061] Step 3: Introduce the fitted oil production error coefficient parameter, and calculate the corrected oil content in combination with the oil content formula obtained in step S2 under the condition that the liquid production remains unchanged;

[0062] Step 4: Determine the corrected water content through the corrected oil content obtained in step S3, and establish a water phase permeability correction equation of the phase permeability curve in combination with Darcy's law;

[0063] Step 5: Prepare a VBA phase permeability adjustment program based on Excel VBA technology, and use the adjustment program to program the quantitative processing of the data calculated by the correction equation obtained in step S4.

[0064] In this embodiment, 10 phase permeability curves are provided in the early stage of numerical simulation. On the basis of the established static model, the dynamic model is loaded, and the numerical simulation history fitting is performed. The trial calculation of the cumulative oil production error is -5.9%. According to the quantitative phase permeability curve processing method, 5.9% is used as the error coefficient. It takes 9 minutes to realize the data loading, phase permeability data import, phase permeability data quantitative processing, phase permeability data detection and phase permeability processing result data export of the phase permeability curve; the numerical simulation calculation is performed again according to the phase permeability curve after quantitative processing, and the calculation result is that the cumulative oil production error is -1.14%; the second phase permeability quantitative processing is performed on the basis of the most original phase permeability, and 10% is used as the error coefficient. The phase permeability processing and numerical simulation trial calculation are performed again. The second numerical simulation calculation result is that the cumulative oil production error is 1.4%; the third phase permeability quantitative processing, the linear difference method is used to determine the optimal adjustment range of 7.74% (Formula 1), and the cumulative oil production error of the third numerical simulation calculation result is -0.12%.

[0065]

[0066] In practical applications, the quantitative processing method of phase permeability curves is used to achieve high-quality and efficient numerical simulation, reduce the number of numerical simulations, and improve the accuracy of numerical simulations.

[0067] The quantitative processing method of the phase permeability curve provided by the present invention quantitatively adjusts the phase permeability block by applying the method, the phase permeability curve processing time is short, the phase permeability curve shape is smooth, the number of numerical simulation times is reduced, the historical fitting conformity is high, and the difficulties in the traditional numerical simulation historical fitting work are solved. In practical applications, the phase permeability curve processing time is reduced, and the programmed processing of the phase permeability curve only takes a few minutes, which greatly improves the work efficiency compared with the traditional processing time of 4 hours; the number of historical fittings is reduced, the traditional reservoir numerical simulation history number is more than 10 times, and the number of historical fittings corresponding to the quantitative processing of the phase permeability curve is only 3 times, which improves the historical fitting work efficiency by more than 3 times; the historical fitting accuracy is improved, the traditional historical fitting accuracy is within 2%, and the historical fitting accuracy of the quantitative programmed processing of the phase permeability curve can be controlled within 0.3%.

[0068] The above implementation modes are only used to illustrate the present invention, but not to limit the present invention. Ordinary technicians in the relevant technical field can make various changes and modifications without departing from the spirit and scope of the present invention. Therefore, all equivalent technical solutions also belong to the scope of the present invention. The patent protection scope of the present invention should be defined by the claims.

Claims

1. A quantitative processing method for phase permeability curve, characterized in that: include: Step S1: derive the water content formula according to the phase permeability curve and Darcy's law; Step S2: Obtain an oil content formula according to the water content formula obtained in step S1; Step S3: introducing the fitted oil production error coefficient parameter, and calculating the corrected oil content in combination with the oil content formula obtained in step S2 under the condition that the liquid production remains unchanged; Step S4: determining the corrected water content through the corrected oil content obtained in step S3, and establishing a water phase permeability correction equation of the phase permeability curve in combination with Darcy's law; Step S5: compile a VBA phase permeability adjustment program based on Excel VBA technology, and program the data calculated by the correction equation obtained in step S4 to be quantitatively processed through the adjustment program.

2. The quantitative processing method of phase permeability curve according to claim 1 is characterized in that: In step S1, the calculation formula of the water content is: Among them, f w , containing water; K rw , water phase permeability; K row , oil phase permeability; μ o , oil phase viscosity; μ w , water phase viscosity.

3. The quantitative processing method of phase permeability curve according to claim 1 is characterized in that: In step S1, the calculation formula of Darcy's law is: Q=KFI Where Q is the seepage rate per unit time, F is the water flow section, I is the hydraulic gradient, and K is the permeability coefficient.

4. The quantitative processing method of phase permeability curve according to claim 1 is characterized in that: In step S2, the water content and the oil content add up to 1, and the calculation formula of the oil content is: f r =1-f w Among them, f r , oil content; f w , containing water.

5. The quantitative processing method of phase permeability curve according to claim 1 is characterized in that: In step S3, the fitted oil production error coefficient parameter is the ratio of the difference between the calculated value of the numerical simulation cumulative oil production and the actual value of the development history cumulative oil production to the development history cumulative oil production.

6. The quantitative processing method of phase permeability curve according to claim 1 is characterized in that: In step S4, the calculation formula of the water phase permeability correction equation of the phase permeability curve is: Among them, K rw , corrected water phase permeability; K rw , water phase permeability; K row , oil phase permeability; α, error coefficient of fitted oil production; μ o , oil phase viscosity; μ w , water phase viscosity.

7. The quantitative processing method of phase permeability curve according to claim 1 is characterized in that: In step S5, the VBA phase permeability adjustment program includes three functional modules, namely, an initialization data module, a quantitative processing module and an auxiliary function module. The quantitative processing program process is to pass the acquired data through the initialization data module, the quantitative processing module and the auxiliary function module in sequence.

8. The quantitative processing method of phase permeability curve according to claim 7 is characterized in that: The initialization data module inputs the data of oil phase viscosity and water phase viscosity; The quantitative processing module fills the fitted oil production error coefficient value into the corresponding cell according to the numerical simulation results, and performs the quantitative processing function to obtain the specific values ​​of the oil phase and water phase permeability of the modified quantitative processing phase permeability curve; The auxiliary function module is to detect the monotonicity of the processed data and convert the results into a txt file.

Images