Shale pore effectiveness evaluation processing method and device

Abstract

The invention provides a shale pore effectiveness evaluation processing method and device, and relates to the technical field of oil and gas field development. The method comprises the following steps: acquiring a pore structure complex index, a cumulative ink bottle pore volume ratio, an improved pore structure complex index and an improved cumulative ink bottle pore volume ratio; calculating to obtain a pore effective index according to the pore structure complex index and the cumulative ink bottle pore volume ratio, and calculating to obtain an improved pore effective index according to the improved pore structure complex index and the improved cumulative ink bottle pore volume ratio; performing pore effectiveness evaluation on the shale sample according to the pore effective index and the improved pore effective index to obtain a shale pore effectiveness evaluation result. The device executes the method. According to the shale pore effectiveness evaluation processing method and device provided by the embodiment of the invention, the shale pore effectiveness can be accurately and quantitatively evaluated.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas field development, and particularly relates to a method and device for evaluating the effectiveness of shale pores. Background Art

[0002] With the continuous growth of the demand for oil and gas resources, the development of unconventional oil and gas resources such as shale gas has attracted increasing attention. Shale, as an important unconventional oil and gas reservoir, has a complex and diverse pore structure. Accurately evaluating the effectiveness of shale pores is crucial for the exploration and development of shale oil and gas resources. At present, there are many challenges in the evaluation of shale pore effectiveness. Traditional pore evaluation methods often have difficulty in comprehensively considering various factors to accurately distinguish effective pores and ineffective pores, and cannot fully reflect the actual contribution of shale pores to oil and gas storage and migration. In addition, the existing evaluation technologies have deficiencies in quantitative analysis and are difficult to give accurate pore effectiveness values, which brings difficulties to the efficient exploration and development of shale oil and gas resources. Summary of the Invention

[0003] In view of the problems in the prior art, embodiments of the present invention provide a method and device for evaluating the effectiveness of shale pores, which can at least partially solve the problems existing in the prior art.

[0004] On the one hand, the present invention proposes a method for evaluating the effectiveness of shale pores, including:

[0005] Performing a pore effectiveness test on a shale sample to obtain a pore structure complexity index and a cumulative volume ratio of ink-bottle pores;

[0006] Performing temperature and pressure improvement on the shale sample, and performing a pore effectiveness test on the improved shale sample again to obtain an improved pore structure complexity index and an improved cumulative volume ratio of ink-bottle pores;

[0007] Calculating a pore effective index according to the pore structure complexity index and the cumulative volume ratio of ink-bottle pores, and calculating an improved pore effective index according to the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores;

[0008] Evaluating the pore effectiveness of the shale sample according to the pore effective index and the improved pore effective index to obtain a shale pore effectiveness evaluation result.

[0009] Wherein, performing a pore effectiveness test on a shale sample to obtain a pore structure complexity index includes:

[0010] Using a digital core test sample to obtain the total pore number, the pore number with a pore diameter smaller than the average pore diameter, the pore number with a pore diameter larger than the average pore diameter, and the connected pore number;

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

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

[0013] Using artificial cores to separately change the number of pores with pore diameters smaller than the average pore diameter, the number of pores with pore diameters larger than the average pore diameter, and the number of connected pores multiple times;

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

[0015] Based on the numerical interval where the first correlation coefficient is located, determine the first preset weight corresponding to the first ratio;

[0016] Based on the number of pores with pore diameters larger than the average pore diameter after each change, the average value of the number of pores with pore diameters larger than the average pore diameter after multiple changes, and the second changed pore structure complexity index calculated corresponding to the number of pores with pore diameters larger than the average pore diameter after each change, and the average value of the second changed pore structure complexity index after multiple changes, calculate the second correlation coefficient between the number of pores with pore diameters larger than the average pore diameter and the pore structure complexity index;

[0017] Based on the numerical interval where the second correlation coefficient is located, determine the second preset weight corresponding to the second ratio;

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

[0019] Based on the numerical interval where the third correlation coefficient is located, determine the third preset weight corresponding to the third ratio.

[0020] Among them, performing a pore effectiveness test on the shale sample to obtain the cumulative ink bottle pore volume ratio includes:

[0021] Obtain the volume of the ink bottle pores corresponding to each mercury injection cycle test respectively, and take the ratio of the sum of the volumes of the ink bottle pores for all times to the total mercury injection volume of the first mercury injection cycle test as the cumulative ink bottle pore volume ratio.

[0022] Among them, the obtaining of the volume of the ink bottle pores corresponding to each mercury injection cycle test respectively includes:

[0023] For the first mercury injection cycle test, the volume of the ink bottle pores is calculated according to the following expression:

[0024]

[0025] Among them, E(r e )1 represents the mercury withdrawal volume corresponding to the pore throat radius r obtained without correction in the first cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i )1 represents the mercury injection volume corresponding to the pore throat radius r obtained after correction in the first cycle with the corrected pore throat radius r during mercury injection as the independent variable;

[0026] For non-first mercury injection cycle tests, the volume of the ink bottle pores is calculated according to the following expression:

[0027]

[0028] Among them, E′(r e ) i represents the mercury withdrawal volume corresponding to the pore throat radius r obtained after correction in the i-th cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i ) i represents the mercury injection volume corresponding to the pore throat radius r obtained after correction in the i-th cycle with the corrected pore throat radius r during mercury injection as the independent variable, r imin is the minimum value of the pore throat radius r in the mercury injection curve, and r emax is the maximum value of the pore throat radius r in the mercury withdrawal curve.

[0029] Among them, the evaluating the pore effectiveness of the shale sample according to the pore effective index and the improved pore effective index to obtain the shale pore effectiveness evaluation result includes:

[0030] Conduct a hierarchical evaluation of the pore effectiveness of the shale sample according to the first comparison result of the pore effective index with the preset interval value and the second comparison result of the pore effective index with the improved pore effective index to obtain the shale pore effectiveness hierarchical evaluation result.

[0031] On the one hand, the present invention proposes a shale pore effectiveness evaluation processing device, including:

[0032] An acquisition unit for performing a pore effectiveness test on a shale sample to obtain a pore structure complexity index and a cumulative volume ratio of ink bottle pores;

[0033] A test unit for performing temperature and pressure improvement on the shale sample and performing a pore effectiveness test on the improved shale sample again to obtain an improved pore structure complexity index and an improved cumulative volume ratio of ink bottle pores;

[0034] A calculation unit for calculating a pore effective index based on the pore structure complexity index and the cumulative volume ratio of ink bottle pores, and calculating an improved pore effective index based on the improved pore structure complexity index and the improved cumulative volume ratio of ink bottle pores;

[0035] An evaluation unit for evaluating the pore effectiveness of the shale sample based on the pore effective index and the improved pore effective index to obtain a shale pore effectiveness evaluation result.

[0036] On the other hand, an embodiment of the present invention provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the following method is implemented:

[0037] Perform a pore effectiveness test on a shale sample to obtain a pore structure complexity index and a cumulative volume ratio of ink bottle pores;

[0038] Perform temperature and pressure improvement on the shale sample and perform a pore effectiveness test on the improved shale sample again to obtain an improved pore structure complexity index and an improved cumulative volume ratio of ink bottle pores;

[0039] Calculate a pore effective index based on the pore structure complexity index and the cumulative volume ratio of ink bottle pores, and calculate an improved pore effective index based on the improved pore structure complexity index and the improved cumulative volume ratio of ink bottle pores;

[0040] Evaluate the pore effectiveness of the shale sample based on the pore effective index and the improved pore effective index to obtain a shale pore effectiveness evaluation result.

[0041] An embodiment of the present invention provides a computer-readable storage medium, including:

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

[0043] Perform a pore effectiveness test on a shale sample to obtain a pore structure complexity index and a cumulative volume ratio of ink bottle pores;

[0044] Perform temperature and pressure improvement on the shale sample, and conduct pore effectiveness testing on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores;

[0045] Calculate the pore effectiveness index based on the pore structure complexity index and the cumulative volume ratio of ink-bottle pores, and calculate the improved pore effectiveness index based on the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores;

[0046] Evaluate the pore effectiveness of the shale sample based on the pore effectiveness index and the improved pore effectiveness index to obtain the shale pore effectiveness evaluation result.

[0047] An embodiment of the present invention also provides a computer program product, which includes a computer program. When the computer program is executed by a processor, the following method is implemented:

[0048] Conduct pore effectiveness testing on a shale sample to obtain the pore structure complexity index and the cumulative volume ratio of ink-bottle pores;

[0049] Perform temperature and pressure improvement on the shale sample, and conduct pore effectiveness testing on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores;

[0050] Calculate the pore effectiveness index based on the pore structure complexity index and the cumulative volume ratio of ink-bottle pores, and calculate the improved pore effectiveness index based on the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores;

[0051] Evaluate the pore effectiveness of the shale sample based on the pore effectiveness index and the improved pore effectiveness index to obtain the shale pore effectiveness evaluation result.

[0052] The shale pore effectiveness evaluation processing method and device provided by the embodiment of the present invention conduct pore effectiveness testing on a shale sample to obtain the pore structure complexity index and the cumulative volume ratio of ink-bottle pores; perform temperature and pressure improvement on the shale sample, and conduct pore effectiveness testing on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores; calculate the pore effectiveness index based on the pore structure complexity index and the cumulative volume ratio of ink-bottle pores, and calculate the improved pore effectiveness index based on the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores; evaluate the pore effectiveness of the shale sample based on the pore effectiveness index and the improved pore effectiveness index to obtain the shale pore effectiveness evaluation result, and can accurately and quantitatively evaluate the shale pore effectiveness. Description of the Drawings

[0053] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings. In the drawings:

[0054] Figure 1 It is a schematic structural diagram of a shale pore effectiveness evaluation processing system provided by an embodiment of the present invention.

[0055] Figure 2 It is an enlarged partial structural view of a sample placement table provided by an embodiment of the present invention.

[0056] Figure 3 It is a schematic flowchart of a shale pore effectiveness evaluation processing method provided by an embodiment of the present invention.

[0057] Figure 4 It is a schematic structural diagram of a shale pore effectiveness evaluation processing device provided by an embodiment of the present invention.

[0058] Figure 5 It is a schematic diagram of the physical structure of a computer device provided by an embodiment of the present invention. Detailed implementation manners

[0059] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer and more understandable, the following further elaborates on the embodiments of the present invention in conjunction with the accompanying drawings. Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but not to limit the present invention. It should be noted that, without conflict, the embodiments and features in the present application can be arbitrarily combined with each other.

[0060] Figure 1 It is a schematic structural diagram of a shale pore effectiveness evaluation processing system provided by an embodiment of the present invention. The shale pore effectiveness evaluation processing of the present invention is applied in the above system. As Figure 1 shown, compared with the existing system, the system of the present invention has improved the sample placement table. As Figure 2 shown, specifically, a heating coil, a temperature sensor, and a pressure sensor are added to simulate the heating and pressurization of the test environment.

[0061] The above system is described as follows according to its functions:

[0062] Sample preparation function:

[0063] Through the sample fixing platform (fixator), infrared locator, telescopic drill bit ( Figure 1 not shown), sample collection box ( Figure 1It is implemented by (not shown). The sample fixing platform is used to fully fix samples of any shape on the platform to ensure the precise specifications of the drilled samples; the infrared locator is used to analyze and identify the suitable parts for drilling on the surface of the original sample to determine the drilling position of the sample; the telescopic drill bit is used for drilling the sample. 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. Finally, the testable sample obtained by drilling is a columnar 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 function:

[0065] It is implemented through a computer and the cables connecting the functional components. It is used to collect and store digital core data, mercury injection data, temperature and pressure data, etc. in a timely manner, process the data according to the set processing rules, and evaluate based on the final processing results.

[0066] Artificial core function:

[0067] It is implemented through a micro-CT scanner, a nano-CT scanner, etc. The shale sample is scanned by micro-nano CT, the scanned data is input into a computer for processing to obtain a digital core model, and then the digital core model is input into the artificial core device to produce an artificial core by mixing artificial materials with different proportions and specifications.

[0068] Cyclic mercury injection function:

[0069] It is implemented through a mercury injection instrument ( Figure 1 not shown), a sealed chamber ( Figure 1 not shown), a mercury content monitoring device, a vacuum pump, a mercury recovery tank ( Figure 1 not shown), and a robotic arm. The mercury injection instrument is used to perform mercury injection tests on the sample. By performing multiple cyclic mercury injections, the cumulative ink bottle pore volume ratio is obtained as one of the evaluation parameters for pore effectiveness; the sealed chamber is used to ensure the safety of the test process and prevent mercury leakage during the test; the mercury content monitoring device is set outside the sealed chamber to monitor the mercury content in the air during the test. When a certain threshold is reached, the monitoring device gives an alarm; the vacuum pump is used to evacuate the sample after the cyclic mercury injection is completed to completely discharge the mercury in the sample. One end of the vacuum pump is connected to the mercury recovery tank to recycle the mercury.

[0070] Temperature and pressure control function:

[0071] It is used to transform the sample after the mercury injection test and evacuation under adjustable temperature and pressure, can withstand high temperature and high pressure conditions and has good sealing performance. The transformed sample is sent back to the cyclic mercury injection for testing to obtain the relevant parameters of the core after transformation.

[0072] Place the original sample of any shape on the fixing table, and the infrared locator automatically identifies and delimits the drilling points, and the drill bit drills to obtain a shale sample of a specific specification for testing, that is, a columnar sample with a diameter of 1.5 cm and a length of 1 cm. Then, tests on the pore structure complexity and the volume of ink-bottle pores are carried out respectively.

[0073] Figure 3 is a schematic flow chart of the shale pore effectiveness evaluation processing method provided by an embodiment of the present invention. As Figure 3 shown, the shale pore effectiveness evaluation processing method provided by the embodiment of the present invention includes:

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

[0075] Step S2: Perform temperature and pressure improvement on the shale sample, and perform a pore 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.

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

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

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

[0079] Performing a pore effectiveness test on the shale sample to obtain a pore structure complexity index includes:

[0080] Using a digital core test sample to obtain the total pore number, the pore number with a pore diameter smaller than the average pore diameter, the pore number with a pore diameter larger than the average pore diameter, and the connected pore number; counting the total pore number (N total ) in the digital core test sample and calculating the average pore diameter, and respectively counting the pore number (N small ) with a pore diameter smaller than the average pore diameter and the pore number (N large ) of pores with a pore diameter larger than the average pore diameter. The connected pore number (N connected ) is statistically obtained through image analysis. The determination method of connected pores is as follows:

[0081] Find the shortest path between two pores 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 to avoid connecting tiny pores that may not have the actual function of fluid channels), then these two pores are considered connected.

[0082] According to the first ratio of the number of pores with pore sizes smaller than the average pore size to the total number of pores, the second ratio of the number of pores with pore sizes larger than the average pore size to the total number of pores, the third ratio of the number of connected pores to the total number of pores, and the preset weights corresponding to each ratio respectively, calculate the pore structure complexity index. Define the pore structure complexity index (PSC): Analyze factors such as the size distribution and connectivity of pores based on the digital core model, and calculate the pore structure complexity index. Comprehensively evaluate by calculating indicators such as the proportion of pores of different sizes and the proportion of connected pores. The calculation formula is as follows:

[0083]

[0084] Among them, N small is the number of pores with pore sizes smaller than the average pore size, N total is the total number of pores, N connected is the number of connected pores, N large is the number of pores with pore sizes larger than the average pore size, and a, b, and c are preset weights obtained through artificial core tests.

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

[0086] Use the artificial core to change the number of pores with pore sizes smaller than the average pore size, the number of pores with pore sizes larger than the average pore size, and the number of connected pores separately multiple times;

[0087] According to the number of pores with pore sizes smaller than the average pore size after each change, the average value of the number of pores with pore sizes smaller than the average pore size after multiple changes, and the first changed pore structure complexity index calculated corresponding to the number of pores with pore sizes smaller than the average pore size after each change, and the average value of the first changed pore structure complexity index after multiple changes, calculate the first correlation coefficient between the number of pores with pore sizes smaller than the average pore size and the pore structure complexity index;

[0088] According to the numerical interval where the first correlation coefficient is located, determine the first preset weight corresponding to the first ratio;

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

[0090] Based on the numerical interval where the second correlation coefficient is located, determine the second preset weight corresponding to the second ratio;

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

[0092] Based on the numerical interval where the third correlation coefficient is located, determine the third preset weight corresponding to the third ratio.

[0093] By specifically changing one of the factors that affect the pore structure complexity index (PSC) in the artificial core and observing the change of the PSC value, thereby determining the values of the corresponding preset weights a, b, and c of N small 、N large 、N connected The calculation method is as follows:

[0094] Scan the shale sample by micro-nano CT and transfer the data to a computer. The computer processes the data to obtain a digital core model, copy the digital core model into 9, and divide them into three groups A, B, and C, with 3 digital cores in each group.

[0095] Input the digital core data into an artificial core device to make artificial cores, and perform specific operations on each group. The following operations are performed on group A: keeping other conditions unchanged, gradually change the number of pores with a smaller pore diameter by artificially adjusting the particle size of the materials in the manufacturing process.

[0096] The following operations are performed on group B: keeping other conditions unchanged, gradually change the number of pores with a larger pore diameter by artificially adjusting the particle size of the materials in the manufacturing process.

[0097] The following operations are performed on group C: keeping other conditions unchanged, affect the connectivity between pores by changing the amount of binder in the process of making the artificial core, and gradually change the connectivity between pores.

[0098] All artificial cores in each group were scanned using micro-nano CT. According to the above method, when the preset weights a, b, and c were all set to 0.33, the respective pore structure complexity indices (PSC) were calculated. The specific values of each preset weight were obtained by analyzing using the Pearson correlation coefficient formula, and the calculation method is as follows:

[0099] The number of pores N with pore diameters smaller than the average pore diameter was calculated through the following formula small The first correlation coefficient S1 with the pore structure complexity index PSC was determined, and the value of the first preset weight a was determined according to the magnitude of S1.

[0100]

[0101] where S1 is the first correlation coefficient of N small with PSC, and N smalli is the number of pores with pore diameters smaller than the average pore diameter at the i-th change. Taking i equal to 3 as an example, is the average value of the number of pores with pore diameters smaller than the average pore diameter in 3 changes, and PSCa i is the pore structure complexity index (the first changed pore structure complexity index) at the i-th change in group A, is the average value of the pore structure complexity indices after 3 changes in group A (the average value of the first changed pore structure complexity indices).

[0102] If S1 is close to 1, it indicates that N small is highly positively correlated with PSC, that is, this factor has a greater impact on the PSC value; if S1 is close to 0, it indicates a weak correlation, and this factor has a smaller 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 pore diameters larger than the average pore diameter was calculated through the following formula large The second correlation coefficient S2 with the pore structure complexity index PSC was determined, and the value of the second preset weight b was determined according to the magnitude of S2.

[0104]

[0105] where S2 is the second correlation coefficient of N large with PSC, and N largei is the number of pores with pore diameters larger than the average pore diameter at the i-th change, is the average value of the number of pores with pore diameters larger than the average pore diameter in 3 changes, and PSC bi is the pore structure complexity index (the second changed pore structure complexity index) at the i-th change in group B, It is the average value of the pore structure complexity index after 3 changes in Group B (the average value of the second changed pore structure complexity index).

[0106] If S2 is close to 1, it means that N large Is highly positively correlated with PSC, that is, this factor has a greater impact on the PSC value; if S2 is close to 0, it means that the correlation is weak and this factor has a smaller 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] Calculate the number of connected pores N through the following formula connected The third correlation coefficient S3 with the pore structure complexity index PSC, and determine the value of the third preset weight c according to the size of S3.

[0108]

[0109] Among them, S3 is N connected The third correlation coefficient with PSC, N connectedi Is the number of connected pores at the i-th change, Is the average value of the number of connected pores in 3 changes, PSC ci Is the pore structure complexity index at the i-th change in Group C (the third changed pore structure complexity index), Is the average value of the pore structure complexity index after 3 changes in Group C (the average value of the third changed pore structure complexity index).

[0110] If S3 is close to 1, it means that N connected Is highly positively correlated with PSC, that is, this factor has a greater impact on the PSC value; if S3 is close to 0, it means that the correlation is weak and this factor has a smaller 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] Perform pore effectiveness testing on the shale sample to obtain the cumulative inkbottle pore volume ratio, including:[[]]

[0112] Obtain the inkbottle pore volume corresponding to each mercury injection cycle test respectively, and take the ratio of the sum of the inkbottle pore volumes of all times to the total mercury injection volume of the first mercury injection cycle test as the cumulative inkbottle pore volume ratio.

[0113] The obtaining of the inkbottle pore volume corresponding to each mercury injection cycle test respectively includes:[[]]

[0114] For the first mercury injection cycle test, calculate the inkbottle pore volume according to the following expression:[[]]

[0115]

[0116] Among them, E(r e )1 represents the mercury withdrawal volume corresponding to the pore throat radius r obtained without correction in the first cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i )1 represents the mercury injection volume corresponding to the pore throat radius r obtained with correction in the first cycle with the corrected pore throat radius r during mercury injection as the independent variable;

[0117] For non-first mercury injection cycle tests, the ink bottle pore volume is calculated according to the following expression:

[0118]

[0119] Among them, E′(r e ) i represents the mercury withdrawal volume corresponding to the pore throat radius r obtained with correction in the i-th cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i ) i represents the mercury injection volume corresponding to the pore throat radius r obtained with correction in the i-th cycle with the corrected pore throat radius r during mercury injection as the independent variable, r imin is the minimum value of the pore throat radius r in the mercury injection curve, and r emax is the maximum value of the pore throat radius r in the mercury withdrawal curve.

[0120] The ink bottle pore volume in the sample is regarded as the ineffective pore volume. The higher the proportion of the ink bottle pore volume, the worse the pore effectiveness. The test method for the ink bottle pore volume is as follows:

[0121] The test can set n mercury injection cycles as needed. Each time mercury is injected, the advancing and retreating mercury pressures and the mercury injection volume during the process are recorded. The advancing and retreating mercury pressures P and the pore throat radius r are converted through the Washburn equation to obtain the distribution curve of the mercury injection volume with respect to the pore throat radius r. The expression of the Washburn equation is as follows:

[0122]

[0123] Among them, the hysteresis phenomenon between the advancing and retreating mercury curves mainly considers the influence of three aspects: the ink bottle pore volume in the shale, the contact angle θ, and the surface tension γ. The mercury withdrawal curve of the first cycle is corrected for the converted pore throat radius r from two aspects of the contact angle θ and the surface tension γ through the Kloubek equation, and it is considered that the remaining hysteresis phenomenon after correction is mainly affected by the ink bottle pore volume.

[0124] The expression of the Kloubek equation is as follows:

[0125]

[0126] where r i is the corrected pore throat radius during mercury injection, and r e is the corrected pore throat radius during mercury withdrawal. Specifically, it can refer to the literature "Using multi-cycle mercury intrusion porosimetry to investigate hysteresis phenomenon of different porous media".

[0127] Subtract the mercury injection curve of the first cycle from the corrected mercury withdrawal curve and then integrate r in the interval [r imin , r emax (where r imin is the minimum value of the pore throat radius r in the mercury injection curve, and r emax is the maximum value of the pore throat radius in the mercury withdrawal curve. Both of these values can be directly read from the distribution curve of mercury injection volume versus pore throat radius r mentioned above) to obtain the ink bottle pore volume V1 of the first cycle. The expression is as follows:

[0128]

[0129] In subsequent cycles, the advancing and retreating mercury pressures P and the pore throat radius r are still converted through the Washburn equation. Different from the first cycle, in the distribution curve of mercury injection volume versus pore throat radius r obtained in subsequent cycles, both the mercury injection and withdrawal curves need to be corrected by the Kloubek equation, and then the corrected advancing and retreating mercury curves are subtracted and integrated to obtain the ink bottle pore volume of subsequent cycles. The expression is as follows:

[0130]

[0131] Define the cumulative ink bottle pore volume ratio (CVR): The number of cycles in the test process is set to n times. Each cycle will obtain an ink bottle pore volume for that cycle, denoted as V1, V2,..., V n ,. The ratio of the sum of the ink bottle pore volumes obtained in n cycles to the total mercury injection volume V total of the first cycle is defined as the cumulative ink bottle pore volume ratio (CVR) for quantitatively evaluating the effectiveness of shale pores. The calculation method is as follows:

[0132]

[0133] After the sample completes the cyclic mercury injection test, use a vacuum pump to evacuate the sample to discharge all the residual mercury in the sample, and collect the discharged mercury into the mercury recovery tank.

[0134] In the above step S2, the device improves the temperature and pressure of the shale sample, and conducts a pore 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. The temperature and pressure control of the sample from which mercury has been discharged is improved, the required temperature and pressure parameters are set by itself, and the sample is improved under these conditions. After the improvement is completed, a pore effectiveness test is conducted on the sample to obtain the improved cumulative ink-bottle pore volume ratio (CVRr) and the improved pore structure complexity index (PSCr).

[0135] In the above step S3, the device calculates the pore effective index based on the pore structure complexity index and the cumulative ink-bottle pore volume ratio, and calculates the improved pore effective 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 pores of the shale sample. The larger the CVR value, the more ineffective space occupied by ink-bottle pores in the pore space, and the lower the pore effectiveness. PSC measures the complexity of the pore structure in terms of the size distribution and connectivity of the pores. The higher the PSC value, the more complex the pore structure, which will also lead to a decrease in pore effectiveness.

[0137] Define the pore effective index (PEI): The pore effective index (PEI) is a parameter that comprehensively considers the cumulative ink-bottle pore volume ratio (CVR) and the pore structure complexity index (PSC), and is used to quantitatively characterize the pore effectiveness of shale samples. The calculation method is as follows:

[0138]

[0139] Define the pore effectiveness stability index (PEIr), that is, the improved pore effective index: The pore effectiveness stability index (PEIr) is used to measure the change degree of the pore effectiveness of shale samples under specific temperature and pressure modification conditions. When PEIr is greater than PEI, it indicates that this temperature and pressure condition has an improving effect on the pore effectiveness of shale. If PEIr is less than PEI, the pore effectiveness becomes worse. The calculation method is as follows:

[0140]

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

[0142] Perform a hierarchical evaluation of the pore effectiveness of the shale sample based on the first comparison result between the pore effective index and the preset interval value, and the second comparison result between the pore effective index and the improved pore effective index, to obtain the hierarchical evaluation result of the shale pore effectiveness. The pore effectiveness of the shale sample is hierarchically evaluated through the pore effective index (PEI) and the pore effectiveness stability index (PEIr). The evaluation method is shown in Table 1:

[0143] Table 1

[0144]

[0145]

[0146] The shale pore effectiveness evaluation method provided by the embodiments of the present invention has the following beneficial technical effects:

[0147] A variety of technical means are comprehensively used, such as cyclic mercury injection testing, digital core, artificial core, and temperature and pressure transformation testing, etc. The pore effectiveness of shale is evaluated from multiple perspectives, and the pore effectiveness value can be quantitatively given, making up for the deficiencies of the existing technology in quantitative analysis. By comprehensively considering multiple factors, such as the cumulative inkbottle pore volume ratio (CVR) obtained from cyclic mercury injection testing and the pore structure complexity index (PSC) obtained from digital core and artificial core analysis testing, as well as the results of temperature and pressure transformation testing, the pore effective index (PEI) and the pore effectiveness stability index (PEIr) are calculated, improving the accuracy of pore effectiveness evaluation. At the same time, the method of the present invention further improves the shale pore evaluation technical system, provides new methods and ideas for related research and practice, and helps to promote the development of shale pore evaluation technology.

[0148] The shale pore effectiveness evaluation method provided by the embodiments of the present invention performs a pore effectiveness test on the shale sample to obtain the pore structure complexity index and the cumulative inkbottle pore volume ratio; performs temperature and pressure improvement on the shale sample, and performs a pore effectiveness test on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative inkbottle pore volume ratio; calculates the pore effective index according to the pore structure complexity index and the cumulative inkbottle pore volume ratio, and calculates the improved pore effective index according to the improved pore structure complexity index and the improved cumulative inkbottle pore volume ratio; performs a pore effectiveness evaluation on the shale sample according to the pore effective index and the improved pore effective index to obtain the shale pore effectiveness evaluation result, and can accurately and quantitatively evaluate the shale pore effectiveness.

[0149] Further, performing a pore effectiveness test on the shale sample to obtain the pore structure complexity index includes:

[0150] Obtain the total number of pores, the number of pores with pore sizes smaller than the average pore size, the number of pores with pore sizes larger than the average pore size, and the number of connected pores by using the digital core test sample; it can be described with reference to the above embodiments and will not be elaborated here.

[0151] Calculate the pore structure complexity index according to the first ratio of the number of pores with pore sizes smaller than the average pore size to the total number of pores, the second ratio of the number of pores with pore sizes larger than the average pore size to the total number of pores, the third ratio of the number of connected pores to the total number of pores, and the preset weights corresponding to each ratio respectively; it can be described with reference to the above embodiments and will not be elaborated here.

[0152] Furthermore, determining the preset weights corresponding to each ratio respectively includes:

[0153] Use artificial cores to change the number of pores with pore sizes smaller than the average pore size, the number of pores with pore sizes larger than the average pore size, and the number of connected pores respectively for multiple times; it can be described with reference to the above embodiments and will not be elaborated here.

[0154] According to the number of pores with pore sizes smaller than the average pore size after each change, the average value of the number of pores with pore sizes smaller than the average pore size after multiple changes, and the first changed pore structure complexity index calculated corresponding to the number of pores with pore sizes smaller than the average pore size after each change, and the average value of the first changed pore structure complexity index after multiple changes, calculate the first correlation coefficient between the number of pores with pore sizes smaller than the average pore size and the pore structure complexity index; it can be described with reference to the above embodiments and will not be elaborated here.

[0155] Determine the first preset weight corresponding to the first ratio according to the numerical interval where the first correlation coefficient is located; it can be described with reference to the above embodiments and will not be elaborated here.

[0156] According to the number of pores with pore sizes larger than the average pore size after each change, the average value of the number of pores with pore sizes larger than the average pore size after multiple changes, and the second changed pore structure complexity index calculated corresponding to the number of pores with pore sizes larger than the average pore size after each change, and the average value of the second changed pore structure complexity index after multiple changes, calculate the second correlation coefficient between the number of pores with pore sizes larger than the average pore size and the pore structure complexity index; it can be described with reference to the above embodiments and will not be elaborated here.

[0157] Determine the second preset weight corresponding to the second ratio according to the numerical interval where the second correlation coefficient is located; it can be described with reference to the above embodiments and will not be elaborated here.

[0158] According to the number of connected pores after each change, the average value of the number of connected pores after multiple changes, and the third changed pore structure complexity index calculated corresponding to the number of connected pores after each change, the average value of the third changed 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; it can be described with reference to the above embodiments and will not be elaborated here.

[0159] According to the numerical interval where the third correlation coefficient is located, a third preset weight corresponding to the third ratio is determined. It can be described with reference to the above embodiments and will not be elaborated here.

[0160] Further, a pore effectiveness test is performed on the shale sample to obtain the cumulative inkbottle pore volume ratio, including:

[0161] The inkbottle pore volume corresponding to each mercury injection cycle test is obtained, and the ratio of the sum of the inkbottle pore volumes of all times to the total mercury injection volume of the first mercury injection cycle test is used as the cumulative inkbottle pore volume ratio. It can be described with reference to the above embodiments and will not be elaborated here.

[0162] Further, the obtaining of the inkbottle pore volume corresponding to each mercury injection cycle test includes:

[0163] For the first mercury injection cycle test, the inkbottle pore volume is calculated according to the following expression:

[0164]

[0165] where E(r e )1 represents the mercury withdrawal amount corresponding to the pore throat radius r obtained without correction in the first cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i )1 represents the mercury injection amount corresponding to the pore throat radius r obtained with correction in the first cycle with the corrected pore throat radius r during mercury injection as the independent variable; it can be described with reference to the above embodiments and will not be elaborated here.

[0166] For non - first mercury injection cycle tests, the inkbottle pore volume is calculated according to the following expression:

[0167]

[0168] where E′(r e ) i represents the mercury withdrawal amount corresponding to the pore throat radius r obtained with correction in the i - th cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, I′(r i ) i represents the mercury injection amount corresponding to the pore throat radius r obtained with correction in the i - th cycle with the corrected pore throat radius r during mercury injection as the independent variable, r iminis the minimum value of the pore throat radius r in the mercury intrusion curve, r emax is the maximum value of the pore throat radius r in the mercury extrusion curve. It can be described with reference to the above embodiments and will not be elaborated here.

[0169] Further, the pore effectiveness evaluation of the shale sample is performed according to the pore effective index and the improved pore effective index, and the shale pore effectiveness evaluation result is obtained, including:

[0170] The pore effectiveness grading evaluation of the shale sample is performed according to the first comparison result between the pore effective index and the preset interval value, and the second comparison result between the pore effective index and the improved pore effective index, and the shale pore effectiveness grading evaluation result is obtained. It can be described with reference to the above embodiments and will not be elaborated here.

[0171] Figure 4 is a schematic structural diagram of a shale pore effectiveness evaluation processing device provided by an embodiment of the present invention. As Figure 4 shown, the shale pore effectiveness evaluation processing device provided by the embodiment of the present invention includes an acquisition unit 401, a test unit 402, a calculation unit 403, and an evaluation unit 404, where:

[0172] The acquisition unit 401 is used to perform pore effectiveness tests on the shale sample to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio; the test unit 402 is used to perform temperature and pressure improvement on the shale sample and perform pore effectiveness tests on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; the calculation unit 403 is used to calculate the pore effective index according to the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculate the improved pore effective index according to the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; the evaluation unit 404 is used to perform pore effectiveness evaluation on the shale sample according to the pore effective index and the improved pore effective index to obtain the shale pore effectiveness evaluation result.

[0173] Specifically, the acquisition unit 401 in the device is used to perform pore effectiveness tests on shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; the test unit 402 is used to perform temperature and pressure improvement on the shale samples, and perform pore effectiveness tests on the improved shale samples 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 pore effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculate an improved pore 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 pore effectiveness of the shale samples based on the pore effectiveness index and the improved pore effectiveness index to obtain a shale pore effectiveness evaluation result.

[0174] The shale pore effectiveness evaluation processing device provided in the embodiment of the present invention performs pore effectiveness tests on shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; performs temperature and pressure improvement on the shale samples, and performs pore effectiveness tests on the improved shale samples again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; calculates a pore effectiveness index based on the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculates an improved pore effectiveness index based on the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; evaluates the pore effectiveness of the shale samples based on the pore effectiveness index and the improved pore effectiveness index to obtain a shale pore effectiveness evaluation result, and can accurately and quantitatively evaluate the pore effectiveness of shale.

[0175] The embodiments of the shale pore effectiveness evaluation processing device provided in the embodiments of the present invention can specifically be used to execute the processing procedures of the above method embodiments, and their functions will not be described in detail here. Reference can be made to the detailed descriptions of the above method embodiments.

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

[0177] Perform pore effectiveness tests on shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio;

[0178] Perform temperature and pressure improvement on the shale samples, and perform pore effectiveness tests on the improved shale samples again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio;

[0179] The pore effective 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 pore effectiveness of the shale sample is evaluated according to the pore effective index and the improved pore effective index, and the shale pore effectiveness evaluation result is obtained.

[0181] This embodiment discloses a computer program product, which includes a computer program. When the computer program is executed by a processor, the following method is implemented:

[0182] The pore effectiveness test is carried out on the shale sample to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0183] The temperature and pressure of the shale sample are improved, and the pore effectiveness test is carried out on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio;

[0184] The pore effective 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 pore effectiveness of the shale sample is evaluated according to the pore effective index and the improved pore effective index, and the shale pore effectiveness evaluation result is obtained.

[0186] This embodiment provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the following method is implemented:

[0187] The pore effectiveness test is carried out on the shale sample to obtain the pore structure complexity index and the cumulative ink bottle pore volume ratio;

[0188] The temperature and pressure of the shale sample are improved, and the pore effectiveness test is carried out on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio;

[0189] The pore effective 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 pore effectiveness of the shale sample is evaluated according to the pore effective index and the improved pore effective index, and the shale pore effectiveness evaluation result is obtained.

[0191] Compared with the technical solutions in the prior art, the shale pore effectiveness evaluation method provided in the embodiments of the present invention performs pore effectiveness tests on shale samples to obtain a pore structure complexity index and a cumulative ink bottle pore volume ratio; performs temperature and pressure improvement on the shale samples, and performs pore effectiveness tests on the improved shale samples again to obtain an improved pore structure complexity index and an improved cumulative ink bottle pore volume ratio; calculates a pore effectiveness index according to the pore structure complexity index and the cumulative ink bottle pore volume ratio, and calculates an improved pore effectiveness index according to the improved pore structure complexity index and the improved cumulative ink bottle pore volume ratio; evaluates the pore effectiveness of the shale samples according to the pore effectiveness index and the improved pore effectiveness index to obtain a shale pore effectiveness evaluation result, which can accurately and quantitatively evaluate the shale pore effectiveness.

[0192] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. 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. Moreover, the present invention can take the form of a computer program product implemented 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] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to the embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram, and the combination of flows and / or blocks in the flowchart and / or block diagram, can be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate means for realizing the functions specified in Figure 1 one or more of the flows Figure 1 or multiple flows and / or blocks

[0194] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including instruction means, and the instruction means realizes the functions specified in Figure 1 one or more of the flows Figure 1 or multiple flows and / or blocks

[0195] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus, so that a series of operation steps are performed on the computer or other programmable apparatus to produce a computer-implemented process, thereby providing instructions for implementing the functions specified in one process Figure 1 one process or more processes and / or blocks Figure 1 steps of the functions specified in one block or more blocks.

[0196] In the description of the present specification, the description with reference to the terms "one embodiment", "a specific embodiment", "some embodiments", "for example", "example", "specific example", or "some examples", etc. means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0197] The above specific embodiments have further elaborated on the objectives, technical solutions and beneficial effects of the present invention. It should be understood that the above are only specific embodiments of the present invention and are not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for evaluating the effectiveness of shale pores, characterized in that, Including: Carry out pore effectiveness tests on shale samples to obtain a pore structure complexity index and a cumulative volume ratio of ink-bottle pores; Carry out temperature and pressure improvement on the shale samples, and carry out pore effectiveness tests on the improved shale samples again to obtain an improved pore structure complexity index and an improved cumulative volume ratio of ink-bottle pores; Calculate a pore effectiveness index based on the pore structure complexity index and the cumulative volume ratio of ink-bottle pores, and calculate an improved pore effectiveness index based on the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores; Carry out pore effectiveness evaluation on the shale samples according to the pore effectiveness index and the improved pore effectiveness index to obtain a shale pore effectiveness evaluation result.

2. The shale pore effectiveness evaluation processing method according to claim 1, characterized in that Carry out pore effectiveness tests on shale samples to obtain a pore structure complexity index, including: Use digital core test samples to obtain the total number of pores, the number of pores with pore diameter less than the average pore diameter, the number of pores with pore diameter greater than the average pore diameter, and the number of connected pores; Calculate the pore structure complexity index according to the first ratio of the number of pores with pore diameter less than the average pore diameter to the total number of pores, the second ratio of the number of pores with pore diameter greater than the average pore diameter to the total number of pores, the third ratio of the number of connected pores to the total number of pores, and the preset weights corresponding to each ratio respectively.

3. The shale pore effectiveness evaluation processing method according to claim 2, wherein Determine the preset weights corresponding to each ratio respectively, including: Use artificial cores to change the number of pores with pore diameter less than the average pore diameter, the number of pores with pore diameter greater than the average pore diameter, and the number of connected pores respectively multiple times; According to the number of pores with pore diameter less than the average pore diameter after each change, the average value of the number of pores with pore diameter less than the average pore diameter after multiple changes, and the first changed pore structure complexity index calculated corresponding to the number of pores with pore diameter less than the average pore diameter after each change, and the average value of the first changed pore structure complexity index after multiple changes, calculate the first correlation coefficient between the number of pores with pore diameter less than the average pore diameter and the pore structure complexity index; Determine the first preset weight corresponding to the first ratio according to the numerical interval where the first correlation coefficient is located; According to the number of pores with pore diameter greater than the average pore diameter after each change, the average value of the number of pores with pore diameter greater than the average pore diameter after multiple changes, and the second changed pore structure complexity index calculated corresponding to the number of pores with pore diameter greater than the average pore diameter after each change, and the average value of the second changed pore structure complexity index after multiple changes, calculate the second correlation coefficient between the number of pores with pore diameter greater than the average pore diameter and the pore structure complexity index; Determine the second preset weight corresponding to the second ratio according to the numerical interval where the second correlation coefficient is located; According to the number of connected pores after each change, the average value of the number of connected pores after multiple changes, and the third changed pore structure complexity index calculated corresponding to the number of connected pores after each change, and the average value of the third changed pore structure complexity index after multiple changes, calculate the third correlation coefficient between the number of connected pores and the pore structure complexity index; Determine the third preset weight corresponding to the third ratio according to the numerical interval where the third correlation coefficient is located.

4. The shale pore effectiveness evaluation processing method according to claim 1, wherein Perform pore effectiveness tests on shale samples to obtain the cumulative volume ratio of ink-bottle pores, including: Obtain the volume of ink-bottle pores corresponding to each mercury intrusion cycle test, and use the ratio of the sum of the volumes of ink-bottle pores for all cycles to the total mercury intrusion volume in the first mercury intrusion cycle test as the cumulative volume ratio of ink-bottle pores.

5. The shale pore effectiveness evaluation processing method according to claim 4, characterized in that The obtaining of the volume of ink-bottle pores corresponding to each mercury intrusion cycle test includes: For the first mercury intrusion cycle test, calculate the volume of ink-bottle pores according to the following expression: Among them, E(r e )1 represents the mercury withdrawal volume corresponding to the pore throat radius r obtained without correction in the first cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i )1 represents the mercury injection volume corresponding to the pore throat radius r obtained with correction in the first cycle with the corrected pore throat radius r during mercury injection as the independent variable; For non-first mercury intrusion cycle tests, calculate the volume of ink-bottle pores according to the following expression: where E′(r e ) i represents the mercury withdrawal volume corresponding to the corrected pore throat radius r obtained in the i-th cycle with the corrected pore throat radius r during mercury withdrawal as the independent variable, and I′(r i ) i represents the mercury injection volume corresponding to the corrected pore throat radius r obtained in the i-th cycle with the corrected pore throat radius r during mercury injection as the independent variable, r imin is the minimum value of the pore throat radius r in the mercury injection curve, and r emax is the maximum value of the pore throat radius r in the mercury withdrawal curve.

6. The shale pore effectiveness evaluation and treatment method according to any one of claims 1 to 5, characterized in that, The evaluation of the pore effectiveness of the shale sample according to the pore effective index and the improved pore effective index to obtain the shale pore effectiveness evaluation result includes: Perform a hierarchical evaluation of the pore effectiveness of the shale sample according to the first comparison result between the pore effective index and the preset interval value, and the second comparison result between the pore effective index and the improved pore effective index, to obtain the shale pore effectiveness hierarchical evaluation result.

7. A device for evaluating the effectiveness of shale pores, characterized in that, Includes: An acquisition unit for performing pore effectiveness tests on shale samples to obtain the pore structure complexity index and the cumulative volume ratio of ink-bottle pores; A test unit for performing temperature and pressure improvement on the shale sample and then performing pore effectiveness tests on the improved shale sample again to obtain the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores; A calculation unit for calculating the pore effective index according to the pore structure complexity index and the cumulative volume ratio of ink-bottle pores, and calculating the improved pore effective index according to the improved pore structure complexity index and the improved cumulative volume ratio of ink-bottle pores; An evaluation unit for evaluating the pore effectiveness of the shale sample according to the pore effective index and the improved pore effective index to obtain the shale pore effectiveness evaluation result.

8. A computer device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method according to any one of claims 1 to 6.

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

10. A computer program product, characterized in that, The computer program product includes a computer program, and when the computer program is executed by the processor, it implements the method according to any one of claims 1 to 6.

Images