Full-aperture pore structure characterization method and device based on multidimensional data driving

Abstract

The present invention provides a full-aperture pore structure characterization method and device driven by multidimensional data, the method comprising: obtaining a target tight reservoir rock sample; subjecting the target tight reservoir rock sample to brine saturation treatment and a first nuclear magnetic resonance measurement to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation; subjecting the brine-saturated rock sample to high-pressure mercury injection and nitrogen adsorption to obtain high-pressure mercury injection and gas adsorption experimental data; calibrating the transverse relaxation time under brine saturation based on a full-aperture cumulative pore volume frequency distribution curve to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; subjecting the target tight reservoir rock sample to constant-rate mercury injection treatment to obtain a pore size distribution curve based on constant-rate mercury injection; determining the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection, and using the overlapping position as a joining point; and determining the full-aperture distribution range of the target tight reservoir.

Description

Technical Field

[0001] The present invention relates to the field of oil and gas exploration technology, and in particular to a full-aperture pore structure characterization method and device based on multidimensional data driving. Background Art

[0002] With the development of unconventional oil and gas exploration and development technologies, the proportion of tight sandstone oil and gas has increased annually. Compared with conventional sandstone reservoirs, tight reservoirs have a more complex microscopic pore structure, generally characterized by small pore throats, uneven pore space distribution, and relatively complex pore throat configurations.

[0003] Tight sandstone reservoirs have a multi-scale pore structure, characterized by both macro- and mesopores as well as micro- and nanopores. Conventional cast thin-section methods and field-emission scanning electron microscopy, due to their limited measurement range, can only characterize pores in the macro- and mesopores range, making it difficult to characterize microscopic pore types at the micron level.

[0004] It can be seen from this that the related art methods for characterizing the pore structure of tight reservoirs have the technical problem of being unable to achieve the full pore size distribution characterization of tight reservoir rocks. Summary of the Invention

[0005] The present invention provides a full-aperture pore structure characterization method and device based on multidimensional data driving, which is used to solve the defect of the tight reservoir pore structure characterization method in the existing technology that it is difficult to achieve full-aperture distribution characterization of tight reservoir rocks, and to establish a complete full-aperture cumulative distribution curve, thereby characterizing the full-aperture pore structure of tight reservoir rocks.

[0006] The present invention provides a full-aperture pore structure characterization method based on multidimensional data driving, which includes the following steps. Obtain a target tight reservoir rock sample; perform brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain a brine-saturated rock sample and a transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; perform high-pressure mercury injection and nitrogen adsorption on the brine-saturated rock sample to obtain high-pressure mercury injection and gas adsorption experimental data; draw a full-pore cumulative pore volume frequency distribution curve based on the high-pressure mercury injection and gas adsorption experimental data and pore throat size distribution; calibrate the transverse relaxation time under brine saturation based on the full-pore cumulative pore volume frequency distribution curve to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; perform constant-rate mercury injection treatment on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection; determine the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection, and use the overlapping position as a junction point; and determine the full-pore size distribution range of the target tight reservoir based on the junction point.

[0007] According to a full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention, after the target tight reservoir rock sample is subjected to brine saturation treatment and a first nuclear magnetic resonance measurement to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, the method further includes: centrifuging the brine-saturated rock sample and performing a second nuclear magnetic resonance measurement to obtain the transverse relaxation time under each centrifugal force, wherein the centrifugation treatment and the second nuclear magnetic resonance measurement are performed simultaneously.

[0008] According to a full-pore structure characterization method based on multidimensional data driving provided by the present invention, the full-pore size distribution range of the target tight reservoir is determined based on the junction point, including: determining the small pore throat data on the left side of the junction point and the large pore throat data on the right side of the junction point; and using the small pore throat data and the large pore throat data as the full-pore size distribution range of the target tight reservoir.

[0009] According to a full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention, before obtaining the target tight reservoir rock sample, the method further includes: pre-processing the first tight reservoir rock sample and the second tight reservoir rock sample to obtain a standard cylindrical core, wherein the diameter of the standard cylindrical core is the target diameter and the length of the standard cylindrical core is within the target length range; cleaning and drying the standard cylindrical core to obtain the target tight reservoir rock sample, wherein the constant temperature of the target tight reservoir rock sample is the target temperature.

[0010] According to a full-pore structure characterization method based on multidimensional data driving provided by the present invention, after drawing a full-pore cumulative pore volume frequency distribution curve based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution, the method further includes: determining the mercury saturation corresponding to the maximum mercury saturation and the displacement pressure; determining the pore throat radius at the mercury saturation and the transverse relaxation time corresponding to the pore throat radius.

[0011] According to a full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention, after performing constant-rate mercury injection on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection, the method further includes: determining a lower limit value of the minimum measurable space based on the pore size distribution curve based on constant-rate mercury injection.

[0012] The present invention also provides a full-aperture pore structure characterization device driven by multidimensional data, comprising the following modules: an acquisition module for acquiring a target tight reservoir rock sample; a first processing module for performing a brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain a brine-saturated rock sample and a transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; a second processing module for performing high-pressure mercury injection and nitrogen adsorption on the brine-saturated rock sample to obtain high-pressure mercury injection and gas adsorption experimental data; a first drawing module for drawing a pore structure characterization device based on the high-pressure mercury injection and gas adsorption experimental data and pore throat size distribution. A full-pore cumulative pore volume frequency distribution curve is obtained; a second drawing module is used to calibrate the transverse relaxation time under brine saturation based on the full-pore cumulative pore volume frequency distribution curve to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; a third processing module is used to perform constant-rate mercury injection processing on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection; a first determination module is used to determine the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection, and use the overlapping position as a combination point; a second determination module is used to determine the full-pore size distribution range of the target tight reservoir based on the combination point.

[0013] The present invention also provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the full-aperture pore structure characterization method based on multidimensional data driving as described above is implemented.

[0014] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements any of the above-described full-aperture pore structure characterization methods driven by multidimensional data.

[0015] The present invention also provides a computer program product, comprising a computer program, which, when executed by a processor, implements any of the above-described full-aperture pore structure characterization methods based on multidimensional data driving.

[0016] The present invention provides a full-aperture pore structure characterization method and device based on multidimensional data drive. The method and device obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation by subjecting the target tight reservoir rock sample to brine saturation treatment and first nuclear magnetic resonance measurement. Subsequently, the brine-saturated rock sample is subjected to high-pressure mercury injection and nitrogen adsorption to draw a full-aperture cumulative pore volume frequency distribution curve. Based on the full-aperture cumulative pore volume frequency distribution curve, the transverse relaxation time under brine saturation is calibrated to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance. The target tight reservoir rock sample is subjected to constant-rate mercury injection treatment to obtain a pore size distribution curve based on constant-rate mercury injection. The overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection is determined as a combination point to determine the full-aperture distribution range of the target tight reservoir. Thus, the data and information of multiple experiments such as high-pressure mercury injection, constant-rate mercury injection, nuclear magnetic resonance and gas adsorption are effectively integrated to well reflect the reservoir quality characteristics and the distribution of nanopores to macropores. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced one by one below. 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 any creative work.

[0018] Figure 1 It is a schematic flow chart of the full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention.

[0019] Figure 2 This is the relaxation time spectrum of the rock samples 1 and 2 tested by nuclear magnetic resonance provided by the present invention.

[0020] Figure 3 These are three experimental cumulative distribution frequency diagrams of rock sample 1 provided by the present invention.

[0021] Figure 4 It is a schematic diagram of the three-stage fitting of the rock sample 1 provided by the present invention.

[0022] Figure 5 It is the pore size distribution curve of the transverse relaxation time of rock sample 1 and rock sample 2 provided by the present invention.

[0023] Figure 6 It is a schematic diagram of the full pore size distribution range of the tight reservoir provided by the present invention.

[0024] Figure 7 This is a conceptual diagram of the full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention.

[0025] Figure 8 It is a schematic diagram of the nitrogen adsorption instrument and nitrogen adsorption data provided by the present invention.

[0026] Figure 9 It is a schematic diagram of the nuclear magnetic resonance instrument and T2 spectrum provided by the present invention.

[0027] Figure 10 It is a schematic diagram of the automatic pore structure tester (constant rate mercury intrusion tester) provided by the present invention.

[0028] Figure 11 It is a structural schematic diagram of a full-aperture pore structure table device driven by multidimensional data provided by the present invention.

[0029] Figure 12 It is a schematic diagram of the physical structure of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0030] To make the objectives, technical solutions, and advantages of the present invention more clear, the technical solutions of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the embodiments described are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0031] With the advancement of unconventional oil and gas exploration and development technologies, the proportion of tight sandstone oil and gas has increased annually. Compared with conventional sandstone reservoirs, tight reservoirs have a more complex micropore structure, generally characterized by small pore throats, uneven pore distribution, and relatively complex pore throat configurations. Tight reservoir micropore structure directly influences reservoir storage capacity, seepage capacity, and fluid distribution. Therefore, clarifying the impact of full-pore-size pore structure characterization on tight reservoir characteristics is beneficial for revealing the formation mechanism of tight sandstone reservoir densification and better controlling, optimizing, and evaluating oil and gas reservoirs.

[0032] Tight sandstone reservoirs have both large and medium pores and are rich in micron and nanoscale pores, resulting in a multi-scale pore structure. Conventional cast thin section methods and field emission scanning electron microscopy methods can only characterize pores in the large and medium pore range due to the limited measurement range, making it difficult to characterize microscopic pore types at the micron level. High-pressure mercury injection destroys the pore structure in rock samples, making it impossible to detect the reservoir space between clay layers where oil, gas, and water can enter. Gas adsorption methods mainly focus on submicron-nanoscale pore throat structures, but their response to large pore throats is weak and cannot be accurately measured. Existing nuclear magnetic resonance experiments often use T2 spectra of tight sandstone in a single oil- or water-saturated state to characterize the pore size distribution of tight sandstone. However, due to the complex and heterogeneous mixed wettability characteristics of tight sandstone, a single fluid cannot saturate the pores within the reservoir. The T2 spectra obtained in this state cannot characterize unsaturated nanopores, making it difficult to characterize the full pore size distribution of shale.

[0033] With the deepening of exploration and development, the pore structure characterization by a single method can no longer meet the current development needs. Therefore, it is necessary to combine experimental methods at different scales to carry out detailed characterization of the full-aperture pore structure.

[0034] refer to Figure 1 , Figure 1 This is a flow chart of the full-aperture pore structure characterization method based on multidimensional data drive provided by the present invention, such as Figure 1 As shown, the method includes the following steps.

[0035] Step 101: Obtain target tight reservoir rock samples.

[0036] Among them, tight sandstone reservoirs are usually sandstone reservoirs with low reservoir permeability.

[0037] In the embodiment of the present invention, two tight reservoir rock samples were collected from the 8th section of the Lower Shihezi Formation in the Jiaxian area in the northeast of the Yishan Slope in the Ordos Basin. The lithology of both samples is sandstone, of which sample 1 has a burial depth of 2078.3m and a porosity of 5.22%; sample 2 has a burial depth of 1976.7m and a porosity of 3.76%.

[0038] Step 102 : performing brine saturation treatment and nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously.

[0039] In the embodiment of the present invention, the obtained tight reservoir rock samples were sequentially saturated with brine and centrifuged, and a nuclear magnetic resonance experiment was performed simultaneously. The nuclear magnetic resonance experiment was measured using a Meso MR12-040 nuclear magnetic resonance spectrometer. Using a low magnetic field intensity of (0.5±0.08T), a waiting time of 1500ms, and 128 scans as measurement parameters, an amplitude curve of the transverse relaxation time T2 from 0.01ms to 10000ms with 60 logarithmic intervals was obtained, and the T2 spectrum under brine saturation and the T2 spectrum under each centrifugal force were obtained in turn.

[0040] refer to Figure 2 , Figure 2 This is the relaxation time spectrum of the nuclear magnetic resonance test of rock samples 1 and 2 provided by the present invention, wherein the abscissa is the T2 relaxation time (ms) and the ordinate is the amplitude; including the saturated state and after centrifugation.

[0041] Specifically, rock samples 1 and 2 were uniformly processed into standard cylindrical cores with a diameter of 25 mm and a length of 35-50 mm, washed to remove residual oil and residual salt, and dried and cooled to a constant temperature of 65°C.

[0042] The rock samples 1 and 2 treated with (1) were saturated with pressurized brine using simulated formation water. During this process, nuclear magnetic resonance measurements were performed simultaneously to obtain the T2 spectrum under brine saturation (for details, please refer to Figure 2 ).

[0043] The rock samples 1 and 2 saturated with salt water were centrifuged and centrifuged at high speed for 3 hours. During this process, nuclear magnetic resonance measurements were performed to obtain T2 spectra (for details, please refer to Figure 2 ).

[0044] Step 103 , performing high-pressure mercury injection and nitrogen adsorption on the rock sample saturated with brine to obtain high-pressure mercury injection and gas adsorption experimental data.

[0045] In the embodiment of the present invention, high-pressure mercury injection experiments and nitrogen adsorption experiments were performed on the rock samples 1 and 2 saturated with brine, respectively, to obtain high-pressure mercury injection and gas adsorption experimental data.

[0046] It should be noted that, here, the rock samples 1 and 2 (i.e., without brine saturation treatment) can also be directly subjected to high-pressure mercury injection and nitrogen adsorption treatments, respectively.

[0047] Step 104 : Based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution, a full-pore-diameter cumulative pore volume frequency distribution curve is plotted.

[0048] In an embodiment of the present invention, the pore throat size distribution determined by high pressure mercury injection testing, nitrogen adsorption and nuclear magnetic resonance testing is plotted as a cumulative distribution frequency diagram (ie, a full pore size cumulative pore volume frequency distribution curve).

[0049] The above pore throat size distribution is determined based on nuclear magnetic resonance testing.

[0050] refer to Figure 3 , Figure 3 The three experimental cumulative distribution frequency diagrams of rock sample 1 provided by the present invention are drawn, wherein the abscissa is the cumulative frequency and the ordinate is the aperture ( ), transverse relaxation time (T2 / ms), including nuclear magnetic resonance, high-pressure mercury injection, and nitrogen adsorption.

[0051] refer to Figure 4 , Figure 4 It is a schematic diagram of the three-stage fitting of the rock sample 1 provided by the present invention.

[0052] In the embodiment of the present invention, a "three-stage" fitting method is used to obtain constants C and n using the least squares method. For details, please refer to Figure 4 .

[0053] Step 105 , based on the full pore size cumulative pore volume frequency distribution curve, calibrate the transverse relaxation time under brine saturation to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance.

[0054] refer to Figure 5 , Figure 5 is the pore size distribution curve of the transverse relaxation time of rock sample 1 and rock sample 2 provided by the present invention, wherein, Figure 5 (a) Pore size distribution curve of transverse relaxation time of rock sample 1, Figure 5 (b) is the pore size distribution curve of transverse relaxation time of rock sample 2. The horizontal axis is the pore throat radius ( ), the vertical axis is the growth frequency.

[0055] Step 106: Perform constant-rate mercury injection on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection.

[0056] In an embodiment of the present invention, constant-rate mercury injection experiments are performed on the rock samples 1 and 2 saturated with brine to obtain constant-rate mercury injection experiment data. Based on the constant-rate mercury injection experiment data, a pore size distribution curve based on the constant-rate mercury injection experiment is determined.

[0057] It should be noted that, here, the constant-rate mercury injection experiments can also be directly performed on rock samples 1 and 2 (i.e., without brine saturation treatment).

[0058] Step 107 : determining the overlapping position of the pore size distribution curve based on the transverse relaxation time of nuclear magnetic resonance and the pore size distribution curve based on constant rate mercury injection, and using the overlapping position as a joining point.

[0059] In the embodiment of the present invention, the lower limit of the measurable minimum pore diameter is determined based on the constant-rate mercury injection pore throat distribution data, and the overlapping position (10 μm) of the nuclear magnetic resonance pore throat radius distribution curve and the constant-rate mercury injection pore throat radius distribution curve is used as the combination point.

[0060] Step 108: Determine the full pore size distribution range of the target tight reservoir based on the binding points.

[0061] refer to Figure 6 , Figure 6 is a schematic diagram of the full pore size distribution range of the tight reservoir provided by the present invention, wherein the abscissa is the pore throat distribution ( ), the vertical axis is frequency ( / %), specifically including the partial NMR images and constant-rate mercury injection images of rock samples 1 and 2.

[0062] In the embodiment of the present invention, the data reflecting small pore throats on the left side of the junction point and the data reflecting large pore throats on the right side are combined to obtain the full pore size distribution range of the tight reservoir (for details, please refer to Figure 6 ).

[0063] Through the present invention, nuclear magnetic resonance (NMR) is combined with high-pressure mercury injection and gas adsorption experiments to calibrate the pore size distribution of tight reservoirs using the NMR T2 spectrum. Based on the overlap between constant-rate mercury injection and the NMR T2 spectrum, the full-aperture pore size distribution of tight sandstone reservoirs, from nanopores to macropores, can be determined. This method is highly operational and is suitable for application in various exploration and research areas. The proposed method and the related research during its development contribute to the quantitative study of the full-aperture pore structure of tight reservoirs.

[0064] The target tight reservoir rock sample is obtained through the above steps provided by the embodiment of the present invention; the target tight reservoir rock sample is subjected to brine saturation treatment and nuclear magnetic resonance measurement to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, wherein the brine saturation treatment and the nuclear magnetic resonance measurement are performed simultaneously; the brine-saturated rock sample is subjected to high-pressure mercury injection and nitrogen adsorption to obtain high-pressure mercury injection and gas adsorption experimental data; based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution, a full-pore cumulative pore volume frequency distribution curve is plotted; based on the full-pore cumulative pore volume frequency distribution curve, the transverse relaxation time under brine saturation is calibrated to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; the target tight reservoir rock sample is subjected to constant-rate mercury injection treatment to obtain a pore size distribution curve based on constant-rate mercury injection; the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection is determined, and the overlapping position is used as a junction point; based on the junction point, the full-pore size distribution range of the target tight reservoir is determined. Therefore, the effective integration of data and information from multiple experiments such as high-pressure mercury injection, constant-rate mercury injection, nuclear magnetic resonance, and gas adsorption can well reflect the reservoir quality characteristics and the distribution of nanopores to macropores.

[0065] According to a full-aperture pore structure characterization method based on multidimensional data drive provided by the present invention, after subjecting a target tight reservoir rock sample to brine saturation treatment and nuclear magnetic resonance measurement to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, the method further comprises:

[0066] After performing brine saturation treatment and nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, the method further includes:

[0067] The rock sample saturated with salt water is centrifuged and subjected to nuclear magnetic resonance measurement to obtain the transverse relaxation time under each centrifugal force, wherein the centrifugation and the second nuclear magnetic resonance measurement are performed simultaneously.

[0068] In the embodiment of the present invention, rock samples 1 and 2 are uniformly processed into standard cylindrical cores with a diameter of 25 mm and a length of 35-50 mm, washed to remove residual oil and residual salt, and dried and cooled to a constant temperature of 65°C.

[0069] The rock samples 1 and 2 treated in (1) were saturated with pressurized brine using simulated formation water. During this process, nuclear magnetic resonance measurements were performed to obtain T2 spectra under brine saturation.

[0070] The rock samples 1 and 2 saturated with salt water were centrifuged and centrifuged at high speed for 3 hours. During this process, nuclear magnetic resonance measurements were performed to obtain T2 spectra (for details, please refer to Figure 2 ).

[0071] According to the present invention, a full-pore size pore structure characterization method based on multidimensional data drive is provided, which determines the full-pore size distribution range of the target tight reservoir based on the binding points, including:

[0072] Determine the small pore throat data on the left side of the junction and the large pore throat data on the right side of the junction;

[0073] The small pore throat data and large pore throat data are taken as the full pore size distribution range of the target tight reservoir.

[0074] In the embodiment of the present invention, the overlapping position (10 μm) of the nuclear magnetic resonance pore throat radius distribution curve and the constant rate mercury injection pore throat radius distribution curve is used as the combination point, and the data reflecting small pore throats on the left side of the combination point and the data reflecting large pore throats on the right side of the combination point are combined to obtain the full pore size distribution range of the tight reservoir. For details, please refer to Figure 6 .

[0075] According to a full-aperture pore structure characterization method based on multidimensional data drive provided by the present invention, before obtaining a target tight reservoir rock sample, the method further includes:

[0076] Preprocessing the first tight reservoir rock sample and the second tight reservoir rock sample to obtain a standard cylindrical core, wherein the diameter of the standard cylindrical core is a target diameter and the length of the standard cylindrical core is within a target length range;

[0077] The standard cylindrical core is cleaned and dried to obtain a target tight reservoir rock sample, wherein the constant temperature of the target tight reservoir rock sample is the target temperature.

[0078] In the embodiment of the present invention, rock samples 1 and 2 are uniformly processed into standard cylindrical cores with a diameter of 25 mm and a length of 35-50 mm, washed to remove residual oil and residual salt, and dried and cooled to a constant temperature of 65°C.

[0079] According to a full-pore-size pore structure characterization method based on multidimensional data drive provided by the present invention, after obtaining a full-pore-size cumulative pore volume frequency distribution curve based on high-pressure mercury injection and gas adsorption experimental data and pore throat size distribution, the method further includes:

[0080] Determine the mercury saturation corresponding to the maximum mercury saturation and displacement pressure;

[0081] Determine the pore throat radius at mercury saturation and the lateral relaxation time corresponding to the pore throat radius.

[0082] In the embodiment of the present invention, the mercury saturation (S(i)) corresponding to the maximum mercury saturation (Smax) and the displacement pressure (Spe) is obtained by interpolation method. At the same saturation (S(i)), the pore throat radius r(i) and the corresponding transverse relaxation time T2(i) are obtained (for details, please refer to Figure 3 ).

[0083] According to a full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention, after performing constant-rate mercury injection on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection, the method further includes:

[0084] Based on the pore size distribution curve based on constant-rate mercury intrusion, the lower limit of the minimum measurable space is determined.

[0085] In the embodiment of the present invention, the core of the full-pore experiment is to reflect the reservoir characteristics and its heterogeneity through microscopic observation. The experiments involved in this method include: gas adsorption method, nuclear magnetic resonance method, high-pressure mercury injection method, and constant-rate mercury injection method. The gas adsorption method can quantitatively analyze the distribution of small pore sizes, but it is not accurate in characterizing larger pore sizes; high-pressure mercury injection is a destructive experiment that will destroy the sample pore size during the experiment, and its effect on small pore size characterization is worse than that of gas adsorption; nuclear magnetic resonance can balance the above two methods, but two methods are required for pore size conversion. Constant-rate mercury injection mainly characterizes the distribution of large pore sizes. Combining the advantages of the above four experiments not only helps to evaluate the potential production capacity of the reservoir, but also provides a solid geological basis for the formulation, adjustment and optimization of oilfield development plans.

[0086] refer to Figure 7 , Figure 7 This is a conceptual diagram of the full-aperture pore structure characterization method based on multidimensional data driving provided by the present invention.

[0087] refer to Figure 7 On the left, in an embodiment of the present invention, a target tight reservoir rock sample is subjected to gas adsorption, nuclear magnetic resonance, high-pressure mercury injection, and other processing processes to determine the cumulative distribution frequency and T2 spectrum conversion (i.e., transverse relaxation time); a pore throat distribution curve is determined based on the cumulative distribution frequency and T2 spectrum conversion; and based on the pore throat distribution curve, a quantitative characterization of the small pore structure is determined.

[0088] refer to Figure 7 On the right, in an embodiment of the present invention, a target tight reservoir rock sample is subjected to constant-rate mercury injection, and the constant-rate mercury injection data is transformed to obtain a pore-throat distribution curve; and based on the pore-throat distribution curve, a quantitative characterization of the macropore structure is determined.

[0089] Subsequently, based on the quantitative characterization of small pore structure and large pore structure, the full pore size cumulative distribution curve is determined to obtain the full pore size pore characterization.

[0090] refer to Figure 8 , Figure 8 It is a schematic diagram of the nitrogen adsorption instrument and nitrogen adsorption data provided by the present invention.

[0091] The gas adsorption method implemented in this invention uses an ASAP2020 physical adsorption instrument with a diameter measurement range of 0.35–500.00 nm. Before the experiment, 1–3 g of the sample was ground to less than 250 μm. The sample was then dried in a vacuum at 150°C for 3 hours until all moisture and volatile gases were removed, and then placed in the instrument for analysis. The experimental temperature was 77 K (−196°C), the pressure was 97.3–127.0 MPa, the powder particle size was 0.28–0.18 mm, and the nitrogen adsorption concentration was 99.99%.

[0092] In gas adsorption experiments, the adsorption and desorption isotherms of nitrogen are measured at different relative pressures (P / P0). When P / P0 is between 0.05 and 0.35, nitrogen is assumed to exist only on the pore surface. The nitrogen adsorption capacity and P / P0 conform to the Brunauer-Emmet-Teller (BET) equation, allowing the specific surface area (SBET) of the sample to be calculated. As P / P0 increases, nitrogen begins to solidify in the micropores.

[0093] Assuming that all pores in the sample are filled with liquid adsorbent at the maximum relative pressure (P / P0=1), the total amount of nitrogen adsorbed is the total pore volume (Vads) of the sample.

[0094] refer to Figure 9 , Figure 9 It is a schematic diagram of the nuclear magnetic resonance instrument and T2 spectrum provided by the present invention.

[0095] In the nuclear magnetic resonance implemented in the present invention, the experiment consists of a saturation experiment and a centrifugation experiment. First, the dense sandstone sample is removed from impurities and hydrocarbons in an oven at 105°C for 72 hours to remove the remaining water. Then, the mass and size of the dried sample are measured. A CaCl2 solution with a salinity of 25,000 ppm is used as the nuclear magnetic resonance experimental liquid. The sample is then saturated with the experimental liquid for 48 hours at a saturation pressure of 20 MPa, and the mass of the sample under saturated water conditions is measured. After sample treatment, the transverse relaxation time (T2) spectrum of the sample is measured. After measuring the T2 spectrum, the sample is centrifuged at 413 PSI. Previous studies have shown that dense sandstone may lose almost all movable water at a centrifugal pressure of 413 psi. The T2 spectrum of the sample is measured after centrifugation.

[0096] refer to Figure 10 , Figure 10 It is a schematic diagram of the automatic pore structure tester (constant rate mercury intrusion tester) provided by the present invention.

[0097] First, the sample was placed in a closed dilatometer for vacuum testing. To ensure that the pore space was not destroyed by mercury, the maximum mercury pressure was set to 200.62 MPa, corresponding to a pore throat radius of approximately 3.7 nm. According to the Washburn equation, the relationship between pore throat radius and mercury pressure is:

[0098] (1)

[0099] in, is the mercury capillary pressure (MPa); is the pore throat radius (μm); θ is the wetting angle (140°); is the interfacial tension (0.48 N / m).

[0100] In a constant-rate mercury injection test, mercury was injected into the rock at a constant rate of 5×10⁻⁵ mL / min, resulting in a maximum intrusion pressure of 6.2 MPa and a throat radius of approximately 0.12 μm. The overall upward trend in intrusion pressure exhibited some fluctuations: the pressure decreased as mercury entered the pores through narrow throats and increased as the mercury filled the pores. The pressure fluctuations of the intruding mercury can be used to distinguish between pores and throats.

[0101] In the examples of the present invention, the distribution characteristics of the T2 spectrum are related to the pore throat structure and can directly reflect the pore throat size distribution. Therefore, it is necessary to establish a reasonable formula to convert the transverse relaxation time into the pore throat radius. The pore throat radius obtained by high-pressure mercury injection experiments on the same sample provides a solution to this problem. According to the basic theory of nuclear magnetic resonance, the transverse relaxation time is expressed as:

[0102] (2)

[0103] in, is the transverse relaxation time, 、 、 are surface relaxation time, diffusion-induced relaxation time, and bulk relaxation time, respectively.

[0104] The above formula (2) can be converted to:

[0105] (3)

[0106] in, is the transverse relaxation time, is the surface relaxation time, is the transverse surface relaxation strength, represents the total surface area of ​​the object, Indicates the volume of an object, represents the pore shape factor (cylindrical pore Fs=2, spherical pore Fs=3), Represents the pore throat radius.

[0107] The above formula (3) can be converted to:

[0108] (4)

[0109] in, Indicates the permeability or flowability of porous media. is the transverse surface relaxation strength, represents the pore shape factor (Fs=2 for cylindrical pores and Fs=3 for spherical pores).

[0110] The relationship between the 2 distribution and the pore throat radius is a power function, which can be obtained:

[0111] (5)

[0112] in, is the transverse relaxation time, represents the pore throat radius, represents the power exponent, is the transverse surface relaxation strength, represents the pore shape factor.

[0113] The above formula (5) can be converted to: (6)

[0114] in, represents the pore throat radius, Indicates the permeability or flowability of porous media. is the transverse relaxation time, Represents the power exponent.

[0115] In this embodiment of the present invention, the mercury injection rate must be low enough to negligibly affect the system pressure. During this period, the injection pressure and injected volume are accurately recorded. The recorded data is then analyzed and processed to obtain a core macropore throat distribution curve.

[0116] In an embodiment of the present invention, on the basis of respectively obtaining the macro- and micro-pore distribution curves of each sample, the data of the nuclear magnetic resonance pore throat size distribution curve reflecting the distribution characteristics of small pore throats is combined with the data of the constant-rate mercury injection pore size distribution curve reflecting the distribution characteristics of large pore throats. Each sample has its own splicing point to obtain the full pore size distribution, and the full pore size distribution curve obtained can more accurately reflect the characteristics of the full pore size pore structure of each sample.

[0117] The present invention provides a full-aperture pore structure characterization method based on multidimensional data driving, the method is as follows:

[0118] Step 1: Select the research target and obtain tight reservoir rock samples;

[0119] Step 2: The tight reservoir rock samples obtained in step 1 are sequentially saturated with brine and centrifuged, and nuclear magnetic resonance experiments are performed simultaneously to obtain T2 spectra under brine saturation;

[0120] Step 3: Conduct high-pressure mercury injection experiments and N2 adsorption experiments on the tight reservoir rock samples obtained in step 1 to obtain high-pressure mercury injection and gas adsorption experimental data;

[0121] Step 4: Using high-pressure mercury injection experiments and N2 adsorption experiments to obtain a full-pore cumulative pore volume frequency distribution curve, the T2 spectrum obtained in step 2 under brine saturation is calibrated and converted into pore size to obtain a pore size distribution curve based on the nuclear magnetic resonance T2 spectrum;

[0122] Step 5: Perform a constant-rate mercury injection experiment on the tight reservoir rock sample obtained in step 1 to obtain a pore size distribution curve based on the constant-rate mercury injection experiment;

[0123] Step 6: Determine the lower limit of the minimum measurable pore diameter based on the constant-rate mercury injection pore throat distribution data. Use the overlapping position of the nuclear magnetic resonance pore throat radius distribution curve and the constant-rate mercury injection pore throat radius distribution curve as the junction point. Combine the data reflecting small pore throats on the right side of the nuclear magnetic resonance junction point with the data reflecting large pore throats on the left side of the constant-rate mercury injection junction point to obtain the full pore diameter distribution range of the tight reservoir.

[0124] In this embodiment, data and information from multiple experiments, including high-pressure mercury injection, constant-rate mercury injection, nuclear magnetic resonance, and gas adsorption, are effectively integrated to form a comprehensive full-aperture pore distribution curve. Unlike previous methods, this embodiment, through a multi-experiment collaborative fitting approach, can effectively reflect reservoir quality characteristics and the distribution of nanopores to macropores.

[0125] The full-aperture pore structure characterization method provided by the embodiment of the present invention combines the advantages of four microscopic experiments to complement the pore size distribution characterized by different experiments. This not only helps to evaluate the potential production capacity of the reservoir, but also provides a solid geological basis for the formulation, adjustment and optimization of oilfield development plans.

[0126] The full-aperture pore structure table device based on multidimensional data driving provided by the present invention is described below. The full-aperture pore structure table device based on multidimensional data driving described below and the full-aperture pore structure table method based on multidimensional data driving described above can be referenced to each other.

[0127] refer to Figure 11 , Figure 11Schematic diagram of the structure of the multi-dimensional data-driven full-aperture pore structure table device provided by the present invention. The device includes an acquisition module 1101, a first processing module 1102, a second processing module 1103, a first drawing module 1104, a second drawing module 1105, a third processing module 1106, a first determination module 1107, and a second determination module 1108.

[0128] The acquisition module 1101 is used to obtain a target tight reservoir rock sample; the first processing module 1102 is used to perform a brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the rock sample saturated with brine and the transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; the second processing module 1103 is used to perform high-pressure mercury injection and nitrogen adsorption on the rock sample saturated with brine to obtain high-pressure mercury injection and gas adsorption experimental data; the first drawing module 1104 is used to draw a full-pore cumulative pore volume frequency distribution curve based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution. ; The second drawing module 1105 is used to calibrate the transverse relaxation time under brine saturation based on the full-pore cumulative pore volume frequency distribution curve, and obtain the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; the third processing module 1106 is used to perform constant-rate mercury injection processing on the target tight reservoir rock sample to obtain the pore size distribution curve based on constant-rate mercury injection; the first determination module 1107 is used to determine the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection, and use the overlapping position as the combination point; the second determination module 1108 is used to determine the full-pore size distribution range of the target tight reservoir based on the combination point.

[0129] Specifically, the above-mentioned full-aperture pore structure characterization device based on multi-dimensional data driving provided by the present invention can implement all the method steps implemented by the above-mentioned full-aperture pore structure characterization method embodiment based on multi-dimensional data driving, and can achieve the same technical effect. The parts and beneficial effects of this embodiment that are the same as those of the method embodiment will not be described in detail here.

[0130] Figure 12 This is a schematic diagram of the physical structure of the electronic device provided by the present invention, such as Figure 12As shown, the electronic device may include: a processor 1210, a communications interface 1220, a memory 1230 and a communications bus 1240, wherein the processor 1210, the communications interface 1220 and the memory 1230 communicate with each other via the communications bus 1240. The processor 1210 may call the logic instructions in the memory 1230 to execute a full-aperture pore structure characterization method driven by multi-dimensional data, the method comprising: obtaining a target tight reservoir rock sample; performing a brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain a rock sample saturated with brine and a transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; performing high-pressure mercury injection and nitrogen adsorption on the rock sample saturated with brine to obtain high-pressure mercury injection and gas adsorption experimental data; based on the high-pressure mercury injection and gas adsorption experimental data According to the pore throat size distribution, a full-pore cumulative pore volume frequency distribution curve is drawn; based on the full-pore cumulative pore volume frequency distribution curve, the transverse relaxation time under brine saturation is calibrated to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; the target tight reservoir rock sample is subjected to constant-rate mercury injection treatment to obtain a pore size distribution curve based on constant-rate mercury injection; the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection is determined, and the overlapping position is used as the combining point; based on the combining point, the full-pore size distribution range of the target tight reservoir is determined.

[0131] Furthermore, the logic instructions in the aforementioned memory 1230 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product, stored in a storage medium, includes instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a mobile hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0132] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the full-aperture pore structure characterization method based on multidimensional data driving provided by the above methods, the method comprising: obtaining a target tight reservoir rock sample; performing brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the rock sample after brine saturation and the transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; performing high-pressure mercury injection and nitrogen adsorption on the rock sample after brine saturation Attached, high-pressure mercury injection and gas adsorption experimental data are obtained; based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution, the full-pore cumulative pore volume frequency distribution curve is drawn; based on the full-pore cumulative pore volume frequency distribution curve, the transverse relaxation time under brine saturation is calibrated to obtain the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; the target tight reservoir rock sample is subjected to constant-rate mercury injection treatment to obtain the pore size distribution curve based on constant-rate mercury injection; the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection is determined, and the overlapping position is used as the combination point; based on the combination point, the full-pore size distribution range of the target tight reservoir is determined.

[0133] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the full-aperture pore structure characterization method based on multidimensional data drive provided by the above methods, the method comprising: obtaining a target tight reservoir rock sample; performing a brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain a brine-saturated rock sample and a transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; performing high-pressure mercury injection and nitrogen adsorption on the brine-saturated rock sample to obtain a high-pressure mercury injection and gas adsorption experiment. experimental data; based on high-pressure mercury injection and gas adsorption experimental data and pore throat size distribution, a full-pore cumulative pore volume frequency distribution curve is drawn; based on the full-pore cumulative pore volume frequency distribution curve, the transverse relaxation time under brine saturation is calibrated to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; constant-rate mercury injection is performed on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection; the overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant-rate mercury injection is determined, and the overlapping position is used as the joining point; based on the joining point, the full-pore size distribution range of the target tight reservoir is determined.

[0134] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they may be located in one location or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of the present embodiment. Persons of ordinary skill in the art will be able to understand and implement the present invention without inventive effort.

[0135] Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented using software plus a necessary general-purpose hardware platform, or of course, hardware. Based on this understanding, the essence of the above technical solution, or the portion that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, or an optical disk, and includes a number of instructions for causing a computer device (such as a personal computer, server, or network device) to execute the methods described in each embodiment or certain portions of the embodiments.

[0136] 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 various embodiments of the present invention.

Claims

1. A full-aperture pore structure characterization method based on multidimensional data driving, characterized in that: include: Obtain target tight reservoir rock samples; performing a brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; Performing high-pressure mercury injection and nitrogen adsorption on the rock sample saturated with brine to obtain high-pressure mercury injection and gas adsorption experimental data; Based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution, a full-pore-diameter cumulative pore volume frequency distribution curve is drawn; Based on the full-pore cumulative pore volume frequency distribution curve, the transverse relaxation time under brine saturation is calibrated to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; Performing a constant-rate mercury injection process on the target tight reservoir rock sample to obtain a pore size distribution curve based on the constant-rate mercury injection process; Determining an overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant rate mercury injection, and using the overlapping position as a joining point; Based on the combination points, the full pore size distribution range of the target tight reservoir is determined.

2. The full-aperture pore structure characterization method based on multidimensional data driving according to claim 1, characterized in that: After performing brine saturation treatment and first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, the method further includes: The rock sample saturated with brine is centrifuged and subjected to a second nuclear magnetic resonance measurement to obtain a transverse relaxation time under each centrifugal force, wherein the centrifugation and the second nuclear magnetic resonance measurement are performed simultaneously.

3. The full-aperture pore structure characterization method based on multidimensional data driving according to claim 1, characterized in that: Determining the full pore size distribution range of the target tight reservoir based on the combined points includes: Determine the small pore throat data on the left side of the junction and the large pore throat data on the right side of the junction; The small pore throat data and the large pore throat data are used as the full pore size distribution range of the target tight reservoir.

4. The full-aperture pore structure characterization method based on multidimensional data driving according to claim 1, characterized in that: Before obtaining the target tight reservoir rock sample, the method further includes: Preprocessing the first tight reservoir rock sample and the second tight reservoir rock sample to obtain a standard cylindrical core, wherein the diameter of the standard cylindrical core is a target diameter and the length of the standard cylindrical core is within a target length range; The standard cylindrical core is cleaned and dried to obtain a target tight reservoir rock sample, wherein the constant temperature of the target tight reservoir rock sample is a target temperature.

5. The full-aperture pore structure characterization method based on multidimensional data driving according to claim 1, characterized in that: After obtaining a full-pore-diameter cumulative pore volume frequency distribution curve based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution, the method further includes: Determine the mercury saturation corresponding to the maximum mercury saturation and displacement pressure; The pore throat radius at the mercury saturation and the transverse relaxation time corresponding to the pore throat radius are determined.

6. The full-aperture pore structure characterization method based on multidimensional data driving according to claim 1, characterized in that: After performing constant-rate mercury injection on the target tight reservoir rock sample to obtain a pore size distribution curve based on constant-rate mercury injection, the method further includes: Based on the pore size distribution curve based on constant-rate mercury intrusion, the lower limit of the minimum measurable space is determined.

7. A full-aperture pore structure characterization device based on multidimensional data drive, characterized in that: include: An acquisition module is used to obtain target tight reservoir rock samples; a first processing module, configured to perform a brine saturation treatment and a first nuclear magnetic resonance measurement on the target tight reservoir rock sample to obtain the brine-saturated rock sample and the transverse relaxation time under brine saturation, wherein the brine saturation treatment and the first nuclear magnetic resonance measurement are performed simultaneously; The second processing module is used to perform high-pressure mercury injection and nitrogen adsorption on the rock sample saturated with brine to obtain high-pressure mercury injection and gas adsorption experimental data; A first drawing module is used to draw a full pore size cumulative pore volume frequency distribution curve based on the high-pressure mercury injection and gas adsorption experimental data and the pore throat size distribution; A second drawing module is used to calibrate the transverse relaxation time under brine saturation based on the full-pore cumulative pore volume frequency distribution curve to obtain a transverse relaxation time pore size distribution curve based on nuclear magnetic resonance; A third processing module is used to perform a constant-rate mercury injection process on the target tight reservoir rock sample to obtain a pore size distribution curve based on the constant-rate mercury injection process; A first determining module is configured to determine an overlapping position of the transverse relaxation time pore size distribution curve based on nuclear magnetic resonance and the pore size distribution curve based on constant rate mercury injection, and use the overlapping position as a joining point; The second determining module is configured to determine the full pore size distribution range of the target tight reservoir based on the combination point.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the full-aperture pore structure characterization method based on multidimensional data driving as described in any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the full-aperture pore structure characterization method based on multidimensional data driving as claimed in any one of claims 1 to 6 is implemented.

10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the full-aperture pore structure characterization method based on multidimensional data driving as claimed in any one of claims 1 to 6 is implemented.

Images