Correction method, device and equipment for rock high-pressure mercury injection data, medium and product

Abstract

The invention discloses a rock high-pressure mercury injection data correction method, device, equipment, medium and product, and relates to the field of oil and gas reservoir characterization, the method comprises the following steps: under the same high-pressure mercury injection experiment condition, obtaining target rock to-be-corrected data and high-pressure mercury injection data of a preset sample; calculating an elastic deformation quantity based on a preset sample volume modulus, and deducting the elastic deformation quantity to obtain blank error data; preliminarily correcting the target rock data by using the blank error data; determining deformation data of the target rock by calculating a volume change rate between adjacent pressures and screening a minimum value; deducting the deformation data of the target rock to obtain corrected high-pressure mercury injection data. According to the method, through two-stage correction combining the blank error data and the deformation data of the target rock, the interference between the blank effect system error of the instrument and the compression of the target rock is synchronously eliminated, the problem of incomplete correction of an existing method is solved, and the precision and reliability of the high-pressure mercury injection data are remarkably improved.

Description

Technical Field

[0001] This application relates to the field of oil and gas reservoir characterization, and in particular to a method, apparatus, equipment, medium, and product for correcting high-pressure mercury intrusion data of rocks. Background Technology

[0002] Oil and gas reservoirs generally exhibit strong heterogeneity and multi-scale pore structures, and the size, distribution, and connectivity of their pore throats directly affect fluid occurrence and seepage characteristics. High-pressure mercury intrusion porosimetry (HIP) is widely used for quantitative characterization of the pore structure of sandstone, carbonate, shale, and coal reservoirs because it can cover the pore throat range from micrometers to nanometers.

[0003] However, under high pressure conditions of 200 MPa–400 MPa, both rock samples and experimental apparatus undergo compression, making it difficult to distinguish the compressed volume from the volume of mercury entering the pore throat, leading to misinterpretations of the results. Related techniques for correcting mercury intrusion porosimetry data often utilize blank experiments or constant compressibility coefficients for simplified correction, but these methods have limitations. Blank experiments cannot correct for sample compression and do not consider the influence of the sample's own volume. Compressibility coefficient correction does not account for the instrument blank effect, which is particularly significant in dense samples. These issues limit the accuracy and reliability of mercury intrusion porosimetry results.

[0004] Therefore, there is an urgent need for a correction method for high-pressure mercury intrusion data in rocks to solve the technical problem of accurately separating the sample compression effect while deducting systematic errors, thereby improving the accuracy and reliability of quantitative characterization of reservoir pore structure. Summary of the Invention

[0005] The purpose of this application is to provide a method, apparatus, equipment, medium, and product for correcting high-pressure mercury intrusion data of rocks. It can obtain the true mercury intrusion and retreat curve data of rocks by correcting instrument blank errors and sample compression deformation, which solves the problems of insufficient and incomplete existing correction methods, provides a new perspective for mercury intrusion data correction, and effectively improves the accuracy and reliability of the data.

[0006] To achieve the above objectives, this application provides the following solution: In a first aspect, this application provides a method for correcting high-pressure mercury intrusion data in rocks, including: Under the same high-pressure mercury intrusion test conditions, high-pressure mercury intrusion data to be corrected for the target rock and high-pressure mercury intrusion data for a preset sample are obtained; the high-pressure mercury intrusion data includes the cumulative mercury intrusion amount under different pressures; the preset sample is a solid with a smooth surface, no interconnected pores, a known bulk modulus, and only undergoing elastic deformation; Based on the bulk modulus of the preset sample, the elastic deformation of the preset sample under different pressures is calculated, and the elastic deformation under the corresponding pressure is subtracted from the high pressure mercury intrusion data of the preset sample to obtain blank error data under different pressures. Based on the high-pressure mercury intrusion data to be corrected for the target rock and the blank error data under different pressures, preliminary corrected high-pressure mercury intrusion data for the target rock are obtained. Based on the preliminary corrected high-pressure mercury intrusion data of the target rock, the volume change rate between adjacent pressures is calculated, and the deformation data of the target rock under different pressures is determined based on the minimum volume change rate. Subtract the deformation data of the target rock under the corresponding pressure from the preliminary corrected high-pressure mercury intrusion data of the target rock to obtain the corrected high-pressure mercury intrusion data of the target rock.

[0007] Secondly, this application provides a correction device for high-pressure mercury intrusion data in rocks, comprising: The data acquisition module is used to acquire, under the same high-pressure mercury intrusion test conditions, the high-pressure mercury intrusion data to be corrected for the target rock and the high-pressure mercury intrusion data for a preset sample; the high-pressure mercury intrusion data includes the cumulative mercury intrusion amount under different pressures; the preset sample is a solid with a smooth surface, no interconnected pores, a known bulk modulus, and only undergoing elastic deformation; The blank error calculation module is used to calculate the elastic deformation of the preset sample under different pressures based on the bulk modulus of the preset sample, and to subtract the elastic deformation under the corresponding pressure from the high pressure mercury intrusion data of the preset sample to obtain blank error data under different pressures. The preliminary correction module is used to obtain preliminary corrected high-pressure mercury intrusion data of the target rock based on the high-pressure mercury intrusion data to be corrected and blank error data under different pressures. The sample deformation acquisition module is used to calculate the volume change rate between adjacent pressures based on the preliminary corrected high-pressure mercury intrusion data of the target rock, and to determine the deformation data of the target rock under different pressures based on the minimum volume change rate. The final correction module is used to subtract the deformation data of the target rock under the corresponding pressure from the preliminary corrected high-pressure mercury intrusion data of the target rock to obtain the corrected high-pressure mercury intrusion data of the target rock.

[0008] Thirdly, this application provides a computer device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the method for correcting rock high-pressure mercury intrusion data as described above.

[0009] Fourthly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method for correcting high-pressure mercury intrusion data in rocks as described above.

[0010] Fifthly, this application provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the method for correcting rock high-pressure mercury intrusion data as described above.

[0011] According to the specific embodiments provided in this application, this application has the following technical effects: This application provides a method, apparatus, equipment, medium, and product for correcting high-pressure mercury intrusion (HPI) data of rocks. By acquiring the HPI data of the target rock to be corrected and the HPI data of a preset sample under the same HPI experimental conditions, it solves the problem of lacking benchmark data for simultaneous calibration of instrument blank effect and sample deformation in existing technologies, and establishes a foundation for separating experimental system errors from sample interference. By calculating the elastic deformation of the preset sample under different pressures based on its bulk modulus, and subtracting the elastic deformation from the HPI data of the preset sample, blank error data is obtained. This solves the technical deficiency of traditional blank experiments in quantifying the elastic compression volume of the instrument under high pressure, and achieves accurate quantification and subtraction of the instrument blank effect. Furthermore, by utilizing the blank error data… Preliminary corrections were made to the target rock data to be corrected, resulting in preliminarily corrected high-pressure mercury intrusion (HPI) data. This solved the problem of instrument system errors contaminating the target rock data and achieved the first correction of HPI data. By calculating the volume change rate between adjacent pressures based on the preliminary correction data and selecting the minimum value to determine the deformation data of the target rock under different pressures, the technical challenge of dynamically capturing the nonlinear deformation of the sample under high pressure was solved, achieving accurate identification and quantification of the compression effect of the target rock. By subtracting the deformation data of the target rock from the preliminary correction data, the final corrected HPI data was obtained. This solved the limitation of existing methods that cannot simultaneously eliminate two types of errors, achieving high-precision restoration of HPI data and providing reliable data support for the characterization of the pore structure of tight reservoirs. Attached Figure Description

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

[0013] Figure 1 This is an application environment diagram of a method for correcting high-pressure mercury intrusion data in rocks according to an embodiment of this application; Figure 2 A flowchart illustrating a method for correcting high-pressure mercury intrusion data in rocks, provided as an embodiment of this application; Figure 3 A flowchart illustrating a method for correcting high-pressure mercury intrusion data in rocks, provided as another embodiment of this application; Figure 4 The various oil and gas reservoir rocks provided in one embodiment of this application have been subjected to... Figure 3 A comparison of the mercury ingress and egress curves before and after each step of the correction process; among them, Figure 4 (a) is a comparison diagram of the advance and retreat of mercury curves in sandstone; Figure 4 (b) is a comparison diagram of the advance and retreat of mercury curves in carbonate rocks; Figure 4 (c) is a comparison diagram of the advance and retreat curves of oil shale; Figure 4 (d) is a comparison diagram of the mercury advance and retreat curves of gas shale; 4(e) is a comparison diagram of the mercury advance and retreat curves of coal and rock. Figure 5 A schematic diagram of the blank error curve of a high-pressure mercury intrusion test on a quartz sample provided in an embodiment of this application; Figure 6 A schematic diagram of the functional modules of a device for correcting high-pressure mercury intrusion data in rocks, provided in an embodiment of this application; Figure 7 This is a schematic diagram of the structure of a computer device provided in an embodiment of this application. Detailed Implementation

[0014] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0015] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0016] The method for correcting high-pressure mercury intrusion data in rocks provided in this application can be applied to, for example... Figure 1The application environment shown is illustrated. Terminal 101 communicates with server 102 via a network. A data storage system can store the data that server 102 needs to process. The data storage system can be set up independently, integrated into server 102, or placed in the cloud or on another server. Terminal 101 can send the high-pressure mercury intrusion (HPI) data of the target rock to be corrected and the HPI data of a preset sample under the same high-pressure mercury intrusion experimental conditions to server 102. Upon receiving the HPI data of the target rock to be corrected and the HPI data of the preset sample under the same high-pressure mercury intrusion experimental conditions, server 102 calculates the elastic deformation of the preset sample under different pressures based on the bulk modulus of the preset sample, and subtracts the elastic deformation under the corresponding pressure from the HPI data of the preset sample to obtain blank error data under different pressures. Based on the HPI data of the target rock to be corrected and the blank error data under different pressures, preliminary corrected HPI data of the target rock is obtained. Based on the preliminary corrected HPI data of the target rock, the volume change rate between adjacent pressures is calculated, and the deformation data of the target rock under different pressures is determined based on the minimum volume change rate. The preliminary corrected HPI data of the target rock is subtracted from the deformation data of the target rock under the corresponding pressure to obtain the corrected HPI data of the target rock. Server 102 can then feed back the obtained corrected HPI data of the target rock to terminal 101. Furthermore, in some embodiments, the method for correcting high-pressure mercury intrusion data of rocks can also be implemented separately by the server 102 or the terminal 101. For example, the terminal 101 can directly perform correction processing on the high-pressure mercury intrusion data of the target rock to be corrected and the high-pressure mercury intrusion data of the preset sample under the same high-pressure mercury intrusion experimental conditions. Alternatively, the server 102 can obtain the high-pressure mercury intrusion data of the target rock to be corrected and the high-pressure mercury intrusion data of the preset sample under the same high-pressure mercury intrusion experimental conditions from the data storage system, and perform correction processing on the high-pressure mercury intrusion data of the target rock to be corrected and the high-pressure mercury intrusion data of the preset sample under the same high-pressure mercury intrusion experimental conditions.

[0017] The terminal 101 can be, but is not limited to, various desktop computers, laptops, smartphones, tablets, IoT devices, and portable wearable devices. IoT devices can include smart speakers, smart TVs, smart air conditioners, and smart in-vehicle devices. Portable wearable devices can include smartwatches, smart bracelets, and head-mounted devices. The server 102 can be implemented using a standalone server or a server cluster composed of multiple servers, or it can be a cloud server.

[0018] In one exemplary embodiment, such as Figure 2As shown, a method for correcting high-pressure mercury intrusion data in rocks is provided. This method is executed by a computer device, specifically a terminal or server, or both. In this embodiment, the method is applied to... Figure 1 Taking server 102 as an example, the explanation includes the following steps 201 to 205. Wherein: Step 201: Under the same high-pressure mercury intrusion test conditions, obtain the high-pressure mercury intrusion data to be corrected for the target rock and the high-pressure mercury intrusion data for the preset sample; the high-pressure mercury intrusion data includes the cumulative mercury intrusion amount under different pressures; the preset sample is a solid with a smooth surface, no interconnected pores, a known bulk modulus, and only undergoing elastic deformation.

[0019] Step 202: Based on the bulk modulus of the preset sample, calculate the elastic deformation of the preset sample under different pressures, and subtract the elastic deformation under the corresponding pressure from the high pressure mercury intrusion data of the preset sample to obtain blank error data under different pressures.

[0020] Step 203: Based on the high-pressure mercury intrusion data of the target rock to be corrected and the blank error data under different pressures, the preliminary corrected high-pressure mercury intrusion data of the target rock are obtained.

[0021] Step 204: Based on the preliminary corrected high-pressure mercury intrusion data of the target rock, calculate the volume change rate between adjacent pressures, and determine the deformation data of the target rock under different pressures based on the minimum volume change rate.

[0022] Step 205: Subtract the deformation data of the target rock under the corresponding pressure from the preliminary corrected high-pressure mercury intrusion data of the target rock to obtain the corrected high-pressure mercury intrusion data of the target rock.

[0023] By implementing steps 201 to 205 above, this application can simultaneously eliminate the coupling interference of two error sources in high-pressure mercury intrusion porosimetry experiments: instrument blank effect and sample nonlinear deformation. By quantifying the instrument system error by pre-setting the sample elastic deformation, and combining the principle of minimizing the volume change rate, the compression effect of the target rock can be accurately identified. Ultimately, high-precision correction of high-pressure mercury intrusion porosimetry data can be achieved, providing a reliable data basis for the quantitative characterization of pore structure parameters (such as porosity, permeability, and pore size distribution) of tight reservoirs, and effectively improving the accuracy of oil and gas reservoir evaluation.

[0024] In another exemplary embodiment of this application, the elastic deformation of the preset sample under different pressures is calculated based on the bulk modulus of the preset sample, specifically including: .

[0025] in, This represents the elastic deformation of the preset sample, in cm. 3 ; This indicates the initial volume of the preset sample, in cm³. 3 K represents the pre-defined bulk modulus of the sample, in GPa. This indicates the change in pressure, in MPa.

[0026] In another exemplary embodiment of this application, preliminary corrected high-pressure mercury intrusion data of the target rock are obtained based on the high-pressure mercury intrusion data to be corrected and blank error data under different pressures, specifically including: If the volume of the target rock is equal to the volume of the preset sample, then the blank error data at the corresponding pressure is subtracted from the high-pressure mercury intrusion data of the target rock to be corrected, so as to obtain the preliminary corrected high-pressure mercury intrusion data of the target rock.

[0027] If the volume of the target rock is not equal to the volume of the preset sample, a preliminary correction is made using the following formula: .

[0028] in, This indicates the preliminary corrected high-pressure mercury intrusion data for the target rock; This represents the high-pressure mercury intrusion data of the target rock that needs correction. Indicates blank error data; Indicates the volume of the preset sample; The volume of the target rock is represented by P; the pressure is represented by MPa. This indicates the bulk modulus of mercury, measured in gigabytes of water (GPa).

[0029] In another exemplary embodiment of this application, the method for calculating the bulk modulus of mercury specifically includes: Under the same high-pressure mercury intrusion test conditions as the target rock, a blank high-pressure mercury intrusion test was conducted in the air dilatometer without placing any solid sample inside, and blank high-pressure mercury intrusion data were obtained.

[0030] Subtract the blank error data at the corresponding pressure from the blank high-pressure mercury intrusion data to obtain the mercury deformation data.

[0031] Based on the deformation data of mercury, the bulk modulus of mercury is calculated using the following formula: .

[0032] in, This represents the rate of change of mercury pressure-volume, calculated based on mercury deformation data.

[0033] In another exemplary embodiment of this application, before the high-pressure mercury intrusion test, the target rock and the preset sample are dried according to preset drying conditions; the preset drying conditions include drying at a temperature of 70°C for 24 hours.

[0034] In another exemplary embodiment of this application, the preset sample includes quartz, high-purity aluminum, or high-purity copper. The preset sample preferably uses a standard substance with a publicly known and definitive bulk modulus. Specifically, materials including, but not limited to, α-quartz (bulk modulus approximately 37 GPa), high-purity aluminum, or high-purity copper can be used. Among these, α-quartz is considered the most suitable preferred option due to its recognized and stable bulk modulus and its mineralogical correlation with rock samples. High-purity aluminum and high-purity copper, when meeting purity requirements (e.g., purity higher than 99.99%), also have referable and definitive bulk modulus values ​​and can be used as applicable standards.

[0035] This application provides a method for correcting high-pressure mercury intrusion data of rocks, comprising the following steps: performing high-pressure mercury intrusion experiments on dried rock samples to obtain mercury intrusion and withdrawal curves of the samples; performing high-pressure mercury intrusion experiments on smooth and non-porous quartz samples, with all experimental conditions consistent with those of the rock samples; obtaining the change in quartz deformation with pressure under the experimental conditions based on the bulk modulus of quartz (37 GPa), subtracting the quartz deformation from the mercury intrusion and withdrawal curves of the quartz samples to obtain a blank error curve under the experimental conditions; subtracting the blank error curve from the pressure-cumulative mercury intrusion curve of the samples to complete the preliminary correction, where the mercury intrusion includes the sample deformation and the amount of actual intrusion into the pores; calculating the minimum value of the rate of change in the preliminary corrected curves and assuming that only deformation occurs here, subtracting the deformation curve from this curve to obtain the true mercury intrusion and withdrawal curves.

[0036] The following example illustrates this application using the correction process for specific high-pressure mercury intrusion data from rocks.

[0037] like Figure 3 As shown, a method for correcting high-pressure mercury intrusion data in rocks is provided, including the following steps: S1. The dried rock sample was subjected to a high-pressure mercury intrusion test to obtain the mercury intrusion and retraction curves of the sample.

[0038] A sample of a certain volume of the target rock is cut, such as the commonly used 1cm sample. 3 For volumetric samples, place them in an oven and dry them at 70℃ for at least 24 hours. The time and temperature can be adjusted appropriately according to the size of the target rock sample.

[0039] After drying, high-pressure mercury intrusion porosimetry (HIP) was performed on the samples using a Micromeritics AutoPore V 9620 mercury intrusion porosimeter. The appropriate dilatometer was selected based on the sample type or porosity of the target rock; for dense samples, a dilatometer with a finer capillary tube was preferred to improve resolution. For relatively dense samples (low porosity), such as shale samples, a 0.096 cm⁻¹ dilatometer could be used. 3 For stem volume dilatometers, 0.412 cm⁻¹ can be used for samples with high porosity, such as sandstone samples.3 The purpose of choosing a stem volume dilatometer is to improve resolution.

[0040] 0.096cm 3 The initial pressure of the dilatometer for stem volume measurement is typically set to 5 psi, 0.412 cm. 3 The initial pressure of the dilatometer for stem volume is typically set at 0.5 psi, with a maximum pressure of 60,000 psi. Regarding the choice of pressure gradient, according to the Washburn equation, pressure and pore size have an inverse functional relationship. Therefore, in low-pressure regions, smaller pressure changes correspond to larger pore size changes, while in high-pressure regions, even a significant increase in pressure will only result in a slight decrease in pore size. Based on this, to ensure a relatively uniform pore size distribution in experimental measurements, a smaller pressure gradient is typically chosen in low-pressure regions, while a larger pressure gradient is typically chosen in high-pressure regions.

[0041] The equilibration time is usually chosen to be 10-15 seconds to ensure that mercury can fully penetrate the sample.

[0042] After the experiment, the data was exported, and the pressure and cumulative mercury injection volumes were extracted. In the original data (the high-pressure mercury intrusion data of the target rock to be corrected), the pressure unit was psi, and the cumulative mercury injection volume was mL / g. The pressure was divided by 145 to convert the unit to MPa, and the cumulative mercury injection volume was multiplied by the sample mass to convert the unit to mL. These two columns of unit-converted data served as the original mercury injection and withdrawal data for the sample. The resulting curve is shown in the figure below. Figure 4 As shown.

[0043] S2. Perform high-pressure mercury intrusion tests on smooth and non-porous quartz samples, with all experimental conditions consistent with those of the rock samples.

[0044] Select a quartz standard sample, also at a depth of 1 cm. 3 For example, the standard sample must have a smooth surface and be free of pores. This requirement is to avoid any mercury ingress other than that caused by deformation. Additional mercury ingress caused by surface effects or inherent porosity will amplify the blank error, thus interfering with the calibration results.

[0045] High-pressure mercury intrusion porosimetry was also performed on the same sample, with all experimental conditions, including the selection of the dilatometer, the setting of the pressure gradient, and the equilibration time, kept consistent with those of the target rock sample. This was done to facilitate point-by-point calibration and to prevent additional errors introduced by inconsistent conditions.

[0046] Similarly, after the experiment, the data was exported and the columns for pressure and cumulative mercury injection were extracted. Unit conversion was also performed, and the pressure data points could be matched point-by-point with the pressure points of the target rock sample. Figure 5 A schematic diagram of the blank error curve for high-pressure mercury intrusion testing of quartz is shown.

[0047] S3. Based on the bulk modulus of quartz (37 GPa), the change in quartz deformation with pressure under these experimental conditions is obtained. The quartz deformation is subtracted from the mercury advance / retreat curve to obtain the blank error curve under these experimental conditions. Specifically, this includes: Extract the pressure data column from the original quartz or rock sample data, and calculate the quartz deformation corresponding to each pressure point according to the formula: .

[0048] in, This represents the elastic deformation of the preset sample, in cm. 3 ; This indicates the initial volume of the preset sample, in cm³. 3 K represents the pre-defined bulk modulus of the sample, in GPa; for quartz, it is 37 GPa. This indicates the change in pressure, in MPa.

[0049] By subtracting the calculated quartz deformation from the cumulative mercury ingress obtained by quartz mercury intrusion at each point, the blank error corresponding to each pressure point under the experimental conditions can be obtained. Figure 5 The blank error curve after subtracting the deformation is shown.

[0050] S4. Subtract the blank error curve from the pressure-cumulative mercury ingress curve of the sample to complete the initial correction. At this time, the mercury ingress includes the amount of sample deformation and the actual amount of mercury intrusion into the pores.

[0051] In this step, it is important to ensure that the volume of the sample used is consistent with that of the quartz standard. If they are consistent, the initial correction can be completed by subtracting the blank error curve from the pressure-cumulative mercury injection curve of the sample.

[0052] When the sample volume differs significantly from the quartz standard volume, further correction is required, based on the following formula: .

[0053] in, This indicates the preliminary corrected high-pressure mercury intrusion data for the target rock; This represents the high-pressure mercury intrusion data of the target rock that needs correction. Indicates blank error data; Indicates the volume of the preset sample, in mL; The volume of the target rock is expressed in mL; P represents pressure, in MPa. This indicates the bulk modulus of mercury, measured in gigabytes of water (GPa).

[0054] There is no known data available for the bulk modulus of mercury, so it needs to be calculated. This can be done by performing a blank experiment (without any sample in the dilatometer). The deformation curve of mercury can be obtained by subtracting the blank error from the cumulative mercury amount in the blank experiment point by point. The modulus can then be calculated using the following formula: .

[0055] in, V t For the volume of quartz, This represents the mercury pressure-volume change rate calculated based on mercury deformation data. It can be calculated from the mercury deformation curve, and the calculated bulk modulus of mercury is approximately 32 GPa.

[0056] Based on the above formula, substituting the parameters corresponding to different pressures yields a more accurate preliminary correction for the mercury injection volume. When the volume of quartz is equal to the volume of the sample, that is... V t = V s In this case, correction can be achieved simply by subtracting the blank error from the cumulative mercury amount in the sample point by point. Figure 4 The preliminarily corrected curves for five different types of rocks are shown. Among them, Figure 4 (a) is a comparison diagram of the advance and retreat of mercury curves in sandstone; Figure 4 (b) is a comparison diagram of the advance and retreat of mercury curves in carbonate rocks; Figure 4 (c) is a comparison diagram of the advance and retreat curves of oil shale; Figure 4 (d) is a comparison diagram of the mercury advance and retreat curves of gas shale; 4(e) is a comparison diagram of the mercury advance and retreat curves of coal rock.

[0057] S5. In the initially corrected curve, calculate the minimum rate of change for the entire pressure range (including mercury ingress and regress) and assume that only deformation occurs here, i.e., d V 1 / d P The minimum value of the curve can be obtained by calculating the rate of change between adjacent points since the data is discrete. Subtracting the deformation curve from the curve yields the true advance and retreat curve.

[0058] As an optional implementation, this correction step is based on certain assumptions: Since the cumulative mercury ingress includes the amount of mercury entering and exiting the pores and the sample deformation of the target rock, the cumulative mercury ingress is always greater than or equal to the sample deformation of the target rock, that is: .

[0059] in, , These represent the amount of mercury entering and exiting the pores and the sample deformation of the target rock, respectively. The minimum value obtained here is the upper limit of the sample deformation of the target rock as a function of pressure. Among the dozens or hundreds of pressure points set, there exists a certain pressure point where d... V p =0, meaning that no mercury entered or exited the pores at this pressure point, only the target rock sample deformed. Therefore, the following assumption is made: .

[0060] A deformation curve of the target rock can be calculated using this minimum value as the volume change rate across the entire pressure range. Subtracting the deformation corresponding to the target rock's deformation curve point by point from the initially corrected cumulative mercury ingress amount yields the true mercury ingress / regress curve of the target rock sample. The final corrected true mercury ingress / regress curve of the target rock sample is shown below. Figure 4 As shown.

[0061] The rationale for this hypothesis lies in the fact that in the high-pressure zone of mercury intrusion porosimetry (MIP) experiments, since capillary pressure is inversely proportional to pore size, the pore size difference corresponding to a small change in pressure is often extremely limited. When the pressure reaches hundreds of megapascals, the pore size interval between adjacent pressure points can be less than 0.1 nm, while the pore throat distribution of natural rocks cannot exist continuously at such a dense scale. Therefore, within the subdivided pressure intervals of the high-pressure stage, there must be several pressure points that do not correspond to any actually accessible pore throats. At this point, the amount of mercury intrusion is extremely small or even approaches zero, and the volume change mainly comes from the elastic compression of the rock skeleton under external high pressure. As the density of pressure sampling points increases, these intervals with "no corresponding pore size" are more easily captured, and their volume-pressure slope gradually approaches the intrinsic response of pure matrix deformation. Therefore, the d in the high-pressure section... V 1 / d P The minimum value is considered as the section dominated by matrix compression, and the volume compressibility coefficient of the rock is estimated accordingly. This is a reasonable approximation based on the pore size probability distribution characteristics and capillary mechanism.

[0062] It is important to note that calculating d V 1 / d P When the value is at its minimum, all data is used, including both mercury ingress and mercury exgress data. This is because the mercury exgress process is also a process of sample deformation recovery and mercury exiting the pores. When mercury cannot exit, the change in the cumulative amount of mercury ingress only reflects the deformation recovery.

[0063] Comparison before and after correction for five different types of rocks, including sandstone, carbonate rocks, oil / gas shale, and coal. Figure 4As shown, the most significant differences before and after correction are observed in carbonate rocks and gas shale, both of which are dense rocks. This is because the instrument's blank error interferes with the signal in the high-pressure zone, offsetting the lower mercury ingress in dense rocks. Therefore, the original data for dense rocks has a large error, making correction necessary. For high-porosity rocks, the instrument's own error has a smaller impact; the instrument error is insufficient to significantly interfere with the original results. In the final corrected curve, only coal shows a lower mercury ingress than the original data. This is because the coal matrix is ​​relatively soft, and compression contributes significantly to its mercury ingress during mercury intrusion intrusion (MII). This phenomenon is consistent with the general understanding in academia that coal matrices are soft and have a significant compression effect under high pressure. This result indirectly verifies that the method in this application can effectively identify and remove the contribution of sample compression to the mercury ingress signal.

[0064] In summary, this application provides a method for correcting high-pressure mercury intrusion data in rocks: First, samples dried at 70℃ for more than 24 hours are subjected to high-pressure mercury intrusion experiments. After the experiment, the raw data is exported, and the pressure and cumulative mercury intrusion data are extracted. The pressure unit (psi) in the raw data is divided by 145 to convert it to MPa, and the cumulative mercury intrusion unit (mL / g) is multiplied by the sample mass to convert it to mL. Further, high-pressure mercury intrusion experiments are performed on smooth, non-porous quartz samples, and all experimental conditions and data processing are consistent with the method used for the rock samples. Based on the bulk modulus of quartz (37 GPa), the change in quartz deformation with pressure under these experimental conditions is calculated. The cumulative mercury intrusion obtained from quartz mercury intrusion is subtracted point by point from the calculated quartz deformation to obtain the blank error corresponding to each pressure point under these experimental conditions. Further, the blank error is subtracted point by point from the cumulative mercury intrusion of the sample to complete the preliminary correction. The cumulative mercury intrusion at this point includes the sample deformation and the actual amount intruded into the pores. If the volume difference between the quartz and the sample is large, the interference caused by mercury compression needs to be considered and subtracted using a formula. Finally, in the preliminarily corrected curve, the minimum rate of change is calculated, that is, the d between adjacent points of the discrete data is calculated. V / d P The minimum value of the volume change rate can be used to calculate a deformation curve. By subtracting the deformation corresponding to the deformation curve from the initially corrected cumulative mercury intrusion amount point by point, the true mercury intrusion and retraction curve of the sample can be obtained. The above-mentioned method for correcting high-pressure mercury intrusion data in rocks solves the problems of insufficient and incomplete existing correction methods. It is more detailed in data processing, has higher reliability in principle, and can effectively improve the accuracy and reliability of the data.

[0065] This application also provides an application scenario in which the above-mentioned method for correcting high-pressure mercury intrusion porosimetry (HSI) data in rocks is applied. Specifically, the method for correcting HSI data in rocks provided in this embodiment can be applied in oil and gas reservoir research scenarios. Oil and gas reservoir research scenarios include field sample collection, laboratory experimental analysis, and comprehensive reservoir evaluation. Rock samples enter the laboratory experimental analysis stage from the field sample collection stage, where they undergo HSI experiments to obtain raw mercury intrusion data and other relevant information, and then proceed to subsequent steps involving data correction and interpretation. The method for correcting HSI data in rocks provided in this embodiment belongs to the data optimization processing step within the laboratory experimental analysis stage.

[0066] Specifically, in the laboratory experimental analysis stage, raw data is obtained after conducting high-pressure mercury intrusion porosimetry (HSI) experiments on collected rock samples. Due to interference from factors such as instrument blank effect and sample deformation during the experiment, the raw data contains errors. To correct this, the proposed method for high-pressure mercury intrusion porosimetry data is used. First, under the same HSI experimental conditions, the target rock (i.e., the collected oil and gas reservoir rock sample) is given HSI data to be corrected, along with the data for a pre-defined sample. Based on the bulk modulus of the pre-defined sample, its elastic deformation under different pressures is calculated. The elastic deformation under the corresponding pressure is then subtracted from the HSI data of the pre-defined sample to obtain blank error data. The blank error data is used to initially correct the target rock's HSI data to be corrected, resulting in initially corrected HSI data. Next, based on the initially corrected HSI data, the volume change rate between adjacent pressures is calculated, and the deformation data of the target rock sample under different pressures is determined based on the minimum volume change rate. Finally, the sample deformation data under the corresponding pressure is subtracted from the initially corrected HSI data to obtain the corrected HSI data. These corrected and accurate data enable a more precise analysis of the pore structure characteristics of oil and gas reservoirs, providing reliable data support for subsequent comprehensive reservoir evaluation.

[0067] Based on the same inventive concept, this application also provides a rock high-pressure mercury intrusion data correction device for implementing the above-mentioned method for correcting rock high-pressure mercury intrusion data. The solution provided by this device is similar to the solution described in the above-described method. Therefore, the specific limitations of one or more rock high-pressure mercury intrusion data correction device embodiments provided below can be found in the limitations of the rock high-pressure mercury intrusion data correction method described above, and will not be repeated here.

[0068] In one exemplary embodiment, such as Figure 6 As shown, a correction device for high-pressure mercury intrusion data of rocks is provided, comprising: The data acquisition module 301 is used to acquire the high-pressure mercury intrusion data to be corrected for the target rock and the high-pressure mercury intrusion data of the preset sample under the same high-pressure mercury intrusion test conditions; the high-pressure mercury intrusion data includes the cumulative mercury intrusion amount under different pressures; the preset sample is a solid with a smooth surface, no interconnected pores, a known bulk modulus, and only undergoing elastic deformation.

[0069] Blank error calculation module 302 is used to calculate the elastic deformation of the preset sample under different pressures based on the bulk modulus of the preset sample, and subtract the elastic deformation under the corresponding pressure from the high pressure mercury intrusion data of the preset sample to obtain blank error data under different pressures.

[0070] The preliminary correction module 303 is used to obtain preliminary corrected high-pressure mercury intrusion data of the target rock based on the high-pressure mercury intrusion data to be corrected and blank error data under different pressures.

[0071] The sample deformation acquisition module 304 is used to calculate the volume change rate between adjacent pressures based on the preliminary corrected high-pressure mercury intrusion data of the target rock, and to determine the deformation data of the target rock under different pressures based on the minimum volume change rate.

[0072] The final correction module 305 is used to subtract the deformation data of the target rock under the corresponding pressure from the preliminary corrected high-pressure mercury intrusion data of the target rock to obtain the corrected high-pressure mercury intrusion data of the target rock.

[0073] In one exemplary embodiment, a computer device is provided, which may be a server or a terminal, and its internal structure diagram may be as follows. Figure 7 As shown, the computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides the environment for the operation of the operating system and computer programs in the non-volatile storage media. The database stores correction processing data for high-pressure mercury intrusion porosimetry (HSI) data in rocks. The I / O interfaces are used for information exchange between the processor and external devices. The communication interface is used for communication with external terminals via a network connection. When the computer program is executed by the processor, it implements a method for correcting high-pressure mercury intrusion porosimetry (HSI) data in rocks.

[0074] Those skilled in the art will understand that Figure 7The structures shown are merely block diagrams of some structures related to the present application and do not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than shown in the figures, or combine certain components, or have different component arrangements. In an exemplary embodiment, a computer device is provided, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0075] In one exemplary embodiment, a computer-readable storage medium is provided storing a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0076] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0077] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0078] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0079] The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0080] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0081] This document uses specific examples to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. Furthermore, those skilled in the art will recognize that, based on the ideas of this application, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for correcting high-pressure mercury intrusion data in rocks, characterized in that, The correction method for the high-pressure mercury intrusion data of the rock includes: Under the same high-pressure mercury intrusion test conditions, high-pressure mercury intrusion data to be corrected for the target rock and high-pressure mercury intrusion data for a preset sample are obtained; the high-pressure mercury intrusion data includes the cumulative mercury intrusion amount under different pressures; the preset sample is a solid with a smooth surface, no interconnected pores, a known bulk modulus, and only undergoing elastic deformation; Based on the bulk modulus of the preset sample, the elastic deformation of the preset sample under different pressures is calculated, and the elastic deformation under the corresponding pressure is subtracted from the high pressure mercury intrusion data of the preset sample to obtain blank error data under different pressures. Based on the high-pressure mercury intrusion data to be corrected for the target rock and the blank error data under different pressures, preliminary corrected high-pressure mercury intrusion data for the target rock are obtained. Based on the preliminary corrected high-pressure mercury intrusion data of the target rock, the volume change rate between adjacent pressures is calculated, and the deformation data of the target rock under different pressures is determined based on the minimum volume change rate. Subtract the deformation data of the target rock under the corresponding pressure from the preliminary corrected high-pressure mercury intrusion data of the target rock to obtain the corrected high-pressure mercury intrusion data of the target rock.

2. The method for correcting high-pressure mercury intrusion data in rocks according to claim 1, characterized in that, Based on the bulk modulus of the preset sample, the elastic deformation of the preset sample under different pressures is calculated, specifically including: ； in, This represents the elastic deformation of the preset sample, in cm. 3 ; This indicates the initial volume of the preset sample, in cm³. 3 K represents the pre-defined bulk modulus of the sample, in GPa. This indicates the change in pressure, in MPa.

3. The method for correcting high-pressure mercury intrusion data in rocks according to claim 1, characterized in that, Based on the high-pressure mercury intrusion data to be corrected for the target rock and the blank error data under different pressures, preliminary corrected high-pressure mercury intrusion data for the target rock are obtained, specifically including: If the volume of the target rock is equal to the volume of the preset sample, the blank error data under the corresponding pressure is subtracted from the high pressure mercury intrusion data of the target rock to be corrected, so as to obtain the preliminary corrected high pressure mercury intrusion data of the target rock. If the volume of the target rock is not equal to the volume of the preset sample, a preliminary correction is made using the following formula: ； in, This indicates the preliminary corrected high-pressure mercury intrusion data for the target rock; This represents the high-pressure mercury intrusion data of the target rock that needs correction. Indicates blank error data; Indicates the volume of the preset sample; The volume of the target rock is represented by P; the pressure is represented by MPa. This indicates the bulk modulus of mercury, measured in gigabytes of water (GPa).

4. The method for correcting high-pressure mercury intrusion data in rocks according to claim 3, characterized in that, The methods for calculating the bulk modulus of mercury include: Under the same high-pressure mercury intrusion test conditions as the target rock, no solid sample was placed into the air dilatometer to conduct a blank high-pressure mercury intrusion test and obtain blank high-pressure mercury intrusion data. Subtract the blank error data at the corresponding pressure from the blank high-pressure mercury intrusion data to obtain the mercury deformation data; Based on the deformation data of mercury, the bulk modulus of mercury is calculated using the following formula: ； in, This represents the rate of change of mercury pressure-volume, calculated based on mercury deformation data.

5. The method for correcting high-pressure mercury intrusion data in rocks according to claim 1, characterized in that, Before the high-pressure mercury intrusion test, the target rock and the preset sample were dried according to the preset drying conditions, which included drying at 70°C for 24 hours.

6. The method for correcting high-pressure mercury intrusion data in rocks according to claim 1, characterized in that, The preset samples include quartz, high-purity aluminum, or high-purity copper.

7. A device for correcting high-pressure mercury intrusion data in rocks, characterized in that, The rock high-pressure mercury intrusion data correction device applies the rock high-pressure mercury intrusion data correction method according to any one of claims 1-6, and the rock high-pressure mercury intrusion data correction device comprises: The data acquisition module is used to acquire, under the same high-pressure mercury intrusion test conditions, the high-pressure mercury intrusion data to be corrected for the target rock and the high-pressure mercury intrusion data for a preset sample; the high-pressure mercury intrusion data includes the cumulative mercury intrusion amount under different pressures; the preset sample is a solid with a smooth surface, no interconnected pores, a known bulk modulus, and only undergoing elastic deformation; The blank error calculation module is used to calculate the elastic deformation of the preset sample under different pressures based on the bulk modulus of the preset sample, and to subtract the elastic deformation under the corresponding pressure from the high pressure mercury intrusion data of the preset sample to obtain blank error data under different pressures. The preliminary correction module is used to obtain preliminary corrected high-pressure mercury intrusion data of the target rock based on the high-pressure mercury intrusion data to be corrected and blank error data under different pressures. The sample deformation acquisition module is used to calculate the volume change rate between adjacent pressures based on the preliminary corrected high-pressure mercury intrusion data of the target rock, and to determine the deformation data of the target rock under different pressures based on the minimum volume change rate. The final correction module is used to subtract the deformation data of the target rock under the corresponding pressure from the preliminary corrected high-pressure mercury intrusion data of the target rock to obtain the corrected high-pressure mercury intrusion data of the target rock.

8. A computer device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the method for correcting rock high-pressure mercury intrusion data according to any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the method for correcting high-pressure mercury intrusion data in rocks as described in any one of claims 1-6.

10. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the method for correcting high-pressure mercury intrusion data in rocks as described in any one of claims 1-6.

Images