Pore throat radius distribution spectrum construction method, construction device and construction system

Abstract

The application provides a method, device and system for constructing a pore throat radius distribution spectrum, the method comprising: obtaining core mercury injection experiment data of at least three rock samples to obtain mercury injection pressure and corresponding mercury injection saturation, the sampling depths of the at least three rock samples being different; calculating a predetermined pore throat radius according to the mercury injection pressure; calculating amplitude values corresponding to the predetermined pore throat radius according to the corresponding mercury injection pressure and the corresponding mercury injection saturation of the at least three rock samples, at least three amplitude values corresponding to one predetermined pore throat radius; obtaining logging data, the logging data comprising porosity and target depth in one-to-one correspondence; determining the porosity corresponding to the sampling depth according to the correspondence between the porosity and the target depth to obtain sampling porosity; determining a fitting formula according to the at least three sampling porosities and the at least three amplitude values corresponding to one predetermined pore throat radius; and constructing the pore throat radius distribution spectrum according to the fitting formula.

Description

Technical Field

[0001] This application relates to the field of oil and gas reservoir evaluation technology, and more specifically, to a method, apparatus, computer-readable storage medium, processor, and system for constructing a pore throat radius distribution spectrum. Background Technology

[0002] Pore ​​structure is defined as the size and connectivity of pores in a rock. It includes two elements: porosity and connectivity. Generally, reservoirs with high porosity and good connectivity between pores are called reservoirs with good pore structure, while reservoirs with low porosity and poor pore connectivity, or high porosity but poor pore connectivity, or low porosity but good pore connectivity are all called reservoirs with poor pore structure.

[0003] As oil and gas exploration deepens, high-quality reservoirs with good pore structure and strong permeability have been largely exhausted. Instead, low-permeability, tight reservoirs with relatively complex pore structures are increasingly becoming the primary targets for oil and gas exploration. These reservoirs exhibit strong heterogeneity, and their permeability does not satisfy the classic Darcy's law, posing significant challenges to reservoir effectiveness evaluation, parameter calculation, and fluid property identification. For this type of reservoir, pore structure is a crucial factor determining its effectiveness. Generally, reservoirs with relatively good pore structure have strong permeability, and the more saturated the oil and gas reservoir, the easier it is to form a high-yield oil and gas reservoir. Conversely, reservoirs with relatively poor pore structure have incomplete oil and gas filling, or even cannot enter the rock pore space to be stored; these reservoirs are considered poor oil and gas layers or non-reservoirs. Therefore, quantitatively evaluating rock pore structure is a major way to improve the exploration efficiency of complex reservoirs and reduce development risks.

[0004] Capillary pressure curves are effective data for quantitatively evaluating reservoir pore structure. By processing capillary pressure curves, the pore throat radius distribution spectrum can be obtained. Using this spectrum, the distribution range and dominant pore throat radius intervals of the reservoir rock can be determined. Furthermore, parameters such as the average pore throat radius, maximum pore throat radius, and median radius can be calculated to classify reservoir types and determine reservoir effectiveness. However, when using capillary pressure curves to obtain pore throat radius distribution spectra for quantitative characterization of reservoir pore structure, there is a problem: the number of capillary pressure curves is very limited, making it impossible to continuously obtain pore throat radius distribution spectra for continuous quantitative characterization of reservoir pore structure. This results in the inability to quantitatively characterize the rock pore structure in well sections where coring and mercury intrusion porosimetry (MIP) experiments are not performed.

[0005] The information disclosed above in the background section is only intended to enhance the understanding of the background art of the art described herein. Therefore, the background art may contain certain information that does not constitute prior art known to those skilled in the art in this country. Summary of the Invention

[0006] The main objective of this application is to provide a method, apparatus, computer-readable storage medium, processor, and system for constructing a pore throat radius distribution spectrum, in order to solve the problem that the pore throat radius distribution spectrum of a reservoir cannot be obtained in the prior art.

[0007] According to one aspect of the present invention, a method for constructing a pore throat radius distribution spectrum is provided, comprising: acquiring core mercury injection test data of at least three rock samples to obtain mercury injection pressure and corresponding mercury injection saturation, wherein the sampling depths of the at least three rock samples are different; calculating a predetermined pore throat radius based on the mercury injection pressure, wherein the predetermined pore throat radius corresponds one-to-one with the mercury injection pressure; calculating an amplitude value corresponding to the predetermined pore throat radius based on the mercury injection pressure and corresponding mercury injection saturation of the at least three rock samples, wherein one predetermined pore throat radius corresponds to at least three amplitude values; acquiring well logging data, wherein the well logging data includes one-to-one porosity and target depth; determining the porosity corresponding to the sampling depth based on the correspondence between the porosity and the target depth, thereby obtaining the sampled porosity; determining a fitting formula based on the at least three sampled porosities and at least three amplitude values ​​corresponding to one predetermined pore throat radius, wherein the fitting formula is a relationship between the amplitude value and the porosity, and the fitting formula corresponds one-to-one with the predetermined pore throat radius; and constructing a pore throat radius distribution spectrum based on the fitting formula.

[0008] Optionally, the amplitude value corresponding to the predetermined pore throat radius is calculated based on the mercury ingress pressure and the corresponding mercury ingress saturation for at least three of the rock samples, including: determining the corresponding mercury ingress pressure based on the predetermined pore throat radius; determining the corresponding mercury ingress saturation based on the mercury ingress pressure; and calculating the amplitude value corresponding to each of the predetermined pore throat radii based on the mercury ingress saturation.

[0009] Optionally, a fitting formula is determined based on at least three sampling porosities and at least three amplitude values ​​corresponding to a predetermined pore throat radius, including: constructing multiple quadratic functions, where the independent variable of the quadratic functions is the porosity and the dependent variable is the amplitude value, and the quadratic functions correspond one-to-one with the predetermined pore throat radii; inputting the amplitude value and porosity corresponding to each rock sample into the corresponding quadratic function, and calculating the coefficients of each quadratic function; and determining the fitting formula based on the quadratic functions and the corresponding coefficients.

[0010] Optionally, constructing multiple univariate quadratic functions includes: constructing a univariate quadratic function of the amplitude value and the porosity corresponding to the predetermined throat radius to obtain an initial univariate quadratic function; and iterating the initial univariate quadratic function successively to obtain multiple univariate quadratic functions.

[0011] Optionally, constructing the distribution spectrum of the pore throat radius according to the fitting formula includes: constructing the distribution spectrum of the pore throat radius by using the logarithm of the predetermined pore throat radius as the abscissa and the amplitude value as the ordinate.

[0012] Optionally, the maximum mercury inlet pressure is greater than a predetermined pressure.

[0013] According to another aspect of the present invention, an apparatus for constructing a pore throat radius distribution spectrum is also provided, comprising: a first acquisition unit, configured to acquire core mercury intrusion porosimetry (CIMS) data of at least three rock samples, obtaining mercury intrusion pressure and corresponding mercury intrusion saturation, wherein the sampling depths of the at least three rock samples are all different; a first calculation unit, configured to calculate a predetermined pore throat radius based on the mercury intrusion pressure, wherein the predetermined pore throat radius corresponds one-to-one with the mercury intrusion pressure; and a second calculation unit, configured to calculate an amplitude value corresponding to the predetermined pore throat radius based on the mercury intrusion pressure and corresponding mercury intrusion saturation of the at least three rock samples, wherein one predetermined pore throat radius corresponds to at least three The system comprises: a first amplitude value; a second acquisition unit for acquiring logging data, the logging data including a one-to-one correspondence between porosity and target depth; a first determination unit for determining the porosity corresponding to the sampling depth based on the correspondence between the porosity and the target depth, thereby obtaining the sampled porosity; a second determination unit for determining a fitting formula based on at least three sampled porosities and at least three amplitude values ​​corresponding to a predetermined pore throat radius, the fitting formula being a relationship between the amplitude value and the porosity, the fitting formula corresponding one-to-one with the predetermined pore throat radius; and a construction unit for constructing a pore throat radius distribution spectrum based on the fitting formula.

[0014] According to another aspect of the present invention, a computer-readable storage medium is also provided, the computer-readable storage medium including a stored program, wherein the program executes any one of the methods described.

[0015] According to another aspect of the present invention, a processor is also provided, the processor being configured to run a program, wherein the program, when running, executes any of the methods described.

[0016] According to another aspect of the present invention, a system for constructing a pore throat radius distribution spectrum is also provided, including a device for constructing a pore throat radius distribution spectrum, the device being used to perform any of the methods described.

[0017] In this embodiment of the invention, the method for constructing the pore throat radius distribution spectrum involves the following steps: First, multiple predetermined pore throat radii are determined according to a sampling point layout method; experimental data from core mercury intrusion porosimetry (CIMS) experiments are obtained from at least three rock samples; then, amplitude values ​​corresponding to each predetermined pore throat radius are calculated based on the experimental data, with at least three rock samples having different sampling depths; subsequently, a fitting formula is determined based on porosity and each amplitude value, where the fitting formula is a relationship between the amplitude value and porosity, and the fitting formula corresponds one-to-one with each predetermined pore throat radius, and the porosity corresponds one-to-one with each sampling depth; finally, the pore throat radius distribution spectrum is constructed based on the fitting formula. This construction method determines the amplitude value corresponding to each predetermined pore throat radius through core mercury intrusion porosimetry. At least three of the rock samples are sampled at different depths, meaning at least three of the rock samples correspond to three different porosities. By fitting the porosity and each of the aforementioned amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude value corresponding to the predetermined pore throat radius at any depth can be calculated. Thus, a pore throat radius distribution spectrum is constructed based on the aforementioned fitting formula, so as to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the prior art. Attached Figure Description

[0018] The accompanying drawings, which form part of this application, are used to provide a further understanding of this application. The illustrative embodiments and descriptions of this application are used to explain this application and do not constitute an undue limitation of this application. In the drawings:

[0019] Figure 1 A flowchart illustrating a method for constructing a pore throat radius distribution spectrum according to an embodiment of this application is shown;

[0020] Figure 2 A schematic diagram of a continuous pore throat radius distribution spectrum according to an embodiment of this application is shown;

[0021] Figure 3 A schematic diagram of an apparatus for constructing a pore throat radius distribution spectrum according to an embodiment of this application is shown;

[0022] Figure 4 The correlation curve between the porosity of 90 core samples and the amplitude value amp(12) of the 12th pore throat radius according to one embodiment of this application is shown;

[0023] Figure 5 A morphological comparison diagram is shown between the pore throat radius distribution spectrum obtained by the method of Example 1 for a core sample according to one embodiment of this application and the pore throat radius distribution spectrum obtained by the core mercury intrusion capillary pressure experiment. Detailed Implementation

[0024] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0025] To enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present application, and not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative effort should fall within the scope of protection of the present application.

[0026] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate for the embodiments of this application described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0027] It should be understood that when an element (such as a layer, film, region, or substrate) is described as being "on" another element, the element may be directly on the other element, or there may be an intermediate element present. Furthermore, in the specification and claims, when an element is described as being "connected" to another element, the element may be "directly connected" to the other element, or "connected" to the other element via a third element.

[0028] For ease of description, the following explains some of the nouns or terms used in the embodiments of this application:

[0029] A pore throat refers to a narrow passage connecting pores in a rock or soil mass.

[0030] As mentioned in the background section, the existing technology cannot obtain the pore throat radius distribution spectrum of the reservoir. In order to solve the above problem, in a typical embodiment of this application, a method, apparatus, computer-readable storage medium, processor and system for constructing the pore throat radius distribution spectrum are provided.

[0031] According to an embodiment of this application, a method for constructing a pore throat radius distribution spectrum is provided.

[0032] Figure 1This is a flowchart of a method for constructing a pore throat radius distribution spectrum according to an embodiment of this application. Figure 1 As shown, the method includes the following steps:

[0033] Step S101: Obtain core mercury intrusion test data of at least three rock samples to obtain mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are different.

[0034] Step S102: Calculate the predetermined throat radius based on the mercury inlet pressure, where the predetermined throat radius corresponds one-to-one with the mercury inlet pressure.

[0035] Step S103: Based on the mercury ingress pressure and mercury ingress saturation corresponding to at least three of the above-mentioned rock samples, calculate the amplitude value corresponding to the predetermined pore throat radius, wherein one predetermined pore throat radius corresponds to at least three of the above-mentioned amplitude values.

[0036] Step S104: Obtain well logging data, which includes one-to-one corresponding porosity and target depth;

[0037] Step S105: Based on the correspondence between the porosity and the target depth, determine the porosity corresponding to the sampling depth to obtain the sampling porosity.

[0038] Step S106: Based on at least three of the above-mentioned sampling porosities and at least three of the above-mentioned amplitude values ​​corresponding to the above-mentioned predetermined throat radius, a fitting formula is determined. The fitting formula is the relationship between the above-mentioned amplitude value and the above-mentioned porosity. The fitting formula corresponds one-to-one with the above-mentioned predetermined throat radius.

[0039] Step S107: Construct the pore throat radius distribution spectrum according to the above fitting formula.

[0040] In the above method for constructing the pore throat radius distribution spectrum, firstly, core mercury injection test data of at least three rock samples are obtained to obtain the mercury injection pressure and the corresponding mercury injection saturation. The sampling depths of the at least three rock samples are all different. Then, a predetermined pore throat radius is calculated based on the mercury injection pressure, and the predetermined pore throat radius corresponds one-to-one with the mercury injection pressure. Next, based on the mercury injection pressure and the corresponding mercury injection saturation of the at least three rock samples, the amplitude value corresponding to the predetermined pore throat radius is calculated, and one predetermined pore throat radius corresponds to at least three amplitude values. Then, well logging data is obtained, and the well logging data includes one-to-one porosity and target depth. Next, based on the correspondence between the porosity and the target depth, the porosity corresponding to the sampling depth is determined to obtain the sampled porosity. Next, based on the at least three sampled porosities and at least three amplitude values ​​corresponding to the predetermined pore throat radius, a fitting formula is determined. The fitting formula is the relationship between the amplitude value and the porosity, and the fitting formula corresponds one-to-one with the predetermined pore throat radius. Finally, the pore throat radius distribution spectrum is constructed based on the fitting formula. This construction method determines multiple predetermined pore throat radii through core mercury intrusion porosimetry. The amplitude values ​​corresponding to each predetermined pore throat radius are determined by sampling at least three rock samples at different depths, i.e., at least three rock samples correspond to three different porosities. By fitting the porosities and the amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude values ​​corresponding to predetermined pore throat radii at any depth can be calculated. Thus, a pore throat radius distribution spectrum can be constructed based on the fitting formula to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the existing technology.

[0041] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0042] It should also be noted that 13 mercury ingress pressures were recorded in the core mercury intrusion test, therefore there are 13 predetermined pore throat radii. The formula for calculating the predetermined pore throat radius is R. i =73.5 / 2 i-2 i = 1, 2…13, R i Let be the radius of the i-th predetermined pore throat, in μm.

[0043] In one embodiment of this application, the amplitude value corresponding to the predetermined pore throat radius is calculated based on the mercury injection pressure and the corresponding mercury injection saturation for at least three of the aforementioned rock samples. This includes: determining the corresponding mercury injection pressure based on the predetermined pore throat radius; determining the corresponding mercury injection saturation based on the mercury injection pressure; and calculating the amplitude value corresponding to each of the predetermined pore throat radii based on the mercury injection saturation. Specifically, the predetermined pore throat radius corresponds one-to-one with the mercury injection pressure. The relationship curve between mercury injection pressure and mercury injection saturation is obtained from the core mercury injection experiment. The predetermined pore throat radius determines the corresponding mercury injection pressure, thereby determining the corresponding mercury injection saturation. Then, the amplitude value corresponding to the predetermined pore throat radius is calculated based on the mercury injection saturation. The formula for calculating the amplitude value amp is amp(1)=SHg(1), amp(i)=SHg(i)-SHg(i-1), i=2,3…13, where SHg(i) is the mercury injection saturation corresponding to the i-th mercury injection pressure, %.

[0044] In one embodiment of this application, a fitting formula is determined based on at least three sampled porosities and at least three amplitude values ​​corresponding to a predetermined pore throat radius. This includes: constructing multiple quadratic functions, where the independent variable of each quadratic function is the porosity, and the dependent variable is the amplitude value; each quadratic function corresponds one-to-one with a predetermined pore throat radius; inputting the amplitude value and porosity corresponding to each rock sample into the corresponding quadratic function to calculate the coefficients of each quadratic function; and determining the fitting formula based on the quadratic function and its corresponding coefficients. Specifically, the sampled porosities corresponding to at least three rock samples are determined based on well logging data, and the corresponding amplitude values ​​are calculated using mercury injection pressure and mercury saturation obtained through core mercury intrusion testing. The at least three sampled porosities and their corresponding amplitude values ​​for the 12th predetermined pore throat radius are then substituted into the... That is, a can be calculated. 12 b 12 and c 12 Similarly, substitute the other three sampling porosities and corresponding amplitude values ​​of the predetermined throat radius mentioned above into amp(i) = a i · amp(i+1) 2 +b i ﹒ amp(i+1)+c i Then the other a can be calculated. i b i and c i Thus, the above fitting formula is obtained.

[0045] In one embodiment of this application, constructing multiple univariate quadratic functions includes: constructing a univariate quadratic function relating the amplitude value and porosity to the predetermined throat radius, obtaining an initial univariate quadratic function; and iterating the initial univariate quadratic function sequentially to obtain multiple univariate quadratic functions. Specifically, constructing a univariate quadratic function relating the amplitude value and porosity to the predetermined throat radius, obtaining an initial univariate quadratic function, and the 12th predetermined throat radius R... i It is 0.072 μm, R 12 The corresponding amplitude value as a function of porosity is an initial quadratic function, i.e. By iterating through the initial quadratic function, multiple quadratic functions are obtained, namely amp(i) = a. i · amp(i+1) 2 +b i ﹒ amp(i+1)+c i , i=12, 11...1, where, a i b i and c i These are coefficients to be determined.

[0046] In one embodiment of this application, constructing the distribution spectrum of the orifice throat radius according to the above-mentioned fitting formula includes: constructing the orifice throat radius distribution spectrum by using the logarithm of the predetermined orifice throat radius as the abscissa and the amplitude value as the ordinate. Specifically, according to the fitting formula, constructing the orifice throat radius distribution spectrum by using the logarithm of the predetermined orifice throat radius as the abscissa and the amplitude value as the ordinate, a continuous orifice throat radius distribution spectrum is obtained, such as... Figure 2 As shown.

[0047] In one embodiment of this application, the maximum mercury ingress pressure is greater than a predetermined pressure. Specifically, the mercury ingress pressure in the core mercury intrusion test increases sequentially, and the predetermined pore throat radius corresponds one-to-one with the mercury ingress pressure. The maximum mercury ingress pressure being greater than the predetermined pressure ensures that the number of predetermined pore throat radii is sufficient to characterize the reservoir pore structure through a pore throat radius distribution spectrum. Furthermore, the predetermined pressure can be adjusted according to actual conditions.

[0048] This application also provides an apparatus for constructing a pore throat radius distribution spectrum. It should be noted that this apparatus can be used to execute the method for constructing a pore throat radius distribution spectrum provided in this application. The apparatus for constructing a pore throat radius distribution spectrum provided in this application is described below.

[0049] Figure 3 This is a schematic diagram of an apparatus for constructing a pore throat radius distribution spectrum according to an embodiment of this application. Figure 3 As shown, the device includes:

[0050] The first acquisition unit 10 is used to acquire core mercury intrusion test data of at least three rock samples to obtain mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are different.

[0051] The first calculation unit 20 is used to calculate the predetermined orifice throat radius based on the mercury inlet pressure, and the predetermined orifice throat radius corresponds one-to-one with the mercury inlet pressure.

[0052] The second calculation unit 30 is used to calculate the amplitude value corresponding to the predetermined pore throat radius based on the mercury ingress pressure and the corresponding mercury ingress saturation of at least three of the above-mentioned rock samples, wherein one predetermined pore throat radius corresponds to at least three of the above-mentioned amplitude values.

[0053] The second acquisition unit 40 is used to acquire logging data, which includes a one-to-one correspondence of porosity and target depth.

[0054] The first determining unit 50 is used to determine the porosity corresponding to the sampling depth based on the correspondence between the porosity and the target depth, and to obtain the sampling porosity.

[0055] The second determining unit 60 is used to determine a fitting formula based on at least three of the above-mentioned sampling porosities and at least three of the above-mentioned amplitude values ​​corresponding to the above-mentioned predetermined throat radius. The fitting formula is a relationship between the above-mentioned amplitude values ​​and the above-mentioned porosity. The fitting formula corresponds one-to-one with the above-mentioned predetermined throat radius.

[0056] Construction unit 70 is used to construct the pore throat radius distribution spectrum according to the above fitting formula.

[0057] In the aforementioned apparatus for constructing the pore throat radius distribution spectrum, the first acquisition unit acquires core mercury intrusion porosimetry data from at least three rock samples, obtaining the mercury intrusion pressure and the corresponding mercury intrusion saturation, wherein the sampling depths of the at least three rock samples are all different; the first calculation unit calculates a predetermined pore throat radius based on the mercury intrusion pressure, wherein the predetermined pore throat radius corresponds one-to-one with the mercury intrusion pressure; the second calculation unit calculates the amplitude value corresponding to the predetermined pore throat radius based on the mercury intrusion pressure and the corresponding mercury intrusion saturation of the at least three rock samples, wherein one predetermined pore throat radius corresponds to at least three amplitude values; the third... The second acquisition unit acquires well logging data, which includes a one-to-one correspondence between porosity and target depth. The first determination unit determines the porosity corresponding to the sampling depth based on the correspondence between porosity and target depth, thus obtaining the sampled porosity. The second determination unit determines a fitting formula based on at least three sampled porosities and at least three amplitude values ​​corresponding to a predetermined pore throat radius. The fitting formula is a relationship between the amplitude value and the porosity, and the fitting formula corresponds one-to-one with the predetermined pore throat radius. The construction unit constructs a pore throat radius distribution spectrum based on the fitting formula. This constructing device determines multiple predetermined pore throat radii through core mercury intrusion porosimetry. The amplitude values ​​corresponding to each predetermined pore throat radius are determined by sampling at least three rock samples at different depths, i.e., at least three rock samples correspond to three different porosities. By fitting the porosity and each of the aforementioned amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude values ​​corresponding to predetermined pore throat radii at any depth can be calculated. Thus, a pore throat radius distribution spectrum can be constructed based on the aforementioned fitting formula, so as to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the prior art.

[0058] It should also be noted that 13 mercury ingress pressures were recorded in the core mercury intrusion test, therefore there are 13 predetermined pore throat radii. The formula for calculating the predetermined pore throat radius is R. i =73.5 / 2 i-2 i = 1, 2…13, R i Let be the radius of the i-th predetermined pore throat, in μm.

[0059] In one embodiment of this application, the second calculation unit includes a first determining module, a first determining module, and a first calculation module. The first determining module is used to determine the corresponding mercury inlet pressure based on the predetermined pore throat radius; the second determining module is used to determine the corresponding mercury inlet saturation based on the mercury inlet pressure; and the first calculation module is used to calculate the amplitude value corresponding to each predetermined pore throat radius based on the mercury inlet saturation. Specifically, the predetermined pore throat radius corresponds one-to-one with the mercury inlet pressure. A core mercury intrusion porosimetry experiment yields a curve showing the relationship between the mercury inlet pressure and the mercury inlet saturation. The predetermined pore throat radius determines the corresponding mercury inlet pressure, thereby determining the corresponding mercury inlet saturation. Then, the amplitude value corresponding to the predetermined pore throat radius is calculated based on the mercury inlet saturation. The formula for calculating the amplitude value amp is amp(1) = SHg(1), amp(i) = SHg(i) - SHg(i-1), i = 2, 3…13, where SHg(i) is the mercury inlet saturation corresponding to the i-th mercury inlet pressure, %.

[0060] In one embodiment of this application, the second determining unit includes a first construction module, a second calculation module, and a third determining module. The first construction module constructs multiple quadratic functions, where the independent variable of each quadratic function is the porosity, and the dependent variable is the amplitude value. Each quadratic function corresponds one-to-one with a predetermined pore throat radius. The second calculation module inputs the amplitude value and porosity corresponding to each rock sample into the corresponding quadratic function to calculate the coefficients of each quadratic function. The determining module determines the fitting formula based on the quadratic function and its corresponding coefficients. Specifically, the sampling porosity corresponding to at least three rock samples is determined based on well logging data, and the corresponding amplitude value is calculated using mercury injection pressure and mercury saturation obtained through core mercury intrusion testing. The at least three sampling porosities and their corresponding amplitude values ​​for the 12th predetermined pore throat radius are then substituted into the formula. That is, a can be calculated. 12 b 12 and c 12 Similarly, substitute the other three sampling porosities and corresponding amplitude values ​​of the predetermined throat radius mentioned above into amp(i) = a i · amp(i+1) 2 +b i ﹒ amp(i+1)+c i Then the other a can be calculated. i b i and c i Thus, the above fitting formula is obtained.

[0061] In one embodiment of this application, the first construction module includes a construction submodule and a processing submodule. The construction submodule is used to construct a univariate quadratic function relating the amplitude value and porosity to the predetermined throat radius, obtaining an initial univariate quadratic function. The processing submodule is used to iterate the initial univariate quadratic function sequentially to obtain multiple univariate quadratic functions. Specifically, a univariate quadratic function relating the amplitude value and porosity to the predetermined throat radius is constructed to obtain an initial univariate quadratic function, and the 12th predetermined throat radius R... i It is 0.072 μm, R 12 The corresponding amplitude value as a function of porosity is an initial quadratic function, i.e. By iterating through the initial quadratic function, multiple quadratic functions are obtained, namely amp(i) = a. i · amp(i+1) 2 +b i ﹒ amp(i+1)+c i , i=12, 11...1, where, a i b i and c i These are coefficients to be determined.

[0062] In one embodiment of this application, the aforementioned constructing unit includes a second constructing module. This second constructing module is used to construct a distribution spectrum of the orifice throat radius using the logarithm of the predetermined orifice throat radius as the abscissa and the amplitude value as the ordinate. Specifically, according to the fitting formula, a continuous orifice throat radius distribution spectrum is obtained by constructing a distribution spectrum using the logarithm of the predetermined orifice throat radius as the abscissa and the amplitude value as the ordinate. Figure 2 As shown.

[0063] In one embodiment of this application, the maximum mercury ingress pressure is greater than a predetermined pressure. Specifically, the mercury ingress pressure in the core mercury intrusion test increases sequentially, and the predetermined pore throat radius corresponds one-to-one with the mercury ingress pressure. The maximum mercury ingress pressure being greater than the predetermined pressure ensures that the number of predetermined pore throat radii is sufficient to characterize the reservoir pore structure through a pore throat radius distribution spectrum. Furthermore, the predetermined pressure can be adjusted according to actual conditions.

[0064] This application also provides a system for constructing a pore throat radius distribution spectrum, including a device for constructing a pore throat radius distribution spectrum, wherein the device is used to perform any of the methods described above.

[0065] The aforementioned pore throat radius distribution spectrum construction system includes a pore throat radius distribution spectrum construction device. A first acquisition unit acquires core mercury intrusion porosimetry (CIMS) data from at least three rock samples, obtaining the mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are all different. A first calculation unit calculates a predetermined pore throat radius based on the mercury intrusion pressure, with each predetermined pore throat radius corresponding to a specific mercury intrusion pressure. A second calculation unit calculates the amplitude value corresponding to the predetermined pore throat radius based on the mercury intrusion pressure and corresponding mercury intrusion saturation of the at least three rock samples. One predetermined pore throat radius corresponds to at least three... The above amplitude value; the second acquisition unit acquires logging data, the logging data including a one-to-one correspondence between porosity and target depth; the first determination unit determines the porosity corresponding to the sampling depth based on the correspondence between the porosity and the target depth, and obtains the sampled porosity; the second determination unit determines a fitting formula based on at least three of the sampled porosities and at least three of the amplitude values ​​corresponding to the predetermined pore throat radius, the fitting formula being the relationship between the amplitude value and the porosity, the fitting formula corresponding one-to-one with the predetermined pore throat radius; the construction unit constructs a pore throat radius distribution spectrum based on the fitting formula. This constructing device determines multiple predetermined pore throat radii through core mercury intrusion porosimetry. The amplitude values ​​corresponding to each predetermined pore throat radius are determined by sampling at least three rock samples at different depths, i.e., at least three rock samples correspond to three different porosities. By fitting the porosity and each of the aforementioned amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude values ​​corresponding to predetermined pore throat radii at any depth can be calculated. Thus, a pore throat radius distribution spectrum can be constructed based on the aforementioned fitting formula, so as to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the prior art.

[0066] In order to enable those skilled in the art to better understand the technical solution of this application, the technical solution of this application will be described below in conjunction with specific embodiments.

[0067] Example 1

[0068] The method for constructing the above-mentioned pore throat radius distribution spectrum in this embodiment includes the following steps:

[0069] Step 1: Conduct conventional physical properties and mercury intrusion porosimetry experiments on 90 representative core samples to obtain experimental data on core porosity, capillary pressure, and corresponding pore throat radius distribution spectrum.

[0070] Step 2: Based on the above in Step 1, calculate the values ​​of the 13 orifice throat radii Rc(i) and the amplitude values ​​corresponding to different orifice throat radii using the following formula:

[0071] R i =73.5 / 2i-2 i = 1, 2, ..., 13

[0072] amp(1) = SHg(1)

[0073] amp(i)=SHg(i)-SHg(i-1), i=1, 2…13

[0074] In the formula, R i Let P be the pressure of the i-th mercury inlet. i The corresponding orifice throat radius, μm; SHg(i) is the i-th mercury inlet pressure P. i The corresponding mercury saturation, %; amp(i) is the radius R of the i-th pore throat. i The corresponding amplitude, %.

[0075] Step 3: Analyze the different pore throat radii R of the 90 core samples. i The corresponding amplitude amp(i) and porosity The correlation between them was found to be that the amplitude amp(12) of the 12th known throat radius Rc(i) = 0.072 μm is related to porosity. The relationship between them exhibits a very strong quadratic function; as porosity increases, the amplitude value increases significantly, and the correlation is as follows: Figure 4 As shown. Based on the quadratic function relationship between the two, the value of amp(12) is calculated using the following formula: In the formula, a 12 b 12 and c 12 For the coefficient to be determined, porosity Substituting amp(12) obtained from the above mercury porosimetry experimental data into the above formula, we can calculate a. 12 =0.0306, b 12 =0.1189, c 12 =2.9047.

[0076] Step four involves analyzing the correlation between the amplitude values ​​amp(i) corresponding to different pore throat radii Rc(i) in the 90 core samples. The results are shown in Table 1. The table shows that for the 90 core samples, the correlation between the amplitudes corresponding to two adjacent pore throat radii is the strongest. As the difference in pore throat radii increases, the correlation between the corresponding amplitudes also decreases. Therefore, the porosity used in step three... After calculating the value of amp(12), the amplitude values ​​corresponding to other aperture throat radii are calculated by iteratively based on the correlation between the amplitudes corresponding to adjacent aperture throat radii.

[0077] It should be noted that the core samples used in this embodiment were taken from a sandstone-conglomerate reservoir. The rock is relatively dense, and the main pore throat radius distribution ranges from 0.036 to 18.375 μm. Therefore, Table 1 only analyzes the correlation between the amplitudes corresponding to different pore throat radii when the pore throat radius is less than 18.375 μm. Furthermore, when calculating the amplitudes corresponding to different pore throat radii using the iterative method, the following formula is used to calculate only the amplitude values ​​for pore throat radii less than 18.375 μm: amp(i) = a i · amp(i+1) 2 +b i ﹒ amp(i+1)+c i , i=12,11...4, where, a i b i and c i These are coefficients to be determined.

[0078] Table 1

[0079]

[0080] Step 5: Plot a graph with the known pore throat radius Rc(i) as the logarithmic abscissa and the calculated amplitude amp(i) corresponding to different pore throat radii Rc(i) as the linear ordinate. This will allow you to construct the pore throat radius distribution spectrum of the rock using porosity. Figure 5 This is a morphological comparison diagram of the pore throat radius distribution spectrum obtained using the method described in this embodiment and the pore throat radius distribution spectrum obtained using actual core mercury intrusion porosimetry experiments. From... Figure 4 As can be seen, the two are basically the same in terms of their main distribution range and magnitude.

[0081] The aforementioned apparatus for constructing the pore throat radius distribution spectrum includes a processor and a memory. The first acquisition unit, the first calculation unit, the second calculation unit, the second acquisition unit, the first determination unit, the second determination unit, and the construction unit are all stored in the memory as program units. The processor executes the aforementioned program units stored in the memory to achieve the corresponding functions.

[0082] The processor contains a kernel, which retrieves the corresponding program units from memory. One or more kernels can be configured, and adjusting kernel parameters can address the problem of obtaining the pore throat radius distribution spectrum of reservoirs, a problem not found in existing technologies.

[0083] The memory may include non-permanent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0084] This invention provides a computer-readable storage medium having a program stored thereon, which, when executed by a processor, implements the above-described method.

[0085] This invention provides a processor for running a program, wherein the program executes the method described above when it runs.

[0086] This invention provides a device including a processor, a memory, and a program stored in the memory and executable on the processor. When the processor executes the program, it performs at least the following steps:

[0087] Step S101: Obtain core mercury intrusion test data of at least three rock samples to obtain mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are different.

[0088] Step S102: Calculate the predetermined throat radius based on the mercury inlet pressure, where the predetermined throat radius corresponds one-to-one with the mercury inlet pressure.

[0089] Step S103: Based on the mercury ingress pressure and mercury ingress saturation corresponding to at least three of the above-mentioned rock samples, calculate the amplitude value corresponding to the predetermined pore throat radius, wherein one predetermined pore throat radius corresponds to at least three of the above-mentioned amplitude values.

[0090] Step S104: Obtain well logging data, which includes one-to-one corresponding porosity and target depth;

[0091] Step S105: Based on the correspondence between the porosity and the target depth, determine the porosity corresponding to the sampling depth to obtain the sampling porosity.

[0092] Step S106: Based on at least three of the above-mentioned sampling porosities and at least three of the above-mentioned amplitude values ​​corresponding to the above-mentioned predetermined throat radius, a fitting formula is determined. The fitting formula is the relationship between the above-mentioned amplitude value and the above-mentioned porosity. The fitting formula corresponds one-to-one with the above-mentioned predetermined throat radius.

[0093] Step S107: Construct the pore throat radius distribution spectrum according to the above fitting formula.

[0094] The devices mentioned in this article can be servers, PCs, tablets, mobile phones, etc.

[0095] This application also provides a computer program product, which, when executed on a data processing device, is suitable for executing an initialization program having at least the following method steps:

[0096] Step S101: Obtain core mercury intrusion test data of at least three rock samples to obtain mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are different.

[0097] Step S102: Calculate the predetermined throat radius based on the mercury inlet pressure. The predetermined throat radius corresponds one-to-one with the mercury inlet pressure.

[0098] Step S103: Based on the mercury ingress pressure and mercury ingress saturation corresponding to at least three of the above-mentioned rock samples, calculate the amplitude value corresponding to the predetermined pore throat radius, wherein one predetermined pore throat radius corresponds to at least three of the above-mentioned amplitude values.

[0099] Step S104: Obtain well logging data, which includes one-to-one corresponding porosity and target depth;

[0100] Step S105: Based on the correspondence between the porosity and the target depth, determine the porosity corresponding to the sampling depth to obtain the sampling porosity.

[0101] Step S106: Based on at least three of the above-mentioned sampling porosities and at least three of the above-mentioned amplitude values ​​corresponding to the above-mentioned predetermined throat radius, a fitting formula is determined. The fitting formula is the relationship between the above-mentioned amplitude value and the above-mentioned porosity. The fitting formula corresponds one-to-one with the above-mentioned predetermined throat radius.

[0102] Step S107: Construct the pore throat radius distribution spectrum according to the above fitting formula.

[0103] In the above embodiments of the present invention, the descriptions of each embodiment have different focuses. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0104] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. The device embodiments described above are merely illustrative; for example, the division of units described above can be a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between units or modules may be electrical or other forms.

[0105] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0106] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0107] If the aforementioned integrated units are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a computer-readable storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned computer-readable storage medium includes various media capable of storing program code, such as USB flash drives, read-only memory (ROM), random access memory (RAM), portable hard drives, magnetic disks, or optical disks.

[0108] As can be seen from the above description, the embodiments of this application achieve the following technical effects:

[0109] 1) In the method for constructing the pore throat radius distribution spectrum of this application, firstly, core mercury intrusion porosimetry data of at least three rock samples are obtained to obtain the mercury intrusion pressure and the corresponding mercury intrusion saturation, and the sampling depths of the at least three rock samples are all different; then, a predetermined pore throat radius is calculated based on the mercury intrusion pressure, and the predetermined pore throat radius corresponds one-to-one with the mercury intrusion pressure; subsequently, based on the mercury intrusion pressure and the corresponding mercury intrusion saturation of the at least three rock samples, the amplitude value corresponding to the predetermined pore throat radius is calculated, and one predetermined pore throat radius corresponds to at least three amplitude values. The process involves several steps: First, acquiring well logging data, which includes a one-to-one correspondence between porosity and target depth. Second, determining the porosity corresponding to the sampling depth based on the correspondence between porosity and target depth, thus obtaining the sampled porosity. Third, determining a fitting formula based on at least three sampled porosities and at least three amplitude values ​​corresponding to a predetermined pore throat radius, where the fitting formula is a relationship between the amplitude value and porosity, and the fitting formula corresponds one-to-one with the predetermined pore throat radius. Finally, constructing a pore throat radius distribution spectrum based on the fitting formula. This construction method determines multiple predetermined pore throat radii through core mercury intrusion porosimetry. The amplitude values ​​corresponding to each predetermined pore throat radius are determined by sampling at least three rock samples at different depths, i.e., at least three rock samples correspond to three different porosities. By fitting the porosities and the amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude values ​​corresponding to predetermined pore throat radii at any depth can be calculated. Thus, a pore throat radius distribution spectrum can be constructed based on the fitting formula to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the existing technology.

[0110] 2) In the pore throat radius distribution spectrum construction device of this application, the first acquisition unit acquires core mercury injection test data of at least three rock samples to obtain mercury injection pressure and corresponding mercury injection saturation, and the sampling depths of the at least three rock samples are all different; the first calculation unit calculates a predetermined pore throat radius based on the mercury injection pressure, and the predetermined pore throat radius corresponds one-to-one with the mercury injection pressure; the second calculation unit calculates the amplitude value corresponding to the predetermined pore throat radius based on the mercury injection pressure and corresponding mercury injection saturation of the at least three rock samples, and one predetermined pore throat radius corresponds to at least three amplitude values. The second acquisition unit acquires well logging data, which includes a one-to-one correspondence between porosity and target depth. The first determination unit determines the porosity corresponding to the sampling depth based on the correspondence between porosity and target depth, thus obtaining the sampled porosity. The second determination unit determines a fitting formula based on at least three sampled porosities and at least three amplitude values ​​corresponding to a predetermined pore throat radius. The fitting formula is a relationship between the amplitude value and the porosity, and the fitting formula corresponds one-to-one with the predetermined pore throat radius. The construction unit constructs a pore throat radius distribution spectrum based on the fitting formula. This constructing device determines multiple predetermined pore throat radii through core mercury intrusion porosimetry. The amplitude values ​​corresponding to each predetermined pore throat radius are determined by sampling at least three rock samples at different depths, i.e., at least three rock samples correspond to three different porosities. By fitting the porosity and each of the aforementioned amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude values ​​corresponding to predetermined pore throat radii at any depth can be calculated. Thus, a pore throat radius distribution spectrum can be constructed based on the aforementioned fitting formula, so as to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the prior art.

[0111] 3) The pore throat radius distribution spectrum construction system of this application includes a pore throat radius distribution spectrum construction device. A first acquisition unit acquires core mercury intrusion porosimetry (MIP) data from at least three rock samples, obtaining the mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are all different. A first calculation unit calculates a predetermined pore throat radius based on the mercury intrusion pressure, and the predetermined pore throat radius corresponds one-to-one with the mercury intrusion pressure. A second calculation unit calculates the amplitude value corresponding to the predetermined pore throat radius based on the mercury intrusion pressure and corresponding mercury intrusion saturation of the at least three rock samples. One predetermined pore throat radius corresponds to... The first unit acquires at least three of the aforementioned amplitude values; the second acquisition unit acquires logging data, which includes a one-to-one correspondence between porosity and target depth; the first determination unit determines the porosity corresponding to the sampling depth based on the correspondence between the aforementioned porosity and the aforementioned target depth, thereby obtaining the sampled porosity; the second determination unit determines a fitting formula based on at least three of the aforementioned sampled porosities and at least three of the aforementioned amplitude values ​​corresponding to the aforementioned predetermined pore throat radius, wherein the fitting formula is a relationship between the aforementioned amplitude values ​​and the aforementioned porosity, and the fitting formula corresponds one-to-one with the aforementioned predetermined pore throat radius; the construction unit constructs a pore throat radius distribution spectrum based on the aforementioned fitting formula. This constructing device determines multiple predetermined pore throat radii through core mercury intrusion porosimetry. The amplitude values ​​corresponding to each predetermined pore throat radius are determined by sampling at least three rock samples at different depths, i.e., at least three rock samples correspond to three different porosities. By fitting the porosity and each of the aforementioned amplitude values, a fitting formula can be obtained. Based on the fitting formula, the amplitude values ​​corresponding to predetermined pore throat radii at any depth can be calculated. Thus, a pore throat radius distribution spectrum can be constructed based on the aforementioned fitting formula, so as to achieve the purpose of continuously and quantitatively characterizing the reservoir pore structure and evaluating the reservoir effectiveness. This solves the problem that the pore throat radius distribution spectrum of the reservoir cannot be obtained in the prior art.

[0112] The above description is merely a preferred embodiment of this application and is not intended to limit this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method for constructing a pore throat radius distribution spectrum, characterized in that, include: Obtain core mercury intrusion test data from at least three rock samples to obtain mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are different. The predetermined throat radius is calculated based on the mercury inlet pressure, and the predetermined throat radius corresponds one-to-one with the mercury inlet pressure. Based on the mercury ingress pressure and mercury ingress saturation corresponding to at least three of the rock samples, the amplitude value corresponding to the predetermined pore throat radius is calculated, and one predetermined pore throat radius corresponds to at least three amplitude values. Acquire well logging data, which includes a one-to-one correspondence of porosity and target depth; Based on the correspondence between the porosity and the target depth, the porosity corresponding to the sampling depth is determined, and the sampling porosity is obtained. Based on at least three amplitude values ​​corresponding to at least three sampling porosities and one predetermined throat radius, a fitting formula is determined, wherein the fitting formula is a relationship between the amplitude value and the porosity, and the fitting formula corresponds one-to-one with the predetermined throat radius. The pore throat radius distribution spectrum is constructed based on the fitting formula.

2. The method according to claim 1, characterized in that, Based on the mercury ingress pressure and mercury ingress saturation corresponding to at least three of the aforementioned rock samples, the amplitude value corresponding to the predetermined pore throat radius is calculated, including: The corresponding mercury inlet pressure is determined based on the predetermined orifice radius; The corresponding mercury saturation is determined based on the mercury inlet pressure. The amplitude value corresponding to each predetermined throat radius is calculated based on the mercury saturation.

3. The method according to claim 1, characterized in that, A fitting formula is determined based on at least three sampling porosities and at least three amplitude values ​​corresponding to a predetermined throat radius, including: Multiple quadratic functions are constructed, wherein the independent variable of the quadratic function is the porosity, the dependent variable of the quadratic function is the amplitude value, and the quadratic function corresponds one-to-one with the predetermined pore throat radius; The amplitude value and porosity corresponding to each rock sample are input into the corresponding quadratic function to calculate the coefficients of each quadratic function. The fitting formula is determined based on the quadratic function and its corresponding coefficients.

4. The method according to claim 3, characterized in that, Construct multiple quadratic functions in one variable, including: Construct a univariate quadratic function of the amplitude value and the porosity corresponding to the predetermined throat radius to obtain the initial univariate quadratic function; The initial quadratic function is iterated successively to obtain multiple quadratic functions.

5. The method according to claim 1, characterized in that, Constructing the distribution spectrum of the pore throat radius according to the fitting formula includes: constructing the distribution spectrum of the pore throat radius by using the logarithm of the predetermined pore throat radius as the abscissa and the amplitude value as the ordinate.

6. The method according to claim 1, characterized in that, The maximum mercury inlet pressure is greater than the predetermined pressure.

7. A device for constructing a pore throat radius distribution spectrum, characterized in that, include: The first acquisition unit is used to acquire core mercury intrusion test data of at least three rock samples, and obtain mercury intrusion pressure and corresponding mercury intrusion saturation. The sampling depths of the at least three rock samples are different. The first calculation unit is used to calculate the predetermined throat radius based on the mercury inlet pressure, wherein the predetermined throat radius corresponds one-to-one with the mercury inlet pressure. The second calculation unit is used to calculate the amplitude value corresponding to the predetermined pore throat radius based on the mercury ingress pressure and the corresponding mercury ingress saturation of at least three of the rock samples, wherein one predetermined pore throat radius corresponds to at least three amplitude values. The second acquisition unit is used to acquire well logging data, which includes a one-to-one correspondence of porosity and target depth. The first determining unit is used to determine the porosity corresponding to the sampling depth based on the correspondence between the porosity and the target depth, and to obtain the sampling porosity. The second determining unit is used to determine a fitting formula based on at least three amplitude values ​​corresponding to at least three sampling porosities and one predetermined throat radius, wherein the fitting formula is a relationship between the amplitude value and the porosity, and the fitting formula corresponds one-to-one with the predetermined throat radius. A construction unit is used to construct the pore throat radius distribution spectrum according to the fitting formula.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium includes a stored program, wherein the program performs the method according to any one of claims 1 to 6.

9. A processor, characterized in that, The processor is used to run a program, wherein the program executes the method according to any one of claims 1 to 6 when it runs.

10. A system for constructing a pore throat radius distribution spectrum, comprising an apparatus for constructing a pore throat radius distribution spectrum, characterized in that, The construction apparatus is used to perform the method according to any one of claims 1 to 6.

Images