Disclosure of Invention
The invention provides a method for quantitatively evaluating the gas diffusion capacity in rock, which is used for solving the problem of inaccurate gas diffusion coefficient measured by using an isobaric diffusion experiment in the prior art, and can quantitatively characterize the gas diffusion capacity in rock more objectively and scientifically.
The invention provides a method for quantitatively evaluating the diffusion capability of gas in rock, which is characterized by comprising the following steps:
Under the conditions of target temperature and pressure, introducing a preset molar quantity of gas a into a first diffusion chamber positioned at one side of a rock sample, and introducing a preset molar quantity of gas b into a second diffusion chamber positioned at the other side of the rock sample, wherein at least one of the gas a and the gas b is a gas to be detected;
Detecting the molar amount of the gas b diffused into the first diffusion chamber and the molar amount of the gas a diffused into the second diffusion chamber at set time intervals;
And determining the diffusion index I of the gas to be detected based on a pre-constructed diffusion index model and the molar concentration change condition of the gas a and the gas b between the first diffusion chamber and the second diffusion chamber.
According to the method for quantitatively evaluating the gas diffusion capacity in the rock, the diffusion index model is constructed based on the mass flux and the molar concentration gradient of substances of the gas to be tested and combined with the Fick first law.
According to the method for quantitatively evaluating the gas diffusion capacity in the rock, which is provided by the invention, the diffusion index model is as follows:
;
Wherein, The molar concentration difference of the gas to be measured at the initial moment is the molar concentration difference of the gas to be measured in the first diffusion chamber and the second diffusion chamber; the first diffusion chamber and the second diffusion chamber are arranged in The molar concentration difference of the gas to be measured at the moment; a volume for the first diffusion chamber; a volume of the second diffusion chamber; S is the cross-sectional area of the rock sample to be measured, and H is the length of the rock sample to be measured.
According to the method for quantitatively evaluating the gas diffusion capacity in the rock, in the diffusion index model, the volume of the first diffusion chamberAnd the volume of the second diffusion chamberIs a fixed value.
According to the method for quantitatively evaluating the gas diffusion capacity in the rock, provided by the invention, the initial time is set to be 0, the molar quantity of the gas to be measured at the initial time in the diffusion chamber at the other side is set to be 0, and then the diffusion index model is as follows:
;
Wherein, An initial molar quantity of the gas to be measured; For the gas to be measured Molar mass of moment diffusion.
According to the method for quantitatively evaluating the gas diffusion capacity in the rock, in the diffusion index model, the volume of a first diffusion chamberAnd the volume of the second diffusion chamberEqual.
According to the method for quantitatively evaluating the gas diffusion capacity in the rock, which is provided by the invention, the diffusion index model is as follows:
;
Wherein V is the volume of the first diffusion chamber and the second diffusion chamber.
According to the method for quantitatively evaluating the diffusion capacity of the gas in the rock, the method for determining the diffusion index of the gas to be measured based on a pre-constructed diffusion index model and the molar concentration change condition of the gas a and the gas b in the first diffusion chamber and the second diffusion chamber comprises the following steps:
determining the linear relation between C and t based on the diffusion molar quantity and the diffusion index model of the gas to be detected at different moments;
And obtaining a final diffusion index I through fitting based on the linear relation between C and t.
According to the method for quantitatively evaluating the diffusion capability of the gas in the rock, the molar quantity of the gas to be measured is obtained by the volume fraction of the gas to be measured by a gas chromatograph.
According to the method for quantitatively evaluating the gas diffusion capability in the rock, under the condition of physical gas, the volume fraction and the molar quantity meet the following relation:
;
Wherein, Is thatThe volume fraction of the gas a in the first diffusion chamber at the moment; Is that The volume fraction of the gas a in the second diffusion chamber at the moment; The molar quantity of the gas a in the first diffusion chamber is the initial moment; Is that A molar amount of gas a diffused into the second diffusion chamber at the moment; The molar quantity of the gas b in the second diffusion chamber is the initial moment; Is that The molar quantity of gas b diffused into the first diffusion chamber at the moment.
Compared with the volume concentration, the molar concentration is taken as a dimensional physical quantity, the method for quantitatively evaluating the gas diffusion capacity in the rock can reflect the gas diffusion flux in the rock more intuitively and accurately, so that the gas diffusion capacity in the rock can be quantitatively represented more objectively and scientifically, reliable geological theory data are provided for the closure of a cover layer in a resource gas reservoir and the evaluation of a gas storage medium, and the method has important significance for researching the mechanism of the resource gas reservoir.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In order to facilitate understanding of the method for quantitatively evaluating the gas diffusion capability in the rock, the application background is introduced first, and for gas resources such as natural gas, natural hydrogen, helium and the like which are hidden in the stratum, the quantitative characterization of the gas diffusion capability in the stratum is an important parameter for evaluating the plugging performance of a gas reservoir cover layer, the gas preservation condition and the potential of the gas resources.
In the related technology, the diffusion capacity of gas in rock is evaluated at home and abroad, the diffusion coefficient of gas is obtained by utilizing the volume content change of gas diffusion in an isobaric diffusion experiment (specific reference can be made to the existing industry standard SY/T6129-2016, the method for measuring the diffusion coefficient of hydrocarbon gas in rock).
In the experiment, a diffusion cavity is arranged in the isobaric gas device, and the rock to be tested divides the diffusion cavity into two gas chambersAndAt a specific temperature and pressure, toIntroducing a gas a such as methane, hydrogen, etc. into the gas chamberIntroducing gas b, such as nitrogen or argon, with stable chemical properties, collecting and diffusing into the container at set time intervals by gas chromatographThe volume concentration of the medium gas b and diffusion toThe volume concentration of the medium gas a is based on the volume concentration change and the diffusion coefficient of the gas to be measured is obtained by combining the Fick first law.
However, the inventors found in practical experiments that by the above-described test method, the gas diffusion coefficient obtained by the method was equal regardless of which two gases were used for isobaric diffusion, and the gas chamber was set at the initial timeThe volume content of the medium gas a is 100%, and the air chamberThe volume content of the gas b is 100%, and at a certain moment, the gas chamber is measured by using a chromatographThe content of the gas b in the gas chamber is x percentThe content of the gas a in the gas chamber is y percentThe content of the gas a in the gas chamber is (100-x)%, and the gas chamberThe content of the gas b is (100-y)%, the concentration difference of the gas a in the two air chambers is calculated to be (100-x-y)%, the concentration difference of the gas b is calculated to be (100-x-y)%, and the concentration differences of the gas a and the gas b between the two air chambers are equal.
Specifically, taking methane and nitrogen as examples, the detection results of the method for measuring the diffusion coefficient of hydrocarbon gas in rock by using SY/T6129-2016 are shown in the following table:
Table 1, one of the diffusion coefficient detection reports;
table 2, second diffusion coefficient detection report;
It is clear from the above two sets of detection results that the concentration differences of the two gases in the two diffusion chambers are equal at any time, so that the diffusion coefficients of the two gases calculated according to SY/T6129-2016, the method for measuring the diffusion coefficient of hydrocarbon gas in rock are equal, and the two diffusion chambers are two completely different gases, the diffusion coefficients of any two completely different gases in a rock sample are obviously unreasonable, in other words, the experimental method is accurate only when the diffusion coefficients of the two gases in the rock sample are equal or very close, but obviously, the diffusion coefficients of the two gases in the rock sample are unknown.
In addition, in the actual test, since the gas to be tested (such as methane, hydrogen, helium, etc.) is the resource gas actually buried in the ground, and the evaluation of the rock covering these gases is the main objective of the experiment, researchers generally only pay attention to the concentration difference and diffusion coefficient of the gas to be tested, so as to measure the diffusion capacity of the resource gas in the cover rock layer, while another gas is used as an aid, the relationship between the concentration difference, diffusion coefficient and the diffusion coefficient thereof and the diffusion coefficient of the gas to be tested is often ignored, so that the calculation result cannot truly reflect the diffusion capacity of the gas, and further influence the evaluation result.
Therefore, how to objectively and scientifically quantitatively characterize the gas diffusion capacity is an important topic to be solved currently.
Based on the technical problems and findings, the invention provides a method for quantitatively evaluating the gas diffusion capability in the rock, which can quantitatively characterize the gas diffusion capability in the rock more objectively and scientifically, provides reliable geological theory data for the closure of a cover layer in a resource gas reservoir and the evaluation of a gas storage medium, and has important significance for researching the mechanism of resource gas reservoir formation.
The method of quantitatively evaluating the ability of a gas to diffuse in rock according to the present invention is described below with reference to fig. 1-2.
Referring to fig. 1, a method for quantitatively evaluating the ability of a gas to diffuse in rock, comprising the steps of:
And 101, introducing a preset molar quantity of gas a into a first diffusion chamber positioned at one side of the rock sample under the condition of target temperature and pressure, and introducing a preset molar quantity of gas b into a second diffusion chamber positioned at the other side, wherein at least one of the gas a and the gas b is the gas to be detected.
Step 102, detecting the molar amount of the gas b diffused into the first diffusion chamber and the molar amount of the gas a diffused into the second diffusion chamber at predetermined time intervals.
And step 103, determining the diffusion index of the gas to be measured based on a pre-constructed diffusion index model and the molar concentration change condition of the gas a and the gas b between the first diffusion chamber and the second diffusion chamber.
Specifically, when the test is actually performed, a rock sample is placed in a diffusion chamber of an isobaric gas device, the rock sample divides the diffusion chamber of the isobaric gas device into the first diffusion chamber and the second diffusion chamber, a gas a, for example, a resource gas such as hydrogen, methane, helium and the like is introduced into the first diffusion chamber, a gas b, for example, a chemically stable gas such as nitrogen, argon and the like is introduced into the diffusion chamber at the other side, the molar amount of the gas b diffused into the first diffusion chamber and the molar amount of the gas a diffused into the second diffusion chamber are detected at preset intervals, and a diffusion index is obtained according to a molar concentration gradient of the gas a and the gas b between the first diffusion chamber and the second diffusion chamber and a diffusion index model constructed in advance.
At the initial time, the molar concentration of the gas a in the first diffusion chamber isWherein, the method comprises the steps of, wherein,In terms of the molar amount of the gas a,The molar concentration of the gas b in the second diffusion chamber is equal to the volume of the first diffusion chamberWherein, the method comprises the steps of, wherein,In terms of the molar amount of the gas b,Is the volume of the second diffusion chamber.
At a certain time, when the molar amount of the gas a diffused into the second diffusion from the first diffusion chamber is x and the molar amount of the gas b diffused into the first diffusion from the second diffusion chamber is y, the molar concentration of the gas a in the first diffusion chamber isThe molar concentration of the gas a in the second diffusion chamber isThe molar concentration of the gas b in the second diffusion chamber isThe molar concentration of the gas b in the first diffusion chamber isThe molar concentration difference of the gas a between the first diffusion chamber and the second diffusion chamber is-And the molar concentration difference of the gas b between the first diffusion chamber and the second diffusion chamber isThe diffusion index is obtained according to the molar concentration variation and the diffusion index model.
Compared with the volume concentration, the molar concentration is taken as a dimensional physical quantity, the diffusion flux of the gas in the rock can be reflected more intuitively and accurately, so that the diffusion capacity of the gas in the rock can be quantitatively represented more objectively and scientifically, reliable geological theory data are provided for the closure of a cover layer in a resource gas reservoir and the evaluation of a gas storage medium, and the method has important significance for researching the mechanism of the reservoir formation of the resource gas.
In one embodiment of the invention, the diffusion index model is constructed based on the mass flux and molar concentration gradient of the substance of the gas to be measured in combination with the Fick first law.
Specifically, the construction of the diffusion index model includes the following steps.
The rock sample is cut into a cylinder with the bottom surface section area S and the height H, the gas diffusion quantity passing through the rock sample in a certain time is measured in two diffusion chambers, and then the gas diffusion index is obtained according to the Fick first law by the measured values. According to Fick's first law, the flux from the high concentration region to the low concentration stream is proportional to the concentration gradient, and for one-dimensional diffusion, the flux of a substance in the x (i.e. the height direction of the rock to be measured) direction can be considered to be proportional to the concentration gradient there, and the specific formula is as follows (it is to be noted that the units in the formula are all determined according to the experimental practical conditions):
Wherein N is the quantity of a substance diffusing gas, mol, S is the bottom surface sectional area of a rock sample, cm 3, t is the duration of a diffusion test, S, I is the gas diffusion index, cm 2/S; h is the height of the rock sample, i.e. the diffusion distance of the diffusion gas, cm.
Setting the volume of the first diffusion chamber asThe volume of the second diffusion chamber isAfter t time diffusion, the molar concentrations of the gases to be detected in the first diffusion chamber and the second diffusion chamber are respectivelyAndNamely there is。
As is known from the principle of conservation of mass,Substitution simplification can be obtained:
Setting initial time, wherein the molar concentration difference of the gas to be measured between the first diffusion chamber and the second diffusion chamber is as follows Through the process ofAfter time diffusion, the molar concentration difference of the gas to be measured between the first diffusion chamber and the second diffusion chamber isThe integral of the above formula can be obtained:
Then 。
The relation between the molar concentration variation of the gas to be measured between two diffusion chambers and the test time can be measured through the diffusion index model.
In one embodiment of the present invention, in the diffusion test, the pressure between the first diffusion chamber and the second diffusion chamber is equal in an ideal state, however, due to the difference of the diffusion capacities of the gas a and the gas b, a tiny pressure difference is necessarily generated in the first diffusion chamber and the second diffusion chamber along with the progress of the experiment, so that the volumes of the first diffusion chamber and the second diffusion chamber can be adjusted along with the progress of the experiment in an ideal state, and the diffusion chambers on two sides are ensured to be in an isobaric state.
However, the above method can certainly increase the difficulty of the experiment, considering that the pressure difference between the two diffusion chambers is small and the effect on the result of the diffusion experiment is limited, thus, in the diffusion index model, the volume of the first diffusion chamberAnd the volume of the second diffusion chamberIs a fixed value.
In one embodiment of the present invention, the initial time is set to 0, and the molar quantity of the gas to be measured in the diffusion chamber at the other side at the initial time is set to 0, and then the diffusion index model is as follows:
。
Wherein, For the initial molar amount of the gas to be measured,For the gas to be measuredMolar mass of moment diffusion.
Through the diffusion index model, the relation between the change of the molar quantity of the gas to be detected between the two diffusion chambers and the test time can be measured.
In one embodiment of the present invention, in order to further improve the convenience of the experiment, in the diffusion index model, the volume of the first diffusion chamberVolume with the second diffusion chamberEqual, the diffusion index model is:
wherein V is the volume of the first diffusion chamber and the second diffusion chamber, and C is 。
By measuringThe molar quantity of gas diffusion to be measured at any time can be obtainedThe value of the gas C to be measured at the moment.
In one embodiment of the present invention, step 103 includes the steps of:
And S1, determining the linear relation between C and t based on the diffusion molar quantity and the diffusion index model of the gas to be detected at different moments.
And S2, fitting based on the linear relation between C and t to obtain a final diffusion index.
Specifically, by measuring differencesThe molar quantity of the gas diffusion to be measured at the moment can be obtained correspondingly to differentAnd fitting the value of the gas C to be measured at the moment by using a least square method to obtain the linear relation between C and t, wherein the slope is the final diffusion index.
In the following, two gases of hydrogen and argon are taken as specific examples, an isobaric diffusion experiment is carried out under the conditions of 3MPa and 25 ℃, the hydrogen diffusion index is calculated as an example, and the acquired data are substituted into a formula to obtain the following table through calculation.
Table 3, diffusion index detection report for hydrogen and argon.
Referring to FIG. 2, a scatter diagram is formed from times t and C, and a slope is obtained by linear fitting using a least square method= 3.78683E-05 (cm 2/s), coefficient of determination=0.9771。
In one embodiment of the present invention, the molar amount of the gas to be measured is obtained by the volume fraction of the gas to be measured by a gas chromatograph.
Setting the initial time, the mole quantity of the gas a in the first diffusion chamber asThe molar amount of the gas b in the second diffusion chamber isPass throughAfter diffusion of time, the molar amount of the gas a entering the second diffusion chamber from the first diffusion chamber isThe molar amount of the gas b entering the first diffusion chamber from the second diffusion chamber isThe volume fraction of the gas a in the first diffusion chamber was measured by gas chromatographThe volume fraction of the gas a in the second diffusion chamber isIn an ideal case, the molar quantity and the volume fraction satisfy the following relationship:
Finally, the molar amount of diffusion of the gas a is obtained as follows: ;
the molar amount of diffusion of gas b is: 。
by adopting the scheme, the molar quantity of the gas diffusion to be detected can be measured by using a gas chromatograph.
It is to be understood that the various embodiments or examples described in this specification and the features of the various embodiments or examples may be combined and combined by persons skilled in the art without contradiction.
By the method for quantitatively evaluating the gas diffusion capability in the rock, provided by the embodiment of the invention, the molar concentration of the gas a in the first diffusion chamber is as follows at the initial momentWherein, the method comprises the steps of, wherein,In terms of the molar amount of the gas a,The molar concentration of the gas b in the second diffusion chamber is equal to the volume of the first diffusion chamberWherein, the method comprises the steps of, wherein,In terms of the molar amount of the gas b,Is the volume of the second diffusion chamber.
At a certain time, when the molar amount of the gas a diffused into the second diffusion from the first diffusion chamber is x and the molar amount of the gas b diffused into the first diffusion from the second diffusion chamber is y, the molar concentration of the gas a in the first diffusion chamber isThe molar concentration of the gas a in the second diffusion chamber isThe molar concentration of the gas b in the second diffusion chamber isThe molar concentration of the gas b in the first diffusion chamber isThe molar concentration difference of the gas a between the first diffusion chamber and the second diffusion chamber is-And the molar concentration difference of the gas b between the first diffusion chamber and the second diffusion chamber isThe diffusion index is obtained according to the molar concentration variation and the diffusion index model, compared with the volume concentration, the molar concentration is taken as a dimensional physical quantity, and the diffusion flux of the gas in the rock can be reflected more intuitively and accurately, so that the diffusion capacity of the gas in the rock can be quantitatively represented more objectively and scientifically, reliable geological theory data is provided for the closure of a cover layer in a resource gas reservoir and the evaluation of a gas storage medium, and the method has important significance for researching the mechanism of the resource gas reservoir.
It should be noted that the above-mentioned embodiments are merely for illustrating the technical solution of the present invention, and not for limiting the same, and although the present invention has been described in detail with reference to the above-mentioned embodiments, it should be understood by those skilled in the art that the technical solution described in the above-mentioned embodiments may be modified or some technical features may be equivalently replaced, and these modifications or substitutions do not make the essence of the corresponding technical solution deviate from the spirit and scope of the technical solution of the embodiments of the present invention.