Method for calculating surface relaxation rate of shale

Abstract

A method for calculating a surface relaxation rate of a shale includes: a relaxation time T distribution curve and a pore throat radius r distribution curve are obtained through experiments; abscissas of the two distribution curves are standardized, and the abscissa of the relaxation time T distribution curve is expanded or shrunk to ensure an abscissa value corresponding to a maximum ordinate value in the transformed relaxation time T distribution curve is same as an abscissa value corresponding to a maximum ordinate value in the pore throat radius r distribution curve; straight lines with a number of N parallel to a y-axis of a combined curve graph including the two distribution curves are drawn and a value corresponding to each straight line is calculated; and value with the number of N are processed to obtain a final surface relaxation rate .

Claims

1. A method for calculating a surface relaxation rate of a shale, comprising: S1, selecting a rock sample from a target shale reservoir corresponding to the shale to perform a nuclear magnetic resonance (NMR) experiment under saturated formation water conditions to obtain a relaxation time T distribution curve; S2, selecting a target rock sample from the rock sample to perform a constant-rate high-pressure mercury injection experiment to obtain a pore throat radius r distribution curve; S3, standardizing abscissas of the relaxation time T distribution curve and the pore throat radius r distribution curve to obtain a standard relaxation time T distribution curve and a standard pore throat radius r distribution curve with a same abscissa starting value and a same abscissa ending value on a logarithmic scale; S4, determining an abscissa value corresponding to a maximum ordinate value on the standard relaxation time T distribution curve, denoted as a, determining an abscissa value corresponding to a maximum ordinate value on the standard pore throat radius r distribution curve, denoted as b, calculating a value of b/, expanding, in response to the value of b/>1, the abscissa of the standard relaxation time T distribution curve by b/units to obtain a transformed relaxation time T distribution curve, or shrinking, in response to the value of b/<1, the abscissa of the standard relaxation time T distribution curve by a/b units to obtain a transformed relaxation time T distribution curve, to make an abscissa value corresponding to a maximum ordinate value in the transformed relaxation time T distribution curve be same as the abscissa value corresponding to the maximum ordinate value in the standard pore throat radius r distribution curve, and combining the transformed relaxation time T distribution curve and the standard pore throat radius r distribution curve in a same coordinate system to obtain a combined curve graph with both the transformed relaxation time T distribution curve and the standard pore throat radius r distribution curve; S5, drawing N straight lines, wherein N is a positive integer, parallel to a y-axis of the combined curve graph with each straight line passing through both the transformed relaxation time T distribution curve and the standard pore throat radius r distribution curve, determining abscissa values of the straight lines as high-pressure mercury injection pore throat radius r values and transformed NMR relaxation time calculated values, multiplying the transformed NMR relaxation time calculated values by b/ or /b or dividing the transformed NMR relaxation time calculated values by b/ or /b to obtain pre-transformation NMR relaxation time calculated values, and obtaining real relaxation time T values which are closest to the pre-transformation NMR relaxation time calculated values from original NMR experimental data obtained in the step S1; S6, calculating a surface relaxation rate based on the high-pressure mercury injection pore throat radius r value and the real relaxation time T value for each straight line to thereby obtain N of surface relaxation rates respectively corresponding to the N of straight lines, wherein a calculation formula of the surface relaxation rate is as follows: $=\frac{R}{F_{S}T},$ where r represents the high-pressure mercury injection pore throat radius value, with a unit in nanometers (nm), T represents the real relaxation time value, and F.sub.s represents a shape factor, determined by a pore throat structure shape of the target shale reservoir; and S7, performing data processing to the surface relaxation rates to obtain a final surface relaxation rate of the target rock sample, wherein the step S7 specifically comprises: S71, sorting N of the surface relaxation rates in an ascending order to obtain sorted surface relaxation rates, wherein a minimum value of the sorted surface relaxation rates is .sub.min and a maximum value of the sorted surface relaxation rates is .sub.max; S72, calculating a difference between every two adjacent values of the sorted surface relaxation rates to obtain difference values d, and calculating an average value D of the difference values d; S73, dividing an interval [.sub.min, .sub.max] formed by the minimum value .sub.min and the maximum value .sub.max into a plurality of secondary intervals with the average value D as an interval, counting values and a number of the values in each secondary interval, and calculating an average value of the values in each secondary interval; S74, calculating a weight coefficient of the average value in each secondary interval based on the number of the values in each secondary interval; and S75, multiplying the average value in each secondary interval with the corresponding weight coefficient to thereby obtain products corresponding to the plurality of secondary intervals, and summing the products to obtain the final surface relaxation rate of the target rock sample.

2. The method for calculating the surface relaxation rate of the shale as claimed in claim 1, wherein in the step S74, the weight coefficient of the average value in each secondary interval is equal to a quotient of the number of the values in each secondary interval divided by a total number of values in the interval [.sub.min, .sub.max].

3. The method for calculating the surface relaxation rate of the shale as claimed in claim 1, wherein in the step S5, the number of the straight lines parallel to the y-axis is N, and N15.

4. The method for calculating the surface relaxation rate of the shale as claimed in claim 1, wherein in the step S6, in response to the target shale reservoir having a spherical pore throat structure, the shape factor F.sub.s=3, in response to the target shale reservoir having a cylindrical pore throat structure, the shape factor F.sub.s=2.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 illustrates a relaxation time distribution curve diagram of a NMR experiment.

(2) FIG. 2 illustrates a pore throat radius distribution curve diagram of a high-pressure mercury injection experiment.

(3) FIG. 3 illustrates a comparison diagram of transverse pore throat structure distribution curves between the NMR experiment and the high-pressure mercury injection experiment.

(4) FIG. 4 illustrates a schematic diagram of enlarged transverse pore throat structure distribution curves between the NMR experiment and the high-pressure mercury injection experiment.

(5) FIG. 5 illustrates a schematic diagram of multiple intersection points of the transverse pore throat structure distribution curves between the NMR experiment and the high-pressure mercury injection experiment with straight lines parallel to a y-axis.

DETAILED DESCRIPTION OF EMBODIMENTS

(6) Preferred embodiments of the disclosure are described below in conjunction with the accompanying drawings. It should be understood that the preferred embodiments described herein are only for illustrating and explaining the disclosure, and are not intended to limit the disclosure.

(7) The method for calculating a surface relaxation rate of a shale includes following steps.

(8) Step 1, full-size (also referred to as full diameter) core samples are selected from a target shale reservoir of the shale to respectively perform a NMR experiment under saturated formation water conditions and to obtain a relaxation time T distribution curve under an influence-eliminated saturated state, i.e., transverse relaxation time distribution of a saturated water sample pore fluid; and a small rock sample with a diameter of 25 millimeters (mm) and a height of 50 mm is selected from the full-size core samples to perform a constant-rate high-pressure mercury injection experiment to obtain a pore throat radius r distribution curve, i.e., transverse relaxation time of the saturated water sample pore fluid. The relaxation time T distribution curve obtained in the embodiment is shown in FIG. 1, and the pore throat radius r distribution curve obtained in the embodiment is shown in FIG. 2. It should be noted that x-axes in FIG. 1 and FIG. 2 should have a same abscissa starting value and a same abscissa ending value which are on a logarithmic scale.

(9) According to known formulas, a relationship between a pore throat radius (r) and relaxation time (T) is as follows:
r=TF.sub.s(1) where r represents the pore throat radius, with a unit in nanometers (nm), T represents the relaxation time, represents a surface relaxation rate of a rock sample, and F.sub.s represents a shape factor, in response to the target shale reservoir having a spherical pore throat structure, the shape factor F.sub.s=3, in response to the target shale reservoir having a cylindrical pore throat structure, the shape factor F.sub.s=2. In the embodiment of the disclosure, F.sub.s=2.

(10) A formula is obtained from the formula (1):
r=T2(2)

(11) An expression formula of the can be obtained by deforming the formula (2):

(12) $\begin{array}{cc}=\frac{r}{2T}& \left(3\right)\end{array}$

(13) It can be seen from the formula (3) that the surface relaxation rate depends on the relaxation time T and the pore throat radius r. In the transverse relaxation time distribution diagram of the saturated water sample pore fluid, under the NMR experiment, the relaxation time serves as an abscissa, and a signal intensity serves as an ordinate. Under the high-pressure mercury injection experiment, the pore throat radius serves as the abscissa, and a pore volume serves as the ordinate. When the ordinates correspond to each other, there is a relationship between the abscissas as stated in the formula (3), and the value obtained by the disclosure is the relationship between the relaxation time T under the NMR experiment and the pore throat radius r under the high-pressure mercury injection experiment.

(14) Step 2, a longitudinal comparison and averaging method is used in the disclosure. It is found in experiments that the distribution curves obtained from the NMR experiment and the high-pressure mercury injection experiment are similar in shape. When the abscissas are expanded, there are situations that peaks correspond to peaks and troughs correspond to troughs on the ordinates.

(15) Firstly, abscissas of the NMR experiment relaxation time distribution curve (i.e., the relaxation time T distribution curve) and the high-pressure mercury injection experiment pore throat radius distribution curve (i.e., the pore throat radius r distribution curve) are standardized to obtain a standard relaxation time distribution curve and a standard pore throat radius distribution curve with the same abscissa starting value and the same abscissa ending value on a logarithmic scale. FIG. 1 and FIG. 2 can be combined into a same coordinate system to obtain a comparison diagram as shown in FIG. 3. It can be seen that an abscissa value corresponding to a maximum ordinate value of the standard relaxation time distribution curve is not the same as an abscissa value corresponding to a maximum ordinate value of the standard pore throat radius distribution curve, and the difference is significant.

(16) In order to make the abscissa value corresponding to the maximum ordinate value in the NMR experiment relaxation time distribution curve same as the abscissa value corresponding to the maximum ordinate value read from the high-pressure mercury injection experiment pore throat radius distribution curve, the following processing is carried out by the disclosure: the abscissa value corresponding to the maximum ordinate value on the standard relaxation time distribution curve is read and denoted as a, the abscissa value corresponding to the maximum ordinate value on the standard pore throat radius distribution curve is read and denoted as b, and a value of b/ is calculated. In response to the value of b/>1, the abscissa of the standard relaxation time distribution curve is expanded by b/units to obtain a transformed relaxation time distribution curve; or in response to the value of b/<1, the abscissa of the standard relaxation time distribution curve is shrunk by /b units to obtain a transformed relaxation time distribution curve, to make an abscissa value corresponding to a maximum ordinate value in the transformed relaxation time distribution curve be same as the abscissa value corresponding to the maximum ordinate value in the standard pore throat radius distribution curve, as shown in FIG. 4.

(17) Then, straight lines with a number of N parallel to a y-axis passing through original experimental data points of the high-pressure mercury injection experiment are arbitrarily drawn in a distribution curve area in FIG. 4, as shown in FIG. 5. Each straight line has a unique intersection point with each of the transformed relaxation time distribution curve and the standard pore throat radius distribution curve. An abscissa value of the unique intersection point is read as a standard pore throat radius value and a transformed NMR relaxation time calculated value. The abscissa values of the unique intersection points (i.e., the transformed NMR relaxation time calculated values) are multiplied by b/ or /b, or divided by b/ or /b to obtain pre-transformation relaxation time calculated values. real relaxation time values which are closest to the pre-transformation relaxation time calculated values are obtained from original NMR experimental data. Each straight line corresponds to a group of data. The pore throat radius (nm) and the corresponding relaxation time (millisecond, ms) are read from the intersection points in FIG. 5. The relaxation time (ms) before expansion corresponding to the pore throat radius (nm) is calculated, and the real relaxation time values most similar to the NMR relaxation calculated time values before expansion are found from the original data, as shown in Table 1.

(18) TABLE-US-00001 TABLE 1 pore throat radius and its corresponding relaxation time real relaxation time pre- value closest to transformed transformation pre-transformation pore throat NMR relaxation relaxation relaxation time radius time calculated time calculated calculated (nm) value (ms) value (ms) value T (ms) 2.409488668 2.409488668 0.053028 0.052 3.75160677 3.75160677 0.082565 0.090 4.546619668 4.546619668 0.100061 0.117 5.635665338 5.635665338 0.124029 0.131 9.95563586 9.95563586 0.219102 0.227 12.15250112 12.15250112 0.26745 0.258 15.49449685 15.49449685 0.341 0.341 21.12414516 21.12414516 0.464896 0.475 30.89845885 30.89845885 0.680008 0.709 40.26639709 40.26639709 0.886175 0.921 62.52133179 62.52133179 1.375958 1.351 95.44831543 95.44831543 2.100609 2.196 120.8368652 120.8368652 2.659355 2.583 150.97677 150.97677 3.322669 3.193 183.019812 183.019812 4.027866 4.199

(19) Step 3, the pore throat radius and the real relaxation time value of each group are substituted to the formula (3) =r/2T to obtain a value of the intersection point of each group, as shown in Table 2.

(20) TABLE-US-00002 TABLE 2 the value corresponding to the intersection point of each group intersection point number 1 23.16816027 2 20.84225983 3 19.42999858 4 21.51017304 5 21.92871335 6 23.55135877 7 22.71920359 8 22.23594227 9 21.79016844 10 21.86015043 11 23.13890888 12 21.73231226 13 23.39079854 14 23.64183683 15 21.79326173

(21) Step 4, the values are processed as follows, to obtain a final surface relaxation rate of the small rock sample. Specifically, the step 4 includes the following five sub-steps. 1, the values are sorted in ascending order to obtain sorted surface relaxation rates: 19.42999858, 20.84225983, 21.51017304, 21.73231226, 21.79016844, 21.79326173, 21.86015043, 21.92871335, 22.23594227, 22.71920359, 23.13890888, 23.16816027, 23.39079854, 23.55135877, 23.64183683.

(22) It can be seen that the values are distributed in an interval [19.42999858, 23.64183683]. 2, a difference d between every two adjacent values of the sorted surface relaxation rates is calculated to obtain difference values d, and an average value D of the difference values d are calculated.

(23) d.sub.1=1.412261249, d.sub.2=0.667913214, d.sub.3=0.222139211, d.sub.4=0.057856186, d.sub.5=0.003093288, d.sub.6=0.066888701, d.sub.7=0.068562917, d.sub.8=0.307228922, d.sub.9=0.483261323, d.sub.10=0.419705284, d.sub.11=0.029251392, d.sub.12=0.222638266, d.sub.13=0.160560232, d.sub.14=0.090478066.

(24) $D=\frac{\begin{array}{c}{d}_1+d_2+d_3+d_4+d_5+d_6+\\ d_7+d_8+d_9+d_{10}+d_{11}+d_{12}+d_{13}+d_{14}\end{array}}{14}=0.300845589.$ 3, the interval [19.42999858, 23.64183683] are divided into multiple secondary intervals based on the D=0.300845589 as an interval, the values and a number of the p values in each secondary interval are counted, and an average value of the values in each secondary interval is calculated, shown in Table 3.

(25) TABLE-US-00003 TABLE 3 the value corresponding to each secondary interval Number Average Corresponding of value of Secondary interval value values the values [19.42999858, 19.73084417) 19.42999858 1 19.42999858 (19.73084417, 20.03168976) / 0 / (20.03168976, 20.33253535) / 0 / (20.33253535, 20.63338094) / 0 / (20.63338094, 20.93422653) 20.84225983 1 20.84225983 (20.93422653, 21.23507212) / 0 / (21.23507212, 21.53591771) 21.51017304 1 21.51017304 (21.53591771, 21.83676329) 21.73231226 3 21.77191414 21.79016844 21.79326173 (21.83676329, 22.13760888) 21.86015043 2 21.89443189 21.92871335 (22.13760888, 22.43845447) 22.23594227 1 22.23594227 (22.43845447, 22.73930006) 22.71920359 1 22.71920359 (22.73930006, 23.04014565) / 0 / (23.04014565, 23.34099124) 23.13890888 2 23.15353458 23.16816027 (23.34099124, 23.64183683] 23.39079854 3 23.52799805 23.55135877 23.64183683 4, a weight coefficient of the average value in each secondary interval is calculated based on the number of the values in each secondary interval: when the number of the p values is 0, p=0/15=0; when the number of the values is 1, p=1/15; when the number of the values is 2, p=2/15; and when the number of the values is 3, p=3/15. 5, the average value in each secondary intervals are multiplied by the corresponding weight coefficient to obtain products corresponding to the secondary intervals and the products are summed up to obtain the final surface relaxation rate of the small rock sample:

(26) ${}^{}=19.42999858\frac{1}{15}+20.84225983\frac{1}{15}+21.51017304\frac{1}{15}+21.77191414\frac{3}{15}+21.89443189\frac{2}{15}+22.23594227\frac{1}{15}+22.71920359\frac{1}{15}+23.15353458\frac{2}{15}+23.52799805\frac{3}{15}=22.18221645,$ i.e., the final surface relaxation rate of the small rock sample is 22.18221645.

(27) The above description is only the preferred embodiment of the disclosure and does not limit the disclosure in any form. Although the disclosure has been disclosed in the preferred embodiment, it is not intended to limit the disclosure. Those skilled in the art can use the disclosed technical content to make slight changes or modifications to equivalent embodiments without departing from the scope of the technical solution of the disclosure. Any simple modifications, equivalent changes, and modifications made to the above embodiments based on the technical essence of the disclosure without departing from the technical solution of the disclosure still belong to the scope of the technical solution of the disclosure.