Microscopic pore throat characteristic data determination method, device, equipment and medium

Abstract

The invention relates to the field of oil and gas reservoir evaluation, and discloses a microscopic pore throat characteristic data determination method, device and equipment and a medium. A group of specific standard pore throat radiuses can be set, and experimental data obtained in high-pressure mercury injection experiments in different laboratories are converted to be under the group of specific standard pore throat radiuses, so that unification of the experimental data of the different laboratories is realized, and multiple standard experimental data are obtained; the microcosmic pore throat characteristic data is determined according to multiple standard experiment data, the difference of the microcosmic pore throat characteristic data generated under the experiment operation difference of different laboratories is effectively reduced, the comparability of the microcosmic pore throat characteristic data is improved, and errors caused by experiment operation of different laboratories are reduced.

Description

Technical Field

[0001] The invention relates to the field of oil and gas reservoir evaluation, and in particular to a method, device, equipment and medium for determining microscopic pore throat characteristic data. Background Art

[0002] In recent years, oil and gas exploration and development, both onshore and offshore, have gradually moved towards medium and deep layers, which are the direction of future oil and gas development. However, as the depth increases, the oil and gas reservoirs inevitably deteriorate, so medium and deep oil and gas reservoirs are mostly low-permeability or even tight oil and gas reservoirs.

[0003] The main characteristics of low-permeability and tight oil and gas reservoirs include poor reservoir properties (low porosity and permeability) and complex pore-permeability relationships. At present, laboratories can conduct high-pressure mercury injection experiments on reservoirs to obtain corresponding experimental data, determine the corresponding microscopic pore throat characteristic data based on the experimental data, and evaluate and develop low-permeability and tight oil and gas reservoirs based on the microscopic pore throat characteristic data of the reservoir.

[0004] However, there are differences in the operations performed by different laboratories when conducting high-pressure mercury injection experiments, such as differences in the values ​​of experimental parameters, which may lead to the incomparability of the microscopic pore throat characteristic data calculated by different laboratories for the same reservoir. For example, the relationship between two microscopic pore throat characteristic parameters calculated by different laboratories for the same reservoir may be completely opposite. Summary of the invention

[0005] The present invention provides a method, device, equipment and medium for determining microscopic pore throat characteristic data, which are used to solve the defect that the microscopic pore throat characteristic data calculated by different laboratories for the same reservoir are not comparable due to the differences in the operations performed by different laboratories in the related technology when performing high-pressure mercury injection experiments, and the difference in microscopic pore throat characteristic data generated by the differences in experimental operations of different laboratories is effectively reduced, the comparability of microscopic pore throat characteristic data is increased, and the errors caused by experimental operations of different laboratories are reduced.

[0006] In a first aspect, the present invention provides a method for determining microscopic pore throat characteristic data, comprising: Acquire a plurality of experimental data obtained by performing high-pressure mercury injection experiments on rock samples of oil and gas reservoirs, each of the experimental data including a corresponding pore throat radius and a cumulative mercury injection saturation; and acquire a plurality of standard pore throat radii; According to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation, the plurality of experimental data are converted into standard experimental data under each of the standard pore throat radii; Based on the standard experimental data at each standard pore throat radius, the microscopic pore throat characteristic data of the rock sample is determined.

[0007] Optionally, converting the plurality of experimental data into standard experimental data under each of the standard pore throat radii according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation comprises: Determine the standard cumulative mercury injection saturation corresponding to each standard pore throat radius according to the multiple standard pore throat radii, the pore throat radii in the multiple experimental data, and the cumulative mercury injection saturation; For any of the standard pore throat radius, the standard pore throat radius and the corresponding standard cumulative mercury injection saturation are taken as a whole as standard experimental data under the standard pore throat radius; The method of determining the microscopic pore throat characteristic data of the rock sample based on the standard experimental data under each standard pore throat radius includes: Based on the standard pore throat radius and the standard cumulative mercury saturation in each of the standard experimental data, the microscopic pore throat characteristic data of the rock sample are determined.

[0008] Optionally, determining the standard cumulative mercury injection saturation corresponding to each standard pore throat radius according to the multiple standard pore throat radii, the pore throat radii in the multiple experimental data, and the cumulative mercury injection saturation includes: Determining a minimum first pore throat radius among the pore throat radii of the plurality of experimental data; and determining a second pore throat radius corresponding to the displacement pressure among the pore throat radii of the plurality of experimental data; Among the plurality of standard pore throat radii, select a plurality of first standard pore throat radii that are not less than the first pore throat radius and not greater than the second pore throat radius; Using each of the first standard pore throat radii as an interpolated pore throat radius, interpolating the pore throat radii and cumulative mercury injection saturations of the plurality of experimental data based on each of the interpolated pore throat radii to obtain the interpolated cumulative mercury injection saturation corresponding to each of the first standard pore throat radii; According to each of the first standard pore throat radii and the corresponding interpolated cumulative mercury injection saturation, the standard cumulative mercury injection saturation corresponding to each of the standard pore throat radii is determined.

[0009] Optionally, determining a second pore throat radius corresponding to the displacement pressure from the pore throat radii of the plurality of experimental data includes: Sorting the plurality of experimental data in ascending order of the cumulative mercury saturation in the experimental data to obtain a corresponding sorting array; Traversing each of the experimental data in the sorted array in order from small to large, for any of the traversed experimental data, determining a second-order variation corresponding to the cumulative mercury intrusion saturation in the experimental data, and judging whether the second-order variation is greater than 1; If it is determined that the second-order variation is less than 1, continue to traverse the next experimental data; If it is determined that the second-order variation is not less than 1, the traversal is stopped, and the pore throat radius in the previous experimental data of the experimental data is determined as the second pore throat radius.

[0010] Optionally, determining the second-order variation corresponding to the cumulative mercury intrusion saturation in the experimental data includes: The second-order variation corresponding to the cumulative mercury saturation in the experimental data is calculated based on a preset formula; wherein the preset formula is: ; ; in, , and the order of arrangement of the experimental arrays in the sorted array; The order of arrangement is The second-order variation corresponding to the cumulative mercury saturation in the experimental data; The order of arrangement is A change in the cumulative mercury saturation in the experimental data compared to the cumulative mercury saturation in the previous experimental data; The order of arrangement is The cumulative mercury saturation in the experimental data is for , , or ,when hour, .

[0011] Optionally, determining the standard cumulative mercury injection saturation corresponding to each of the first standard pore throat radii and the corresponding interpolated cumulative mercury injection saturation comprises: For any of the first standard pore throat radius, the interpolated cumulative mercury injection saturation corresponding to the first standard pore throat radius is directly used as the standard cumulative mercury injection saturation corresponding to the first standard pore throat radius; Determining a maximum radius and a minimum radius among the plurality of first standard pore throat radii, and determining the interpolated cumulative mercury injection saturation corresponding to the minimum radius; For any second standard pore throat radius, if it is determined that the second standard pore throat radius is greater than the maximum radius, 0 is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius; if it is determined that the second standard pore throat radius is less than the minimum radius, the interpolated cumulative mercury injection saturation corresponding to the minimum radius is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius.

[0012] Optionally, the microscopic pore throat characteristic data includes at least one of a permeability contribution value, a pore throat radius mean, an average pore throat radius and a sorting coefficient.

[0013] In a second aspect, the present invention provides a device for determining microscopic pore throat characteristic data, comprising: A first acquisition unit is used to acquire a plurality of experimental data obtained by performing a high-pressure mercury injection experiment on a rock sample of an oil and gas reservoir, each of the experimental data including a corresponding pore throat radius and a cumulative mercury injection saturation; The second acquisition unit is used to acquire a plurality of standard pore throat radii; a conversion unit, configured to convert the plurality of experimental data into standard experimental data under each of the standard pore throat radii according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation; A determination unit is used to determine the microscopic pore throat characteristic data of the rock sample based on the standard experimental data under each standard pore throat radius.

[0014] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for determining microscopic pore throat characteristic data of the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0015] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for determining microscopic pore throat characteristic data of the first aspect or any corresponding embodiment thereof.

[0016] The method, device, equipment and medium for determining microscopic pore throat characteristic data provided by the present invention can obtain multiple experimental data obtained by performing high-pressure mercury injection experiments on rock samples of oil and gas reservoirs, and obtain multiple standard pore throat radii. Wherein, each experimental data includes a corresponding pore throat radius and a cumulative mercury injection saturation. According to multiple standard pore throat radii, the pore throat radii and the cumulative mercury injection saturation in multiple experimental data, the multiple experimental data are converted into standard experimental data under each standard pore throat radius. Based on the standard experimental data under each standard pore throat radius, the microscopic pore throat characteristic data of the rock sample is determined. The present invention can set a group of specific standard pore throat radii, convert the experimental data obtained by different laboratories in high-pressure mercury injection experiments to the specific group of standard pore throat radii, thereby realizing the unification of experimental data of different laboratories, obtaining multiple standard experimental data, and determining microscopic pore throat characteristic data according to multiple standard experimental data, effectively reducing the difference in microscopic pore throat characteristic data generated by experimental operation differences of different laboratories, increasing the comparability of microscopic pore throat characteristic data, and reducing the error caused by experimental operation of different laboratories. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0018] Figure 1 A flow chart of a method for determining microscopic pore throat characteristic data provided by an embodiment of the present invention; Figure 2 A schematic diagram of the relationship between the radius mean and the permeability calculated by a related technology provided in an embodiment of the present invention; Figure 3 A schematic diagram of the relationship between the sorting coefficient and the permeability calculated by a related technology provided in an embodiment of the present invention; Figure 4 A schematic table of experimental data provided by an embodiment of the present invention; Figure 5 A schematic table of second-order variation corresponding to the cumulative mercury intrusion saturation provided in an embodiment of the present invention; Figure 6 A schematic diagram of comparison between interpolated experimental data and original experimental data provided in an embodiment of the present invention; Figure 7 A schematic diagram of standard experimental data provided by an embodiment of the present invention; Figure 8A schematic diagram of the relationship between the radius mean and the permeability calculated based on standard experimental data provided in an embodiment of the present invention; Fig. 9 A schematic diagram of the relationship between a sorting coefficient and a permeability calculated based on standard experimental data provided in an embodiment of the present invention; Fig.10 A schematic structural diagram of a device for determining microscopic pore throat characteristic data provided by an embodiment of the present invention; Fig.11 A schematic diagram of the structure of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0019] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0020] Combine the following Figure 1-Figure 9 The method for determining microscopic pore throat characteristic data of the present invention is described.

[0021] like Figure 1 As shown, this embodiment proposes a first method for determining microscopic pore throat characteristic data, which may include the following steps: S101, obtaining a plurality of experimental data obtained by performing high-pressure mercury injection experiments on rock samples of oil and gas reservoirs, wherein each experimental data includes a corresponding pore throat radius and cumulative mercury injection saturation.

[0022] The oil and gas reservoir may be a reservoir at a certain layer in a certain oil field.

[0023] Specifically, this embodiment can perform high-pressure mercury injection experiments on rock samples and record experimental data. The above multiple experimental data can be recorded at different times during the experiment or under different experimental parameters (such as mercury injection pressure and pore throat radius).

[0024] It should be noted that the basic principle of the high-pressure mercury injection experiment is to make mercury enter the pores of porous materials by applying high pressure. The inventors of the present invention have found that under different mercury injection pressures, the pore sizes of the pores that the mercury fluid can enter are different. The greater the pressure, the smaller the pore size of the pore that the mercury can enter. There is a linear correlation between the mercury injection pressure and the pore throat radius, specifically an inverse proportional relationship of r=0.735 / Pc, where r is the pore throat radius of the pore and Pc is the experimental mercury injection pressure.

[0025] Optionally, each of the above experimental data may also include mercury injection pressure. In this case, each of the experimental data may include the corresponding mercury injection pressure, pore throat radius and cumulative mercury injection saturation.

[0026] S102. Obtain multiple standard pore throat radii.

[0027] Specifically, the standard pore throat radius may be a specific pore throat radius set by a technician according to actual conditions. The above-mentioned multiple standard pore throat radii may be a group of specific pore throat radii set by a technician according to actual conditions.

[0028] It should be noted that in the relevant technology, different laboratories generally have different operations when conducting high-pressure mercury injection experiments. For example, the difference in the value of the pore throat radius may lead to the incomparability of the microscopic pore throat characteristic data calculated by different laboratories for the same oil field and the same layer of oil and gas reservoirs. Figure 2 As shown in the related art, the relationship between the radius mean and the permeability calculated by Laboratory 1 and Laboratory 2 for the same oil field and the same layer of oil and gas reservoirs is completely opposite, as shown in Figure 3 As shown, in the related art, there are obvious differences in the relationship between the sorting coefficient and the permeability calculated by Laboratory 1 and Laboratory 2 for the oil and gas reservoirs in the same oil field and the same layer, and they are not comparable.

[0029] This embodiment can set a set of specific standard pore throat radii, and convert the experimental data obtained by different laboratories in high-pressure mercury injection experiments to the set of specific standard pore throat radii, so as to unify the experimental data of different laboratories. Afterwards, this embodiment can use the unified experimental data to determine the microscopic pore throat characteristic data, effectively reduce the difference in microscopic pore throat characteristic data caused by the experimental operation differences of different laboratories, increase the comparability of microscopic pore throat characteristic data, and reduce the errors caused by experimental operations of different laboratories.

[0030] Optionally, the plurality of standard pore throat radii may include 630.0000, 400.0000, 250.0000, 160.0000, 100.0000, 63.0000, 40.0000, 25.0000, 16.0000, 10.0000, 6.3000, 4.0000, 2.5000, 1.6000, 1.0000, 0.6300, 0.4000, 0.2500, 0.1600, 0.1000, 0.0630, 0.0400, 0.0250, 0.0160, 0.0100, 0.0063, 0.0040, 0.0025, 0.0016, and 0.0010. This embodiment may take this group of standard pore throat radii as an example to illustrate the subsequent execution process.

[0031] It can be understood that, in this embodiment, based on the inverse proportional relationship between the pore throat radius and the mercury injection pressure (mercury injection pressure=0.735 / pore throat radius), the standard mercury injection pressure corresponding to each standard pore throat radius can be calculated.

[0032] S103, converting the plurality of experimental data into standard experimental data under each standard pore throat radius according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation.

[0033] Specifically, this embodiment can convert the plurality of experimental data into standard experimental data under the standard pore throat radius according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation.

[0034] S104. Determine microscopic pore throat characteristic data of the rock sample based on standard experimental data at each standard pore throat radius.

[0035] Specifically, this embodiment can calculate the microscopic pore throat characteristic data of the rock sample according to the standard experimental data under each standard pore throat radius.

[0036] Optionally, in the second method for determining microscopic pore throat characteristic data proposed in this embodiment, the above step S103 may include: According to a plurality of standard pore throat radii, the pore throat radii in a plurality of experimental data and the cumulative mercury injection saturation, the standard cumulative mercury injection saturation corresponding to each standard pore throat radius is determined; For any standard pore throat radius, the standard pore throat radius and the corresponding standard cumulative mercury saturation are taken as the standard experimental data under the standard pore throat radius. At this time, the above step S104 may include: The microscopic pore throat characteristic data of the rock sample are determined based on the standard pore throat radius and standard cumulative mercury saturation in each standard experimental data.

[0037] Optionally, each standard experimental data may also include a standard mercury injection pressure corresponding to a standard pore throat radius. In this case, each standard experimental data may include a corresponding standard mercury injection pressure, a standard pore throat radius, and a standard cumulative mercury injection saturation.

[0038] The method for determining microscopic pore throat characteristic data proposed in this embodiment can obtain multiple experimental data obtained by performing high-pressure mercury injection experiments on rock samples of oil and gas reservoirs, and obtain multiple standard pore throat radii. Among them, each experimental data includes a corresponding pore throat radius and cumulative mercury injection saturation. According to multiple standard pore throat radii, the pore throat radii in multiple experimental data and the cumulative mercury injection saturation, the multiple experimental data are converted into standard experimental data under each standard pore throat radius. Based on the standard experimental data under each standard pore throat radius, the microscopic pore throat characteristic data of the rock sample is determined. In this embodiment, a group of specific standard pore throat radii can be set, and the experimental data obtained by different laboratories in high-pressure mercury injection experiments can be converted to the specific standard pore throat radius, so as to unify the experimental data of different laboratories, obtain multiple standard experimental data, and determine the microscopic pore throat characteristic data according to the multiple standard experimental data, effectively reducing the difference in microscopic pore throat characteristic data generated by the experimental operation differences of different laboratories, increasing the comparability of microscopic pore throat characteristic data, and reducing the errors caused by experimental operations in different laboratories.

[0039] based on Figure 1 This embodiment proposes a second method for determining microscopic pore throat characteristic data. In this method, step S103 may include: Determine a minimum first pore throat radius from among the pore throat radii of the multiple experimental data; and determine a second pore throat radius corresponding to the displacement pressure from among the pore throat radii of the multiple experimental data; Among the multiple standard pore throat radii, select multiple first standard pore throat radii that are not less than the first pore throat radius and not greater than the second pore throat radius; Taking each first standard pore throat radius as an interpolated pore throat radius, interpolating the pore throat radii and cumulative mercury injection saturations of multiple experimental data based on each interpolated pore throat radius to obtain the interpolated cumulative mercury injection saturation corresponding to each first standard pore throat radius; According to each first standard pore throat radius and the corresponding interpolated cumulative mercury injection saturation, the standard cumulative mercury injection saturation corresponding to each standard pore throat radius is determined.

[0040] The first pore throat radius is the minimum pore throat radius among the pore throat radii of the above-mentioned multiple experimental data.

[0041] The second pore throat radius is the pore throat radius corresponding to the displacement pressure of this high-pressure mercury injection experiment among the pore throat radii of the above-mentioned multiple experimental data.

[0042] Optionally, determining the second pore throat radius corresponding to the displacement pressure from the pore throat radii of the plurality of experimental data includes: According to the order of the cumulative mercury saturation in the experimental data from small to large, a plurality of experimental data are sorted to obtain a corresponding sorting array; According to the order of arrangement from small to large, each experimental data in the sorted array is traversed in turn, and for any experimental data traversed, the second-order variation corresponding to the cumulative mercury saturation in the experimental data is determined, and it is judged whether the second-order variation is greater than 1; If it is determined that the second-order variation is less than 1, continue to traverse the next experimental data; If it is determined that the second-order variation is not less than 1, the traversal is stopped, and the pore throat radius in the previous experimental data of the experimental data is determined as the second pore throat radius.

[0043] In order to better illustrate the processing of the test process data, this embodiment proposes and combines Figure 4 The original mercury injection experimental data are shown in the figure. Figure 4 Each row of data from number 1 to 30 is an experimental data. At this time, the experimental data includes the corresponding mercury injection pressure, pore throat radius and cumulative mercury injection saturation. The pore throat radius can be calculated based on the mercury injection pressure. Figure 4 P, R and Shg represent the mercury injection pressure, pore throat radius and cumulative mercury injection saturation, respectively. Figure 4 The experimental data with sequence numbers 1 to 30 can be regarded as a sorted array obtained by sorting all the experimental data in the order of cumulative mercury injection saturation from small to large in this embodiment. Figure 4 In the experimental data, the mercury injection pressure is inversely proportional to the pore throat radius, the mercury injection pressure is positively correlated with the mercury injection saturation, and the pore throat radius is negatively correlated with the cumulative mercury injection saturation.

[0044] Optionally, the second-order variation corresponding to the cumulative mercury saturation in the experimental data determined above includes: The second-order variation corresponding to the cumulative mercury saturation in the experimental data is calculated based on the preset formula; wherein the preset formula is: ; ; in, , and the order of the experimental array in the sorted array; The order of arrangement is The second-order variation corresponding to the cumulative mercury saturation in the experimental data; The order of arrangement is The change in the cumulative mercury saturation in the experimental data compared with the cumulative mercury saturation in the previous experimental data; The order of arrangement is The cumulative mercury saturation in the experimental data is for , , or Among them, when hour, .

[0045] Specifically, this embodiment can be based on the above preset formula Figure 4 The experimental data shown in the figure are processed to obtain the corresponding and ,like Figure 5 shown.

[0046] Specifically, Figure 5 As shown, this embodiment can determine that the first pore throat radius is 0.0037. This embodiment can determine that the first pore throat radius is not less than 1 Corresponding =9, then the second pore throat radius can be determined to be R8=2.7974. It can be understood that the displacement pressure is P8=0.2627, and the second pore throat radius corresponding to the displacement pressure is the maximum pore throat radius of the pore.

[0047] Among them, the first standard pore throat radius is a standard pore throat radius that is not less than the first pore throat radius and not greater than the second pore throat radius.

[0048] Specifically, this embodiment can select all first standard pore throat radii from the above-mentioned multiple standard pore throat radii. At this time, this embodiment can select all first standard pore throat radii from the standard pore throat radii of 630.0000, 400.0000, 250.0000, 160.0000, 100.0000, 63.0000, 40.0000, 25.0000, 16.0000, 10.0000, 6.3000, 4.0000, 2.5000, 1.6000, 1.0000, 0.6300, 0.4000, 0.2500, 0.1600, 0.1000, 0.0630, 0.0400, 0.0250, 0.01 60, 0.0100, 0.0063, 0.0040, 0.0025, 0.0016 and 0.0010, all the first standard pore throat radii of 2.5000, 1.6000, 1.0000, 0.6300, 0.4000, 0.2500, 0.1600, 0.1000, 0.0630, 0.0400, 0.0250, 0.0160, 0.0100, 0.0063 and 0.0040, which are not less than the first pore throat radius 0.0037 and not greater than the second pore throat radius 2.7974, are selected.

[0049] Afterwards, this embodiment can use each first standard pore throat radius as an interpolated pore throat radius, interpolate the pore throat radius and cumulative mercury injection saturation of the above multiple experimental data, and obtain the interpolated cumulative mercury injection saturation corresponding to each first standard pore throat radius. Figure 6The interpolated experimental data are shown compared with the original experimental data. Figure 6 The original experimental data only reflect the corresponding relationship between the mercury injection pressure and the cumulative mercury injection saturation in the experimental data, and the interpolated experimental data also only reflect the corresponding relationship between the interpolated mercury injection pressure and the cumulative mercury injection saturation.

[0050] Optionally, each first standard pore throat radius is used as an interpolated pore throat radius, and based on each interpolated pore throat radius, interpolation is performed among the pore throat radii and cumulative mercury injection saturations of multiple experimental data to obtain the interpolated cumulative mercury injection saturation corresponding to each first standard pore throat radius, including: The pore throat radii in the multiple experimental data are taken as the original X data, the cumulative mercury injection saturation in the multiple experimental data are taken as the original Y data, and all the first standard pore throat radii are taken as the interpolation X data. The original X data, the original Y data and the interpolation X data are input into the quadratic B-spline interpolation tool for interpolation to obtain the interpolated cumulative mercury injection saturation corresponding to each first standard pore throat radius.

[0051] Optionally, the method of determining the standard cumulative mercury injection saturation corresponding to each standard pore throat radius according to each first standard pore throat radius and the corresponding interpolated cumulative mercury injection saturation includes: For any first standard pore throat radius, the interpolated cumulative mercury injection saturation corresponding to the first standard pore throat radius is directly used as the standard cumulative mercury injection saturation corresponding to the first standard pore throat radius; Determining a maximum radius and a minimum radius among a plurality of first standard pore throat radii, and determining an interpolated cumulative mercury injection saturation corresponding to the minimum radius; For any second standard pore throat radius, if it is determined that the second standard pore throat radius is greater than the maximum radius, 0 is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius; if it is determined that the second standard pore throat radius is less than the minimum radius, the interpolated cumulative mercury injection saturation corresponding to the minimum radius is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius.

[0052] Specifically, this embodiment can determine the standard cumulative mercury injection saturation corresponding to the standard pore throat radius 630.0000, 400.0000, 250.0000, 160.0000, 100.0000, 63.0000, 40.0000, 25.0000, 16.0000, 10.0000, 6.3000, 4.0000, 2.5000, 1.6000, 1.0000, 0.6300, 0.4000, 0.2500, 0.1600, 0.1000, 0.0630, 0.0400, 0.0250, 0.0160, 0.0100, 0.0063, 0.0040, 0.0025, 0.0016 and 0.0010. Figure 7 The data shown are converted into standard experimental data under standard pore throat radius, where Pz, Rz and Shgz are standard mercury injection pressure, standard pore throat radius and standard cumulative mercury injection saturation, respectively.

[0053] Figure 7 The standard cumulative mercury injection saturation in can be divided into two parts, one part is the interpolated cumulative mercury injection saturation corresponding to the first standard pore throat radius obtained by interpolation, and the other part is the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius. The standard cumulative mercury injection saturation corresponding to the second standard pore throat radius greater than the above maximum radius is set to 0, and the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius less than the above minimum radius is set to the interpolated cumulative mercury injection saturation corresponding to the above minimum radius.

[0054] The method for determining microscopic pore throat characteristic data proposed in this embodiment can convert experimental data into standard experimental data under standard pore throat radius according to the set standard pore throat radius array, effectively ensuring the realization of unified experimental data of different laboratories.

[0055] based on Figure 1 In the third method for determining microscopic pore throat characteristic data proposed in this embodiment, the microscopic pore throat characteristic data includes at least one of a permeability contribution value, a pore throat radius mean, an average pore throat radius and a sorting coefficient.

[0056] Specifically, in this embodiment, all standard experimental data can be sorted in descending order of the standard pore throat radius in the standard experimental data to obtain a target sorting array, such as Figure 7 Standard experimental data shown.

[0057] Among them, regarding the calculation of the permeability contribution value, this embodiment can calculate the permeability contribution value based on the target sorting array and the following formula: , and .

[0058] ; ; .

[0059] in, is the target sorted array, the order is The change in the standard cumulative mercury saturation in the standard experimental data compared to the standard cumulative mercury saturation in the previous standard experimental data. In the target sorted array, the order is The standard cumulative mercury saturation in the standard experimental data is: for or Among them, when hour, . is the target sorted array, the order is The standard mercury injection pressure in the standard experimental data. is the target sorted array, the order is The permeability contribution value corresponding to the standard pore throat radius in the standard experimental data. The total number of standard experimental data in the target sorting array.

[0060] Among them, regarding the calculation of the mean pore throat radius, this embodiment can first calculate Afterwards, the present embodiment can calculate the mean pore throat radius according to the following formula: .

[0061] in, is the mean pore throat radius. In the target sorted array, the order is The standard pore throat radius in the standard experimental data.

[0062] Specifically, regarding the calculation of the average pore throat radius, this embodiment can first calculate Afterwards, the average pore throat radius can be calculated according to the following formula in this embodiment: .

[0063] Specifically, regarding the calculation of the sorting coefficient, this embodiment can be calculated according to the following formulas: , , and : ; ; ; .

[0064] Afterwards, the present embodiment can calculate the sorting coefficient according to the following formula: .

[0065] It should be noted that this embodiment combines the conversion of experimental data with the calculation method of microscopic pore throat characteristic data, which can increase the comparability of microscopic pore throat characteristic data and reduce the errors caused by experimental operations in different laboratories. Figure 8As shown, based on the conversion of experimental data and microscopic pore throat characteristic data, the relationship between the radius mean and the permeability calculated for the oil and gas reservoirs in the same oil field and the same layer is similar, as shown in FIG. Fig. 9 As shown, the relationship between the sorting coefficient and the permeability calculated in this embodiment is similar based on the conversion of experimental data and microscopic pore throat characteristic data.

[0066] The method for determining microscopic pore throat characteristic data proposed in this embodiment proposes a calculation method for different microscopic pore throat parameters in the microscopic pore throat characteristic data, effectively realizes the calculation of microscopic pore throat characteristic data, and improves the calculation accuracy of microscopic pore throat characteristic data. This embodiment combines the conversion of experimental data and the calculation method of microscopic pore throat characteristic data, which can further reduce the differences in microscopic pore throat parameters caused by different laboratory conditions, increase the comparability of microscopic pore throat parameters, and reduce the errors caused by experimental operations in different laboratories. This embodiment can enhance the accuracy and comparability of microscopic pore throat parameters obtained by mercury injection experiments in different laboratories, and provide a solid and accurate foundation for oil and gas reservoir evaluation, analog research, etc.

[0067] like Fig.10 As shown, this embodiment provides a device for determining microscopic pore throat characteristic data, comprising: The first acquisition unit 101 is used to acquire a plurality of experimental data obtained by performing a high-pressure mercury injection experiment on a rock sample of an oil and gas reservoir, each of which includes a corresponding pore throat radius and a cumulative mercury injection saturation; The second acquisition unit 102 is used to acquire a plurality of standard pore throat radii; A conversion unit 103 is used to convert the multiple experimental data into standard experimental data under each standard pore throat radius according to the multiple standard pore throat radii, the pore throat radii in the multiple experimental data, and the cumulative mercury injection saturation; The determination unit 104 is used to determine the microscopic pore throat characteristic data of the rock sample based on the standard experimental data under each standard pore throat radius. It should be noted that the processing of the first acquisition unit 101, the second acquisition unit 102, the conversion unit 103 and the determination unit 104 and the beneficial effects thereof can be respectively referred to in Figure 1 Steps S101 to S104 in the above are not described in detail.

[0068] Optionally, the conversion unit 103 is further configured to: According to a plurality of standard pore throat radii, the pore throat radii in a plurality of experimental data and the cumulative mercury injection saturation, the standard cumulative mercury injection saturation corresponding to each standard pore throat radius is determined; For any standard pore throat radius, the standard pore throat radius and the corresponding standard cumulative mercury saturation are taken as the standard experimental data under the standard pore throat radius. The determining unit 104 is further configured to: The microscopic pore throat characteristic data of the rock sample are determined based on the standard pore throat radius and standard cumulative mercury saturation in each standard experimental data.

[0069] Optionally, the conversion unit 103 is further configured to: Determine a minimum first pore throat radius from among the pore throat radii of the multiple experimental data; and determine a second pore throat radius corresponding to the displacement pressure from among the pore throat radii of the multiple experimental data; Among the multiple standard pore throat radii, select multiple first standard pore throat radii that are not less than the first pore throat radius and not greater than the second pore throat radius; Taking each first standard pore throat radius as an interpolated pore throat radius, interpolating the pore throat radii and cumulative mercury injection saturations of multiple experimental data based on each interpolated pore throat radius to obtain the interpolated cumulative mercury injection saturation corresponding to each first standard pore throat radius; According to each first standard pore throat radius and the corresponding interpolated cumulative mercury injection saturation, the standard cumulative mercury injection saturation corresponding to each standard pore throat radius is determined.

[0070] Optionally, the conversion unit 103 is further configured to: According to the order of the cumulative mercury saturation in the experimental data from small to large, a plurality of experimental data are sorted to obtain a corresponding sorting array; According to the order of arrangement from small to large, each experimental data in the sorted array is traversed in turn, and for any experimental data traversed, the second-order variation corresponding to the cumulative mercury saturation in the experimental data is determined, and it is judged whether the second-order variation is greater than 1; If it is determined that the second-order variation is less than 1, continue to traverse the next experimental data; If it is determined that the second-order variation is not less than 1, the traversal is stopped, and the pore throat radius in the previous experimental data of the experimental data is determined as the second pore throat radius.

[0071] Optionally, the conversion unit 103 is further configured to: The second-order variation corresponding to the cumulative mercury saturation in the experimental data is calculated based on the preset formula; wherein the preset formula is: ; ; in, , and the order of the experimental array in the sorted array; The order of arrangement is The second-order variation corresponding to the cumulative mercury saturation in the experimental data; The order of arrangement is The change in the cumulative mercury saturation in the experimental data compared with the cumulative mercury saturation in the previous experimental data; The order of arrangement is The cumulative mercury saturation in the experimental data is for , , or ,when hour, .

[0072] Optionally, the conversion unit 103 is further configured to: For any first standard pore throat radius, the interpolated cumulative mercury injection saturation corresponding to the first standard pore throat radius is directly used as the standard cumulative mercury injection saturation corresponding to the first standard pore throat radius; Determining a maximum radius and a minimum radius among a plurality of first standard pore throat radii, and determining an interpolated cumulative mercury injection saturation corresponding to the minimum radius; For any second standard pore throat radius, if it is determined that the second standard pore throat radius is greater than the maximum radius, 0 is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius; if it is determined that the second standard pore throat radius is less than the minimum radius, the interpolated cumulative mercury injection saturation corresponding to the minimum radius is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius.

[0073] Optionally, the microscopic pore throat characteristic data includes at least one of a permeability contribution value, a pore throat radius mean, an average pore throat radius and a sorting coefficient.

[0074] The device for determining microscopic pore throat characteristic data proposed in this embodiment can obtain multiple experimental data obtained by performing high-pressure mercury injection experiments on rock samples of oil and gas reservoirs, and obtain multiple standard pore throat radii. Among them, each experimental data includes a corresponding pore throat radius and cumulative mercury injection saturation. According to multiple standard pore throat radii, the pore throat radii and cumulative mercury injection saturation in multiple experimental data, multiple experimental data are converted into standard experimental data under each standard pore throat radius. Based on the standard experimental data under each standard pore throat radius, the microscopic pore throat characteristic data of the rock sample is determined. In this embodiment, a group of specific standard pore throat radii can be set, and the experimental data obtained by different laboratories in high-pressure mercury injection experiments can be converted to the specific standard pore throat radius, so as to unify the experimental data of different laboratories, obtain multiple standard experimental data, and determine the microscopic pore throat characteristic data according to the multiple standard experimental data, effectively reducing the difference in microscopic pore throat characteristic data generated by the experimental operation differences of different laboratories, increasing the comparability of microscopic pore throat characteristic data, and reducing the errors caused by experimental operations in different laboratories.

[0075] The microscopic pore throat characteristic data determination device in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0076] The embodiment of the present invention also provides a computer device having the above Fig.10 The device for determining microscopic pore throat characteristic data is shown.

[0077] See also Fig.11 , a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. The various components are connected to each other using different buses for communication, and can be installed on a common motherboard or installed in other ways as needed. The processor can process instructions executed in the computer device, including instructions stored in or on the memory to display graphical information of the GUI on an external input / output device (such as a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides part of the necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Fig. 9 A processor 10 is taken as an example.

[0078] The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

[0079] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiment.

[0080] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function. The data storage area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage devices. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0081] The memory 20 may include a volatile memory, such as a random access memory. The memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive. The memory 20 may also include a combination of the above-mentioned types of memory.

[0082] The computer device further comprises a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0083] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0084] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining microscopic pore throat characteristic data, characterized in that: include: Acquire a plurality of experimental data obtained by performing high-pressure mercury injection experiments on rock samples of oil and gas reservoirs, each of the experimental data including a corresponding pore throat radius and a cumulative mercury injection saturation; and acquire a plurality of standard pore throat radii; According to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation, the plurality of experimental data are converted into standard experimental data under each of the standard pore throat radii; Based on the standard experimental data at each standard pore throat radius, the microscopic pore throat characteristic data of the rock sample is determined.

2. The method according to claim 1, characterized in that The converting the plurality of experimental data into standard experimental data under each of the standard pore throat radii according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation comprises: Determine the standard cumulative mercury injection saturation corresponding to each standard pore throat radius according to the multiple standard pore throat radii, the pore throat radii in the multiple experimental data, and the cumulative mercury injection saturation; For any of the standard pore throat radius, the standard pore throat radius and the corresponding standard cumulative mercury injection saturation are taken as a whole as standard experimental data under the standard pore throat radius; The method of determining the microscopic pore throat characteristic data of the rock sample based on the standard experimental data under each standard pore throat radius includes: Based on the standard pore throat radius and the standard cumulative mercury saturation in each of the standard experimental data, the microscopic pore throat characteristic data of the rock sample are determined.

3. The method according to claim 2, characterized in that The step of determining the standard cumulative mercury injection saturation corresponding to each standard pore throat radius according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation comprises: Determining a minimum first pore throat radius among the pore throat radii of the plurality of experimental data; and determining a second pore throat radius corresponding to the displacement pressure among the pore throat radii of the plurality of experimental data; Among the plurality of standard pore throat radii, select a plurality of first standard pore throat radii that are not less than the first pore throat radius and not greater than the second pore throat radius; Using each of the first standard pore throat radii as an interpolated pore throat radius, interpolating the pore throat radii and cumulative mercury injection saturations of the plurality of experimental data based on each of the interpolated pore throat radii to obtain the interpolated cumulative mercury injection saturation corresponding to each of the first standard pore throat radii; According to each of the first standard pore throat radii and the corresponding interpolated cumulative mercury injection saturation, the standard cumulative mercury injection saturation corresponding to each of the standard pore throat radii is determined.

4. The method according to claim 3, characterized in that The step of determining a second pore throat radius corresponding to the displacement pressure from the pore throat radii of the plurality of experimental data comprises: Sorting the plurality of experimental data in ascending order of the cumulative mercury saturation in the experimental data to obtain a corresponding sorting array; Traversing each of the experimental data in the sorted array in order from small to large, for any of the traversed experimental data, determining a second-order variation corresponding to the cumulative mercury intrusion saturation in the experimental data, and judging whether the second-order variation is greater than 1; If it is determined that the second-order variation is less than 1, continue to traverse the next experimental data; If it is determined that the second-order variation is not less than 1, the traversal is stopped, and the pore throat radius in the previous experimental data of the experimental data is determined as the second pore throat radius.

5. The method according to claim 4, characterized in that The determining of the second-order variation corresponding to the cumulative mercury intrusion saturation in the experimental data comprises: The second-order variation corresponding to the cumulative mercury saturation in the experimental data is calculated based on a preset formula; wherein the preset formula is: ； ； in, , and the order of arrangement of the experimental arrays in the sorted array; The order of arrangement is The second-order variation corresponding to the cumulative mercury saturation in the experimental data; The order of arrangement is A change in the cumulative mercury saturation in the experimental data compared to the cumulative mercury saturation in the previous experimental data; The order of arrangement is The cumulative mercury saturation in the experimental data is for , , or ,when hour, .

6. The method according to claim 3, characterized in that The step of determining the standard cumulative mercury injection saturation corresponding to each of the first standard pore throat radii and the corresponding interpolated cumulative mercury injection saturation comprises: For any of the first standard pore throat radius, the interpolated cumulative mercury injection saturation corresponding to the first standard pore throat radius is directly used as the standard cumulative mercury injection saturation corresponding to the first standard pore throat radius; Determining a maximum radius and a minimum radius among the plurality of first standard pore throat radii, and determining the interpolated cumulative mercury injection saturation corresponding to the minimum radius; For any second standard pore throat radius, if it is determined that the second standard pore throat radius is greater than the maximum radius, 0 is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius; if it is determined that the second standard pore throat radius is less than the minimum radius, the interpolated cumulative mercury injection saturation corresponding to the minimum radius is used as the standard cumulative mercury injection saturation corresponding to the second standard pore throat radius.

7. The method according to any one of claims 1 to 6, characterized in that The microscopic pore throat characteristic data includes at least one of a permeability contribution value, a pore throat radius mean value, an average pore throat radius and a sorting coefficient.

8. A device for determining microscopic pore throat characteristic data, characterized in that: include: A first acquisition unit is used to acquire a plurality of experimental data obtained by performing a high-pressure mercury injection experiment on a rock sample of an oil and gas reservoir, each of the experimental data including a corresponding pore throat radius and a cumulative mercury injection saturation; The second acquisition unit is used to acquire a plurality of standard pore throat radii; a conversion unit, configured to convert the plurality of experimental data into standard experimental data under each of the standard pore throat radii according to the plurality of standard pore throat radii, the pore throat radii in the plurality of experimental data, and the cumulative mercury injection saturation; A determination unit is used to determine the microscopic pore throat characteristic data of the rock sample based on the standard experimental data under each standard pore throat radius.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for determining microscopic pore throat characteristic data according to any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for determining microscopic pore throat characteristic data according to any one of claims 1 to 7.

Images