Method and device for determining polymer-driven movable unit, medium and processor

Abstract

The embodiment of the invention provides a method and device for determining a polymer flooding movable unit, a medium and a processor, and belongs to the technical field of oil development. The method comprises the following steps: acquiring a plurality of core samples in a working area and relative permeability experiment data of the core samples, and obtaining a fitting function of water saturation parameters and water production rate parameters in the working area according to the relative permeability experiment data; according to the fitting function and a preset water production rate parameter threshold interval, a corresponding water saturation parameter threshold interval is calculated, and then a deep lateral resistivity parameter threshold interval corresponding to a water saturation parameter is solved through an Azig formula; and obtaining parameter curves of the interval transit time and the deep lateral resistivity at different depth positions, and determining a polymer flooding movable unit according to whether the parameters of the deep lateral resistivity and the interval transit time at each depth position fall into corresponding parameter threshold intervals at the same time or not. According to the method for determining the polymer-driven movable unit, the movable unit capable of being driven by the polymer can be found from the used reservoir.

Description

Technical Field

[0001] The invention relates to the technical field of petroleum development, and in particular to a method, device, medium and processor for determining a polymer drive movable unit. Background Art

[0002] Polymer flooding is a tertiary oil recovery technology developed in recent years. It is actually a method of increasing oil recovery by adding water-soluble polymers to the injected water to increase the viscosity of the injected water, improve the mobility ratio, and expand the swept volume. See attached Figure 1 .

[0003] For areas that have been produced once or twice, the situation faced by polymer flooding is that there are many displacement layers, strong heterogeneity, complex injection-production relationship, and high water content in the reservoir. Especially for the situation of comprehensive regional production, the microscopic residual oil after water flooding is extremely scattered, and the residual oil exists in the form of clusters, islands, corners, oil films, and throat residual oil. See attached Figure 2 .

[0004] In order to better explore the potential of remaining oil, restore the vitality of oil fields, and achieve the goals of stable production and increased reserves, it is necessary to study the seepage units of oil fields and fully understand the utilization laws of different flow units in different periods, so as to establish standards for selecting available units and guide polymer flooding development. Summary of the invention

[0005] The purpose of the embodiment of the present invention is to provide a method that can identify the areas that have been mobilized, especially when the areas have been fully mobilized, and the movable units that can be driven by polymers. In order to achieve the above purpose, the embodiment of the present invention provides a method for determining the movable units driven by polymers, comprising:

[0006] Step 1: Obtain phase permeability experimental data of multiple core samples, and obtain fitting functions of water saturation parameters and water production rate parameters based on the phase permeability experimental data. The core samples are obtained by coring in the working area;

[0007] Step 2: Calculate the water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and the preset water production rate parameter threshold interval, and then obtain the deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using the Archie formula;

[0008] Step 3: Acquire the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions. When the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, determine the corresponding depth position as the polymer drive movable position, and the depth interval composed of the polymer drive movable position is the polymer drive movable unit.

[0009] Preferably, the preset water production rate parameter threshold interval includes different sub-intervals, and the corresponding water saturation parameter threshold sub-interval and deep lateral resistivity parameter threshold sub-interval are calculated respectively according to the water production rate parameter threshold sub-interval, and the effective producing level of the polymer drive movable unit is determined according to the deep lateral resistivity parameter falling into different deep lateral resistivity parameter threshold sub-intervals, wherein those identified as the same effective producing level are produced using the polymer formula of the corresponding level.

[0010] Furthermore, the effective production levels include, from low to high: the first level, the second level, the third level, the fourth level, and the fifth level, and from the first level to the fifth level, the upper and lower limits of the deep lateral resistivity parameter threshold sub-interval increase step by step, and the upper and lower limits of the water production rate parameter threshold sub-interval and the water saturation parameter threshold sub-interval decrease step by step.

[0011] Furthermore, the preset lower limits of the water production rate parameter threshold sub-intervals of the first level to the fifth level are: 90%, 80%, 60%, 40%, 10% respectively;

[0012] The lower limit of the threshold sub-interval of the water saturation parameter of each level is calculated by the lower limit of the threshold sub-interval of the water production rate parameter of the corresponding level; and

[0013] The lower limit of the threshold sub-interval of the deep lateral resistivity parameter at each level is calculated by the lower limit of the threshold sub-interval of the water saturation parameter at the corresponding level.

[0014] Preferably, before performing step 1, it is first identified whether each target reservoir has been produced by using a preset determination condition, and then the polymer flooding movable units in the produced reservoir are determined by steps 1-3.

[0015] Furthermore, the preset determination conditions include one or more of a natural potential change method, a gas composition determination method, and an oil content identification method.

[0016] Preferably, the process of obtaining the fitting function in step 1 includes: determining the water saturation and water production rate cross-plot through the phase permeability experimental data, and obtaining the formula of the fitting function based on the water saturation and water production rate cross-plot. The water production rate of each sampling point is calculated by the water relative permeability, oil relative permeability, absolute permeability, cross-sectional area, water viscosity, oil viscosity, and pressure difference of the corresponding sampling point; the water saturation of each sampling point is corrected by the formula of the fitting function; for the point where no sampling is performed, its water saturation and water production rate are obtained by the fitting function.

[0017] Preferably, in step 2, the process of calculating the water saturation parameter threshold interval includes: on the function curve of the fitting function, intercepting a portion of the curve within the upper and lower limits of the preset water production rate parameter threshold interval, and determining the upper and lower limits of the water saturation parameter threshold interval based on the partial curve.

[0018] Preferably, in step 2, Archie's formula is used to determine the corresponding relationship between the water saturation parameter and the deep lateral resistivity, and the expression of Archie's formula is:

[0019] R t =(a×b×R wh )÷(S w n1 ×Φ m1 ), where: S w is water saturation; R t is the deep lateral resistivity value; a and b are the lithology coefficients obtained through lithology experiments; m1 is the cementation index; n1 is the saturation index; Φ is the total porosity of the formation; R wh is the resistivity of the mixed solution.

[0020] Preferably, m1 and n1 are calculated by the mixed solution salinity and the rock electrical parameters obtained by lithology experiments;

[0021] The resistivity of the mixed solution is calculated based on the salinity of the mixed solution and the reservoir temperature of the target layer.

[0022] Preferably, the method for determining the polymer drive movable unit further comprises the following steps:

[0023] According to the sampling data set in the actual production process, the identification plate of the acoustic time difference and resistivity cross-plot is established. The sampling data set includes a set of known sampling points of effective producing level and the measured values ​​of deep lateral resistivity parameters and acoustic time difference parameters of each sampling point;

[0024] According to the distribution of each sampling point in the identification plate, the correctness of the deep lateral resistivity parameter threshold interval obtained in step 2 is verified, and the upper and lower limits of the deep lateral resistivity parameter threshold interval are adjusted and applied according to the verification results to identify the polymer flooding movable units.

[0025] In another aspect, the present invention provides a device for determining a polymer drive movable unit, comprising:

[0026] The acquisition and analysis module obtains the phase permeability experimental data of multiple core samples, and obtains the fitting function of the water saturation parameter and the water production rate parameter based on the phase permeability experimental data. The core samples are obtained by coring in the working area.

[0027] a processing module, which calculates a water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and a preset water production rate parameter threshold interval, and then obtains a deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using the Archie formula; and

[0028] The identification module obtains the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions, and when the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, the corresponding depth position is determined as a polymer drive movable position, and the depth interval composed of the polymer drive movable positions is a polymer drive movable unit.

[0029] Preferably, the water production rate parameter threshold interval preset in the processing module includes different sub-intervals, and the corresponding water saturation parameter threshold sub-interval and deep lateral resistivity parameter threshold sub-interval are calculated respectively according to the water production rate parameter threshold sub-interval, and the effective producing level of the polymer drive movable unit is determined according to the deep lateral resistivity parameter falling into different deep lateral resistivity parameter threshold sub-intervals, wherein those identified as the same effective producing level are produced using the polymer formula of the corresponding level.

[0030] Furthermore, the effective production levels include, from low to high: the first level, the second level, the third level, the fourth level, and the fifth level, and from the first level to the fifth level, the upper and lower limits of the deep lateral resistivity parameter threshold sub-interval increase step by step, and the upper and lower limits of the water production rate parameter threshold sub-interval and the water saturation parameter threshold sub-interval decrease step by step.

[0031] Furthermore, the preset lower limits of the water production rate parameter threshold sub-intervals of the first level to the fifth level are: 90%, 80%, 60%, 40%, 10% respectively;

[0032] The lower limit of the threshold sub-interval of the water saturation parameter of each level is calculated by the lower limit of the threshold sub-interval of the water production rate parameter of the corresponding level; and

[0033] The lower limit of the threshold sub-interval of the deep lateral resistivity parameter at each level is calculated by the lower limit of the threshold sub-interval of the water saturation parameter at the corresponding level.

[0034] Preferably, the device for determining the polymer drive movable unit also includes a preset analysis module. Before obtaining the phase permeability experimental data of the sampling point through the acquisition and analysis module, it is first identified through preset judgment conditions whether each target reservoir has been used, and for the identified reservoir that has been used, the acquisition and analysis module is started.

[0035] Furthermore, the preset determination conditions include one or more of a natural potential change method, a gas composition determination method, and an oil content identification method.

[0036] Preferably, the device for determining the polymer drive movable unit further comprises a verification module, which establishes an identification plate of the acoustic time difference and resistivity cross-plot according to a sampling data set in an actual production process, wherein the sampling data set comprises a set of known sampling points of effective producing level, and the measured values ​​of the deep lateral resistivity parameters and the acoustic time difference parameters of each sampling point;

[0037] According to the distribution of each sampling point in the identification plate, the correctness of the deep lateral resistivity parameter threshold interval obtained in the step processing module is verified, and the upper and lower limits of the deep lateral resistivity parameter threshold interval are adjusted and applied according to the verification result to identify the polymer flooding movable unit.

[0038] On the other hand, the present invention provides a machine-readable storage medium having instructions stored thereon, the instructions being used to enable a machine to execute the method for determining a polymer-driven movable unit of the present application.

[0039] In another aspect, the present invention provides a processor, which is used to run a program, and the program, when run, is used to execute: the method for determining the polymer-driven movable unit of the present application.

[0040] Through the above technical scheme, the phase permeability experimental data of multiple core samples are first obtained, and the fitting function of the water saturation parameter and the water production rate parameter is obtained according to the phase permeability experimental data. Then, according to the obtained fitting function and the preset water production rate parameter threshold interval, the water saturation parameter threshold interval corresponding to the water production rate parameter threshold interval is obtained, and then the corresponding deep lateral resistivity parameter threshold interval is obtained. Finally, according to the parameter curves of the acoustic wave time difference and deep lateral resistivity at different depth positions, it is determined whether the deep lateral resistivity and acoustic wave time difference parameters at each depth position fall into the corresponding parameter threshold interval at the same time, and the polymer drive movable unit is determined. According to the method for determining the polymer drive movable unit of the present application, the movable unit that can be used for polymer drive can be found from the reservoir that has been used.

[0041] Other features and advantages of the embodiments of the present invention will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following specific embodiments, they are used to explain the embodiments of the present invention, but do not constitute a limitation on the embodiments of the present invention. In the accompanying drawings:

[0043] Figure 1 It is a schematic diagram of the redirection effect of polymer flooding in the prior art;

[0044] Figure 2 It is the microscopic distribution diagram of the fully utilized residual oil faced by polymer flooding in the prior art;

[0045] Figure 3 It is a flow chart of an embodiment of a method for determining a polymer drive movable unit of the present application;

[0046] Figure 4A , 4B They are Figure 3 The relative permeability cross-plot of water phase and oil phase obtained in step 1 of the embodiment;

[0047] Figure 5 It is comprehensive Figure 4A , 4B And the phase permeability curve obtained by the water production rate;

[0048] Figure 6 is based on Figure 5 The cross plot of water saturation and water production rate obtained from the phase permeability curve;

[0049] Figure 7 yes Figure 3 A corresponding relationship diagram between the cementation index and the mineralization degree obtained in step 2 of the embodiment;

[0050] Figure 8 yes Figure 3 The corresponding relationship diagram between the saturation index and the mineralization degree obtained in step 2 of the embodiment;

[0051] Fig. 9 is a flow chart of another embodiment of a method for determining a polymer drive movable unit of the present application;

[0052] Fig.10 according to Fig. 9 The acoustic time difference resistivity cross plot obtained in the embodiment;

[0053] Fig.11 To project the parameters of a set of different depth points of a well log onto Fig.10 Schematic diagram in the intersection diagram;

[0054] Fig.12 for Fig.11 The single well comprehensive histogram of the well logging;

[0055] Fig.13 A structural diagram of an embodiment of a device for determining a polymer-driven movable unit of the present application; and

[0056] Fig.14 This is a structural diagram of another embodiment of the device for determining the polymer-driven movable unit of the present application. DETAILED DESCRIPTION

[0057] The specific implementation of the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described here is only used to illustrate and explain the embodiment of the present invention, and is not used to limit the embodiment of the present invention.

[0058] An embodiment of the method for determining the polymer drive movable unit of the present application is applied to the oil production area that has been used, especially for the heterogeneous sandstone reservoir, which is no longer able to produce oil after water drive, and needs to adopt a targeted polymer drive method to implement tertiary oil recovery, but the present invention is not limited to this. In this embodiment, the implementation cost of polymer drive is higher than that of water drive, and its implementation steps are as follows: Figure 3 Shown include;

[0059] Step 1: Obtain phase permeability experimental data of multiple core samples, and obtain fitting functions of water saturation parameters and water production rate parameters based on the phase permeability experimental data. The core samples are obtained by coring in the working area;

[0060] Step 2: Calculate the water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and the preset water production rate parameter threshold interval, and then obtain the deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using the Archie formula;

[0061] Step 3: Acquire the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions. When the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, determine the corresponding depth position as the polymer drive movable position, and the depth interval composed of the polymer drive movable position is the polymer drive movable unit.

[0062] It should be noted that the above preset water production rate parameter threshold interval is selected according to the field conditions of oil production, such as the construction conditions of the prepared polymer flooding technology and the preparation of the corresponding intended formulation, and the corresponding water production rate parameter threshold interval is selected. Different water production rate parameter threshold intervals can be used for different field conditions.

[0063] Furthermore, multiple water production rate parameter threshold sub-intervals can be preset, and different polymer flooding construction technical conditions and corresponding polymer formulas are used for different sub-intervals. Multiple water production rate parameter threshold sub-intervals are preset as shown in Table 1 below. At the same time, different sub-intervals also indicate the effective utilization level of the corresponding seepage unit.

[0064] Table 1 Water production rate parameter threshold sub-interval settings

[0065]

[0066] In some embodiments, it is also determined whether each target reservoir has been produced by using a preset determination condition, and then the polymer drive movable unit in the produced reservoir is determined by steps 1-3. The preset determination condition includes one or more of the natural potential change method, the gas component determination method, and the oil content identification method.

[0067] Furthermore, in some embodiments, the natural potential change method, gas component determination method, and oil content identification method are used in combination to assist in the identification and determination of the polymer drive movable unit, and the identification and determination criteria are shown in Table 2 below.

[0068] Table 2 Qualitative classification evaluation criteria for effective units of polymer flooding

[0069]

[0070]

[0071] It should also be noted that the lower limit of the water saturation parameter threshold sub-interval at each level is calculated through the lower limit of the water production rate parameter threshold sub-interval at the corresponding level; the lower limit of the deep lateral resistivity parameter threshold sub-interval at each level is calculated through the lower limit of the water saturation parameter threshold sub-interval at the corresponding level.

[0072] That is to say, different oil layers, especially different oil fields and work areas, have different storage and previous oil production conditions, and the values ​​of the water production rate parameter threshold sub-interval and the deep lateral resistivity parameter threshold sub-interval should be obtained according to the actual situation. In this embodiment, the obtained water saturation parameter threshold sub-interval is shown in Table 3 below, and the further obtained deep lateral resistivity parameter threshold sub-interval is shown in Table 4 below.

[0073] Table 3 Water saturation parameter threshold sub-intervals corresponding to the 5-level effective producing units

[0074]

[0075] Table 4 The corresponding deep lateral resistivity parameter threshold sub-intervals of the 5-level effective producing units

[0076]

[0077] In some embodiments, the process of obtaining the fitting function in step 1 includes: determining the water saturation and water production rate cross-plot through phase permeability experimental data, and obtaining the formula of the fitting function based on the water saturation and water production rate cross-plot. The water production rate of each sampling point is calculated by the water relative permeability, oil relative permeability, absolute permeability, cross-sectional area, water viscosity, oil viscosity, and pressure difference of the corresponding sampling point; the water saturation of each sampling point is corrected by the formula of the fitting function; for the point where no sampling is performed, its water saturation and water production rate are obtained by the fitting function.

[0078] In some embodiments, in step 2, the process of calculating the water saturation parameter threshold interval includes: on the function curve of the fitting function, intercepting a portion of the curve within the upper and lower limits of the preset water production rate parameter threshold interval, and determining the upper and lower limits of the water saturation parameter threshold interval based on the partial curve.

[0079] Compared with the prior art, the technical advantages of the method for determining the polymer drive movable unit of the present application are:

[0080] (1) For oil layers that have been produced, especially those that have been fully produced, or even those with multiple displacement layers and strong heterogeneity, the scattered movable units can be subdivided to better explore the potential of remaining oil and achieve the goal of stabilizing production and increasing reserves;

[0081] (2) By reasonably setting the water saturation parameter threshold interval, especially the sub-interval, targeted polymer formulations and corresponding polymer flooding technologies can be used to improve the collection efficiency and oil recovery.

[0082] It should be noted that in step 1, the phase permeability experimental data is used to establish Figure 4A The cross-plot of water relative permeability and water saturation shown in Figure 4B The cross-plot of oil relative permeability and water saturation is shown in Figure 1, and the corresponding water production rate of each point in the plot is calculated according to Darcy's formula, thus obtaining Figure 5 The water production rate curve shown in Figure 4A , 4B The cross-plot of water and oil relative permeability and water saturation is as follows: Figure 5 The phase penetration intersection diagram shown in the figure is then sorted out to obtain Figure 6 The cross-plot of water production rate and water saturation is shown, and finally according to Figure 6 The fitting formula of water saturation and water production rate is obtained from the curve.

[0083] It should be noted that the storage conditions in different work areas are different, and the fitting formulas obtained are also different. It is necessary to conduct phase permeability experiments, take coring samples in the work area to obtain multiple core samples, and obtain the fitting functions of water saturation parameters and water production rate parameters based on the phase permeability experimental data of the core samples.

[0084] In this embodiment, the water production rate is calculated using the formula: in, In the formula, Q W is the water production, Q o is the oil production, K rw is the water relative permeability, K ro is the relative permeability of oil, K is the absolute permeability, A is the cross-sectional area, U w is the viscosity of water, U o is the viscosity of the oil, P is the pressure difference, and the above are the phase permeability test data of each core sample.

[0085] In this embodiment, the fitting function of the water production rate parameter and the water saturation parameter is: w =-0.055S w 2 +7.391S w -149.377, of which S w is water saturation, F w is the water production rate.

[0086] It should also be noted that the fitting function obtained by fitting according to the above formula is a reliable formula that discards abnormal data, so it can be used to correct the experimental data of each core sample. At the same time, the fitting function is also used to obtain the water saturation parameters and water production rate parameters of the location points where no core samples were taken.

[0087] In step 2, the effective unit level of polymer flooding is different, the formation mixed fluid salinity value is also different, the parameter value of Archie formula, such as the mixed fluid resistivity R wh , m1, n1 are parameters that change dynamically with the mineralization of the formation mixture. The lithology coefficients a and b are mainly affected by the lithology of the reservoir. The lithology changes slightly after water injection development and are taken according to the actual value of the block.

[0088] The modified formula of Archie's formula is:

[0089] R t =(a×b×R wh )÷(S w n1 ×Φ m1 ),

[0090] Where: S w is water saturation; R tis the deep lateral resistivity value; a and b are the lithology coefficients obtained through lithology experiments; m1 is the cementation index; n1 is the saturation index; Φ is the total porosity of the formation; R wh is the resistivity of the mixed solution.

[0091] In the above formula, m1 and n1 are calculated by the mixed solution salinity and the rock electrical parameters obtained through lithology experiments;

[0092] The resistivity of the mixed solution is calculated based on the salinity of the mixed solution and the reservoir temperature of the target layer.

[0093] Specifically, after the oil layer is produced, the injected water is mixed with the original formation water. Due to the difference in salinity, ion diffusion occurs between the solutions, and the formation water resistivity of the original Archie formula is transformed into the resistivity of the mixed liquid, which is a parameter that changes with the salinity of the mixed liquid and the temperature of the target layer reservoir. The salinity of the mixed liquid and the temperature of the target layer reservoir are taken according to the actual production data.

[0094]

[0095] in:

[0096] P z is the mineralization of the mixed solution, mg / L;

[0097] T is the target reservoir temperature, °C;

[0098] R wh is the resistivity of the mixed solution, Ω.m.

[0099] The cementation index and saturation index are parameters that reflect the tortuosity of current conduction in rock seepage channels, which are mainly determined by the rock's own properties such as pore geometry, microscopic pore structure characteristics and cementation conditions. Therefore, the cementation index and saturation index in the Archie formula rock electrical parameters are closely related to the volume fraction of clay mineralization in the reservoir. As a large amount of clay minerals are washed away, the pore structure of the rock changes. The cementation index m of the original Archie formula is transformed into m1, which changes dynamically with the mineralization of the mixed solution, and the saturation index n of the original Archie formula is transformed into n1, which changes dynamically with the mineralization of the mixed solution.

[0100] In order to make an accurate evaluation, a large number of core samples were taken from the study area to conduct experimental measurements of rock electrical parameters, and the measured data of different mineralizations were statistically regressed to finally obtain the relationship between rock electrical parameters and mineralization.

[0101] m1=5.9×10 -22 P z 5 -3.3×10 -17 P z 4 +7.2×10 -13 Pz 3 -8.0×10 -9 P z 2 +0.0005P z +0.7514,

[0102] Where: m1 is the cementation index that changes dynamically with the mineralization of the mixed solution;

[0103] P z is the mineralization of the mixed solution.

[0104] n1=-2.2×10 -22 P z 5 +1.6×10 -17 P z 4 -4.5×10 -13 P z 3 +6.3×10 -9 P z 2 -0.0004P z +2.5881,

[0105] Where: n1 is the saturation index that changes dynamically with the mineralization of the mixed solution;

[0106] P z is the mineralization of the mixed solution.

[0107] The corresponding relationship between the cementation index and the mineralization degree obtained in this example is as follows: Figure 7 As shown in Figure 2, the corresponding relationship between saturation index and mineralization is as follows: Figure 8 shown.

[0108] For the possible existence of oil layers whose production is unknown, and for further accurately obtaining the corresponding deep lateral resistivity parameter threshold interval as the water saturation parameter changes during oil production, the process of another preferred embodiment of the method for determining the polymer flooding movable unit of the present application is as follows: Fig. 9 As shown, including:

[0109] Pre-identification step S0: First determine whether the corresponding oil layer has been used. If the target reservoir is identified as a reservoir that has been used through the preset judgment conditions, continue to identify the polymer drive movable unit according to the following steps S1-S5. If it has not been used, it means that the reservoir can be used for oil production in other ways, and a more economical way should be adopted according to the actual situation.

[0110] Step S1: acquiring phase permeability experimental data of multiple core samples, and obtaining fitting functions of water saturation parameters and water production rate parameters according to the phase permeability experimental data, wherein the core samples are obtained by coring in the working area;

[0111] Step S2: Calculate the water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and the preset water production rate parameter threshold interval, and then obtain the deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using the Archie formula;

[0112] Step S3: Acquire the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions, and when the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval, and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, determine the corresponding depth position as the polymer drive movable position, and the depth interval composed of the polymer drive movable position is the polymer drive movable unit;

[0113] Step S4: establishing an identification plate of an acoustic time difference and resistivity cross-plot according to a sampling data set in an actual production process, wherein the sampling data set includes a set of known sampling points of effective producing level and measured values ​​of deep lateral resistivity parameters and acoustic time difference parameters of each sampling point;

[0114] Step S5: According to the distribution of each sampling point in the identification plate, verify the correctness of the deep lateral resistivity parameter threshold interval obtained, and adjust and apply the upper and lower limits of the deep lateral resistivity parameter threshold interval according to the verification result to identify the polymer flooding movable unit.

[0115] The above steps S1-S3 are the same as those in the previous embodiment and will not be described again. It should be noted that steps S4-S5 are used to further accurately obtain the corresponding deep lateral resistivity parameter threshold interval as the water saturation parameter changes, and to judge the level of the polymer flooding movable unit, and the polymer flooding construction technical conditions and corresponding polymer formula to be adopted according to the more accurate deep lateral resistivity parameter threshold interval.

[0116] In this embodiment, based on the sampled data set in the actual production process, the acoustic time difference resistivity cross-plot is obtained as follows: Fig.10 As shown, the sampling data obtained under the given lithology parameters and given cementation index and saturation index conditions of this embodiment are marked, including the values ​​of the acoustic wave time difference parameters and the deep resistivity parameters.

[0117] At the same time, for a given well logging curve (number J16-23-42C), if Fig.12 Please focus on the acoustic time difference parameter and deep resistivity parameter curve, and mark the corresponding acoustic time difference and deep resistivity values ​​at the positions marked with dots. Fig.10In the intersection diagram shown, we get Fig.11 ,from Fig.12 It can be seen that compared with the original interpretation conclusion obtained in the prior art, the current interpretation conclusion obtained according to the application technical solution can make a fine division of the original interpretation conclusion. Fig.12 According to the original interpretation conclusion, the area 65 was identified as a first-level movable unit. In the current interpretation conclusion, the area is divided into 65a and 65b and identified as a second-level movable unit and a first-level movable unit respectively.

[0118] Compared with the prior art, the technical advantage of this embodiment is that for oil layers with multiple displacement layers and strong heterogeneity, the scope and level of the movable units can be identified more accurately.

[0119] The present application also provides a device for determining a polymer-driven movable unit, the composition structure of which is as follows: Fig.13 As shown, including:

[0120] The acquisition and analysis module obtains the phase permeability experimental data of multiple core samples, and obtains the fitting function of the water saturation parameter and the water production rate parameter based on the phase permeability experimental data. The core samples are obtained by coring in the working area.

[0121] a processing module, which calculates a water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and a preset water production rate parameter threshold interval, and then obtains a deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using the Archie formula; and

[0122] The identification module obtains the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions, and when the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, the corresponding depth position is determined as a polymer drive movable position, and the depth interval composed of the polymer drive movable positions is a polymer drive movable unit.

[0123] In some embodiments, the water production rate parameter threshold interval preset in the processing module includes different sub-intervals, and the corresponding water saturation parameter threshold sub-interval and deep lateral resistivity parameter threshold sub-interval are calculated respectively according to the water production rate parameter threshold sub-interval, and the effective production level of the polymer drive movable unit is determined according to the deep lateral resistivity parameter falling into different deep lateral resistivity parameter threshold sub-intervals, wherein those identified as the same effective production level are produced using the polymer formula of the corresponding level.

[0124] In some embodiments, the effective utilization levels include, from low to high: the first level, the second level, the third level, the fourth level, and the fifth level, and from the first level to the fifth level, the upper and lower limits of the deep lateral resistivity parameter threshold sub-interval gradually increase, and the upper and lower limits of the water production rate parameter threshold sub-interval and the water saturation parameter threshold sub-interval gradually decrease.

[0125] Furthermore, the preset lower limits of the water production rate parameter threshold sub-intervals of the first level to the fifth level are: 90%, 80%, 60%, 40%, 10% respectively;

[0126] The lower limit of the threshold sub-interval of the water saturation parameter of each level is calculated by the lower limit of the threshold sub-interval of the water production rate parameter of the corresponding level; and

[0127] The lower limit of the threshold sub-interval of the deep lateral resistivity parameter at each level is calculated by the lower limit of the threshold sub-interval of the water saturation parameter at the corresponding level.

[0128] In some embodiments, the polymer drive movable unit determination device further includes a preset analysis module and a verification module, and its composition structure diagram is shown in FIG. Fig.14 shown.

[0129] Among them, the preset analysis module, before obtaining the phase permeability experimental data of the sampling point through the acquisition and analysis module, first identifies whether each target reservoir has been produced through the preset judgment conditions, and starts the acquisition and analysis module for the identified reservoir that has been produced. The preset judgment conditions include one or more of the natural potential change method, the gas component determination method, and the oil content identification method.

[0130] The verification module establishes an identification plate of the acoustic time difference and resistivity cross-plot based on the sampling data set in the actual production process. The sampling data set includes a set of known sampling points of effective producing level and the measured values ​​of deep lateral resistivity parameters and acoustic time difference parameters of each sampling point;

[0131] According to the distribution of each sampling point in the identification plate, the correctness of the deep lateral resistivity parameter threshold interval obtained in the step processing module is verified, and the upper and lower limits of the deep lateral resistivity parameter threshold interval are adjusted and applied according to the verification result to identify the polymer flooding movable unit.

[0132] An embodiment of the present invention provides a machine-readable storage medium, on which instructions are stored, and the instructions are used to enable a machine to execute the method for determining a polymer-driven movable unit of the present application.

[0133] An embodiment of the present invention provides a processor, which is used to run a program. When the program is run, it is used to execute: the method for determining the polymer-driven movable unit of the present application.

[0134] The embodiment of the present invention provides a device, which includes a processor, a memory, and a program stored in the memory and can be run on the processor, and the steps of the method for determining the polymer drive movable unit of the present application are implemented when the processor executes the program. The device in this article can be a server, a PC, a PAD, a mobile phone, etc.

[0135] The present application also provides a computer program product, which, when executed on a data processing device, is suitable for executing the program of initializing the steps of the method for determining a polymer drive movable unit of the present application.

[0136] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0137] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0138] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0139] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1The steps for the functions specified in one or more boxes.

[0140] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0141] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0142] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0143] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0144] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A method for determining a polymer drive movable unit, comprising: Step 1: Obtain phase permeability experimental data of multiple core samples, and obtain fitting functions of water saturation parameters and water production rate parameters according to the phase permeability experimental data, wherein the core samples are obtained by coring in the working area; Step 2: Calculate the water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and the preset water production rate parameter threshold interval, and then obtain the deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using the Archie formula; Step 3: Acquire the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions. When the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, determine the corresponding depth position as a polymer drive movable position, and the depth interval composed of the polymer drive movable position is a polymer drive movable unit.

2. The method according to claim 1, characterized in that The preset water production rate parameter threshold interval includes different sub-intervals, and the corresponding water saturation parameter threshold sub-interval and deep lateral resistivity parameter threshold sub-interval are calculated respectively according to the water production rate parameter threshold sub-interval, and the effective producing level of the polymer drive movable unit is determined according to the deep lateral resistivity parameter falling into different deep lateral resistivity parameter threshold sub-intervals, wherein those identified as the same effective producing level are produced using the polymer formula of the corresponding level.

3. The method according to claim 2, characterized in that The effective production levels include, from low to high: the first level, the second level, the third level, the fourth level, and the fifth level, and from the first level to the fifth level, the upper and lower limits of the deep lateral resistivity parameter threshold sub-interval increase step by step, and the upper and lower limits of the water production rate parameter threshold sub-interval and the water saturation parameter threshold sub-interval decrease step by step.

4. The method according to claim 3, characterized in that The preset lower limits of the water production rate parameter threshold sub-intervals from the first level to the fifth level are: 90%, 80%, 60%, 40%, 10% respectively; The lower limit of the threshold sub-interval of the water saturation parameter of each level is calculated by the lower limit of the threshold sub-interval of the water production rate parameter of the corresponding level; and The lower limit of the threshold sub-interval of the deep lateral resistivity parameter at each level is calculated by the lower limit of the threshold sub-interval of the water saturation parameter at the corresponding level.

5. The method according to claim 1, characterized in that Before performing step 1, it is firstly identified whether each target reservoir has been produced by using preset judgment conditions, and then the polymer flooding movable units in the produced reservoir are determined through steps 1-3.

6. The method according to claim 5, characterized in that The preset determination conditions include one or more of a natural potential change method, a gas component determination method, and an oil content identification method.

7. The method according to claim 1, characterized in that The process of obtaining the fitting function in step 1 includes: The water saturation and water production rate cross-plot is determined by the phase permeability experimental data, and the formula of the fitting function is obtained based on the water saturation and water production rate cross-plot. Among them, the water production rate of each sampling point is calculated by the water relative permeability, oil relative permeability, absolute permeability, cross-sectional area, water viscosity, oil viscosity, and pressure difference of the corresponding sampling point; The water saturation of each sampling point is corrected by the formula of the fitting function; and For the points where no sampling is performed, their water saturation and water production rate are calculated by the fitting function.

8. The method according to claim 1, characterized in that In step 2, the process of calculating the water saturation parameter threshold interval includes: on the function curve of the fitting function, intercepting a portion of the curve within the upper and lower limits of the preset water production rate parameter threshold interval, and determining the upper and lower limits of the water saturation parameter threshold interval based on the portion of the curve.

9. The method according to claim 1, characterized in that: In step 2, the Archie formula is used to determine the corresponding relationship between the water saturation parameter and the deep lateral resistivity. The expression of the Archie formula is: R t =(a×b×R wh )÷(S w n1 ×Φ m1 ), where: S w is water saturation; R t is the deep lateral resistivity value; a and b are lithology coefficients obtained through lithology experiments; m1 is the cementation index; n1 is the saturation index; Φ is the total porosity of the formation; R wh is the resistivity of the mixed solution.

10. The method according to claim 9, characterized in that The m1 and n1 are both calculated through the mixed solution salinity and the rock electrical parameters obtained through lithology experiments; The resistivity of the mixed solution is calculated based on the salinity of the mixed solution and the target layer reservoir temperature.

11. The method according to any one of claims 1 to 10, characterized in that Further including: According to the sampling data set in the actual production process, an identification plate of the acoustic time difference and resistivity cross-plot is established, wherein the sampling data set includes a group of known sampling points of effective producing level and the measured values ​​of the deep lateral resistivity parameters and the acoustic time difference parameters of each sampling point; According to the distribution of each sampling point in the identification plate, the correctness of the deep lateral resistivity parameter threshold interval obtained in step 2 is verified, and the upper and lower limits of the deep lateral resistivity parameter threshold interval are adjusted and applied according to the verification result to identify the polymer flooding movable unit.

12. A device for determining a polymer drive movable unit, comprising: An acquisition and analysis module is used to obtain phase permeability experimental data of the plurality of core samples, and obtain fitting functions of water saturation parameters and water production rate parameters according to the phase permeability experimental data, wherein the core samples are obtained by coring in the working area; A processing module, which calculates a water saturation parameter threshold interval corresponding to the water production rate parameter threshold according to the fitting function and a preset water production rate parameter threshold interval, and then obtains a deep lateral resistivity parameter threshold interval corresponding to the water saturation parameter threshold interval by using an Archie formula; and The identification module obtains the acoustic wave time difference parameter curve and the deep lateral resistivity parameter curve at different depth positions, and when the deep lateral resistivity parameter falls into the deep lateral resistivity parameter threshold interval and the acoustic wave time difference parameter falls into the preset acoustic wave time difference parameter threshold interval, the corresponding depth position is determined as a polymer drive movable position, and the depth interval composed of the polymer drive movable position is a polymer drive movable unit.

13. The device according to claim 12, characterized in that The preset water production rate parameter threshold interval in the processing module includes different sub-intervals, and the corresponding water saturation parameter threshold sub-interval and deep lateral resistivity parameter threshold sub-interval are calculated respectively according to the water production rate parameter threshold sub-interval, and the effective producing level of the polymer drive movable unit is determined according to the deep lateral resistivity parameter falling into different deep lateral resistivity parameter threshold sub-intervals, wherein those identified as the same effective producing level are produced using the polymer formula of the corresponding level.

14. The device according to claim 13, characterized in that The effective production levels include, from low to high: the first level, the second level, the third level, the fourth level, and the fifth level, and from the first level to the fifth level, the upper and lower limits of the deep lateral resistivity parameter threshold sub-interval increase step by step, and the upper and lower limits of the water production rate parameter threshold sub-interval and the water saturation parameter threshold sub-interval decrease step by step.

15. The device according to claim 14, characterized in that The preset lower limits of the water production rate parameter threshold sub-intervals from the first level to the fifth level are: 90%, 80%, 60%, 40%, 10% respectively; The lower limit of the threshold sub-interval of the water saturation parameter of each level is calculated by the lower limit of the threshold sub-interval of the water production rate parameter of the corresponding level; and The lower limit of the threshold sub-interval of the deep lateral resistivity parameter at each level is calculated by the lower limit of the threshold sub-interval of the water saturation parameter at the corresponding level.

16. The device according to claim 12, characterized in that It also includes a preset analysis module. Before obtaining the phase permeability experimental data of the sampling point through the acquisition and analysis module, it first identifies whether each target reservoir has been used through preset judgment conditions, and starts the acquisition and analysis module for the identified reservoir that has been used.

17. The device according to claim 16, characterized in that The preset determination conditions include one or more of a natural potential change method, a gas component determination method, and an oil content identification method.

18. The device according to claims 12-17, characterized in that It also includes a verification module, which establishes an identification plate of the acoustic time difference and resistivity cross-plot according to a sampling data set in an actual production process, wherein the sampling data set includes a set of known sampling points of effective producing level and measured values ​​of deep lateral resistivity parameters and acoustic time difference parameters of each sampling point; According to the distribution of each sampling point in the identification plate, the correctness of the deep lateral resistivity parameter threshold interval obtained in the processing module is verified, and the upper and lower limits of the deep lateral resistivity parameter threshold interval are adjusted and applied according to the verification result to identify the polymer flooding movable unit.

19. A machine-readable storage medium having instructions stored thereon, the instructions being used to cause a machine to execute: the method for determining a polymer-driven movable unit as described in any one of claims 1 to 11.

20. A processor, characterized in that: Used to run a program, wherein the program, when run, is used to execute: a method for determining a polymer-driven movable unit as described in any one of claims 1-11.

Images