A two-parameter phase infiltration interpolation method based on normalization and kriging interpolation

By employing a method based on normalization and Kriging interpolation, and selecting an appropriate combination of temperature and interfacial tension for interpolation of the phase permeation curve, this method solves the technical problem in existing technologies that fail to simultaneously consider temperature and interfacial tension. It addresses the technical issues related to temperature, phase point, and the shape and endpoints of the phase permeation curve, thus resolving a technical challenge unresolved in existing technologies and demonstrating its practical contribution to solving technical problems.

CN119887506BActive Publication Date: 2026-01-06CHINA PETROLEUM & CHEMICAL CORP +1
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202311320055.6
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2023-10-12
Publication Date
2026-01-06
Estimated Expiration
2043-10-12

AI Technical Summary

Technical Problem

In existing interpolation methods for relative permeation curves, the effects of temperature and interfacial tension on the relative permeation curve cannot be considered simultaneously, resulting in large errors. Furthermore, the influence of temperature on the curve shape and endpoints is not fully reflected.

Method used

A two-parameter method based on normalization and Kriging interpolation was adopted. Multiple relative permeation curves under different temperatures and interfacial tensions were obtained through experimental testing. The closest combination of temperature and interfacial tension was selected for normalization and logarithmic interpolation. Finally, the relative permeation endpoint value was calculated using Kriging interpolation.

Benefits of technology

It achieves accurate interpolation of the morphology and endpoint values ​​of the phase permeation curve under the combined effects of temperature and interfacial tension, reducing errors and improving interpolation accuracy.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119887506B_ABST
    Figure CN119887506B_ABST
Patent Text Reader

Abstract

The application provides a two-parameter phase permeability interpolation method based on normalization and Kriging interpolation, which comprises the following steps: according to the relationship between the temperature and the interfacial tension of to-be-interpolated phase permeability and the temperature and the interfacial tension of a plurality of measured phase permeability curves, four phase permeability curves meeting the requirements are selected from the measured phase permeability curves, and the selected curves are normalized; then according to the relationship between the to-be-interpolated temperature and the temperature of the selected phase permeability curve and the relationship between the to-be-interpolated interfacial tension and the interfacial tension of the selected phase permeability curve, logarithmic interpolation is performed on the normalized phase permeability curves respectively, so that the normalized phase permeability curve under the condition of the to-be-interpolated temperature and the interfacial tension is obtained; then the end point values of all the measured phase permeability curves are taken as data bodies, and the end point values of the phase permeability under the condition of the to-be-interpolated temperature and the interfacial tension are calculated by using Kriging interpolation; finally, the interpolated normalized phase permeability curve and the end point values of the phase permeability curve are combined, so that the actual phase permeability curve after interpolation is obtained, thereby realizing the common interpolation of the shape and the end point of the phase permeability curve. The application simultaneously considers the influence of temperature and interfacial tension on the shape and the end point value of the phase permeability curve, and provides a method for considering the phase permeability curve interpolation under the common action of two parameters.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of oilfield development technology, and in particular to a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation. Background Technology

[0002] Currently, in the interpolation method for relative permeability curves, the influence of temperature and interfacial tension on relative permeability is achieved through stepwise interpolation. The specific steps include: ① Defining a basic relative permeability curve at a certain temperature, and the relative permeability endpoint values ​​corresponding to different temperatures, and characterizing the influence of temperature on relative permeability by the change of relative permeability endpoint values ​​at different temperatures; ② Interpolating the relative permeability endpoints at different temperatures to obtain the endpoint values ​​of the relative permeability curve at the desired temperature, thus completing the temperature interpolation of the relative permeability curve; ③ Defining two relative permeability curves under water drive and ultra-low interfacial tension conditions respectively; ④ Completing the interfacial tension interpolation of the relative permeability curve based on the relationship between the interfacial tension value to be interpolated and the interfacial tension values ​​of water drive and ultra-low interfacial tension.

[0003] Currently used interpolation methods for relative permeability assume that temperature changes only affect the endpoint values ​​of the relative permeability curve. Therefore, the basic shape of the relative permeability curve obtained through this interpolation method remains unchanged; it is merely a scaling or translation of the basic curve. However, numerous experimental results show that temperature affects both the shape and endpoints of the relative permeability curve. Furthermore, in currently used interpolation methods, the effects of temperature and interfacial tension on the endpoint values ​​of the relative permeability curve cannot be considered simultaneously but must be performed step by step, resulting in significant errors. Summary of the Invention

[0004] In view of the above problems, the present invention is proposed to provide a two-parameter phase infiltration interpolation method based on normalization and kriging interpolation to overcome or at least partially solve the above problems.

[0005] According to one aspect of the present invention, a two-parameter phase permeation interpolation method based on normalization and kriging interpolation is provided, the interpolation method comprising:

[0006] Step 1: Obtain multiple phase permeation curves under different temperatures and interfacial tensions through experimental testing;

[0007] Step 2: From the multiple relative permeability curves, select the temperature value T that differs the least from the interpolation temperature T. a T b The interfacial tension value σ that differs least from the interfacial tension σ to be interpolated. m σ n Temperature and interfacial tension values ​​were selected from the measured phase permeation curves, corresponding to the first phase permeation curve (T). a , σ m ), second phase permeation curve (T) b , σ m ), third phase permeation curve (T)a , σ n ), fourth phase permeation curve (T b , σ n The four relative permeability curves were obtained and normalized.

[0008] Step 3: Based on the first relative permeability curve (T) in Step 2 a , σ m ), the second phase permeation curve (T) b , σ m Temperature value T a T b The relationship between the interfacial tension σ and the interfacial tension σ at temperature T is obtained by logarithmic interpolation of the relative permeability curve. m Normalized phase infiltration curve;

[0009] Step 4: Based on the third phase permeation curve (T) described in Step 2 a , σ n The fourth phase permeation curve (T) b , σ n Temperature value T a T b The relationship between the interfacial tension σ and the interfacial tension σ at temperature T is obtained by logarithmic interpolation of the relative permeability curve. n Normalized phase infiltration curve;

[0010] Step 5: The (T, σ) obtained in step 3 m The normalized phase infiltration curve under the given conditions and the curve obtained in step 4 (T, σ) n Based on the normalized phase permeation interpolation curve under the given conditions, and according to the interfacial tension value σ in the two phase permeation curves... m σ n The relationship between the interfacial tension σ and the interfacial tension σ is used to perform logarithmic interpolation of the phase permeation curve, resulting in a normalized phase permeation interpolation curve corresponding to temperature T and interfacial tension σ.

[0011] Step 6: Using all the measured relative permeability curve endpoint values ​​from Step 1 as the data volume, calculate the relative permeability endpoint values ​​corresponding to temperature T and interfacial tension σ using Kriging interpolation.

[0012] Step 7: Based on the normalized phase permeation interpolation curve obtained in Step 5 under the conditions of temperature T and interfacial tension σ, and the interpolated phase permeation endpoint value obtained in Step 6, the actual phase permeation curve corresponding to temperature T and interfacial tension σ is obtained by reverse calculation.

[0013] Optionally, in step 1, the multiple phase permeation curves obtained through experimental testing include at least two sets of temperature test points and at least two sets of interfacial tension test points.

[0014] Optionally, in step 2, the temperature value T that differs the least from the temperature T to be interpolated is selected. a T b The selection condition satisfies: among all temperature test points T(1,2,...,n), T a T b The square root deviation from the temperature T to be interpolated is the smallest. Its value is the smallest.

[0015] Optionally, in step 2, the temperature value σ that differs least from the interpolated temperature T is... m σ n The selection satisfies the following condition: among all temperature test points σ(1,2,……,n), the square root deviation between σ1, σ2 and the temperature to be interpolated σ is the smallest. Its value is the smallest.

[0016] Optionally, in step 2, the selection of the four measured relative permeability curves is based on the temperature value T selected earlier in this step. a T b With interfacial tension value σ m σ n Pairwise combinations, select those corresponding to (T) a , σ m ), (T b , σ m ), (T a , σ n ), (T b , σ n ) 4 measured relative permeability curves.

[0017] Optionally, in step 3, the two interpenetration curves involved in the interpolation have the same interfacial tension σ. m However, the temperatures are respectively T a T b The normalized phase permeation curve.

[0018] Optionally, in step 4, the two interpenetration curves involved in the interpolation have the same interfacial tension σ. n However, the temperatures are respectively T a T b The normalized phase permeation curve.

[0019] Optionally, in step 5, the two interpenetration curves involved in the interpolation are obtained through interpolation in steps 3 and 4, respectively. The two interpenetration curves have the same temperature T, but their interfacial tension values ​​correspond to σ. m σ n .

[0020] Optionally, in step 6, when interpolating the relative permeation endpoint values ​​corresponding to temperature T and interfacial tension σ, the interpolation method is Kriging interpolation, and the data volume used for interpolation is the endpoint values ​​of all measured relative permeation curves in step 1.

[0021] Optionally, in step 6, the relative permeability endpoint values ​​include bound water saturation, residual oil saturation, oil phase relative permeability, and water phase relative permeability.

[0022] The present invention provides a two-parameter relative permeability interpolation method based on normalization and Kriging interpolation. According to the relationship between the interpolation temperature, interfacial tension, and the temperature and interfacial tension values ​​of multiple measured relative permeability curves, four relative permeability curves that meet the requirements are selected from the measured curves. These selected curves are then normalized. Next, based on the relationship between the interpolation temperature and the selected relative permeability curve temperature values, and the relationship between the interfacial tension and the selected relative permeability curve interfacial tension values, logarithmic interpolation is performed on the normalized relative permeability curves to obtain the normalized relative permeability curves under the interpolation temperature and interfacial tension conditions. Then, using the endpoint values ​​of all measured relative permeability curves as the data volume, Kriging interpolation is used to calculate the relative permeability endpoint values ​​under the interpolation temperature and interfacial tension conditions. Finally, by combining the interpolated normalized relative permeability curves with the endpoint values ​​of the relative permeability curves, the actual relative permeability curve after interpolation can be obtained, thus achieving joint interpolation of the shape and endpoints of the relative permeability curves. This invention takes into account the effects of temperature and interfacial tension on the shape and endpoint values ​​of the relative permeation curve, providing a way to interpolate the relative permeation curve considering the combined effect of two parameters.

[0023] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and in order to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description

[0024] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the following description of the embodiments will be briefly introduced. Obviously, the 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.

[0025] Figure 1 The flowchart illustrates a specific example of a two-parameter phase infiltration interpolation method based on normalization and Kriging interpolation provided in this embodiment of the invention.

[0026] Figure 2 This invention relates to a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation, which provides phase permeation curves corresponding to different temperatures and interfacial tension values ​​obtained through experiments.

[0027] Figure 3 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, four measured phase permeation curves that meet the interpolation requirements are provided.

[0028] Figure 4 This is the normalized result of four measured phase permeation curves that meet the interpolation requirements in a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation according to the present invention.

[0029] Figure 5 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, the normalized phase permeation curve obtained by interpolation under the conditions of a temperature of 95°C and an interface tension of 0.1 Nm / m is obtained.

[0030] Figure 6 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, the normalized phase permeation curve obtained by interpolation under the conditions of a temperature of 95°C and an interface tension of 1.0 Nm / m is obtained.

[0031] Figure 7 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, the normalized phase permeation curve obtained by interpolation under the conditions of a temperature of 95°C and an interface tension of 0.5 Nm / m is obtained.

[0032] Figure 8 This invention relates to a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation, which provides the actual phase permeation curve obtained under the conditions of a temperature of 95°C and an interface tension of 0.5 Nm / m.

[0033] Figure 9 This invention relates to a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation, which provides phase permeation curves corresponding to different temperatures and interfacial tension values ​​obtained through experiments.

[0034] Figure 10 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, four measured phase permeation curves that meet the interpolation requirements are provided.

[0035] Figure 11 This is the normalized result of four measured phase permeation curves that meet the interpolation requirements in a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation according to the present invention.

[0036] Figure 12 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, the normalized phase permeation curve obtained by interpolation under the conditions of a temperature of 95°C and an interface tension of 0.1 Nm / m is obtained.

[0037] Figure 13 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, the normalized phase permeation curve obtained by interpolation under the conditions of a temperature of 95°C and an interface tension of 1.0 Nm / m is obtained.

[0038] Figure 14 In a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation of the present invention, the normalized phase permeation curve obtained by interpolation under the conditions of a temperature of 95°C and an interface tension of 0.5 Nm / m is obtained.

[0039] Figure 15 This invention relates to a two-parameter phase permeation interpolation method based on normalization and Kriging interpolation, which provides the actual phase permeation curve obtained under the conditions of a temperature of 95°C and an interface tension of 0.5 Nm / m. Detailed Implementation

[0040] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0041] The terms "comprising" and "having," and any variations thereof, in the specification, embodiments, claims, and drawings of this invention are intended to cover non-exclusive inclusion, such as including a series of steps or units.

[0042] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments.

[0043] Example 1

[0044] In a specific embodiment 1 of the present invention, the invention is applied. Figure 1 This is a flowchart of a two-parameter phase infiltration interpolation method based on normalization and Kriging interpolation. In Example 1, the temperature to be interpolated is set to 95℃, and the interfacial tension value to be interpolated is 0.5 Nm / m.

[0045] In step 101, multiple relative permeation curves were obtained through experimental testing at temperatures of 65℃, 120℃, and 200℃, and interfacial tensions of 1.0 mN / m, 0.1 mN / m, and 0.01 mN / m, respectively. Figure 2 As shown.

[0046] In step 102, two measured temperature values ​​with the smallest difference from the interpolation temperature of 90℃ are selected: 65℃ and 120℃; two measured interfacial tension values ​​with the smallest difference from the interfacial tension of 0.5 Nm / m are selected: 0.1 mN / m and 1.0 mN / m. The selected measured temperature values ​​and measured interfacial tension values ​​are then combined pairwise to obtain the following... Figure 3 The four measured relative permeability curves shown correspond to the following temperatures and interfacial tension values: (65℃, 0.1 mN / m), (120℃, 0.1 mN / m), (65℃, 1.0 mN / m), and (120℃, 1.0 mN / m). These four measured relative permeability curves were then normalized. The normalized four relative permeability curves are shown below. Figure 4 As shown.

[0047] In step 103, from Figure 4 Two normalized relative permeability curves with the same interfacial tension of 0.1 Nm / m were selected: the relative permeability curves corresponding to (65℃, 0.1 mN / m) and (120℃, 0.1 mN / m). Interpolation was performed based on the relationship between the interpolation temperature of 95℃ and the measured temperatures of 65℃ and 120℃ to obtain the relative permeability curves under the condition of 95℃ and an interfacial tension of 0.1 Nm / m, as shown below. Figure 5 As shown.

[0048] In step 104, from Figure 4 Two normalized relative permeability curves with the same interfacial tension of 1.0 Nm / m were selected: the relative permeability curves corresponding to (65℃, 1.0 mN / m) and (120℃, 1.0 mN / m). Interpolation was performed based on the relationship between the interpolation temperature of 95℃ and the measured temperatures of 65℃ and 120℃ to obtain the relative permeability curves under the condition of 95℃ and an interfacial tension of 1.0 Nm / m, as shown below. Figure 6 As shown.

[0049] In step 105, the two relative permeability curves obtained in step 103 (95℃, 0.1 Nm / m) and step 104 (95℃, 1.0 Nm / m) are interpolated based on the relationship between the interfacial tension to be interpolated (0.5 Nm / m) and the measured interfacial tensions (0.1 Nm / m and 1.0 Nm / m) to obtain the relative permeability curve under the condition of 95℃ and interfacial tension of 0.5 Nm / m. Figure 7 As shown.

[0050] In step 106, using the endpoint values ​​of all nine measured relative permeability curves from step 101 as the data volume, Kriging interpolation is used to obtain the relative permeability endpoint values ​​at 95℃ and 0.5 mN / m: bound water saturation Swc, residual oil saturation Sor, oil phase relative permeability Kro, and water phase relative permeability Krw, as shown in the table below:

[0051]

[0052] In step 107, combining the normalized phase permeation interpolation curve obtained in step 105 under the conditions of 95℃ temperature and 0.5 Nm / m interfacial tension, and the interpolated phase permeation endpoint values ​​obtained in step 106 under the conditions of 95℃ temperature and 0.5 Nm / m interfacial tension, the actual phase permeation curve at 95℃ temperature and 0.5 Nm / m interfacial tension can be obtained by reverse calculation, as follows: Figure 8 As shown.

[0053] Example 2

[0054] In a specific embodiment 2 of the present invention, the invention is applied. Figure 1 This is a flowchart of a two-parameter phase infiltration interpolation method based on normalization and Kriging interpolation. In Example 2, the interpolation temperature is set to 95℃, and the interfacial tension value is set to 0.5 Nm / m.

[0055] In step 101, multiple relative permeation curves were obtained through experimental testing at temperatures of 65℃, 120℃, and 200℃, and interfacial tensions of 1.0 mN / m, 0.1 mN / m, and 0.01 mN / m, respectively. Figure 9 As shown.

[0056] In step 102, two measured temperature values ​​with the smallest difference from the interpolation temperature of 90℃ are selected: 65℃ and 120℃; two measured interfacial tension values ​​with the smallest difference from the interfacial tension of 0.5 Nm / m are selected: 0.1 mN / m and 1.0 mN / m. The selected measured temperature values ​​and measured interfacial tension values ​​are then combined pairwise to obtain the following... Figure 10 The four measured relative permeability curves shown correspond to the following temperatures and interfacial tension values: (65℃, 0.1 mN / m), (120℃, 0.1 mN / m), (65℃, 1.0 mN / m), and (120℃, 1.0 mN / m). These four measured relative permeability curves were then normalized. The normalized four relative permeability curves are shown below. Figure 11 As shown.

[0057] In step 103, from Figure 11 Two normalized relative permeability curves with the same interfacial tension of 0.1 Nm / m were selected: the relative permeability curves corresponding to (65℃, 0.1 mN / m) and (120℃, 0.1 mN / m). Interpolation was performed based on the relationship between the interpolation temperature of 95℃ and the measured temperatures of 65℃ and 120℃ to obtain the relative permeability curves under the condition of 95℃ and an interfacial tension of 0.1 Nm / m, as shown below. Figure 12 As shown.

[0058] In step 104, from Figure 11Two normalized relative permeability curves with the same interfacial tension of 1.0 Nm / m were selected: the relative permeability curves corresponding to (65℃, 1.0 mN / m) and (120℃, 1.0 mN / m). Interpolation was performed based on the relationship between the interpolation temperature of 95℃ and the measured temperatures of 65℃ and 120℃ to obtain the relative permeability curves under the condition of 95℃ and an interfacial tension of 1.0 Nm / m, as shown below. Figure 13 As shown.

[0059] In step 105, the two relative permeability curves obtained in step 103 (95℃, 0.1 Nm / m) and step 104 (95℃, 1.0 Nm / m) are interpolated based on the relationship between the interfacial tension to be interpolated (0.5 Nm / m) and the measured interfacial tensions (0.1 Nm / m and 1.0 Nm / m) to obtain the relative permeability curve under the condition of 95℃ and interfacial tension of 0.5 Nm / m. Figure 14 As shown.

[0060] In step 106, using the endpoint values ​​of all nine measured relative permeability curves from step 101 as the data volume, Kriging interpolation is used to obtain the relative permeability endpoint values ​​at 95℃ and 0.5 mN / m: bound water saturation Swc, residual oil saturation Sor, oil phase relative permeability Kro, and water phase relative permeability Krw, as shown in the table below:

[0061]

[0062] In step 107, combining the normalized phase permeation interpolation curve obtained in step 105 under the conditions of 95℃ temperature and 0.5 Nm / m interfacial tension, and the interpolated phase permeation endpoint values ​​obtained in step 106 under the conditions of 95℃ temperature and 0.5 Nm / m interfacial tension, the actual phase permeation curve at 95℃ temperature and 0.5 Nm / m interfacial tension can be obtained by reverse calculation, as follows: Figure 15 As shown.

[0063] Beneficial Effects: The dual-parameter relative permeability interpolation method based on normalization and Kriging interpolation of this invention selects four relative permeability curves that meet the requirements from the measured relative permeability curves based on the relationship between the temperature to be interpolated, the interfacial tension, and the temperature and interfacial tension values ​​of multiple measured relative permeability curves. These selected curves are then normalized. Logarithmic interpolation is performed on the normalized relative permeability curves based on the relationship between the temperature to be interpolated and the selected relative permeability curve temperature values, and between the interfacial tension to be interpolated and the selected relative permeability curve interfacial tension values, respectively, to obtain the normalized relative permeability curves under the conditions of the temperature and interfacial tension to be interpolated. Then, using the endpoint values ​​of all measured relative permeability curves as the data volume, Kriging interpolation is used to calculate the relative permeability endpoint values ​​under the conditions of the temperature and interfacial tension to be interpolated. Finally, by combining the interpolated normalized relative permeability curves with the endpoint values ​​of the relative permeability curves, the actual relative permeability curve after interpolation can be obtained, thus achieving joint interpolation of the shape and endpoints of the relative permeability curves. This invention takes into account the effects of temperature and interfacial tension on the shape and endpoint values ​​of the relative permeation curve, providing a way to interpolate the relative permeation curve considering the combined effect of two parameters.

[0064] The above specific embodiments further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A two-parameter phase permeability interpolation method based on normalized and Kriging interpolation, characterized in that, The interpolation method comprises: Step 1: a plurality of phase permeability curves under different temperatures and different interfacial tensions are obtained by experimental testing; Step 2: From the multiple relative permeability curves, select the temperature value T that differs the least from the interpolation temperature T. a T b The interfacial tension value σ that differs least from the interfacial tension σ to be interpolated. m σ n Temperature and interfacial tension values ​​were selected from the measured phase permeation curves, corresponding to the first phase permeation curve (T). a , σ m ), second phase permeation curve (T) b , σ m ), third phase permeation curve (T) a , σ n ), fourth phase permeation curve (T b , σ n The four relative permeability curves were obtained and normalized. Step 3: According to the relationship between the temperature values T a , T m of the first phase permeability curve (T b , σ m ) and the second phase permeability curve (T a , σ b ) and the temperature T to be interpolated, logarithmic interpolation of the phase permeability curves is performed to obtain a normalized phase permeability interpolation curve corresponding to the temperature T and the interfacial tension σ m . Step 4: According to the relationship between the temperature values T a , T n of the third phase permeability curve (T b , σ n ) and the fourth phase permeability curve (T a , σ b ) and the temperature T to be interpolated, logarithmic interpolation of the phase permeability curves is performed to obtain a normalized phase permeability interpolation curve corresponding to the temperature T and the interfacial tension σ n ; Step 5: Based on the normalized phase permeability interpolation curve under the condition of (T, σ m ) obtained in Step 3 and the normalized phase permeability interpolation curve under the condition of (T, σ n ) obtained in Step 4, according to the relationship between the interfacial tension values σ m , σ n and the interfacial tension σ to be interpolated in the two phase permeability curves, the phase permeability interpolation curve corresponding to the temperature T and the interfacial tension σ is obtained by logarithmic interpolation. Step 6: the end point value data of all the measured phase permeability curves in step 1 are taken as a data body, and the end point value corresponding to the temperature T and the interfacial tension σ is calculated by using Kriging interpolation; Step 7: according to the normalized phase permeability interpolation curve under the condition of the temperature T and the interfacial tension σ obtained in step 5 and the interpolated end point value obtained in step 6, the actual phase permeability curve corresponding to the temperature T and the interfacial tension σ is obtained.

2. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 1, the plurality of phase permeability curves obtained by experimental testing contain at least two groups of temperature test points and at least two groups of interfacial tension test points.

3. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 2, the temperature value T a , T b is selected to meet the condition that among all the temperature test points T(1, 2, …, n), T a , T b has the minimum square root deviation from the temperature to be interpolated T, and the value of T has the minimum value.

4. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 2, the temperature value σ that is closest to the temperature T to be interpolated m , σ n is selected to satisfy the following condition: among all the temperature test points σ (1, 2, …, n), the square root deviation between σ1, σ2 and the temperature σ to be interpolated is the smallest, the value of the function is minimum.

5. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 2, the selection of the four measured relative permeability curves is based on the temperature value T selected earlier in this step. a T b With interfacial tension value σ m σ n Pairwise combinations, select those corresponding to (T) a , σ m ), (T b , σ m ), (T a , σ n ), (T b , σ n ) 4 measured relative permeability curves.

6. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 3, the two phase curves involved in the interpolation are the normalized phase curves with the same interfacial tension σ m but temperatures T a , T b , respectively.

7. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 4, the two phase curves involved in the interpolation are the normalized phase curves with the same interfacial tension σ n but temperatures T a , T b , respectively.

8. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 5, two interfacial tension curves participating in the interpolation are obtained by step 3, step 4 interpolation respectively, and the two interfacial tension curves have the same temperature T, but the interfacial tension values correspond to σ m , σ n respectively.

9. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 6, when the end point value corresponding to the temperature T and the interfacial tension σ is interpolated and calculated, the interpolation method adopts Kriging interpolation, and the data body used for interpolation is the end point value of all the measured phase permeability curves in step 1.

10. The two-parameter phase permeability interpolation method based on normalization and Kriging interpolation of claim 1, wherein, In step 6, the phase permeability end point value includes irreducible water saturation, residual oil saturation, oil phase relative permeability and water phase relative permeability.

Citation Information

Patent Citations

  • Theoretical water drive curve determination method and device

    CN112883530A

  • Water-containing gas reservoir development effect evaluation method considering relative permeability influence

    CN116245423A