A shale organic-inorganic pore volume determination method, device and storage medium

Abstract

The application discloses a shale organic-inorganic pore volume determination method and device and a storage medium, and relates to the technical field of shale exploitation. The method comprises the following steps: obtaining a nitrogen adsorption test data set and a water vapor adsorption test data set of a target shale; processing the nitrogen adsorption test data set and the water vapor adsorption test data set by using a nitrogen adsorption pore calculation model and a water vapor adsorption pore calculation model, to obtain a first pore radius data set, a first pore volume data set, a second pore radius data set, a second pore volume data set and a target pore radius interval set; then determining the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval; and finally determining the organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set. The determination result of the method can comprehensively reflect the pore characteristics of the sample.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of shale mining, in particular to a shale organic-inorganic pore volume determination method, device and storage medium. BACKGROUND

[0002] Shale reservoirs mainly develop micro-nanopores, and their pore structure parameters are important factors affecting shale gas occurrence and percolation. In addition, shale pores are mainly composed of organic pores and inorganic pores. It is generally believed that organic pores exhibit oil-wetness, while inorganic pores exhibit water-wetness. Therefore, the occurrence mode and percolation behavior of shale gas in organic pores and inorganic pores are different, and it is necessary to quantitatively characterize organic pores and inorganic pores.

[0003] Related technologies use image methods to characterize shale pore volume, that is, scanning electron microscopy and other means are used to obtain pore images of shale, then organic pores and inorganic pores in the images are identified, and finally statistical analysis is carried out. For example, Chinese patent (CN112414917A) discloses "a shale oil reservoir organic pore and inorganic pore division and characterization method", which uses a scanning electron microscope combined with an image segmentation algorithm to identify organic and inorganic pores.

[0004] However, the research field of this image method is small, and it cannot fully reflect the pore characteristics of the sample, so the results are not strongly representative. Moreover, the method requires coring, polishing, and continuous cutting scanning of the target shale oil reservoir section for not less than 700 times during the preparation process, and the test procedure is complex and the test time is long. SUMMARY

[0005] The shale organic-inorganic pore volume determination method provided by the embodiments of the present application solves the problems of weak representativeness and long test time of shale organic-inorganic pore volume results determined by related technologies. The technical solutions are as follows:

[0006] In a first aspect, a shale organic-inorganic pore volume determination method is provided, which includes:

[0007] Obtaining a nitrogen adsorption test data set and a water vapor adsorption test data set of a target shale, the nitrogen adsorption test data set being data of nitrogen adsorption amount corresponding to N relative pressures at 77 degrees Kelvin, and the water vapor adsorption test data set being data of water vapor adsorption amount corresponding to S relative pressures at 298 degrees Kelvin, N and S being positive integers greater than 1;

[0008] Using a nitrogen adsorption pore calculation model to process the nitrogen adsorption test data set to obtain a first pore radius data set and a first pore volume data set;

[0009] processing the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set;

[0010] determining a target pore radius interval set according to the first pore radius data set and the second pore radius data set;

[0011] determining a cumulative nitrogen pore volume and a cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set;

[0012] determining an organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set.

[0013] Optionally, the processing the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set comprises:

[0014] processing the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set,

[0015]

[0016] Δv gn (p ri )=v gn (p ri+1 )-v gn (p ri )

[0017] Δt n (p ri )=t n (p ri+1 )-t n (p ri )

[0018] Δv ln (p ri )=0.0015468Δv gn (p ri )

[0019]

[0020] wherein, the p rifor the i-th relative pressure in the nitrogen adsorption test data set, unit is dimensionless, i = 1, 2, 3, …, N; the t n (p ri ) is the adsorption layer thickness when the relative pressure is p ri , unit is nm; the r kn (p ri ) is the corresponding curvature radius when the relative pressure is p ri , to the average pore radius, unit is nm; the v gn (p ri ) is the corresponding nitrogen adsorption amount when the relative pressure is p ri , unit is cm 3 / g; the Δv gn (p ri ) is the nitrogen desorption amount after the relative pressure decreases from p ri+1 to p ri , unit is cm 3 / g; the Δt n (p ri ) is the change amount of the adsorption layer thickness after the relative pressure decreases from p ri+1 to p ri , unit is nm; the Δv ln (p ri ) is the liquid nitrogen volume desorbed after the relative pressure decreases from p ri+1 to p ri , unit is cm 3 / g; the V is the volume of the pores with the pore radius of r , unit is cm 3 / g.

[0021] Optionally, the water vapor adsorption pore calculation model is used to process the water vapor adsorption test data set to obtain a second pore radius data set and a second pore volume data set, including:

[0022] The water vapor adsorption test data set is processed to obtain a second pore radius data set and a second pore volume data set by the following formula,

[0023] t h (p ri ) = 0.395-0.189xln(-ln(p ri ))

[0024]

[0025] Δv gh (p ri ) = vgh (p ri+1 )-v gh (p ri )

[0026] Δt h (p ri )=t h (p ri+1 )-t h (p ri )

[0027] Δv lh (p ri ) = 8.0357 × 10 -4 Δv gh (p ri )

[0028]

[0029] Wherein, p ri The t represents the i-th relative pressure among the S relative pressures in the water vapor adsorption test dataset, in dimensionless units, i = 1, 2, 3, ..., S; h (p ri When the relative pressure is p ri The thickness of the water adsorption layer at that time, in nm; the r kh (p ri When the relative pressure is p ri The radius of curvature of the water at that time, in nm; When the relative pressure is (p) ri +p ri+1 The radius of curvature of water at ) / 2, in nm; When the relative pressure is (p) ri +p ri+1 The average pore radius detected in water at ) / 2 is expressed in nm; the v gh (p ri The relative pressure is p. ri The amount of water adsorbed under standard conditions, in cm³. 3 / g; the Δv gh (p ri () represents the relative pressure from p ri+1 Drop to p ri The amount of water adsorbed under standard conditions after desorption, expressed in cm³. 3 / g; the Δt n (p ri () represents the relative pressure from p ri+1 Drop to p ri The change in the thickness of the post-water adsorption layer, in nm; the Δvlh (p ri ) is the relative pressure from p ri+1 to p ri , the volume of water desorbed after p 3 , in cm / g; the volume of pores with a radius of r , in cm 3 / g.

[0030] Optionally, the determining the target pore radius interval set according to the first pore radius data set and the second pore radius data set comprises:

[0031] obtaining a maximum value of the pore radius in the first pore radius data set and a maximum value of the pore radius in the second pore radius data set;

[0032] determining the pore radius with a larger value as the upper limit of the pore radius interval from the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set;

[0033] determining the value 0 as the lower limit of the pore radius interval;

[0034] dividing the pore radius interval according to a preset interval length to obtain a plurality of target pore radius intervals;

[0035] determining the plurality of target pore radius intervals as the target pore radius interval set.

[0036] Optionally, the determining the organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval in the target pore radius interval set comprises:

[0037] letting r = 1, comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set;

[0038] if the cumulative nitrogen pore volume of the rth target pore radius interval is greater than the cumulative water vapor pore volume, then determining the difference between the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval as the organic pore volume of the rth target pore radius interval, and determining the cumulative water vapor pore volume of the rth target pore radius interval as the inorganic pore volume of the rth target pore radius interval;

[0039] if the cumulative nitrogen pore volume of the rth target pore radius interval is less than the cumulative water vapor pore volume, then the cumulative water vapor pore volume of the rth target pore radius interval is determined as the inorganic pore volume of the rth target pore radius interval, and the organic pore volume of the rth target pore radius interval is determined as 0;

[0040] determining whether r is equal to M, where M is the number of target pore radius intervals in the target pore radius interval set;

[0041] if r is not equal to M, then r = r + 1, and returning to the step of comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set;

[0042] if r is equal to M, then the organic pore volume and the inorganic pore volume of the M target pore radius intervals are added to obtain the total organic pore volume and the total inorganic pore volume of the target shale.

[0043] Optionally, before the step of obtaining the nitrogen adsorption test data set and the water vapor adsorption test data set of the target shale, the method further comprises:

[0044] obtaining the original nitrogen adsorption relationship and the original water vapor adsorption relationship of the target shale, where the original nitrogen adsorption relationship refers to the linear relationship between the relative pressure and the nitrogen adsorption amount at 77 K, and the original water vapor adsorption relationship refers to the linear relationship between the relative pressure and the water vapor adsorption amount at 298 K;

[0045] performing interpolation processing on the original nitrogen adsorption relationship by using the spline interpolation method to obtain the data of the nitrogen adsorption amount corresponding to N relative pressures with the same interval, and determining the N relative pressures and the data of the nitrogen adsorption amount corresponding to the N relative pressures as the nitrogen adsorption test data set of the target shale;

[0046] performing interpolation processing on the original water vapor adsorption relationship by using the spline interpolation method to obtain the data of the water vapor adsorption amount corresponding to S relative pressures with the same interval, and determining the S relative pressures and the data of the water vapor adsorption amount corresponding to the S relative pressures as the water vapor adsorption test data set of the target shale.

[0047] In a second aspect, a shale organic-inorganic pore volume determination device is provided, and the device comprises:

[0048] The acquisition module is configured to acquire a nitrogen adsorption test data set and a water vapor adsorption test data set of a target shale, the nitrogen adsorption test data set refers to data of nitrogen adsorption amounts corresponding to N relative pressures at 77 degrees Kelvin, and the water vapor adsorption test data set refers to data of water vapor adsorption amounts corresponding to S relative pressures at 298 degrees Kelvin, N and S are positive integers greater than 1;

[0049] The first processing module is configured to process the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set.

[0050] The second processing module is configured to process the water vapor adsorption test data set by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set.

[0051] The first determining module is configured to determine a target pore radius interval set according to the first pore radius data set and the second pore radius data set.

[0052] The second determining module is configured to determine cumulative nitrogen pore volumes and cumulative water vapor pore volumes of each target pore radius interval of the target pore radius interval set according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set, and the second pore volume data set.

[0053] The third determining module is configured to determine an organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volumes and the cumulative water vapor pore volumes of each target pore radius interval of the target pore radius interval set.

[0054] Optionally, the first processing module is specifically configured to:

[0055] The nitrogen adsorption test data set is processed by the following formula to obtain the first pore radius data set and the first pore volume data set,

[0056]

[0057] Δv gn (p ri )=v gn (p ri+1 )-v gn (p ri )

[0058] Δt n (p ri )=t n (p ri+1 )-t n (p ri)

[0059] Δv ln (p ri )=0.0015468Δv gn (p ri )

[0060]

[0061] Wherein, p ri The t represents the i-th relative pressure among N relative pressures in the nitrogen adsorption test dataset, in dimensionless units, where i = 1, 2, 3, ..., N; n (p ri When the relative pressure is p ri The adsorption layer thickness at that time is expressed in nm; the r kn (p ri When the relative pressure is p ri The corresponding radius of curvature, in nm; When the relative pressure is (p) ri +p ri+1 The radius of curvature corresponding to ) / 2 is expressed in nm; When the relative pressure is (p) ri +p ri+1 The average pore radius detected by nitrogen at ) / 2, in nm; the v gn (p ri The relative pressure is p. ri The corresponding nitrogen adsorption amount at that time, in cm³. 3 / g; the Δv gn (p ri () represents the relative pressure from p ri+1 Drop to p ri The amount of nitrogen adsorbed after desorption is expressed in cm³. 3 / g; the Δt n (p ri () represents the relative pressure from p ri+1 Drop to p ri The change in the thickness of the post-adsorption layer, in nm; the Δv ln (p ri () represents the relative pressure from p ri+1 Drop to p ri The volume of liquid nitrogen subsequently desorbed is expressed in cm³. 3 / g; the aforementioned The pore radius detected by nitrogen gas is The volume of the pores, in cm³. 3 / g.

[0062] Optionally, the second processing module is specifically used for:

[0063] The water vapor adsorption test data set is processed by the following formula to obtain a second pore radius data set and a second pore volume data set,

[0064] t h (p ri ) = 0.395 - 0.189 x ln(-ln(p ri ))

[0065]

[0066] Δv gh (p ri ) = v gh (p ri+1 ) - v gh (p ri )

[0067] Δt h (p ri ) = t h (p ri+1 ) - t h (p ri )

[0068] Δv lh (p ri ) = 8.0357 x 10 -4 Δv gh (p ri )

[0069]

[0070] wherein, the p ri is the i-th relative pressure of S relative pressures in the water vapor adsorption test data set, unit is dimensionless, i = 1, 2, 3, …, S; the t h (p ri ) is the adsorbed layer thickness of water when the relative pressure is p ri , unit is nm; the r kh (p ri ) is the corresponding curvature radius of water when the relative pressure is p ri , unit is nm; the is the corresponding curvature radius of water when the relative pressure is (p ri +p ri+1 ) / 2, unit is nm; the is the average pore radius detected by water when the relative pressure is (p ri +p ri+1 ) / 2, unit is nm; the v gh (pri ) is the adsorbed amount of water in the standard state when the relative pressure is p ri ; the Δv 3 (p gh ) is the adsorbed amount of water in the standard state after desorption when the relative pressure decreases from p ri to p ri+1 ; the Δv ri (p 3 ) is the change in the thickness of the water adsorption layer after the relative pressure decreases from p h to p ri ; the Δv ri+1 (p ri ) is the volume of water desorbed after the relative pressure decreases from p lh to p ri ; the Δv ri+1 (p ri ) is the volume of water desorbed after the relative pressure decreases from p 3 to p ; the V is the volume of the pores with a pore radius of r 3 probed by water, in units of cm

[0071] Optionally, the first determining module comprises:

[0072] an acquisition unit, configured to acquire a maximum value of the pore radius in the first pore radius data set and a maximum value of the pore radius in the second pore radius data set;

[0073] a first determining unit, configured to determine the pore radius with a larger value between the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set as the upper limit of the pore radius interval;

[0074] a second determining unit, configured to determine the value 0 as the lower limit of the pore radius interval;

[0075] a division unit, configured to divide the pore radius interval according to a preset interval length to obtain a plurality of target pore radius intervals;

[0076] a third determining unit, configured to determine the plurality of target pore radius intervals as the target pore radius interval set.

[0077] Optionally, the third determining module comprises:

[0078] a comparison unit, configured to let r = 1, and compare the sizes of the cumulative nitrogen pore volume and the cumulative water vapor pore volume in the rth target pore radius interval in the target pore radius interval set;

[0079] the fourth determining unit is configured to determine, if the cumulative nitrogen pore volume of the rth target pore radius interval is greater than the cumulative water vapor pore volume, a difference between the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval as an organic pore volume of the rth target pore radius interval, and determine the cumulative water vapor pore volume of the rth target pore radius interval as an inorganic pore volume of the rth target pore radius interval;

[0080] the fifth determining unit is configured to determine, if the cumulative nitrogen pore volume of the rth target pore radius interval is less than the cumulative water vapor pore volume, the cumulative water vapor pore volume of the rth target pore radius interval as the inorganic pore volume of the rth target pore radius interval, and determine the organic pore volume of the rth target pore radius interval as 0;

[0081] the first judging unit is configured to judge whether r is equal to M, where M is a number of target pore radius intervals in the target pore radius interval set;

[0082] the second judging unit is configured to, if r is not equal to M, set r = r + 1, and return to the step of comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set;

[0083] the accumulating unit is configured to, if r is equal to M, accumulate the organic pore volumes and the inorganic pore volumes of the M target pore radius intervals respectively to obtain a total organic pore volume and a total inorganic pore volume of the target shale.

[0084] Optionally, the shale organic-inorganic pore volume determination apparatus further comprises:

[0085] the second obtaining module is configured to obtain an original nitrogen adsorption relationship and an original water vapor adsorption relationship of the target shale, where the original nitrogen adsorption relationship refers to a linear relationship between a relative pressure and a nitrogen adsorption amount at 77 degrees Kelvin, and the original water vapor adsorption relationship refers to a linear relationship between a relative pressure and a water vapor adsorption amount at 298 degrees Kelvin;

[0086] the third processing module is configured to perform interpolation processing on the original nitrogen adsorption relationship by using a spline interpolation method to obtain data of nitrogen adsorption amounts corresponding to N relative pressures with the same interval, and determine the N relative pressures and the data of nitrogen adsorption amounts corresponding to the N relative pressures as a nitrogen adsorption test data set of the target shale;

[0087] The fourth processing module is configured to perform spline interpolation on the original water vapor adsorption relationship to obtain data of water vapor adsorption amounts corresponding to S relative pressures with equal intervals, and determine the S relative pressures and the data of water vapor adsorption amounts corresponding to the S relative pressures as the water vapor adsorption test data set of the target shale.

[0088] In a third aspect, a shale organic-inorganic pore volume determination device is provided, and the device comprises:

[0089] a processor and a memory for storing processor-executable instructions;

[0090] The processor is configured to execute any of the methods of the first aspect.

[0091] In a fourth aspect, a computer-readable storage medium is provided, and the storage medium stores a computer program which, when executed by a processor, implements any of the methods of the first aspect.

[0092] The technical solutions provided by the embodiments of the present application can bring at least the following beneficial effects:

[0093] In the embodiments of the present application, the nitrogen adsorption test data set and the water vapor adsorption test data set of the target shale can be obtained, the nitrogen adsorption test data set is processed by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set, the water vapor adsorption test data set is processed by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set, the target pore radius interval set is determined according to the first pore radius data set and the second pore radius data set, the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set are determined according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set, and the organic-inorganic pore volume of the target shale is determined according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set. That is, after the nitrogen adsorption test and the water vapor adsorption test of the target shale gas are performed, the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the corresponding interval of different pore radii are obtained by processing and calculating the nitrogen adsorption test and the water vapor adsorption test, and the organic-inorganic pore volume of the target shale is determined by comparing the nitrogen pore volume and the cumulative water vapor pore volume, and the determination result can fully reflect the pore characteristics of the sample. Moreover, the determination result can be directly output by operation after the two adsorption tests, which simplifies the test procedure and shortens the test time. BRIEF DESCRIPTION OF DRAWINGS

[0094] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiments will be briefly introduced. Obviously, the drawings in the following description only constitute some embodiments of the present application, and for those skilled in the art, other drawings can be obtained without creative labor based on these drawings.

[0095] Figure 1 is a flowchart of a shale organic-inorganic pore volume determination method provided by an embodiment of the present application;

[0096] Figure 2 is a flowchart of another shale organic-inorganic pore volume determination method provided by an embodiment of the present application;

[0097] Figure 3 is a possible target shale organic-inorganic pore distribution diagram shown by an embodiment of the present application;

[0098] Figure 4 is a structural diagram of a shale organic-inorganic pore volume determination device provided by an embodiment of the present application;

[0099] Figure 5 is a structural diagram of a terminal 500 provided by an embodiment of the present application. DETAILED DESCRIPTION

[0100] In order to make the purpose, technical solutions and advantages of the present application more clear, the embodiments of the present application will be further described in detail below with reference to the drawings.

[0101] Before the embodiments of the present application are explained in detail, the terms, application scenarios and system architecture involved in the embodiments of the present application will be explained.

[0102] First, the terms involved in the embodiments of the present application will be introduced.

[0103] Nitrogen adsorption test

[0104] Nitrogen adsorption test refers to testing the nitrogen adsorption amount of a target shale sample by changing the relative pressure under the condition of nitrogen as the adsorbate at 77 degrees Kelvin, and the obtained result is a linear relative pressure-adsorption amount curve.

[0105] Water vapour adsorption test

[0106] Water vapor adsorption test refers to testing the water vapor adsorption amount of a target shale sample by changing the relative pressure under the condition of water vapor as the adsorbate at 298 degrees Kelvin, and the obtained result is a linear relative pressure-adsorption amount curve.

[0107] Organic-inorganic pore volume of target shales

[0108] The organic-inorganic pore volume of the target shale refers to the total of the organic pore volume and the inorganic pore volume of the target shale.

[0109] Secondly, the application scenarios related to the embodiments of the present application are introduced.

[0110] Before the development of the shale reservoir, the pore structure of the shale reservoir needs to be described, and the pore volume of the organic pores and the inorganic pores is characterized to understand the shale gas occurrence and seepage conditions and to develop a development plan. In this case, the shale organic-inorganic pore volume determination method provided by the embodiments of the present application can quantitatively characterize the shale organic-inorganic pores in a numerical driving manner to facilitate the engineering personnel to develop a shale reservoir development plan according to the characterization data.

[0111] Finally, the system architecture related to the embodiments of the present application is introduced.

[0112] The shale organic-inorganic pore volume determination method provided by the embodiments of the present application can be applied to a terminal with a data processing function. Specifically, the terminal can be a smart phone, a tablet computer, a notebook computer, a desktop computer or other terminals capable of data processing.

[0113] Figure 1 is a flowchart of a shale organic-inorganic pore volume determination method provided by the embodiments of the present application. Referring to Figure 1 The method comprises the following steps:

[0114] Step 101: Obtain a nitrogen adsorption test data set and a water vapor adsorption test data set of the target shale, the nitrogen adsorption test data set refers to the data of nitrogen adsorption amount corresponding to N relative pressures at 77 degrees Kelvin, and the water vapor adsorption test data set refers to the data of water vapor adsorption amount corresponding to S relative pressures at 298 degrees Kelvin, N and S are both positive integers greater than 1.

[0115] Step 102: Process the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set.

[0116] Step 103: Process the water vapor adsorption test data set by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set.

[0117] Step 104: Determine a target pore radius interval set according to the first pore radius data set and the second pore radius data set.

[0118] Step 105: determining the cumulative nitrogen gas pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set.

[0119] Step 106: determining the organic-inorganic pore volume of the target shale according to the cumulative nitrogen gas pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set.

[0120] In the embodiment of the present application, the nitrogen gas adsorption test data set and the water vapor adsorption test data set of the target shale can be obtained, the nitrogen gas adsorption test data set is processed by using the nitrogen gas adsorption pore calculation model to obtain the first pore radius data set and the first pore volume data set, the water vapor adsorption test data set is processed by using the water vapor adsorption pore calculation model to obtain the second pore radius data set and the second pore volume data set, the target pore radius interval set is determined according to the first pore radius data set and the second pore radius data set, the cumulative nitrogen gas pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set are determined according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set, and the organic-inorganic pore volume of the target shale is determined according to the cumulative nitrogen gas pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set. That is, after the nitrogen gas adsorption test and the water vapor adsorption test of the target shale gas are performed once, the cumulative nitrogen gas pore volume and the cumulative water vapor pore volume of the corresponding interval of different pore radii are obtained by processing and calculating the nitrogen gas adsorption test and the water vapor adsorption test, the organic-inorganic pore volume of the target shale is determined by comparing the nitrogen gas pore volume and the cumulative water vapor pore volume, and the determination result can fully reflect the pore characteristics of the sample. Moreover, the determination result can be directly output by operation after the two adsorption tests, the test procedure is simplified, and the test time is shortened.

[0121] Optionally, the nitrogen gas adsorption test data set is processed by using the nitrogen gas adsorption pore calculation model to obtain the first pore radius data set and the first pore volume data set, including:

[0122] The nitrogen gas adsorption test data set is processed by using the following formula to obtain the first pore radius data set and the first pore volume data set,

[0123]

[0124]

[0125] Δvgn (p ri ) = v gn (p ri+1 )-v gn (p ri )

[0126] Δt n (p ri )=t n (p ri+1 )-t n (p ri )

[0127] Δv ln (p ri )=0.0015468Δv gn (p ri )

[0128]

[0129] Where, p ri Let t be the i-th relative pressure among N relative pressures in the nitrogen adsorption test dataset, in dimensionless units, i = 1, 2, 3, ..., N; n (p ri When the relative pressure is p ri The thickness of the adsorption layer at that time, in units When the relative pressure is (p) ri +p ri+1 The average pore radius detected by nitrogen at ) / 2, in nm; gn (p ri The relative pressure is p. ri The corresponding nitrogen adsorption amount at that time, in cm³. 3 / g;Δv gn (p ri () represents the relative pressure from p ri+1 Drop to p ri The amount of nitrogen adsorbed after desorption is expressed in cm³. 3 / g;Δt n (p ri () represents the relative pressure from p ri+1 Drop to p ri The change in the thickness of the post-adsorption layer, in nm; Δv ln (p ri () represents the relative pressure from p ri+1 Drop to p ri The volume of liquid nitrogen subsequently desorbed is expressed in cm³. 3 / g; The pore radius detected by nitrogen gas is The volume of the pores, in cm³. 3 / g.

[0130] Optionally, the water vapor adsorption test data set is processed by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set, comprising:

[0131] The water vapor adsorption test data set is processed by using the following formula to obtain the second pore radius data set and the second pore volume data set,

[0132] t h (p ri )=0.395-0.189×ln(-ln(p ri ))

[0133]

[0134] Δv gn (p ri )=v gh (p ri+1 )-v gh (p ri )

[0135] Δt h (p ri )=t h (p ri+1 )-t h (p ri )

[0136] Δv lh (p ri )=8.0357×10 -4 Δv gh (p ri )

[0137]

[0138] wherein p ri is the i-th relative pressure of S relative pressures in the water vapor adsorption test data set, unit is dimensionless, i = 1, 2, 3, …, S; t h (p ri ) is the adsorbed layer thickness of water when the relative pressure is p ri , unit is nm; r kh (p ri ) is the corresponding curvature radius of water when the relative pressure is p ri , unit is nm; is the corresponding curvature radius of water when the relative pressure is (p ri +p ri+1 ) / 2, unit is nm; is the corresponding curvature radius of water when the relative pressure is (pri + p ri+1 ) is the average pore radius of water detection, unit: nm; v gh (p ri ) is the adsorption amount of water under standard state when the relative pressure is p ri , unit: cm 3 / g; Δv gh (p ri ) is the adsorption amount of water under standard state after the relative pressure decreases from p ri+1 to p ri , unit: cm 3 / g; Δt h (p ri ) is the change amount of water adsorption layer thickness after the relative pressure decreases from p ri+1 to p ri , unit: nm; Δv lh (p ri ) is the volume of water desorption after the relative pressure decreases from p ri+1 to p ri , unit: cm 3 / g; is the volume of the pore with a pore radius of r , unit: cm 3 / g.

[0139] Optionally, according to the first pore radius data set and the second pore radius data set, a target pore radius interval set is determined, including:

[0140] Obtaining the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set;

[0141] Determining the pore radius with a larger value between the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set as the upper limit of the pore radius interval;

[0142] Determining the value 0 as the lower limit of the pore radius interval;

[0143] Dividing the pore radius interval according to a preset interval length to obtain a plurality of target pore radius intervals;

[0144] Determining the plurality of target pore radius intervals as the target pore radius interval set.

[0145] Optionally, according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set, the organic-inorganic pore volume of the target shale is determined, including:

[0146] Let r = 1, compare the size of the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set;

[0147] If the cumulative nitrogen pore volume of the rth target pore radius interval is greater than the cumulative water vapor pore volume, then the difference between the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval is determined as the organic pore volume of the rth target pore radius interval, and the cumulative water vapor pore volume of the rth target pore radius interval is determined as the inorganic pore volume of the rth target pore radius interval;

[0148] If the cumulative nitrogen pore volume of the rth target pore radius interval is less than the cumulative water vapor pore volume, then the cumulative water vapor pore volume of the rth target pore radius interval is determined as the inorganic pore volume of the rth target pore radius interval, and the organic pore volume of the rth target pore radius interval is determined as 0;

[0149] Determine whether r is equal to M, M being the number of target pore radius intervals in the target pore radius interval set;

[0150] If r is not equal to M, then let r = r + 1, and return to the step of comparing the size of the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set;

[0151] If r is equal to M, the organic pore volumes and the inorganic pore volumes of the M target pore radius intervals are added respectively to obtain the total organic pore volume and the total inorganic pore volume of the target shale.

[0152] Optionally, before obtaining the nitrogen adsorption test data set and the water vapor adsorption test data set of the target shale, the method further comprises:

[0153] Obtaining the original nitrogen adsorption relationship and the original water vapor adsorption relationship of the target shale, the original nitrogen adsorption relationship being a linear relationship between the relative pressure and the nitrogen adsorption amount at 77 degrees Kelvin, and the original water vapor adsorption relationship being a linear relationship between the relative pressure and the water vapor adsorption amount at 298 degrees Kelvin;

[0154] Performing interpolation processing on the original nitrogen adsorption relationship by a spline interpolation method to obtain data of nitrogen adsorption amounts corresponding to N relative pressures with the same interval, and determining the N relative pressures and the data of nitrogen adsorption amounts corresponding to the N relative pressures as the nitrogen adsorption test data set of the target shale;

[0155] Performing interpolation processing on the original water vapor adsorption relationship by a spline interpolation method to obtain data of water vapor adsorption amounts corresponding to S relative pressures with the same interval, and determining the S relative pressures and the data of water vapor adsorption amounts corresponding to the S relative pressures as the water vapor adsorption test data set of the target shale.

[0156] All the optional technical solutions described above can be combined to form optional embodiments of the present application, and the embodiments of the present application will not be described one by one.

[0157] Figure 2 is a flowchart of another shale organic-inorganic pore volume determination method provided by the embodiments of the present application. Referring to Figure 2 , the method comprises the following steps:

[0158] It should be noted that before obtaining the original nitrogen adsorption relationship and the original water vapor adsorption relationship of the target shale, the experimental sample of the target shale needs to be prepared, and the experimental sample needs to be tested for nitrogen adsorption and water vapor adsorption. The shale sample can be in powder form, and the sample weight is not limited. Preferably, 5g of 40 / 60 mesh powder of the target shale can be prepared as the experimental sample.

[0159] Step 201: Obtain the original nitrogen adsorption relationship and the original water vapor adsorption relationship of the target shale, the original nitrogen adsorption relationship refers to the linear relationship between the relative pressure and the nitrogen adsorption amount at 77 degrees Kelvin, and the original water vapor adsorption relationship refers to the linear relationship between the relative pressure and the water vapor adsorption amount at 298 degrees Kelvin.

[0160] The original nitrogen adsorption relationship of the target shale is the linear relationship between the relative pressure and the nitrogen adsorption amount obtained after testing the experimental sample of the target shale for nitrogen adsorption under the experimental temperature environment of 77 degrees Kelvin. This relationship can be input by the user, sent by other equipment, or obtained by the adsorption testing system after recording the experimental results. For example, the adsorption testing system records the values of the nitrogen adsorption amount corresponding to the continuous change of the relative pressure during the experiment, and generates the linear relationship between the relative pressure and the nitrogen adsorption amount.

[0161] The original water vapor adsorption relationship of the target shale is the linear relationship between the relative pressure and the water vapor adsorption amount obtained after testing the experimental sample of the target shale for water vapor adsorption under the experimental temperature environment of 298 degrees Kelvin. This relationship can be input by the user, sent by other equipment, or obtained by the adsorption testing system after recording the experimental results. For example, the adsorption testing system records the values of the water vapor adsorption amount corresponding to the continuous change of the relative pressure during the experiment, and generates the linear relationship between the relative pressure and the water vapor adsorption amount.

[0162] It should be noted that the relative pressure is the pressure represented with atmospheric pressure as the reference, and the value range of the relative pressure in the nitrogen adsorption test and the water vapor adsorption test is [0-1].

[0163] Step 202: The original nitrogen adsorption relationship is processed by the spline interpolation method to obtain data of nitrogen adsorption amounts corresponding to N relative pressures with the same interval, and the N relative pressures and the data of nitrogen adsorption amounts corresponding to the N relative pressures are determined as the nitrogen adsorption test data set of the target shale, and N is a positive integer greater than 1.

[0164] It should be noted that the spline interpolation method is a mathematical method for making a smooth curve passing through a series of points by using a variable spline. Since it is prior art, the embodiments of the present application do not make specific elaboration. By the spline interpolation method, the linear relationship in the original nitrogen adsorption relationship can be processed into a smooth curve, and then N relative pressures with the same interval and corresponding nitrogen adsorption amounts are selected from the smooth curve as the nitrogen adsorption test data set of the target shale.

[0165] It should be further noted that the value of N is determined by the interval. The larger the interval, the smaller N is, and the smaller the interval, the larger N is. The interval can be input by the user or sent by other devices. In a possible implementation, the user can set the interval to 0.01, and then 101 relative pressures and corresponding nitrogen adsorption amounts in the range of relative pressure [0-1] can be determined, and the 101 relative pressures and corresponding nitrogen adsorption amounts are determined as the nitrogen adsorption test data set of the target shale.

[0166] In addition, since the nitrogen adsorption amounts in the two range intervals of relative pressure [0-0.05) and (0.95-1] are generally not included in the test results, when selecting the relative pressure, the values in the range interval of [0.05-0.95] can be selected. For example, if the interval is set to 0.01, the nitrogen adsorption amounts under the relative pressures of 0.05:0.06:0.07....0.94:0.95 are determined as the nitrogen adsorption test data set of the target shale.

[0167] Step 203: The original water vapor adsorption relationship is processed by the spline interpolation method to obtain data of water vapor adsorption amounts corresponding to S relative pressures with the same interval, and the S relative pressures and the data of water vapor adsorption amounts corresponding to the S relative pressures are determined as the water vapor adsorption test data set of the target shale, and S is a positive integer greater than 1.

[0168] Specifically, by the spline interpolation method, the linear relationship in the original water vapor adsorption relationship can be processed into a smooth curve, and then S relative pressures with the same interval and corresponding water vapor adsorption amounts are selected from the smooth curve as the water vapor adsorption test data set of the target shale.

[0169] It should be noted that the value of S is determined by the size of the interval, the larger the interval, the smaller S, the smaller the interval, the larger S. In a possible implementation, the interval is set to 0.01, and the water vapor adsorption amounts at relative pressures of 0.05:0.06:0.07....0.94:0.95 are determined as the water vapor adsorption test data set of the target shale.

[0170] Step 204: processing the nitrogen adsorption test data set by using the nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set.

[0171] It should be noted that there are many types of nitrogen adsorption pore calculation models, such as BET (Brunauer, Emmett and Teller) calculation model, BJH (Barret, Joyner and Halenda) calculation model, etc. Preferably, the nitrogen adsorption pore calculation model in the embodiment of the present application is a calculation model suitable for shale organic-inorganic pore volume calculation obtained by improving the calculation method based on the BJH calculation model. Specifically, the nitrogen adsorption test data set is processed by the following formula to obtain the first pore radius data set and the first pore volume data set:

[0172]

[0173] Δv gn (p ri )=v gn (p ri+1 )-v gn (p ri ) (5)

[0174] Δt n (p ri )=t n (p ri+1 )-t n (p ri ) (6)

[0175] Δv ln (p ri )=0.0015468Δv gn (p ri ) (7)

[0176]

[0177] wherein p ri is the i-th relative pressure of N relative pressures in the nitrogen adsorption test data set, unit is dimensionless, i = 1, 2, 3,..., N; t n (p ri ) is the water vapor adsorption amount when the relative pressure is pri The thickness of the adsorption layer at that time is expressed in nm; r kn (p ri When the relative pressure is p ri The radius of curvature corresponding to the time, in nm; When the relative pressure is (p) ri +p ri+1 The radius of curvature corresponding to ) / 2 is expressed in nm. When the relative pressure is (p) ri +p ri+1 The average pore radius detected by nitrogen at ) / 2, in nm; gn (p ri The relative pressure is p. ri The corresponding nitrogen adsorption amount at that time, in cm³. 3 / g;Δv gn (p ri () represents the relative pressure from p ri+1 Drop to p ri The amount of nitrogen adsorbed after desorption is expressed in cm³. 3 / g;Δt n (p ri () represents the relative pressure from p ri+1 Drop to p ri The change in the thickness of the post-adsorption layer, in nm; Δv ln (p ri () represents the relative pressure from p ri+1 Drop to p ri The volume of liquid nitrogen subsequently desorbed is expressed in cm³. 3 / g; The pore radius detected by nitrogen gas is The volume of the pores, in cm³. 3 / g.

[0178] It should be noted that when applying the above formulas to calculate the first pore radius dataset and the first pore volume dataset, the i-th and (i+1)-th relative pressures are sequentially selected from the N relative pressures in the nitrogen adsorption test dataset and input into the above formulas (1) and (5) to obtain the relative pressures as (p ri +p ri+1 The average pore radius detected by nitrogen at ) / 2 and the pore volume corresponding to that pore radius are calculated. After traversing and calculating N relative pressures, N-1 pore radii and their corresponding pore volumes are obtained. The N-1 pore radii are determined as the first pore radius dataset, and the N-1 pore volumes are determined as the first pore volume dataset.

[0179] Step 205: processing the water vapor adsorption test data set by using the water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set.

[0180] It should be noted that based on the description of the nitrogen adsorption pore calculation model, the water vapor adsorption test data set is the same as the nitrogen adsorption pore calculation model, and can be of various types, such as BET (Bmnauer, Emmett and Teller) calculation model, BJH (Barret, Joyner and Halenda) calculation model, etc. Preferably, the water vapor adsorption pore calculation model in the embodiment of the present application is based on the BJH calculation model, and a calculation model suitable for shale organic-inorganic pore volume calculation is obtained by improving the calculation method. Specifically, the water vapor adsorption test data set is processed by the following formula to obtain the second pore radius data set and the second pore volume data set,

[0181] t h (p ri )=0.395-0.189×ln(-ln(p ri )) (9)

[0182]

[0183] Δv gh (p ri )=v gh (p ri+1 )-v gh (p ri ) (13)

[0184] Δt h (p ri )=t h (p ri+1 )-t h (p ri ) (14)

[0185] Δv lh (p ri )=8.0357×10 -4 Δv gh (p ri ) (15)

[0186]

[0187] wherein p ri is the i-th relative pressure in the S relative pressures in the water vapor adsorption test data set, unit is dimensionless, i = 1, 2, 3, …, S; t h (p ri ) is the relative pressure pri the adsorbed layer thickness of water at the time, in nm; r kh ri ) is the corresponding radius of curvature of water at the time, in nm; ri ri ri+1 ri ri+1 gh ri ri 3 gh ri ri+1 ri 3 h ri ri+1 ri lh ri ri+1 ri 3 3

[0188] It should be noted that when calculating the second pore radius dataset and the second pore volume dataset using the above formula, the i-th and the i+1-th relative pressures are sequentially selected from the S relative pressures in the water vapor adsorption test dataset and input into the above formula (9) and (13) to obtain the average pore radius detected by the water vapor when the relative pressure is (p ri +p ri+1 ) / 2 and the pore volume corresponding to the pore radius. After traversing the calculation of the S relative pressures, S-1 pore radii and corresponding pore volumes are obtained, the S-1 pore radii are determined as the second pore radius dataset, and the S-1 pore volumes are determined as the second pore volume dataset.

[0189] ​​​​​​​​​​​​​​​​​​​​​​​​​​​​​Step 206: determining a target pore radius interval set according to the first pore radius data set and the second pore radius data set.

[0190] It should be noted that after step 205 is completed, the pore radius interval needs to be divided to obtain the target pore radius interval set, so as to improve the accuracy of determining the target shale organic-inorganic pore volume. Specifically, the target pore radius interval set can be determined through steps 2061-2065.

[0191] Step 2061: obtaining the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set.

[0192] It should be noted that when obtaining the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set, a traversal method or a dichotomy method can be used, and the present application does not make a specific limitation thereon, as long as the speed of obtaining is not affected. For example, the pore radius in the first pore radius data set includes (1, 5, 15, 20), and the pore radius in the second pore radius data set includes (1, 5, 14, 19), and the unit is mm. By using the traversal method, the maximum value of the pore radius in the first pore radius data set is 20 mm, and the maximum value of the pore radius in the second pore radius data set is 19 mm.

[0193] Step 2062: determining the pore radius with the larger value between the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set as the upper limit of the pore radius interval.

[0194] It should be noted that the maximum value of the pore radius in the first pore radius data set and the maximum value of the pore radius in the second pore radius data set are compared, and the pore radius with the larger value is determined as the upper limit of the pore radius interval. For example, the maximum value of the pore radius in the first pore radius data set is 20 mm, and the maximum value of the pore radius in the second pore radius data set is 19 mm, and 20 mm is taken as the upper limit of the pore radius interval.

[0195] Step 2063: determining 0 as the lower limit of the pore radius interval.

[0196] That is, the lower limit of the pore radius interval is 0, and the upper limit is determined according to step 2062. After the upper limit and the lower limit of the pore radius interval are determined, the pore radius interval can be determined. For example, the upper limit of the pore radius interval determined in step 2062 is 20 mm, and the lower limit is 0, and the pore radius interval is [0-20 mm].

[0197] Step 2064: dividing the pore radius interval according to a preset interval length to obtain a plurality of target pore radius intervals.

[0198] It should be noted that the preset interval length can be obtained by user input, and can be obtained by other devices. For example, the pore radius interval is [0-20mm], the user can input the preset interval length as 2mm, and then the pore radius interval is divided into 10 target pore radius intervals according to 2mm, which are [0-2mm], [2-4mm], [4-6mm], [6-8mm], [8-10mm], [10-12mm], [12-14mm], [14-16mm], [16-18mm], and [18-20mm].

[0199] Step 2065: Determine the plurality of target pore radius intervals as a target pore radius interval set.

[0200] That is, the plurality of target pore radius intervals divided in step 2064 are further determined as a target pore radius interval set. For example, the target pore radius intervals [0-2mm], [2-4mm], [4-6mm], [6-8mm], [8-10mm], [10-12mm], [12-14mm], [14-16mm], [16-18mm], and [18-20mm] are determined as the target pore radius interval set.

[0201] Step 207: According to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set, and the second pore volume data set, determine the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set.

[0202] It should be noted that the pore radii in the first pore radius data set and the second pore radius data set are distributed in the target pore radius interval set, so for each target pore radius interval, the sum of the first pore volume data corresponding to the first pore radius data distributed in the pore radius interval, and the sum of the second pore volume data corresponding to the second pore radius data distributed in the pore radius interval, can be determined. For example, for a target pore radius interval [10-12mm], there are three first pore radii 10.5mm, 10.8mm, and 11.5mm in the first pore radius data set, and the first pore radius volume corresponding to the three first pore radii is 0.08cm 3 / g, 0.07cm 3 / g, 0.075cm 3 / g, then the cumulative nitrogen pore volume of the target pore radius interval is 0.225cm 3 / g.

[0203] Step 208: determining the organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set.

[0204] It should be noted that for shale, the proportion of inorganic minerals is much larger than that of organic matter. Therefore, the pore volume of inorganic pores is much larger than that of organic pores. Organic pores are oil-wet, while inorganic pores are water-wet. Therefore, water molecules can only detect inorganic pores, while nitrogen molecules can detect both organic and inorganic pores. In addition, nitrogen molecules are larger than water molecules, so water molecules can detect some pores that nitrogen molecules cannot enter, especially ink bottle-shaped pores, which have a small entrance but a large actual pore volume. Based on the above principle, after determining the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set, the organic-inorganic pore volume of each target pore radius interval can be determined by comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval, and then the organic-inorganic pore volume of each target pore radius interval is added to determine the organic-inorganic pore volume of the target shale.

[0205] Specifically, the organic-inorganic pore volume of the target shale can be determined by the following method:

[0206] Let r = 1, compare the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set; if the cumulative nitrogen pore volume of the rth target pore radius interval is greater than the cumulative water vapor pore volume, then the difference between the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval is determined as the organic pore volume of the rth target pore radius interval, and the cumulative water vapor pore volume of the rth target pore radius interval is determined as the inorganic pore volume of the rth target pore radius interval; if the cumulative nitrogen pore volume of the rth target pore radius interval is less than the cumulative water vapor pore volume, then the cumulative water vapor pore volume of the rth target pore radius interval is determined as the inorganic pore volume of the rth target pore radius interval, and the organic pore volume of the rth target pore radius interval is determined as 0;

[0207] Determine whether r is equal to M, M is the number of target pore radius intervals in the target pore radius interval set; if r is not equal to M, then r = r + 1, and return to the step of comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set; if r is equal to M, the organic pore volume and the inorganic pore volume of the M target pore radius intervals are added respectively to obtain the total organic pore volume and the total inorganic pore volume of the target shale.

[0208] It should be noted that, in order to further illustrate the principle of determining the organic-inorganic pore volume of the target shale, the principle of detecting inorganic pores and organic pores by nitrogen and water vapor is described below.

[0209] The pore volume of the inorganic pores that can be detected by both water molecules and nitrogen molecules is V in1 (cm 3 / g), the pore volume of the inorganic pores that can be detected only by water molecules is V in2 (cm 3 / g), the pore volume of the organic pores detected by nitrogen molecules is V on (cm 3 / g), and the total pore volume is V t (cm 3 / g). Therefore, there is the following relationship among V in1 , V in2 and V on : V in1 >>V in2 , V in1 >>V on . In addition, there is the following relationship among V in1 , V in2 , V on , and :

[0210]

[0211] Finally, based on the sizes of V and V , the organic pore volume and the inorganic pore volume in different pore size ranges are determined, and the specific steps are as follows:

[0212] 1) If V , then V on >V in2 , because V on and V in2 are very small in themselves, so V in2 ≈0. Further, it can be obtained that the organic pore volume of the shale in this pore size range is V , the inorganic pore volume is V , and the total pore volume is V

[0213] 2) If V , then V on <V in2 , because V on and V in2 are very small in themselves, so V on ≈0. Further, it can be obtained that the organic pore volume of the shale in this pore size range is 0, the inorganic pore volume is V , and the total pore volume is

[0214] Here That is, V on ≠ V in2 This is because organic matter mainly develops micropores (pore diameter ≤ 2 nm), and inorganic minerals mainly develop mesopores and macropores (2 nm ≤ pore diameter). For micropores, V on > V in2 ; and for mesopores and macropores, V on < V in2 .

[0215] It should be noted that when determining the organic-inorganic pore volume of the target shale, the organic-inorganic pore volume of each target pore radius interval of the target pore radius interval set needs to be judged and determined first, that is, determined one by one in order of the size of the target pore radius interval from small to large. After determination, the organic-inorganic pore volumes of all target pore radius intervals can be summed to determine the organic-inorganic pore volume of the target shale, or the organic pore volume and the inorganic pore volume corresponding to each target pore radius interval can be output in the form of a column chart to more intuitively display the organic-inorganic pore volume of the target shale.

[0216] Figure 3 is a possible target shale organic-inorganic pore distribution chart shown by an embodiment of the present application, see Figure 3 In a possible implementation, the organic pore volume and the inorganic pore volume corresponding to each target pore radius interval are output in the form of a column chart to more intuitively display the organic-inorganic pore volume of the target shale. That is, the target pore radius interval set is [0-2mm], [2-4mm], [4-6mm], [6-8mm], [8-10mm], [10-12mm], [12-14mm], [14-16mm], [16-18mm], [18-20mm], the organic pore volume, the inorganic pore volume and the total pore volume of each target pore radius interval are determined, and output in the form of a column chart to more intuitively display the organic-inorganic pore volume distribution of the target shale.

[0217] In the embodiment of the present application, the nitrogen adsorption test data set and the water vapor adsorption test data set of the target shale can be obtained, the nitrogen adsorption test data set is processed by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set, the water vapor adsorption test data set is processed by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set, the target pore radius interval set is determined according to the first pore radius data set and the second pore radius data set, the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set are determined according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set, and the organic-inorganic pore volume of the target shale is determined according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set. That is, after the nitrogen adsorption test and the water vapor adsorption test of the target shale gas are performed, the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the corresponding interval of different pore radii are obtained by processing and calculating the nitrogen adsorption test and the water vapor adsorption test, the organic-inorganic pore volume of the target shale is determined by comparing the nitrogen pore volume and the cumulative water vapor pore volume, and the determination result can fully reflect the pore characteristics of the sample. Moreover, the determination result can be directly output by operation after the two adsorption tests, the test procedure is simplified, and the test time is shortened.

[0218] Figure 4 is a structural schematic diagram of a shale organic-inorganic pore volume determination device provided by the embodiment of the present application. Referring to Figure 4 , the device can include:

[0219] The acquisition module 401 is configured to acquire a nitrogen adsorption test data set and a water vapor adsorption test data set of a target shale, the nitrogen adsorption test data set refers to data of nitrogen adsorption amounts corresponding to N relative pressures at 77 degrees Kelvin, and the water vapor adsorption test data set refers to data of water vapor adsorption amounts corresponding to S relative pressures at 298 degrees Kelvin, N and S are positive integers greater than 1;

[0220] The first processing module 402 is configured to process the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set;

[0221] The second processing module 403 is configured to process the water vapor adsorption test data set by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set;

[0222] The first determining module 404 is configured to determine a target pore radius interval set according to the first pore radius data set and the second pore radius data set.

[0223] The second determining module 405 is configured to determine a cumulative nitrogen pore volume and a cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set.

[0224] The third determining module 406 is configured to determine an organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set.

[0225] Optionally, the first processing module is specifically configured to:

[0226] The nitrogen adsorption test data set is processed to obtain the first pore radius data set and the first pore volume data set by the following formula,

[0227]

[0228] Δv gn (p ri )=v gn (p ri+1 )-v gn (p ri )

[0229] Δt n (p ri )=t n (p ri+1 )-t n (p ri )

[0230] Δv ln (p ri )=0.0015468Δv gn (p ri )

[0231]

[0232] wherein, the p ri is the i-th relative pressure of N relative pressures in the nitrogen adsorption test data set, unit is dimensionless, i = 1, 2, 3, …, N; the t n (p ri ) is an adsorption layer thickness when the relative pressure is p ri , unit is nm; the r kn(p ri ) is the corresponding curvature radius, in nm, when the relative pressure is p ri ; the (p ri +p ri+1 ) / 2, in nm; the (p ri +p ri+1 ) / 2, in nm; the v gn (p ri ) is the nitrogen adsorption amount corresponding to the relative pressure p ri , in cm 3 / g; the Δv gn (p ri ) is the nitrogen desorption amount after the relative pressure decreases from p ri+1 to p ri , in cm 3 / g; the Δt n (p ri ) is the change amount of the adsorption layer thickness after the relative pressure decreases from p ri+1 to p ri , in nm; the Δv ln (p ri ) is the liquid nitrogen volume desorbed after the relative pressure decreases from p ri+1 to p ri , in cm 3 / g; the is the volume of the pores with the pore radius of , in cm 3 / g.

[0233] Optionally, the second processing module is specifically configured to:

[0234] The water vapor adsorption test data set is processed to obtain a second pore radius data set and a second pore volume data set by the following formula,

[0235] t h (p ri ) = 0.395 - 0.189 * ln(-ln(p ri ))

[0236]

[0237] Δv gh (p ri ) = v gh (p ri+1 ) - v gh (p ri)

[0238] Δt h (p ri )=t h (p ri+1 )-t h (p ri )

[0239] Δv lh (p ri )=8.0357×10 -4 Δv gh (p ri )

[0240]

[0241] wherein, the p ri is the i-th relative pressure of S relative pressures in the water vapor adsorption test data set, unit is dimensionless, i = 1, 2, 3, …, S; the t h (p ri ) is the adsorbed layer thickness of water when the relative pressure is p ri , unit is nm; the r kh (p ri ) is the corresponding curvature radius of water when the relative pressure is p ri , unit is nm; the is the corresponding curvature radius of water when the relative pressure is (p ri +p ri+1 ) / 2, unit is nm; the is the average pore radius detected by water when the relative pressure is (p ri +p ri+1 ) / 2, unit is nm; the v gh (p ri ) is the adsorption amount of water in standard state corresponding to the relative pressure p ri , unit is cm 3 / g; the Δv gh (p ri ) is the desorption amount of water in standard state after the relative pressure decreases from p ri+ 1 to p ri , unit is cm 3 / g; the Δt h (p ri ) is the change amount of water adsorbed layer thickness after the relative pressure decreases from p ri+1 to p ri , unit is nm; the Δv lh (p ri ) is the change amount of average pore radius detected by water after the relative pressure decreases from p ri+1 to pri volume of water desorbed, in cm 3 / g; the volume of pores with a pore radius of cm for water 3 / g.

[0242] Optionally, the first determining module comprises:

[0243] an acquisition unit, configured to acquire a maximum value of pore radii in the first pore radius data set and a maximum value of pore radii in the second pore radius data set;

[0244] a first determining unit, configured to determine, as an upper limit of the pore radius interval, a pore radius with a larger value between the maximum value of pore radii in the first pore radius data set and the maximum value of pore radii in the second pore radius data set;

[0245] a second determining unit, configured to determine 0 as a lower limit of the pore radius interval;

[0246] a division unit, configured to divide the pore radius interval into a plurality of target pore radius intervals according to a preset interval length;

[0247] a third determining unit, configured to determine the plurality of target pore radius intervals as the target pore radius interval set.

[0248] Optionally, the third determining module comprises:

[0249] a comparison unit, configured to let r = 1, and compare sizes of a cumulative nitrogen pore volume and a cumulative water vapor pore volume of an rth target pore radius interval in the target pore radius interval set;

[0250] a fourth determining unit, configured to, if the cumulative nitrogen pore volume of the rth target pore radius interval is greater than the cumulative water vapor pore volume, determine a difference between the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval as an organic pore volume of the rth target pore radius interval, and determine the cumulative water vapor pore volume of the rth target pore radius interval as an inorganic pore volume of the rth target pore radius interval;

[0251] a fifth determining unit, configured to, if the cumulative nitrogen pore volume of the rth target pore radius interval is less than the cumulative water vapor pore volume, determine the cumulative water vapor pore volume of the rth target pore radius interval as the inorganic pore volume of the rth target pore radius interval, and determine an organic pore volume of the rth target pore radius interval as 0;

[0252] The first judging unit is configured to judge whether r is equal to M, where M is the number of target pore radius intervals in the target pore radius interval set;

[0253] The second judging unit is configured to, if r is not equal to M, set r=r+1, and return to the step of comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the rth target pore radius interval in the target pore radius interval set.

[0254] The accumulating unit is configured to, if r is equal to M, accumulate the organic pore volume and the inorganic pore volume of the M target pore radius intervals respectively to obtain the total organic pore volume and the total inorganic pore volume of the target shale.

[0255] Optionally, the shale organic-inorganic pore volume determination device further comprises:

[0256] The second obtaining module is configured to obtain an original nitrogen adsorption relationship and an original water vapor adsorption relationship of the target shale, where the original nitrogen adsorption relationship refers to a linear relationship between a relative pressure and a nitrogen adsorption amount at 77K, and the original water vapor adsorption relationship refers to a linear relationship between a relative pressure and a water vapor adsorption amount at 298K.

[0257] The third processing module is configured to perform interpolation processing on the original nitrogen adsorption relationship by using a spline interpolation method to obtain data of nitrogen adsorption amounts corresponding to N relative pressures with the same interval, and determine the N relative pressures and the data of nitrogen adsorption amounts corresponding to the N relative pressures as a nitrogen adsorption test data set of the target shale.

[0258] The fourth processing module is configured to perform interpolation processing on the original water vapor adsorption relationship by using a spline interpolation method to obtain data of water vapor adsorption amounts corresponding to S relative pressures with the same interval, and determine the S relative pressures and the data of water vapor adsorption amounts corresponding to the S relative pressures as a water vapor adsorption test data set of the target shale.

[0259] In the embodiment of the present application, the shale organic-inorganic pore volume determination device can acquire a nitrogen adsorption test data set and a water vapor adsorption test data set of a target shale, process the nitrogen adsorption test data set by using a nitrogen adsorption pore calculation model to obtain a first pore radius data set and a first pore volume data set, process the water vapor adsorption test data set by using a water vapor adsorption pore calculation model to obtain a second pore radius data set and a second pore volume data set, determine a target pore radius interval set according to the first pore radius data set and the second pore radius data set, determine the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set according to the target pore radius interval set, the first pore radius data set, the second pore radius data set, the first pore volume data set and the second pore volume data set, and determine the organic-inorganic pore volume of the target shale according to the cumulative nitrogen pore volume and the cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set. That is, after one nitrogen adsorption test and one water vapor adsorption test are performed on the target shale, the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the corresponding interval of different pore radii are obtained by processing and calculating the nitrogen adsorption test and the water vapor adsorption test, the organic-inorganic pore volume of the target shale is determined by comparing the nitrogen pore volume and the cumulative water vapor pore volume, and the determination result can fully reflect the pore characteristics of the sample. Moreover, the determination result can be directly output by operation after two adsorption tests, which simplifies the test procedure and shortens the test time.

[0260] It should be noted that the shale organic-inorganic pore volume determination device provided in the above embodiments is only used as an example to illustrate the division of the functional modules in determining the shale organic-inorganic pore volume, and in actual applications, the above functions can be completed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the embodiments of the shale organic-inorganic pore volume determination device and the shale organic-inorganic pore volume determination method provided in the above embodiments belong to the same concept, and the specific implementation process is described in detail in the method embodiments, which will not be repeated here.

[0261] Figure 5Fig. 1 is a structural schematic diagram of a terminal 500 according to an embodiment of the present application. The terminal 500 can be a smart phone, a tablet computer, an MP3 player, an MP4 player, a notebook computer or a desktop computer. The terminal 500 can also be referred to as a user equipment, a portable terminal, a laptop terminal, a desktop terminal, or other names.

[0262] Generally, the terminal 500 includes a processor 501 and a memory 502.

[0263] The processor 501 can include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 501 can be implemented in at least one of a hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), a PLA (Programmable Logic Array). The processor 501 can also include a main processor and a coprocessor. The main processor is a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit). The coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 501 can be integrated with a GPU (Graphics Processing Unit) for rendering and drawing content required to be displayed by a display screen. In some embodiments, the processor 501 can further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

[0264] The memory 502 can include one or more computer-readable storage media, which can be non-transitory. The memory 502 can also include a high-speed random access memory, and a non-volatile memory such as one or more disk storage devices, flash storage devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 502 is used to store at least one instruction for being executed by the processor 501 to implement a shale organic-inorganic pore volume determination method provided by a method embodiment of the present application.

[0265] In some embodiments, terminal 500 can further include a peripheral device interface 503 and at least one peripheral device. The processor 501, the memory 502 and the peripheral device interface 503 can be connected through a bus or a signal line. Each peripheral device can be connected to the peripheral device interface 503 through a bus, a signal line or a circuit board. Specifically, the peripheral device includes at least one of a radio frequency circuit 504, a touch display screen 505, a camera 506, an audio circuit 507, a positioning component 508 and a power supply 509.

[0266] The peripheral device interface 503 can be used to connect at least one peripheral device related to I / O (Input / Output) to the processor 501 and the memory 502. In some embodiments, the processor 501, the memory 502 and the peripheral device interface 503 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 501, the memory 502 and the peripheral device interface 503 can be implemented on a separate chip or circuit board, and the present embodiments are not limited in this regard.

[0267] The radio frequency circuit 504 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 504 communicates with a communication network and other communication devices through electromagnetic signals. The radio frequency circuit 504 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals into electrical signals. Optionally, the radio frequency circuit 504 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and the like. The radio frequency circuit 504 can communicate with other terminals through at least one wireless communication protocol. The wireless communication protocol includes but is not limited to metropolitan area networks, various generations of mobile communication networks (2G, 3G, 4G and 5G), wireless local area networks and / or WiFi (Wireless Fidelity) networks. In some embodiments, the radio frequency circuit 504 can also include NFC (Near Field Communication) related circuitry, which is not limited in the present application.

[0268] The display screen 504 is configured to display a UI (User Interface). The UI can include graphics, text, icons, video, and any combination thereof. When the display screen 505 is a touch display screen, the display screen 505 is further configured to capture touch signals on or above the surface of the display screen 505. The touch signals can be input to the processor 501 as control signals for processing. In this case, the display screen 505 can also be configured to provide virtual buttons and / or virtual keyboard, also known as soft buttons and / or soft keyboard. In some embodiments, the display screen 505 can be one, arranged on the front panel of the terminal 500; in other embodiments, the display screen 505 can be at least two, arranged on different surfaces of the terminal 500 or in a folding design; in still other embodiments, the display screen 505 can be a flexible display screen, arranged on a curved surface or a folding surface of the terminal 500. Even, the display screen 505 can also be arranged in an irregular shape other than a rectangle, i.e., a special-shaped screen. The display screen 505 can be made of materials such as LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), etc.

[0269] The camera assembly 506 is configured to capture images or videos. Optionally, the camera assembly 506 includes a front camera and a rear camera. Typically, the front camera is arranged on the front panel of the terminal, and the rear camera is arranged on the back of the terminal. In some embodiments, the rear camera is at least two, which are any one of a main camera, a depth-of-field camera, a wide-angle camera, and a telephoto camera, to realize the background blur function by fusing the main camera and the depth-of-field camera, the panoramic shooting and VR (Virtual Reality) shooting function by fusing the main camera and the wide-angle camera, or other fusion shooting functions. In some embodiments, the camera assembly 506 can further include a flash. The flash can be a single-color-temperature flash or a dual-color-temperature flash. The dual-color-temperature flash refers to a combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0270] The audio circuit 507 can include a microphone and a speaker. The microphone is used to collect sound waves of a user and an environment, and convert the sound waves into an electrical signal input to the processor 501 for processing, or input to the radio frequency circuit 504 to realize voice communication. For the purpose of stereo sound collection or noise reduction, the microphone can be multiple, and arranged at different parts of the terminal 500. The microphone can also be an array microphone or an omnidirectional collection type microphone. The speaker is used to convert an electrical signal from the processor 501 or the radio frequency circuit 504 into sound waves. The speaker can be a conventional diaphragm speaker, or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only can it convert an electrical signal into a sound wave audible to humans, but it can also convert an electrical signal into an inaudible sound wave to humans for ranging purposes, etc. In some embodiments, the audio circuit 507 can also include a headphone jack.

[0271] The positioning component 508 is used to position the current geographic location of the terminal 500 to realize navigation or LBS (Location Based Service). The positioning component 508 can be a positioning component based on the GPS (Global Positioning System) of the United States, the Beidou system of China, the GLONASS system of Russia, or the Galileo system of the European Union.

[0272] The power supply 509 is used to supply power to various components in the terminal 500. The power supply 509 can be alternating current, direct current, disposable batteries, or rechargeable batteries. When the power supply 509 includes rechargeable batteries, the rechargeable batteries can support wired charging or wireless charging. The rechargeable batteries can also be used to support fast charging technology.

[0273] In some embodiments, the terminal 500 also includes one or more sensors 510. The one or more sensors 510 include, but are not limited to, an acceleration sensor 511, a gyroscope sensor 512, a pressure sensor 513, a fingerprint sensor 514, an optical sensor 515, and a proximity sensor 516.

[0274] The acceleration sensor 511 can detect the acceleration in three coordinate axes of the coordinate system established by the terminal 500. For example, the acceleration sensor 511 can be used to detect the components of gravitational acceleration in three coordinate axes. The processor 501 can control the touch display 505 to display a user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 511. The acceleration sensor 511 can also be used for gaming or collection of user motion data.

[0275] The gyroscope sensor 512 can detect the body direction and rotation angle of the terminal 500, and can collect 3D motions of a user on the terminal 500 in cooperation with the acceleration sensor 511. The processor 501 can implement the following functions according to the data collected by the gyroscope sensor 512: motion sensing (e.g., changing a UI according to a tilt operation of a user), image stabilization during photographing, game control, and inertial navigation.

[0276] The pressure sensor 513 can be disposed on a side frame of the terminal 500 and / or under the touch display 505. When the pressure sensor 513 is disposed on the side frame of the terminal 500, a grip signal of a user on the terminal 500 can be detected, and left / right hand recognition or a shortcut operation can be performed by the processor 501 according to the grip signal collected by the pressure sensor 513. When the pressure sensor 513 is disposed under the touch display 505, a pressure operation of a user on the touch display 505 can be detected by the processor 501, and an operable control on a UI can be controlled according to the pressure operation. The operable control includes at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0277] The fingerprint sensor 514 is used to collect a fingerprint of a user, and the identity of the user can be recognized by the processor 501 according to the fingerprint collected by the fingerprint sensor 514 or by the fingerprint sensor 514. When the identity of the user is recognized as a trusted identity, the processor 501 authorizes the user to perform a related sensitive operation, which includes unlocking a screen, viewing encrypted information, downloading software, payment, and changing a setting. The fingerprint sensor 514 can be disposed on the front, back, or side of the terminal 500. When a physical button or a manufacturer's logo is disposed on the terminal 500, the fingerprint sensor 514 can be integrated with the physical button or the manufacturer's logo.

[0278] The optical sensor 515 is used to collect ambient light intensity. In an embodiment, the processor 501 can control the display brightness of the touch display 505 according to the ambient light intensity collected by the optical sensor 515. Specifically, when the ambient light intensity is high, the display brightness of the touch display 505 is increased, and when the ambient light intensity is low, the display brightness of the touch display 505 is decreased. In another embodiment, the processor 501 can also dynamically adjust the photographing parameters of the camera assembly 506 according to the ambient light intensity collected by the optical sensor 515.

[0279] The proximity sensor 516, also known as a distance sensor, is typically located on the front panel of the terminal 500. The proximity sensor 516 is used to detect the distance between the user and the front of the terminal 500. In one embodiment, when the proximity sensor 516 detects that the distance between the user and the front of the terminal 500 is gradually decreasing, the processor 501 controls the touchscreen display 505 to switch from a screen-on state to a screen-off state; when the proximity sensor 516 detects that the distance between the user and the front of the terminal 500 is gradually increasing, the processor 501 controls the touchscreen display 505 to switch from a screen-off state to a screen-on state.

[0280] That is, embodiments of the present invention not only provide a terminal, including a processor and a memory for storing processor-executable instructions, wherein the processor is configured to execute... Figure 1 or Figure 2 Furthermore, embodiments of the present invention also provide a computer-readable storage medium storing a computer program that, when executed by a processor, can implement... Figure 1 or Figure 2 The method for determining the organic-inorganic pore volume of shale in the illustrated embodiment.

[0281] Those skilled in the art will understand that Figure 5 The structure shown does not constitute a limitation on terminal 500, and may include more or fewer components than shown, or combine certain components, or use different component arrangements.

[0282] Those skilled in the art will understand that all or part of the steps of the above embodiments can be implemented by hardware or by a program instructing related hardware. The program can be stored in a computer-readable storage medium, such as a read-only memory, a disk, or an optical disk.

[0283] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for determining the organic-inorganic pore volume of shale, characterized in that, The method for determining the organic-inorganic pore volume of shale includes: Obtain nitrogen adsorption test dataset and water vapor adsorption test dataset of the target shale. The nitrogen adsorption test dataset refers to the data of nitrogen adsorption amount corresponding to N relative pressures at 77 K degrees, and the water vapor adsorption test dataset refers to the data of water vapor adsorption amount corresponding to S relative pressures at 298 K degrees, where N and S are both positive integers greater than 1. Using a nitrogen adsorption pore calculation model, the nitrogen adsorption test dataset is processed to obtain a first pore radius dataset and a first pore volume dataset. The water vapor adsorption test dataset was processed using a water vapor adsorption pore calculation model to obtain a second pore radius dataset and a second pore volume dataset. Obtain the maximum value of the pore radius in the first pore radius dataset and the maximum value of the pore radius in the second pore radius dataset; The pore radius with the larger value between the maximum pore radius in the first pore radius dataset and the maximum pore radius in the second pore radius dataset is determined as the upper limit of the pore radius range. The value 0 is defined as the lower limit of the pore radius range; The pore radius interval is divided according to a preset interval length to obtain multiple target pore radius intervals; The plurality of target pore radius intervals are defined as the target pore radius interval set; Based on the target pore radius interval set, the first pore radius dataset, the second pore radius dataset, the first pore volume dataset, and the second pore volume dataset, determine the cumulative nitrogen pore volume and the cumulative water vapor pore volume for each target pore radius interval of the target pore radius interval set; The organic-inorganic pore volume of the target shale is determined based on the cumulative nitrogen pore volume and cumulative water vapor pore volume of each target pore radius interval in the target pore radius interval set.

2. The method for determining the organic-inorganic pore volume of shale according to claim 1, characterized in that, The nitrogen adsorption pore calculation model is used to process the nitrogen adsorption test dataset to obtain a first pore radius dataset and a first pore volume dataset, including: The nitrogen adsorption test dataset is processed using the following formulas to obtain the first pore radius dataset and the first pore volume dataset. Δv gn (p ri )=v gn (p ri+1 )-v gn (p ri ) Δt n (p ri )=t n (p ri+1 )-t n (p ri ) Δv ln (p ri )=0.0015468Δv gn (p ri ) Wherein, p ri The t represents the i-th relative pressure among N relative pressures in the nitrogen adsorption test dataset, in dimensionless units, where i = 1, 2, 3, ..., N; n (p ri When the relative pressure is p ri The adsorption layer thickness at that time is expressed in nm; the r kn (p ri When the relative pressure is p ri The corresponding radius of curvature, in nm; When the relative pressure is (p) ri +p ri+1 The radius of curvature corresponding to ) / 2 is expressed in nm; When the relative pressure is (p) ri +p ri+1 The average pore radius detected by nitrogen at ) / 2, in nm; the v gn (p ri The relative pressure is p. ri The corresponding nitrogen adsorption amount at that time, in cm³. 3 / g; the Δv gn (p ri () represents the relative pressure from p ri+1 Drop to p ri The amount of nitrogen adsorbed after desorption is expressed in cm³. 3 / g; the Δt n (p ri () represents the relative pressure from p ri+1 Drop to p ri The change in the thickness of the post-adsorption layer, in nm; the Δv lh (p ri () represents the relative pressure from p ri+1 Drop to p ri The volume of liquid nitrogen subsequently desorbed is expressed in cm³. 3 / g; the aforementioned The pore radius detected by nitrogen gas is The volume of the pores, in cm³. 3 / g.

3. The method for determining the organic-inorganic pore volume of shale according to claim 1, characterized in that, The water vapor adsorption pore calculation model is used to process the water vapor adsorption test dataset to obtain a second pore radius dataset and a second pore volume dataset, including: The water vapor adsorption test dataset is processed using the following formulas to obtain a second pore radius dataset and a second pore volume dataset. t h (p ri )=0.395-0.189×ln(-ln(p ri )) Δv gh (p ri )=v gh (p ri+1 )-v gh (p ri ) Δt h (p ri )=t h (p ri+1 )-t h (p ri ) Δv lh (p ri )=8.0357×10 -4 Δv gh (p ri ) Wherein, p ri The t represents the i-th relative pressure among the S relative pressures in the water vapor adsorption test dataset, in dimensionless units, i = 1, 2, 3, ..., S; h (p ri When the relative pressure is p ri The thickness of the water adsorption layer at that time, in nm; the r kh (p ri When the relative pressure is p ri The radius of curvature of the water at that time, in nm; When the relative pressure is (p) ri +p ri+1 The radius of curvature of water at ) / 2, in nm; When the relative pressure is (p) ri +p ri+1 The average pore radius detected in water at ) / 2 is expressed in nm; the v gh (p ri The relative pressure is p. ri The amount of water adsorbed under standard conditions, in cm³. 3 / g; the Δv gh (p ri () represents the relative pressure from p ri+1 Drop to p ri The amount of water adsorbed under standard conditions after desorption, expressed in cm³. 3 / g; the Δt h (p ri () represents the relative pressure from p ri+1 Drop to p ri The change in the thickness of the post-water adsorption layer, in nm; the Δv lh (p ri () represents the relative pressure from p ri+1 Drop to p ri The volume of water that is subsequently desorbed is expressed in cm³. 3 / g; the aforementioned The pore radius detected for water is The volume of the pores, in cm³. 3 / g.

4. The method for determining the organic-inorganic pore volume of shale according to claim 1, characterized in that, The step of determining the organic-inorganic pore volume of the target shale based on the cumulative nitrogen pore volume and cumulative water vapor pore volume of each target pore radius interval in the target pore radius interval set includes: Let r = 1, and compare the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the r-th target pore radius interval in the target pore radius interval set; If the cumulative nitrogen pore volume in the r-th target pore radius interval is greater than the cumulative water vapor pore volume, then the difference between the cumulative nitrogen pore volume and the cumulative water vapor pore volume in the r-th target pore radius interval is determined as the organic pore volume in the r-th target pore radius interval, and the cumulative water vapor pore volume in the r-th target pore radius interval is determined as the inorganic pore volume in the r-th target pore radius interval. If the cumulative nitrogen pore volume in the r-th target pore radius interval is less than the cumulative water vapor pore volume, then the cumulative water vapor pore volume in the r-th target pore radius interval is determined as the inorganic pore volume in the r-th target pore radius interval, and the organic pore volume in the r-th target pore radius interval is determined as 0. Determine whether r is equal to M, where M is the number of target pore radius intervals in the target pore radius interval set; If r is not equal to M, then let r = r + 1, and return to the step of comparing the cumulative nitrogen pore volume and the cumulative water vapor pore volume of the r-th target pore radius interval in the target pore radius interval set; If r equals M, the organic pore volume and inorganic pore volume of the M target pore radius intervals are summed to obtain the total organic pore volume and total inorganic pore volume of the target shale.

5. The method for determining the organic-inorganic pore volume of shale according to claim 1, characterized in that, Before obtaining the nitrogen adsorption test dataset and water vapor adsorption test dataset of the target shale, the following steps are also included: The original nitrogen adsorption relationship and the original water vapor adsorption relationship of the target shale are obtained. The original nitrogen adsorption relationship refers to the linear relationship between relative pressure and nitrogen adsorption amount at 77 K degrees, and the original water vapor adsorption relationship refers to the linear relationship between relative pressure and water vapor adsorption amount at 298 K degrees. The original nitrogen adsorption relationship was interpolated using spline interpolation to obtain data on nitrogen adsorption amounts corresponding to N relative pressures at equal intervals. The data on the N relative pressures and the nitrogen adsorption amounts corresponding to the N relative pressures were then used as the nitrogen adsorption test dataset for the target shale. The original water vapor adsorption relationship was interpolated using spline interpolation to obtain data on the water vapor adsorption amount corresponding to S relative pressures at equal intervals. The data on the S relative pressures and the water vapor adsorption amount corresponding to the S relative pressures were then determined as the water vapor adsorption test dataset of the target shale.

6. A device for determining the organic-inorganic pore volume of shale, characterized in that, The shale organic-inorganic pore volume determination device includes: The first acquisition module is used to acquire nitrogen adsorption test dataset and water vapor adsorption test dataset of the target shale. The nitrogen adsorption test dataset refers to the data of nitrogen adsorption amount corresponding to N relative pressures at 77 K degrees, and the water vapor adsorption test dataset refers to the data of water vapor adsorption amount corresponding to S relative pressures at 298 K degrees, where N and S are both positive integers greater than 1. The first processing module is used to process the nitrogen adsorption test dataset using a nitrogen adsorption pore calculation model to obtain a first pore radius dataset and a first pore volume dataset. The second processing module is used to process the water vapor adsorption test dataset using a water vapor adsorption pore calculation model to obtain a second pore radius dataset and a second pore volume dataset. The first determining module is used to determine the target pore radius interval set based on the first pore radius dataset and the second pore radius dataset; The second determining module is used to determine the cumulative nitrogen pore volume and cumulative water vapor pore volume of each target pore radius interval of the target pore radius interval set based on the target pore radius interval set, the first pore radius dataset, the second pore radius dataset, the first pore volume dataset, and the second pore volume dataset. The third determining module is used to determine the organic-inorganic pore volume of the target shale based on the cumulative nitrogen pore volume and cumulative water vapor pore volume of each target pore radius interval in the target pore radius interval set. The first determining module includes: The acquisition unit is used to acquire the maximum value of the pore radius in the first pore radius dataset and the maximum value of the pore radius in the second pore radius dataset; The first determining unit is used to determine the pore radius with the larger value between the maximum value of the pore radius in the first pore radius dataset and the maximum value of the pore radius in the second pore radius dataset as the upper limit of the pore radius range; The second determining unit is used to determine the value 0 as the lower limit of the pore radius range; A dividing unit is used to divide the pore radius interval according to a preset interval length to obtain multiple target pore radius intervals; The third determining unit is used to determine the plurality of target pore radius intervals as the target pore radius interval set.

7. A device for determining the organic-inorganic pore volume of shale, characterized in that, The device includes: processor; Memory used to store processor-executable instructions; The processor is configured to perform the method as described in any one of claims 1-5.

8. A computer-readable storage medium, characterized in that, The storage medium stores a computer program, which, when executed by a processor, implements the method described in any one of claims 1-5.

Images