A shale pore effectiveness evaluation processing method and device

Abstract

The present application provides a kind of shale pore availability evaluation processing method and device, it is related to oil and gas field development technical field.The method comprises: obtaining pore structure complex index, cumulative ink bottle pore volume ratio, improved pore structure complex index and improved cumulative ink bottle pore volume ratio;According to the pore structure complex index and the cumulative ink bottle pore volume ratio, pore effective index is calculated, and according to the improved pore structure complex index and the improved cumulative ink bottle pore volume ratio, improved pore effective index is calculated;According to the pore effective index and the improved pore effective index, the pore availability of the shale sample is evaluated, and shale pore availability evaluation result is obtained.The device executes above-mentioned method.The shale pore availability evaluation processing method and device provided in the embodiment of the present application can accurately and quantitatively evaluate shale pore availability.

Description

Technical Field

[0001] This invention relates to the field of oil and gas field development technology, specifically to a method and apparatus for evaluating the porosity of shale. Background Technology

[0002] With the continuous growth in demand for oil and gas resources, the development of unconventional oil and gas resources such as shale gas is receiving increasing attention. Shale, as an important unconventional oil and gas reservoir, has a complex and diverse pore structure. Accurately evaluating the porosity of shale is crucial for the exploration and development of shale oil and gas resources. Currently, there are many challenges in evaluating the porosity of shale. Traditional porosity evaluation methods often struggle to comprehensively consider multiple factors to accurately distinguish between effective and ineffective porosity, failing to fully reflect the actual contribution of shale porosity to oil and gas accumulation and migration. Furthermore, existing evaluation techniques have shortcomings in quantitative analysis, making it difficult to provide precise porosity values, which hinders the efficient exploration and development of shale oil and gas resources. Summary of the Invention

[0003] To address the problems in the prior art, embodiments of the present invention provide a method and apparatus for evaluating the porosity effectiveness of shale, which can at least partially solve the problems existing in the prior art.

[0004] On the one hand, this invention proposes a method for evaluating the porosity effectiveness of shale, comprising:

[0005] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0006] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0007] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0008] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0009] Among these, the porosity effectiveness of shale samples was tested to obtain a pore structure complexity index, including:

[0010] The total number of pores, the number of pores with a diameter smaller than the average pore diameter, the number of pores with a diameter larger than the average pore diameter, and the number of connected pores were obtained using digital core testing samples.

[0011] The pore structure complexity index is calculated based on the first ratio of the number of pores with a diameter smaller than the average pore diameter to the total number of pores, the second ratio of the number of pores with a diameter larger than the average pore diameter to the total number of pores, the third ratio of the number of interconnected pores to the total number of pores, and the preset weights corresponding to each ratio.

[0012] Among them, determining the preset weights corresponding to each ratio includes:

[0013] The number of pores with a diameter smaller than the average pore diameter, the number of pores with a diameter larger than the average pore diameter, and the number of connected pores were changed multiple times using artificial rock cores.

[0014] Based on the number of pores with a diameter smaller than the average pore diameter after each change, the average number of pores with a diameter smaller than the average pore diameter after multiple changes, the first pore structure complexity index calculated corresponding to the number of pores with a diameter smaller than the average pore diameter after each change, and the average of the first pore structure complexity index after multiple changes, the first correlation coefficient between the number of pores with a diameter smaller than the average pore diameter and the pore structure complexity index is calculated.

[0015] Based on the numerical range of the first correlation coefficient, determine the first preset weight corresponding to the first ratio;

[0016] Based on the number of pores with a diameter greater than the average pore diameter after each change, the average number of pores with a diameter greater than the average pore diameter after multiple changes, the second pore structure complexity index calculated corresponding to the number of pores with a diameter greater than the average pore diameter after each change, and the average of the second pore structure complexity index after multiple changes, a second correlation coefficient between the number of pores with a diameter greater than the average pore diameter and the pore structure complexity index is calculated.

[0017] Based on the numerical range of the second correlation coefficient, determine the second preset weight corresponding to the second ratio;

[0018] Based on the number of connected pores after each change, the average number of connected pores after multiple changes, the third pore structure complexity index calculated corresponding to the number of connected pores after each change, and the average of the third pore structure complexity index after multiple changes, the third correlation coefficient between the number of connected pores and the pore structure complexity index is calculated.

[0019] Based on the numerical range of the third correlation coefficient, a third preset weight corresponding to the third ratio is determined.

[0020] Among them, the porosity of shale samples was tested to obtain the cumulative ink bottle pore volume ratio, including:

[0021] Obtain the ink bottle orifice volume corresponding to each mercury intrusion cycle test, and use the ratio of the sum of the ink bottle orifice volumes of all cycles to the total mercury ingress in the first mercury intrusion cycle test as the cumulative ink bottle orifice volume ratio.

[0022] The step of obtaining the ink bottle orifice volume corresponding to each mercury intrusion cycle test includes:

[0023] For the first mercury intrusion cycle test, the ink bottle orifice volume is calculated using the following expression:

[0024]

[0025] Wherein, E(r) e )1 represents the amount of mercury removed in the first cycle, using the corrected pore throat radius r obtained during the mercury removal process as the independent variable, corresponding to the uncorrected pore throat radius r. I′(r) i )1 represents the amount of mercury introduced in the first cycle, with the corrected orifice throat radius r as the independent variable during the mercury introduction process;

[0026] For non-first mercury intrusion cyclic tests, the ink bottle orifice volume is calculated using the following expression:

[0027]

[0028] Where, E′(r e ) i I′(r) represents the amount of mercury removed in the i-th cycle, corresponding to the corrected pore throat radius r obtained during the mercury removal process. i ) i This represents the amount of mercury introduced in the i-th cycle, with the corrected orifice throat radius r as the independent variable during the mercury introduction process. imin The minimum value of the orifice throat radius r in the mercury ingress curve, r emax This represents the maximum value of the pore throat radius r in the mercury removal curve.

[0029] The step of evaluating the porosity of the shale sample based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation result includes:

[0030] The shale sample is evaluated for pore effectiveness based on a first comparison between the pore effectiveness index and a preset range value, and a second comparison between the pore effectiveness index and the improved pore effectiveness index, to obtain the shale pore effectiveness evaluation result.

[0031] On one hand, the present invention proposes a shale porosity effectiveness evaluation and processing device, comprising:

[0032] The acquisition unit is used to test the porosity of shale samples and obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio.

[0033] The testing unit is used to perform temperature and pressure improvement on the shale sample and to test the porosity effectiveness of the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0034] The calculation unit is used to calculate the effective pore index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and to calculate the improved pore effective index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0035] An evaluation unit is used to evaluate the porosity of the shale sample based on the porosity effective index and the improved porosity effective index, and to obtain the shale porosity effective evaluation result.

[0036] In another aspect, embodiments of the present invention provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the following method:

[0037] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0038] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0039] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0040] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0041] This invention provides a computer-readable storage medium, comprising:

[0042] The computer-readable storage medium stores a computer program that, when executed by a processor, implements the following method:

[0043] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0044] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0045] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0046] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0047] This invention also provides a computer program product, which includes a computer program that, when executed by a processor, implements the following method:

[0048] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0049] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0050] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0051] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0052] The shale porosity evaluation method and apparatus provided in this invention test the porosity of shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; perform temperature and pressure modification on the shale samples, and then test the porosity of the modified shale samples again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; calculate a porosity effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculate an improved porosity effectiveness index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; evaluate the porosity of the shale samples based on the porosity effectiveness index and the improved porosity effectiveness index to obtain shale porosity effectiveness evaluation results, which can accurately and quantitatively evaluate the porosity of shale. Attached Figure Description

[0053] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. In the drawings:

[0054] Figure 1 This is a schematic diagram of the structure of a shale porosity effectiveness evaluation and processing system provided in an embodiment of the present invention.

[0055] Figure 2 This is a magnified view of a portion of the sample placement stage provided in an embodiment of the present invention.

[0056] Figure 3 This is a schematic flowchart of a shale porosity effectiveness evaluation method provided in an embodiment of the present invention.

[0057] Figure 4 This is a schematic diagram of the shale porosity effectiveness evaluation and processing device provided in an embodiment of the present invention.

[0058] Figure 5 This is a schematic diagram of the physical structure of a computer device provided in an embodiment of the present invention. Detailed Implementation

[0059] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Here, the illustrative embodiments and descriptions of the present invention are used to explain the present invention, but are not intended to limit the present invention. It should be noted that, unless otherwise specified, the embodiments and features in the embodiments of this application can be arbitrarily combined with each other.

[0060] Figure 1 This is a schematic diagram of a shale porosity effectiveness evaluation and processing system provided in an embodiment of the present invention. The shale porosity effectiveness evaluation and processing system of the present invention is applied in the above system, such as... Figure 1 As shown, compared to existing systems, the system of the present invention has improved the sample placement stage, such as... Figure 2 As shown, a heating coil, a temperature sensor, and a pressure sensor have been added to simulate the heating and pressurization of the test environment.

[0061] The above system is described in terms of its functions as follows:

[0062] Sample preparation function:

[0063] Using a sample fixing platform (fixer), infrared locator, and telescopic drill bit ( Figure 1 (not shown), sample collection box ( Figure 1(Not shown) The sample fixing platform is used to fully fix samples of any shape on the platform to ensure the accuracy of the drilled sample size; the infrared locator is used to analyze and identify suitable drilling locations on the original sample surface to determine the drilling position; the telescopic drill bit is used for drilling the sample, and the telescopic rod above the drill bit can adjust the drill bit to an appropriate height according to the size of the original sample to drill the original sample. The final drilled sample size for testing is a cylindrical sample with a diameter of 1.5 cm and a length of 1 cm. The drill bit is equipped with a weighing sensor to obtain the weight of the sample in its initial state.

[0064] Data processing functions:

[0065] This is achieved through a computer and cables connecting the various functional components. It is used for the timely collection and storage of digital core data, mercury intrusion porosimetry data, temperature and pressure data, etc., and processes the data according to the set processing rules, and evaluates the final processing results.

[0066] Functions of artificial rock cores:

[0067] This is achieved using micro- and nano-CT scanners. Micro- and nano-CT scans are used to scan shale samples, and the acquired data is input into a computer for processing to obtain a digital core model. This digital core model is then input into an artificial core device, where artificial cores are created by mixing artificial materials of different proportions and specifications.

[0068] Circulating mercury injection function:

[0069] Through mercury porosimeter ( Figure 1 (not shown), sealed chamber ( Figure 1 Mercury content monitoring device (not shown), vacuum pump, mercury recovery tank (not shown) Figure 1 (Not shown) and a robotic arm are used to perform mercury porosimetry tests on the samples. Through multiple cycles of mercury porosimetry, the cumulative ink bottle pore volume ratio is obtained as one of the evaluation parameters for porosity effectiveness. The sealed chamber ensures the safety of the testing process and prevents mercury leakage. A mercury content monitoring device is located outside the sealed chamber to monitor the mercury content in the air during the test. When a certain threshold is reached, the monitoring device will sound an alarm. A vacuum pump is used to evacuate the sample after the mercury porosimetry cycle, ensuring complete removal of mercury. One end of the vacuum pump is connected to a mercury recovery tank for mercury recovery and reuse.

[0070] Temperature and pressure control function:

[0071] This device is used to modify samples that have undergone mercury intrusion porosimetry (MIP) testing and vacuuming under adjustable temperature and pressure. It can withstand high temperature and high pressure conditions and has good sealing performance. The modified sample is then sent back for cyclic mercury intrusion testing to obtain relevant parameters of the modified core.

[0072] A raw sample of any shape is placed on the fixed stage. The infrared locator automatically identifies and marks the drilling point, and the drill bit drills to obtain a shale sample of a specific size for testing, namely a columnar sample with a diameter of 1.5 cm and a length of 1 cm. Then, tests are conducted on the pore structure complexity and ink bottle pore volume.

[0073] Figure 3 This is a schematic flowchart of a shale porosity effectiveness evaluation method provided in an embodiment of the present invention, as shown below. Figure 3 As shown, the shale porosity effectiveness evaluation method provided in this embodiment of the invention includes:

[0074] Step S1: Perform a pore effectiveness test on the shale sample to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio.

[0075] Step S2: Perform temperature and pressure improvement on the shale sample, and conduct another pore effectiveness test on the improved shale sample to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0076] Step S3: Calculate the effective pore index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculate the improved effective pore index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0077] Step S4: Evaluate the porosity of the shale sample based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation result.

[0078] In step S1 above, the device performs a porosity effectiveness test on the shale sample to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio. The device can be a computer device, such as a server, that performs the method. The acquisition, storage, use, and processing of data in this application's technical solution all comply with relevant regulations.

[0079] Pore ​​effectiveness testing was performed on shale samples to obtain a pore structure complexity index, including:

[0080] The total number of pores, the number of pores with a diameter smaller than the average pore diameter, the number of pores with a diameter larger than the average pore diameter, and the number of connected pores were obtained using digital core testing samples. The total number of pores (N) in the digital core testing samples was then statistically analyzed. total And calculate the average pore size, and count the number (N) of pores with a pore size smaller than the average pore size. small ) and the number of pores with a pore size larger than the average pore size (N) large The number of connected pores (N) was statistically determined through image analysis. connected The method for determining connected pores is as follows:

[0081] The shortest path between two pores is used to determine whether they are connected. If such a path exists and the pores on the path meet certain geometric and physical conditions (such as the pore size being greater than a certain threshold, avoiding connecting tiny pores that may not have the function of actual fluid channels), then the two pores are considered to be connected.

[0082] The pore structure complexity index is calculated based on three ratios: a first ratio of the number of pores smaller than the average pore diameter to the total number of pores; a second ratio of the number of pores larger than the average pore diameter to the total number of pores; a third ratio of the number of interconnected pores to the total number of pores; and preset weights corresponding to each ratio. The pore structure complexity index (PSC) is defined as follows: It is calculated based on factors such as pore size distribution and connectivity analyzed using a digital core model. The index is comprehensively evaluated by calculating indicators such as the proportion of pores of different sizes and the proportion of interconnected pores. The calculation formula is as follows:

[0083]

[0084] Where, N small N represents the number of pores with a diameter smaller than the average pore diameter. total N represents the total number of pores. connected N represents the number of connected pores. large denoted as the number of pores with a diameter greater than the average pore diameter, and a, b, and c are preset weights obtained through artificial core testing.

[0085] Determine the preset weights corresponding to each ratio, including:

[0086] The number of pores with a diameter smaller than the average pore diameter, the number of pores with a diameter larger than the average pore diameter, and the number of connected pores were changed multiple times using artificial rock cores.

[0087] Based on the number of pores with a diameter smaller than the average pore diameter after each change, the average number of pores with a diameter smaller than the average pore diameter after multiple changes, the first pore structure complexity index calculated corresponding to the number of pores with a diameter smaller than the average pore diameter after each change, and the average of the first pore structure complexity index after multiple changes, the first correlation coefficient between the number of pores with a diameter smaller than the average pore diameter and the pore structure complexity index is calculated.

[0088] Based on the numerical range of the first correlation coefficient, determine the first preset weight corresponding to the first ratio;

[0089] Based on the number of pores with a diameter greater than the average pore diameter after each change, the average number of pores with a diameter greater than the average pore diameter after multiple changes, the second pore structure complexity index calculated corresponding to the number of pores with a diameter greater than the average pore diameter after each change, and the average of the second pore structure complexity index after multiple changes, a second correlation coefficient between the number of pores with a diameter greater than the average pore diameter and the pore structure complexity index is calculated.

[0090] Based on the numerical range of the second correlation coefficient, determine the second preset weight corresponding to the second ratio;

[0091] Based on the number of connected pores after each change, the average number of connected pores after multiple changes, the third pore structure complexity index calculated corresponding to the number of connected pores after each change, and the average of the third pore structure complexity index after multiple changes, the third correlation coefficient between the number of connected pores and the pore structure complexity index is calculated.

[0092] Based on the numerical range of the third correlation coefficient, a third preset weight corresponding to the third ratio is determined.

[0093] By selectively altering a specific factor affecting the pore structure complexity index (PSC) in artificial rock cores and observing the changes in the PSC value, N can be determined. small N large N connected The values ​​of their respective preset weights a, b, and c are calculated as follows:

[0094] Shale samples are scanned using micro-nano CT and the data is transmitted to a computer. The computer processes the data to obtain a digital core model, which is then replicated into nine digital core models, divided into three groups: A, B, and C, with three digital cores in each group.

[0095] Digital core data is input into an artificial core fabrication device to create artificial cores, and specific operations are performed on each group. For group A, the following operations are performed: keeping other conditions constant, the number of smaller pores is gradually changed by manually adjusting the particle size of the material in the fabrication process.

[0096] For group B, the following operation was performed: keeping other conditions unchanged, the number of larger pores was gradually changed by artificially adjusting the particle size of the material in the manufacturing process.

[0097] For group C, the following operation was performed: keeping other conditions unchanged, the connectivity between pores was affected by changing the amount of binder used in the artificial core fabrication process, thereby gradually changing the connectivity between pores.

[0098] Micro-nano CT scans were performed on all artificial core samples from each group. Following the method described above, with preset weights a, b, and c all set to 0.33, the pore structure complexity index (PSC) for each sample was calculated. The specific values ​​of each preset weight were determined using the Pearson correlation coefficient formula, calculated as follows:

[0099] The number of pores N with a pore size smaller than the average pore size is calculated using the following formula. small The first correlation coefficient S1 with the pore structure complexity index PSC is used to determine the value of the first preset weight a based on the magnitude of S1.

[0100]

[0101] Where S1 is N small The first correlation coefficient with PSC, N smalli This represents the number of pores whose diameter is smaller than the average pore diameter at the i-th time of the change. Taking i = 3 as an example... PSCa is the average number of pores with a pore size smaller than the average pore size in three changes. i Let be the pore structure complexity index for the i-th change in group A (the pore structure complexity index for the first change). The average value of the pore structure complexity index after three changes in Group A (the average value of the pore structure complexity index after the first change).

[0102] If S1 is close to 1, it means N small A strong positive correlation with PSC indicates that this factor has a significant impact on the PSC value; if S1 is close to 0, the correlation is weak, and the factor has a small impact on the PSC value. When S1 ≥ 0.8, set a = 0.4; when 0.8 > S1 ≥ 0.5, set a = 0.3; when S1 < 0.5, set a = 0.2.

[0103] The number of pores N with a diameter larger than the average pore diameter is calculated using the following formula. large The second correlation coefficient S2 with the pore structure complexity index PSC is used to determine the value of the second preset weight b.

[0104]

[0105] Where S2 is N large The second correlation coefficient with PSC, N largei This represents the number of pores whose diameter is greater than the average pore diameter at the i-th time. PSC is the average number of pores with a pore size larger than the average pore size during three changes. bi Let be the pore structure complexity index for the i-th change in group B (the pore structure complexity index for the second change). The average value of the pore structure complexity index after three changes in group B (the average value of the pore structure complexity index after the second change).

[0106] If S2 is close to 1, it means N large A strong positive correlation with PSC indicates that this factor has a significant impact on the PSC value; if S2 is close to 0, the correlation is weak, and the factor has a small impact on the PSC value. When S2 ≥ 0.8, set b = 0.4; when 0.8 > S2 ≥ 0.5, set b = 0.3; when S2 < 0.5, set b = 0.2.

[0107] The number of connected pores N is calculated using the following formula. connected The third correlation coefficient S3 with the pore structure complexity index PSC is used to determine the value of the third preset weight c.

[0108]

[0109] Where S3 is N connected The third correlation coefficient with PSC, N connectedi The number of connected pores at the i-th change. PSC is the average number of interconnected pores in three changes. ci Let be the pore structure complexity index for the i-th change in group C (the pore structure complexity index for the third change). The average value of the pore structure complexity index after three changes in group C (the average value of the pore structure complexity index after the third change).

[0110] If S3 is close to 1, it means N connected A strong positive correlation with PSC indicates that this factor has a significant impact on the PSC value; if S3 is close to 0, the correlation is weak, and the factor has a small impact on the PSC value. When S3 ≥ 0.8, set c = 0.4; when 0.8 > S3 ≥ 0.5, set c = 0.3; when S3 < 0.5, set c = 0.2.

[0111] Pore ​​effectiveness tests were performed on shale samples to obtain the cumulative ink bottle pore volume ratio, including:

[0112] Obtain the ink bottle orifice volume corresponding to each mercury intrusion cycle test, and use the ratio of the sum of the ink bottle orifice volumes of all cycles to the total mercury ingress in the first mercury intrusion cycle test as the cumulative ink bottle orifice volume ratio.

[0113] The process of obtaining the ink bottle orifice volume corresponding to each mercury intrusion porosimetry cycle includes:

[0114] For the first mercury intrusion cycle test, the ink bottle orifice volume is calculated using the following expression:

[0115]

[0116] Wherein, E(r) e )1 represents the amount of mercury removed in the first cycle, using the corrected pore throat radius r obtained during the mercury removal process as the independent variable, corresponding to the uncorrected pore throat radius r. I′(r) i )1 represents the amount of mercury introduced in the first cycle, with the corrected orifice throat radius r as the independent variable during the mercury introduction process;

[0117] For non-first mercury intrusion cyclic tests, the ink bottle orifice volume is calculated using the following expression:

[0118]

[0119] Where, E′(r e ) i I′(r) represents the amount of mercury removed in the i-th cycle, corresponding to the corrected pore throat radius r obtained during the mercury removal process. i ) i This represents the amount of mercury introduced in the i-th cycle, with the corrected orifice throat radius r as the independent variable during the mercury introduction process. imin The minimum value of the orifice throat radius r in the mercury ingress curve, r emax This represents the maximum value of the pore throat radius r in the mercury removal curve.

[0120] The volume of ink bottle pores in the sample is considered invalid pore volume. The higher the proportion of ink bottle pore volume, the worse the porosity. The test method for ink bottle pore volume is as follows:

[0121] The test can be set up with n mercury intrusion cycles as needed. During each cycle, the ingress and egress pressures and the amount of mercury introduced are recorded. The Washburn equation is used to convert the ingress and egress pressures P to the orifice throat radius r, yielding a distribution curve of the mercury introduction amount as a function of the orifice throat radius r. The Washburn equation is expressed as follows:

[0122]

[0123] The hysteresis phenomenon between the mercury ingress and egress curves mainly considers the influence of three factors: the pore volume of the ink bottle in the shale, the contact angle θ, and the surface tension γ. The mercury egress curve of the first cycle is corrected for the calculated pore throat radius r from the aspects of contact angle θ and surface tension γ using the Kloubek equation. It is assumed that the hysteresis phenomenon that still exists after correction is mainly affected by the pore volume of the ink bottle.

[0124] The Kloubek equation is expressed as follows:

[0125]

[0126] Where, r i r is the corrected throat radius during mercury injection. e This refers to the corrected pore throat radius during the mercury removal process. For details, please refer to the literature "Using multi-cycle mercury intrusion porosimetry to investigate hysteresis phenomenon of different porous media".

[0127] The difference between the mercury ingress curve of the first cycle and the corrected mercury removal curve is calculated in the interval [r]. imin ,r emax (where r) imin The minimum value of the orifice throat radius r in the mercury ingress curve, r emax The maximum value of the orifice throat radius in the mercury removal curve (both values ​​can be directly read from the above distribution curve of mercury intake with orifice throat radius r) is given by integrating r to obtain the ink bottle orifice volume V1 for the first cycle, as expressed below:

[0128]

[0129] In subsequent cycles, the mercury ingress and egress pressures P and the orifice throat radius r are still converted using the Washburn equation. However, unlike the first cycle, in the subsequent mercury ingress and egress curves, both the ingress and egress curves need to be corrected using the Kloubek equation. The difference between the corrected ingress and egress curves is then integrated to obtain the ink bottle orifice volume for the subsequent cycle, as shown in the following expression:

[0130]

[0131] Define the Cumulative Ink Bottle Orifice Volume Ratio (CVR): The test process is looped n times. In each loop, the ink bottle orifice volume for that loop is obtained and denoted as V1, V2, ..., V... n The accumulated amount of ink in the ink bottle through n cycles is added to the total mercury intake V of the first cycle. total The ratio, defined as the Cumulative Ink Bottle Pore Volume Ratio (CVR), is used to quantitatively evaluate the effectiveness of shale porosity. The calculation method is as follows:

[0132]

[0133] After the sample completes the cyclic mercury intrusion test, a vacuum pump is used to remove all residual mercury from the sample, and the removed mercury is collected in a mercury recovery tank.

[0134] In step S2 above, the device performs temperature and pressure improvement on the shale sample, and then performs a porosity effectiveness test on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio. For the mercury-free sample, temperature and pressure control improvement is performed. The required temperature and pressure parameters are set automatically, and the sample is improved under these conditions. The improved sample is then tested for porosity effectiveness to obtain the improved cumulative ink bottle pore volume ratio (CVRr) and the improved pore structure complexity index (PSCr).

[0135] In step S3 above, the device calculates the effective pore index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculates the improved effective pore index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0136] CVR reflects the proportion of ink bottle pores in the pore space of a shale sample. A higher CVR value indicates that ink bottle pores occupy more ineffective space and the pore effectiveness is lower. PSC measures the complexity of the pore structure in terms of pore size distribution and connectivity. A higher PSC value indicates a more complex pore structure, which also leads to a decrease in pore effectiveness.

[0137] Pore ​​Effectiveness Index (PEI) is defined as follows: The Pore Effectiveness Index (PEI) is a parameter that comprehensively considers the Cumulative Volume Ratio (CVR) and the Pore Complexity Index (PSC) to quantitatively characterize the porosity of shale samples. The calculation method is as follows:

[0138]

[0139] The Pore Effectiveness Stability Index (PEIr) is defined as the Improved Pore Effectiveness Index. The PEIr measures the degree of change in porosity of a shale sample under specific temperature and pressure alteration conditions. When PEIr is greater than PEI, it indicates that the temperature and pressure conditions improve shale porosity; if PEIr is less than PEI, porosity deteriorates. The calculation method is as follows:

[0140]

[0141] In step S4 above, the device evaluates the porosity of the shale sample based on the porosity effective index and the improved porosity effective index to obtain a shale porosity effective evaluation result. The step of evaluating the porosity of the shale sample based on the porosity effective index and the improved porosity effective index to obtain a shale porosity effective evaluation result includes:

[0142] The porosity of the shale sample is graded and evaluated based on a first comparison between the porosity effective index and a preset range value, and a second comparison between the porosity effective index and the improved porosity effective index, resulting in a shale porosity effectiveness grading evaluation. The porosity effectiveness of the shale sample is graded and evaluated using the Porosity Effective Index (PEI) and the Porosity Effectiveness Stability Index (PEIr), and the evaluation method is shown in Table 1.

[0143] Table 1

[0144]

[0145]

[0146] The shale porosity effectiveness evaluation method provided in this invention has the following beneficial technical effects:

[0147] This invention comprehensively utilizes various techniques, such as cyclic mercury intrusion porosimetry, digital core analysis, artificial core analysis, and thermobaric modification testing, to evaluate shale porosity from multiple perspectives and provide quantitative porosity values, thus overcoming the limitations of existing technologies in quantitative analysis. By comprehensively considering multiple factors, such as the cumulative ink bottle volume ratio (CVR) obtained from cyclic mercury intrusion porosimetry, the pore structure complexity index (PSC) obtained from digital and artificial core analysis, and the results of thermobaric modification testing, the Porosity Effectiveness Index (PEI) and Porosity Effectiveness Stability Index (PEIr) are calculated, improving the accuracy of porosity evaluation. Furthermore, this invention further improves the shale porosity evaluation technology system, providing new methods and ideas for related research and practice, and contributing to the development of shale porosity evaluation technology.

[0148] The shale porosity evaluation method provided in this invention involves testing the porosity of a shale sample to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; subjecting the shale sample to temperature and pressure modification, and then testing the porosity of the modified shale sample again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; calculating a porosity effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculating an improved porosity effectiveness index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; and evaluating the porosity effectiveness of the shale sample based on the porosity effectiveness index and the improved porosity effectiveness index to obtain a shale porosity effectiveness evaluation result, which can accurately and quantitatively evaluate the porosity effectiveness of shale.

[0149] Furthermore, pore effectiveness tests were conducted on the shale samples to obtain a pore structure complexity index, including:

[0150] The total number of pores, the number of pores with a diameter smaller than the average pore diameter, the number of pores with a diameter larger than the average pore diameter, and the number of connected pores are obtained by using digital core testing samples; the above embodiments can be referred to for explanation, and will not be repeated here.

[0151] The pore structure complexity index is calculated based on a first ratio of the number of pores smaller than the average pore diameter to the total number of pores, a second ratio of the number of pores larger than the average pore diameter to the total number of pores, a third ratio of the number of interconnected pores to the total number of pores, and preset weights corresponding to each ratio. This can be referred to the above embodiments for explanation, and will not be repeated here.

[0152] Further, the preset weights corresponding to each ratio are determined, including:

[0153] The number of pores with a diameter smaller than the average pore diameter, the number of pores with a diameter larger than the average pore diameter, and the number of interconnected pores were changed multiple times using artificial rock cores; the above embodiments can be referred to for explanation, and will not be repeated here.

[0154] Based on the number of pores with a diameter smaller than the average pore diameter after each change, the average number of pores with a diameter smaller than the average pore diameter after multiple changes, and the first pore structure complexity index calculated corresponding to the number of pores with a diameter smaller than the average pore diameter after each change, and the average of the first pore structure complexity index after multiple changes, a first correlation coefficient between the number of pores with a diameter smaller than the average pore diameter and the pore structure complexity index is calculated; this can be referred to the above embodiment for explanation, and will not be repeated here.

[0155] Based on the numerical range of the first correlation coefficient, a first preset weight corresponding to the first ratio is determined; this can be referred to the above embodiments for explanation, and will not be repeated here.

[0156] Based on the number of pores with a diameter greater than the average pore diameter after each change, the average number of pores with a diameter greater than the average pore diameter after multiple changes, and the second pore structure complexity index calculated corresponding to the number of pores with a diameter greater than the average pore diameter after each change, and the average of the second pore structure complexity index after multiple changes, a second correlation coefficient between the number of pores with a diameter greater than the average pore diameter and the pore structure complexity index is calculated; this can be referred to the above embodiment for explanation, and will not be repeated here.

[0157] Based on the numerical range of the second correlation coefficient, a second preset weight corresponding to the second ratio is determined; this can be referred to the above embodiment for explanation, and will not be repeated here.

[0158] Based on the number of connected pores after each change, the average number of connected pores after multiple changes, the third pore structure complexity index calculated corresponding to the number of connected pores after each change, and the average of the third pore structure complexity index after multiple changes, a third correlation coefficient between the number of connected pores and the pore structure complexity index is calculated; this can be referred to the above embodiments for explanation, and will not be repeated here.

[0159] Based on the numerical range of the third correlation coefficient, a third preset weight corresponding to the third ratio is determined. This can be referred to the above embodiment for further explanation, and will not be repeated here.

[0160] Furthermore, the porosity of the shale samples was tested to obtain the cumulative ink bottle pore volume ratio, including:

[0161] Obtain the ink bottle orifice volume corresponding to each mercury intrusion porosimetry cycle test, and use the ratio of the sum of the ink bottle orifice volumes of all cycles to the total mercury ingress amount in the first mercury intrusion porosimetry cycle test as the cumulative ink bottle orifice volume ratio. Refer to the above embodiment for further details.

[0162] Further, obtaining the ink bottle orifice volume corresponding to each mercury intrusion porosimetry cycle test includes:

[0163] For the first mercury intrusion cycle test, the ink bottle orifice volume is calculated using the following expression:

[0164]

[0165] Wherein, E(r) e )1 represents the amount of mercury removed in the first cycle, using the corrected pore throat radius r obtained during the mercury removal process as the independent variable, corresponding to the uncorrected pore throat radius r. I′(r) i )1 represents the amount of mercury introduced in the first cycle with the corrected orifice throat radius r as the independent variable; this can be referred to the above embodiment for explanation, and will not be repeated here.

[0166] For non-first mercury intrusion cyclic tests, the ink bottle orifice volume is calculated using the following expression:

[0167]

[0168] Where, E′(r e ) i I′(r) represents the amount of mercury removed in the i-th cycle, corresponding to the corrected pore throat radius r obtained during the mercury removal process. i ) i This represents the amount of mercury introduced in the i-th cycle, with the corrected orifice throat radius r as the independent variable during the mercury introduction process. iminThe minimum value of the orifice throat radius r in the mercury ingress curve, r emax This represents the maximum value of the pore throat radius r in the mercury removal curve. Refer to the above examples for further details; no further elaboration will be provided.

[0169] Further, the step of evaluating the porosity of the shale sample based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation result includes:

[0170] The shale sample is graded and evaluated for porosity effectiveness based on a first comparison result between the porosity effective index and a preset range value, and a second comparison result between the porosity effective index and the improved porosity effective index, thus obtaining the shale porosity effectiveness grading evaluation result. This can be referred to the above embodiments for explanation, and will not be repeated here.

[0171] Figure 4 This is a schematic diagram of the shale porosity effectiveness evaluation and treatment device provided in an embodiment of the present invention, as shown below. Figure 4 As shown, the shale porosity effectiveness evaluation and processing device provided in this embodiment of the invention includes an acquisition unit 401, a testing unit 402, a calculation unit 403, and an evaluation unit 404, wherein:

[0172] The acquisition unit 401 is used to perform a porosity effectiveness test on the shale sample to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; the testing unit 402 is used to perform temperature and pressure improvement on the shale sample, and then perform a porosity effectiveness test on the improved shale sample again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; the calculation unit 403 is used to calculate a porosity effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and to calculate an improved porosity effectiveness index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; the evaluation unit 404 is used to evaluate the porosity effectiveness of the shale sample based on the porosity effectiveness index and the improved porosity effectiveness index to obtain a shale porosity effectiveness evaluation result.

[0173] Specifically, the acquisition unit 401 in the device is used to perform a porosity effectiveness test on the shale sample to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; the testing unit 402 is used to perform temperature and pressure improvement on the shale sample, and then perform a porosity effectiveness test on the improved shale sample again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; the calculation unit 403 is used to calculate a porosity effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and to calculate an improved porosity effectiveness index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; the evaluation unit 404 is used to evaluate the porosity effectiveness of the shale sample based on the porosity effectiveness index and the improved porosity effectiveness index to obtain a shale porosity effectiveness evaluation result.

[0174] The shale porosity evaluation and processing device provided in this invention performs porosity effectiveness testing on shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; it then performs temperature and pressure modification on the shale samples and performs porosity effectiveness testing again on the modified shale samples to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, it calculates a porosity effectiveness index, and based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio, it calculates an improved porosity effectiveness index; based on the porosity effectiveness index and the improved porosity effectiveness index, it evaluates the porosity effectiveness of the shale samples to obtain shale porosity effectiveness evaluation results, which can accurately and quantitatively evaluate shale porosity effectiveness.

[0175] The embodiments of the present invention provide a shale porosity effectiveness evaluation and processing device that can be used to execute the processing flow of the above-described method embodiments. Its functions will not be repeated here, but can be referred to the detailed description of the above-described method embodiments.

[0176] Figure 5 This is a schematic diagram of the physical structure of a computer device provided in an embodiment of the present invention, such as... Figure 5 As shown, the computer device includes: a memory 501, a processor 502, and a computer program stored in the memory 501 and executable on the processor 502. When the processor 502 executes the computer program, it implements the following method:

[0177] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0178] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0179] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0180] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0181] This embodiment discloses a computer program product, which includes a computer program that, when executed by a processor, implements the following method:

[0182] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0183] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0184] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0185] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0186] This embodiment provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the following method:

[0187] Pore ​​effectiveness tests were conducted on shale samples to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0188] The shale sample was subjected to temperature and pressure modification, and the porosity of the modified shale sample was tested again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0189] The effective pore index is calculated based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and the improved pore effective index is calculated based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio.

[0190] The porosity of the shale sample is evaluated based on the porosity effective index and the improved porosity effective index to obtain the shale porosity effective evaluation results.

[0191] Compared with existing technologies, the shale porosity evaluation method provided in this invention involves testing the porosity of shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; subjecting the shale samples to temperature and pressure modification, and then testing the porosity of the modified shale samples again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; calculating a porosity effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculating an improved porosity effectiveness index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; and evaluating the porosity effectiveness of the shale samples based on the porosity effectiveness index and the improved porosity effectiveness index to obtain shale porosity effectiveness evaluation results, which can accurately and quantitatively evaluate shale porosity effectiveness.

[0192] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0193] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0194] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0195] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0196] In the description of this specification, the references to terms such as "an embodiment," "a specific embodiment," "some embodiments," "for example," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of the invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples.

[0197] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above descriptions are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for evaluating a shale pore effectiveness treatment, characterized in that, The method comprises the following steps: The shale sample is subjected to a pore effectiveness test to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; The shale sample is subjected to a temperature and pressure improvement, and the improved shale sample is subjected to a pore effectiveness test again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; A pore effectiveness index is calculated according to the pore structure complexity index and the cumulative ink bottle pore volume ratio, and an improved pore effectiveness index is calculated according to the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; The shale sample is subjected to a pore effectiveness evaluation according to the pore effectiveness index and the improved pore effectiveness index to obtain a shale pore effectiveness evaluation result; The pore structure complexity index is calculated according to the following formula: wherein, is the number of pores having a pore diameter smaller than the average pore diameter, is the total number of pores, is the number of connected pores, is the number of pores having a pore diameter larger than the average pore diameter, a, b, c are preset weights derived from artificial core testing; The cumulative ink bottle pore volume ratio is obtained, comprising: The ink bottle pore volume corresponding to each mercury injection cycle test is obtained, and the ratio of the sum of the ink bottle pore volumes of all times to the total mercury injection amount of the first mercury injection cycle test is taken as the cumulative ink bottle pore volume ratio; The improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio are obtained by subjecting the sample from which mercury has been discharged to temperature and pressure control improvement, setting the required temperature and pressure parameters by oneself, and improving the sample under the conditions, and then subjecting the improved sample to a pore effectiveness test; the pore effective exponent is calculated according to the formula: wherein, is the cumulative ink bottle pore volume ratio, is the pore structure complexity index; The improved pore effective exponent is calculated according to the following formula: wherein, is the modified cumulative ink bottle pore volume ratio, is the modified pore structure complexity index; The ink bottle pore volume corresponding to each mercury injection cycle test is obtained, comprising: For the first mercury injection cycle test, the ink bottle pore volume is calculated according to the following expression: wherein, represents the mercury withdrawal amount corresponding to the uncorrected pore throat radius r of the first cycle, with the pore throat radius r corrected in the mercury withdrawal process as the argument, represents the mercury injection amount corresponding to the corrected pore throat radius r of the first cycle, with the pore throat radius r corrected in the mercury injection process as the argument; For the non-first mercury injection cycle test, the ink bottle pore volume is calculated according to the following expression: wherein, represents the mercury withdrawal amount corresponding to the pore throat radius r of the i-th cycle obtained by correction, with the pore throat radius r corrected in the mercury withdrawal process as the argument, represents the mercury injection amount corresponding to the pore throat radius r of the i-th cycle obtained by correction, with the pore throat radius r corrected in the mercury injection process as the argument, is the minimum value of the pore throat radius r in the mercury injection curve, is the maximum value of the pore throat radius r in the mercury withdrawal curve.

2. The method for shale porosity effectiveness evaluation processing of claim 1, wherein, The shale sample is subjected to a pore effectiveness test to obtain a pore structure complexity index, comprising: The total pore quantity, the pore quantity with a pore diameter less than the average pore diameter, the pore quantity with a pore diameter greater than the average pore diameter, and the connected pore quantity are obtained by testing the sample by using a digital core; The pore structure complexity index is calculated according to a first ratio of the pore quantity with a pore diameter less than the average pore diameter to the total pore quantity, a second ratio of the pore quantity with a pore diameter greater than the average pore diameter to the total pore quantity, a third ratio of the connected pore quantity to the total pore quantity, and a preset weight corresponding to each ratio.

3. The method of shale porosity effectiveness evaluation processing of claim 2, wherein, The preset weight corresponding to each ratio is determined, comprising: The pore quantity with a pore diameter less than the average pore diameter, the pore quantity with a pore diameter greater than the average pore diameter, and the connected pore quantity are changed by using an artificial core multiple times; A first correlation coefficient of the pore quantity with a pore diameter less than the average pore diameter to the pore structure complexity index is calculated according to the pore quantity with a pore diameter less than the average pore diameter after each change, the average of the pore quantity with a pore diameter less than the average pore diameter after multiple changes, and a first changed pore structure complexity index corresponding to the pore quantity with a pore diameter less than the average pore diameter after each change; A first preset weight corresponding to the first ratio is determined according to a numerical interval in which the first correlation coefficient is located. According to the number of pores with a pore diameter greater than the average pore diameter after each change, an average of the number of pores with a pore diameter greater than the average pore diameter after multiple changes, and a second change pore structure complexity index calculated corresponding to the number of pores with a pore diameter greater than the average pore diameter after each change, an average of the second change pore structure complexity index after multiple changes, a second correlation coefficient of the number of pores with a pore diameter greater than the average pore diameter and the pore structure complexity index is calculated; According to the numerical interval where the second correlation coefficient is located, a second preset weight corresponding to the second ratio is determined; According to the number of connected pores after each change, an average of the number of connected pores after multiple changes, and a third change pore structure complexity index calculated corresponding to the number of connected pores after each change, an average of the third change pore structure complexity index after multiple changes, a third correlation coefficient of the number of connected pores and the pore structure complexity index is calculated; According to the numerical interval where the third correlation coefficient is located, a third preset weight corresponding to the third ratio is determined.

4. The method for shale porosity effectiveness evaluation processing according to any one of claims 1 to 3, characterized in that, The shale sample is evaluated for pore effectiveness according to the pore effectiveness index and the improved pore effectiveness index, and a shale pore effectiveness evaluation result is obtained. The shale sample is evaluated for pore effectiveness according to the pore effectiveness index and the improved pore effectiveness index, and a shale pore effectiveness evaluation result is obtained.

5. A shale pore effectiveness evaluation processing device characterized by comprising: It comprises: The acquisition unit is used for testing the pore effectiveness of the shale sample to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; The test unit is used for improving the shale sample by temperature and pressure, and testing the improved shale sample again for pore effectiveness to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; The calculation unit is used for calculating a pore effectiveness index according to the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculating an improved pore effectiveness index according to the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; The evaluation unit is used for evaluating the pore effectiveness of the shale sample according to the pore effectiveness index and the improved pore effectiveness index, and obtaining a shale pore effectiveness evaluation result; The pore structure complexity index is calculated according to the following formula: wherein, is the number of pores having a pore diameter smaller than the average pore diameter, is the total number of pores, is the number of connected pores, is the number of pores having a pore diameter larger than the average pore diameter, a, b, c are preset weights derived from artificial core testing; The cumulative ink bottle pore volume ratio is obtained, comprising: The ink bottle pore volume corresponding to each mercury injection cycle test is obtained, and the ratio of the sum of the ink bottle pore volumes of all times to the total mercury injection amount of the first mercury injection cycle test is taken as the cumulative ink bottle pore volume ratio; The improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio are obtained by testing the sample after the mercury has been discharged under the control of temperature and pressure, setting the required temperature and pressure parameters by itself, and improving the sample under the condition, and then testing the improved sample for pore effectiveness; the pore effective exponent is calculated according to the formula: wherein, is the cumulative ink bottle pore volume ratio, is the pore structure complexity index; The improved pore effective index is calculated according to the following formula: wherein, is the modified cumulative ink bottle pore volume ratio, is the modified pore structure complexity index; The ink bottle pore volume corresponding to each mercury injection cycle test is obtained, comprising: For the first mercury injection cycle test, the ink bottle pore volume is calculated according to the following expression: wherein, represents the mercury withdrawal amount corresponding to the uncorrected pore throat radius r of the first cycle, with the pore throat radius r corrected in the mercury withdrawal process as the argument, represents the mercury injection amount corresponding to the corrected pore throat radius r of the first cycle, with the pore throat radius r corrected in the mercury injection process as the argument; For non-first mercury intrusion cycle tests, the inkwell pore volume is calculated according to the following expression: wherein, represents the mercury withdrawal amount corresponding to the pore throat radius r of the i-th cycle obtained by correction, with the pore throat radius r corrected in the mercury withdrawal process as the argument, represents the mercury injection amount corresponding to the pore throat radius r of the i-th cycle obtained by correction, with the pore throat radius r corrected in the mercury injection process as the argument, is the minimum value of the pore throat radius r in the mercury injection curve, is the maximum value of the pore throat radius r in the mercury withdrawal curve.

6. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the method of any one of claims 1 to 4 when executing the computer program.

7. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program, and the computer program implements the method of any one of claims 1 to 4 when executed by a processor.

8. A computer program product, characterised in that, The computer program product comprises a computer program, and the computer program implements the method of any one of claims 1 to 4 when executed by a processor.

Images