CN113484220A - Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment - Google Patents

Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment Download PDF

Info

Publication number
CN113484220A
CN113484220A CN202110843857.XA CN202110843857A CN113484220A CN 113484220 A CN113484220 A CN 113484220A CN 202110843857 A CN202110843857 A CN 202110843857A CN 113484220 A CN113484220 A CN 113484220A
Authority
CN
China
Prior art keywords
size distribution
pore size
distribution curve
water content
sample
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110843857.XA
Other languages
Chinese (zh)
Other versions
CN113484220B (en
Inventor
冯东
吴克柳
李靖
陈掌星
赵文
王鼎涵
李相方
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202110843857.XA priority Critical patent/CN113484220B/en
Publication of CN113484220A publication Critical patent/CN113484220A/en
Application granted granted Critical
Publication of CN113484220B publication Critical patent/CN113484220B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The specification discloses a method, a device and electronic equipment for determining shale organic matter and inorganic matter distribution, wherein the method comprises the following steps: calculating a water saturation scatter diagram of pores in the target clay sample at the second water content according to the pore size distribution curve of the target clay sample at the first water content and the pore size distribution curve of the target clay sample at the second water content; calculating the pore size distribution curve of inorganic substances in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content and the second water content, the pore size distribution curve of the target clay sample under the first water content and a water saturation scatter diagram of pores in the target clay sample; and calculating to obtain the aperture distribution curve of the organic matters in the shale sample under the first water content according to the aperture distribution curve of the shale sample under the first water content and the aperture distribution curve of the inorganic matters in the shale sample under the first water content. According to the scheme, the pore size distribution curves of organic matters and inorganic matters in the shale sample can be distinguished.

Description

Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment
Technical Field
The application relates to the technical field of shale oil-gas exploration and development, in particular to a method, a device and electronic equipment for determining shale organic matter and inorganic matter distribution.
Background
Since the 21 st century, the world energy supply and demand contradiction is aggravated by the increasingly reduced and growing energy demand of conventional oil and gas resources, and the successful development of shale gas has a certain relieving effect on the contradiction. In the actual development process, accurate evaluation of the pore throat characteristics of the reservoir is the basis and the premise for ensuring that the shale gas reservoir can be developed efficiently, accurate pore distribution characteristics can provide key basic data for shale gas reserve prediction (adsorbed gas content and free gas content) and seepage capability evaluation, and the method has important significance for making a subsequent development scheme and predicting the capacity.
Compared with a conventional oil and gas reservoir, the shale reservoir has two characteristics: firstly, the nano-scale pores are developed in a large quantity, so that the pore distribution characteristics in real shale are difficult to accurately reflect through a conventional mercury intrusion test, and the pore distribution of the shale is often measured through a liquid nitrogen adsorption experiment; secondly, the shale has complicated pore types, and the most remarkable characteristic is that the shale contains inorganic pores and organic pores with opposite wettability (to water).
Generally, the organic matter of the high-maturity shale has hydrophobicity, methane gas is mainly generated in pores, one part of the methane gas is adsorbed on the surfaces of the pores of the organic matter in an adsorption state, and the other part of the methane gas exists in the center of the pores in a free state; the seepage of methane in the organic pores is mainly based on bulk phase transmission and surface diffusion. The inorganic pores of the shale have hydrophilicity, the surface of the inorganic pores usually has a water film with a certain thickness or is completely filled with formation water, and methane in the inorganic pores generally exists in a free state and flows through a bulk phase transmission mode. Therefore, the occurrence characteristics and seepage modes of shale organic matter pore and inorganic matter pore fluid have obvious differences, and the differences further influence the evaluation of the gas-containing performance and seepage capability of the actual reservoir. Therefore, quantitative identification of real shale organic matter pores and inorganic matter pores is very important for improving the overall effect of shale gas reservoir development.
However, the traditional liquid nitrogen adsorption can only test the whole pore structure of the shale, and can not effectively distinguish organic matters and inorganic matters in the shale, so that the evaluation of the gas storage capacity and the seepage capacity of a real shale reservoir is greatly influenced.
Disclosure of Invention
An object of the embodiments of the present application is to provide a method, an apparatus and an electronic device for determining the organic and inorganic pore size distribution of shale, so as to distinguish the organic and inorganic pore size distribution curves in a shale sample.
To solve the technical problem, embodiments of the present disclosure provide a method for determining organic matter and inorganic matter distribution of shale, including: acquiring a pore size distribution curve of a target clay sample under a first water content and a pore size distribution curve under a second water content; calculating a water saturation scatter diagram of pores in the target clay sample at a second water content according to the pore size distribution curve of the target clay sample at a first water content and the pore size distribution curve of the target clay sample at a second water content; acquiring a pore size distribution curve of the shale sample under a first water content and a pore size distribution curve under a second water content; calculating a pore size distribution curve of inorganic substances in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content, the pore size distribution curve of the target clay sample under the second water content, and a water saturation scatter diagram of pores in the target clay sample; and calculating to obtain a pore size distribution curve of organic matters in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content and the pore size distribution curve of inorganic matters in the shale sample under the first water content.
Embodiments of the present disclosure also provide an apparatus for determining organic and inorganic substance distribution of shale, including: the first acquisition module is used for acquiring a pore size distribution curve of the target clay sample under a first water content and a pore size distribution curve under a second water content; the first calculation module is used for calculating a water saturation scatter diagram of pores in the target clay sample under the second water content according to the pore size distribution curve of the target clay sample under the first water content and the pore size distribution curve of the pore size distribution curve under the second water content; the second acquisition module is used for acquiring a pore size distribution curve of the shale sample under the first water content and a pore size distribution curve under the second water content; the second calculation module is used for calculating the inorganic substance pore size distribution curve in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content, the pore size distribution curve of the target clay sample under the second water content, and the water saturation scatter diagram of pores in the target clay sample; and the third calculation module is used for calculating the pore size distribution curve of the organic matters in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content and the pore size distribution curve of the inorganic matters in the shale sample under the first water content.
According to the method, the device and the electronic equipment for determining the distribution of organic matters and inorganic matters in shale, which are provided by the embodiment of the specification, according to a quantitative determination theory established and deduced by an inventor, the water saturation of target clay is calculated, and then according to the conclusion that the change of the shale differential pore volume under the water containing condition is mainly caused by the change of the inorganic matter differential pore volume in shale, the water saturation is utilized to calculate the pore size distribution curve of inorganic matters in a shale sample, so that the pore size distribution curve of the organic matters in the shale sample is obtained. The method can obtain the pore diameter distribution curve quantitatively expressed by the pore characteristics of organic matters and inorganic matters in the shale sample.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 illustrates a flow diagram of a method of determining shale organic and inorganic mass distribution in accordance with an embodiment of the present description;
FIG. 2A shows the corresponding liquid nitrogen adsorption curves for clay samples at a first water content and at a second water content;
FIG. 2B illustrates liquid nitrogen adsorption curves corresponding to the shale sample at the first water content and the second water content;
FIG. 3A shows corresponding pore size distribution curves for clay samples at a first water content and at a second water content;
FIG. 3B illustrates corresponding pore size distribution curves for a shale sample at a first water cut and at a second water cut;
FIG. 4 illustrates a water saturation profile of inorganic porosity in a shale sample at a second water cut;
FIG. 5 illustrates a calculated plot of pore size distribution for organic and inorganic matter in a shale sample at a first water cut;
fig. 6 illustrates a flow diagram of another method of determining shale organic and inorganic mass distribution in accordance with an embodiment of the present description;
fig. 7 is a schematic block diagram illustrating an apparatus for determining organic and inorganic material distributions of shale in accordance with an embodiment of the present disclosure;
FIG. 8 shows a functional block diagram of an electronic device according to an embodiment of the invention.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments obtained by a person of ordinary skill in the art without any inventive work based on the embodiments in the present application shall fall within the scope of protection of the present application.
Generally, the organic matter of the high-maturity shale has hydrophobicity, methane gas is mainly generated in pores, one part of the methane gas is adsorbed on the surfaces of the pores of the organic matter in an adsorption state, and the other part of the methane gas exists in the center of the pores in a free state; the seepage of methane in the organic pores is mainly based on bulk phase transmission and surface diffusion. The inorganic pores of the shale have hydrophilicity, the surface of the inorganic pores usually has a water film with a certain thickness or is completely filled with formation water, and methane in the inorganic pores generally exists in a free state and flows through a bulk phase transmission mode. Therefore, the occurrence characteristics and seepage modes of shale organic matter pore and inorganic matter pore fluid have significant differences.
Considering the fluid distribution difference caused by the wettability difference of organic matters and inorganic matters in the shale, the inventor establishes a quantitative determination method for the distribution of organic matter pores and inorganic matter pores in the shale. The theoretical basis of the quantitative determination is described in detail below.
The shale integrally comprises organic and inorganic pores, and the differential pore volume at the first water content can be expressed as:
Figure BDA0003179817560000041
wherein S, or and in represent shale, organic matter and inorganic matter, respectively. D represents the pore diameter in nm; Δ V/Δ D is the differential pore volume, representing the ratio of the volume increase to the difference between the two pore diameters defining the increase, and is given by 10-3cc/g/nm。
The differential pore volume of the shale at the second water cut (the second water cut being greater than the first water cut) may be expressed as:
Figure BDA0003179817560000042
since organic matter in shale is hydrophobic, the pores with hydrophobicity generally do not adsorb water vapor. Assuming that inorganic substances in the shale do not influence organic pores of the shale after absorbing water, the method comprises the following steps
Figure BDA0003179817560000043
The coupling formulae (1), (2) and (3) can be obtained
Figure BDA0003179817560000044
From the equation (4), it can be seen that the change of the differential pore volume in the shale under the first water content and the second water content is mainly caused by the change of the inorganic differential pore volume in the shale. Therefore, the key of the problem is converted into a change rule for describing the pore volume of the inorganic substance at the first water content and at the first water content.
The examples of the present specification use clay having the same composition as the main component of the inorganic substances in the shale instead of the inorganic substances in the shale. By a liquid nitrogen adsorption experiment and a BJH method, pore size distribution curves corresponding to the clay under the first water content and the clay under the first water content can be obtained. Using this curve we can get
Figure BDA0003179817560000045
Wherein D is the diameter of the inorganic pores, Sw-in(D) The water saturation at inorganic pore size D.
By combining formula (4) and formula (5) again, the process can be obtained
Figure BDA0003179817560000046
From equation (5), it can be seen that if the differential pore volume in the shale at the first water content and at the first water content is known, and the water saturation of the inorganic pores is known, the differential pore volume of the clay at the first water content, i.e., the pore size distribution curve corresponding to the clay, can be obtained.
Further, as can be seen from equation (1), if the differential pore volume of the shale and the inorganic matter at the first water content is known, the differential pore volume of the organic matter at the first water content can be obtained.
Based on the theoretical basis, the embodiment of the specification provides a method for determining the distribution of organic matters and inorganic matters in shale, and the differential pore volume of the organic matters and the inorganic matters in the shale can be measured. As shown in fig. 1, the method includes the following steps.
S110: and acquiring a pore size distribution curve of the target clay sample under the first water content and a pore size distribution curve under the second water content.
The pore size distribution curve is used to represent the trend of differential pore volume as a function of pore size. Differential pore volume represents the amount of increase in pore volume by one unit of increase in pore diameter, i.e., Δ V/Δ D as described in the theoretical basis analysis above, calculated as the ratio of the volume increase to the difference between the two pore diameters that determine the increase. Fig. 3A shows a pore size distribution curve for a target clay sample at a first water content (i.e., the target clay in an environment having a humidity RH of 0) and a pore size distribution curve at a second water content (i.e., the target clay in an environment having a humidity RH of 0).
The first water content is the water content when the water content of the target clay sample is less than a preset value. The preset values may be values close to 0, such as 0.3, 0.2, 0.1, 0.05, 0.01, i.e. the target clay sample at the first water content is a dry sample. The second moisture content may be a moisture content of the target clay sample in a higher moisture content state in which the sample moisture content is more stable.
The target clay sample at the first moisture content may be obtained by drying the clay sample, i.e. the second moisture content is greater than the first moisture content.
The target clay sample at the second water content may be obtained by: placing the target clay sample at the first water content in a container with a saturated salt solutionAnd in an environment, measuring the weight of the sample once every preset time period until the weight difference of the sample is less than a preset percentage when the sample is weighed twice. Wherein the saturated salt solution may be K2SO4、CuSO4Or Na2HPO4. These saturated salt solutions correspond to an ambient humidity of 0.98. Alternatively, other moisture saturated salt solutions may be used.
Of course, the target clay sample with the second water content can be obtained by other means, such as obtaining an environment with constant humidity by using other principles of the device, and placing the sample in the environment.
The target clay sample in the examples of the present specification may be any one of several equal small samples separated from the same clay sample.
Of course, in the embodiments of the present description, the first water content and the second water content may be any values, and it is only necessary to ensure that the second water content is greater than the first water content.
The pore size distribution curve of the target clay sample at the first water content can be obtained by the following method: obtaining a liquid nitrogen adsorption curve obtained by putting a target clay sample with a first water content into a specific surface area and porosity analyzer to carry out a liquid nitrogen adsorption experiment; and according to the liquid nitrogen adsorption curve, calculating by using a BJH (BJH) method to obtain a pore size distribution curve of the target clay sample under the first water content. This method may also be used to obtain a pore size distribution curve for the target clay sample at the second water content.
The specific surface area and porosity analyzer will give the values of both the adsorption curve and the desorption curve. Fig. 2A shows a schematic of an adsorption curve and a desorption curve. The method of obtaining the pore size distribution curve by BJH based on the adsorption or desorption curve is a complex calculation process, which is well known to those skilled in the art and will not be described herein again.
As the liquid nitrogen adsorption experiment can give out an adsorption curve and a desorption curve together, the pore size distribution curve can be obtained based on the adsorption curve and also can be obtained based on the desorption curve. For smaller pores, the accuracy of the pore size distribution curve obtained by using the adsorption curve is higher.
S120: and calculating a water saturation scatter diagram of pores in the target clay sample at the second water content according to the pore size distribution curve of the target clay sample at the first water content and the pore size distribution curve of the target clay sample at the second water content.
The water saturation corresponding to pores with different diameters in the target clay sample at the second water content can be calculated according to the formula (5) above, namely the following formula:
Figure BDA0003179817560000061
wherein S isw-in(D) Representing the water saturation of the target clay sample at a pore diameter D,
Figure BDA0003179817560000062
the differential pore volume is expressed in terms of,
Figure BDA0003179817560000063
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content,
Figure BDA0003179817560000064
representing the differential pore volume on the pore size distribution curve of the target clay sample at the second water content.
And then generating a water saturation scatter diagram of the pores in the target clay sample under the second water content according to the water saturation corresponding to the pores with different diameters.
The calculated water saturation scatter plot may be as shown in fig. 4.
S130: and obtaining a pore size distribution curve of the shale sample under the first water content and a pore size distribution curve under the second water content.
Please refer to step S110 for a description of the pore size distribution curve. Fig. 3B shows a pore size distribution curve for a shale sample at a first water cut (RH ═ 0) and a pore size distribution curve at a second water cut (RH ═ 0.98).
The first water content is the water content of the shale sample when the water content is smaller than a preset value. The preset values may be values close to 0, such as 0.3, 0.2, 0.1, 0.05, 0.01, i.e. the shale sample at the first water cut is a dry sample. The second water content may be a water content of the shale sample in a high water content state in which the water content of the sample is relatively stable.
The shale sample with the first water content can be obtained by drying the shale sample, namely the second water content is larger than the first water content.
The shale sample at the second water cut may be obtained by: placing the shale sample with the first water content in a sealed environment with saturated salt solution, and measuring the weight of the sample once every preset time period until the weight difference of the sample is less than a preset percentage when the shale sample is weighed twice. Wherein the saturated salt solution may be K2SO4、CuSO4Or Na2HPO4. The humidity of these saturated salt solutions was 0.98. Alternatively, other moisture saturated salt solutions may be used.
Of course, the shale sample with the second water content may be obtained by other means, for example, by using equipment to obtain an environment with constant humidity by using other principles, and placing the sample in the environment.
The shale sample in the examples of this specification may be any of several small equivalent samples separated from the same shale sample.
Of course, in the embodiments of the present description, the first water content and the second water content may be any values, and it is only necessary to ensure that the second water content is greater than the first water content.
Obtaining a pore size distribution curve of the shale sample under the first water content, wherein the following method can be adopted: obtaining a liquid nitrogen adsorption curve obtained by putting the shale sample with the first water content into a specific surface area and porosity analyzer to carry out a liquid nitrogen adsorption experiment; and according to the liquid nitrogen adsorption curve, calculating by using a BJH (BJH) method to obtain a pore size distribution curve of the shale sample under the first water content. The method may also be used to obtain a pore size distribution curve of the shale sample at the second water cut.
The specific surface area and porosity analyzer will give the values of both the adsorption curve and the desorption curve. Fig. 2B shows a schematic of the adsorption and desorption curves. The method of obtaining the pore size distribution curve by BJH based on the adsorption or desorption curve is a complex calculation process, which is well known to those skilled in the art and will not be described herein again.
As the liquid nitrogen adsorption experiment can give out an adsorption curve and a desorption curve together, the pore size distribution curve can be obtained based on the adsorption curve and also can be obtained based on the desorption curve. For smaller pores, the accuracy of the pore size distribution curve obtained by using the adsorption curve is higher.
The first water content of the shale sample is the same as the first water content of the target clay sample, and the second water content of the shale sample is the same as the second water content of the target clay sample.
S140: and calculating the pore size distribution curve of inorganic substances in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content, the pore size distribution curve of the target clay sample under the second water content, and the water saturation scatter diagram of pores in the target clay sample.
A plurality of differential pore volumes of inorganic material in the shale sample at the first water content may be calculated according to equation (6), i.e., the following equation:
Figure BDA0003179817560000071
wherein S isw-in(D) Representing the water saturation of the target clay sample at a pore diameter D,
Figure BDA0003179817560000072
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content,
Figure BDA0003179817560000073
representing a differential pore volume on a pore size distribution curve of the shale sample at a first water cut,
Figure BDA0003179817560000074
representing a differential pore volume on a pore size distribution curve of the shale sample at a second water cut.
A pore size distribution curve for an inorganic material in the shale sample at the first water content is then generated based on the plurality of differential pore volumes.
S150: and calculating to obtain a pore size distribution curve of organic matters in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content and the pore size distribution curve of inorganic matters in the shale sample under the first water content.
A plurality of differential pore volumes of organic matter in the shale sample at the first water cut may be calculated according to equation (1), i.e., the following equation:
Figure BDA0003179817560000081
wherein the content of the first and second substances,
Figure BDA0003179817560000082
representing a differential pore volume on a pore size distribution curve of the shale sample at a first water cut,
Figure BDA0003179817560000083
representing a differential pore volume on a pore size distribution curve of the shale sample at a second water cut,
Figure BDA0003179817560000084
representing the differential pore volume on the pore size distribution curve of the target clay sample at the first water content.
And then generating a pore size distribution curve of organic matter in the shale sample at the first water content according to the plurality of differential pore volumes.
According to the method, the device and the electronic equipment for determining the distribution of organic matters and inorganic matters in shale, which are provided by the embodiment of the specification, according to a quantitative determination theory established and deduced by an inventor, the water saturation of target clay is calculated, and then according to the theory that the change of the differentiated pore volume of shale is mainly caused by the change of the differentiated pore volume of inorganic matters in shale, the water saturation is utilized to calculate the pore size distribution curve of the inorganic matters in the shale sample, so that the pore size distribution curve of the organic matters in the shale sample is obtained. The method can obtain the pore diameter distribution curve quantitatively expressed by the pore characteristics of organic matters and inorganic matters in the shale sample.
Fig. 6 shows a specific process of the shale assay method performed by the inventors. As shown in fig. 6, the process includes the following steps.
S701: and selecting 100-150 meshes of shale samples and clay samples for drying treatment, wherein the drying temperature is 200 ℃ and the drying time is 48 hours.
S702: placing saturated K in a big beaker2SO4 solution and sealed, thereby leaving the ambient humidity RH in the beaker at 0.98 at all times.
S703: preparing about 20g of each of the dried shale sample and the clay sample, placing the shale sample and the clay sample in a small beaker, placing the small beaker in a large beaker, sealing, measuring the weight once every 6 hours until the mass difference is less than 3% when the shale sample and the clay sample are weighed twice, and thus obtaining the shale sample and the clay sample under the condition of high water content.
In this step, care was taken to wipe off the water bead outside the beaker each time it was weighed.
S704: selecting dried shale samples and clay samples of 0.5-0.8g respectively to carry out liquid nitrogen adsorption experiments to obtain liquid nitrogen adsorption and desorption curves of the samples, and selecting adsorption branches to obtain pore size distribution curves corresponding to the shale samples and the clay samples under the dry condition by using a BJH method.
S705: selecting 0.5-0.8g of each of the shale sample and the clay sample with high water content to carry out a liquid nitrogen adsorption test to obtain liquid nitrogen adsorption and desorption curves of the samples, and selecting an adsorption branch to obtain pore size distribution curves corresponding to the shale sample and the clay sample with high water content by using a BJH method.
The liquid nitrogen adsorption curve corresponding to the dried and high-water-content clay sample obtained in this example is shown in fig. 2A, the liquid nitrogen adsorption curve corresponding to the dried and high-water-content shale sample is shown in fig. 2B, the pore size distribution curve corresponding to the dried and high-water-content clay sample is shown in fig. 3A, and the data are shown in table 1; the pore size distribution curve for the dried and high moisture shale samples is shown in fig. 3B, and the data is shown in table 2.
TABLE 1 data sheet for pore size distribution of clay samples dried and having high water content
Figure BDA0003179817560000091
TABLE 2 pore size distribution data sheet for dry and high moisture content shale samples
Figure BDA0003179817560000101
S706: and (4) calculating to obtain a water saturation scatter diagram of pores in the clay sample by using the formula (5) according to the pore size distribution curve corresponding to the dried and high-water-content clay sample.
The calculated water saturation distribution quantification of the pores in the clay sample is shown in figure 4.
S707: according to the pore size distribution curve corresponding to the dried and high-moisture-content shale sample and the pore size distribution curve corresponding to the dried clay sample, replacing inorganic substances in the shale sample with the clay sample, and calculating by using the formula (6) to obtain the pore size distribution curve corresponding to the inorganic substances in the shale sample.
The calculated pore size distribution curve corresponding to inorganic substances in the shale sample is shown in fig. 5, and the data is shown in table 3.
S708: and (3) calculating to obtain a pore size distribution curve corresponding to organic matters in the shale sample according to the pore size distribution curve corresponding to the dried shale sample, the pore size distribution curve corresponding to inorganic matters in the shale sample and the formula (1).
The calculated pore size distribution curve corresponding to the organic matter in the shale sample is shown in fig. 5, and the data is shown in table 3.
The analysis results of this example show that: for the experimental sample, shale organic matter developed primarily small pores less than 20nm, and shale inorganic matter developed primarily large pores (10-200 nm). The pore volume of the shale organic matter is 4.78 multiplied by 10-3cm3G, pore volume of inorganic substance 16.03X 10-3cm3The organic matter has a pore volume fraction of 22.52%.
It should be noted that due to the difference between the experimental test and the theoretical calculation, the inorganic pore volume corresponding to a certain scale obtained by calculation may be larger than the shale pore volume, that is, the organic pore volume may be smaller than 0, and at this time, the organic pore volume is determined to be equal to 0.
TABLE 3 data table of distribution of organic and inorganic substances in real shale sample
Figure BDA0003179817560000111
The invention provides a method for quantitatively identifying distribution characteristics of organic matter pores and inorganic matter pores in shale, which improves the traditional experimental method on the basis of fully considering wettability difference of organic matter and inorganic matter in shale, simultaneously measures pore distribution curves of shale and clay samples under drying conditions and water-containing conditions, and establishes a method for determining the distribution characteristics of the organic matter pores and the inorganic matter pores in real shale by comparing the distribution characteristics of water in the pores of the shale and the clay. The method has important significance for accurately evaluating the gas-containing performance and the seepage capability of the reservoir.
The embodiment of the present specification also provides an apparatus for determining shale organic matter and inorganic matter distribution, which can be used for executing the method described in fig. 1. As shown in fig. 7, the apparatus includes a first obtaining module 10, a first calculating module 20, a second obtaining module 30, a second calculating module 40, and a third calculating module 50.
The first obtaining module 10 is configured to obtain a pore size distribution curve of the target clay sample at a first water content and a pore size distribution curve of the target clay sample at a second water content.
The first calculation module 20 is configured to calculate a water saturation scatter plot of pores in the target clay sample at the second water content according to the pore size distribution curve of the target clay sample at the first water content and the pore size distribution curve of the pore size distribution curve at the second water content.
The second obtaining module 30 is configured to obtain a pore size distribution curve of the shale sample at the first water content and a pore size distribution curve at the second water content.
The second calculation module 40 is configured to calculate a pore size distribution curve of inorganic substances in the shale sample at the first water content according to the pore size distribution curve of the shale sample at the first water content and the pore size distribution curve at the second water content, the pore size distribution curve of the target clay sample at the first water content, and a water saturation scatter diagram of pores in the target clay sample.
The third calculating module 50 is configured to calculate a pore size distribution curve of organic matter in the shale sample at the first water content according to the pore size distribution curve of the shale sample at the first water content and a pore size distribution curve of inorganic matter in the shale sample at the first water content.
In some embodiments, the second percentage of water content is greater than the first percentage of water content.
In some embodiments, the first computation module 20 includes a first computation submodule and a first generation submodule.
The first calculation submodule is used for calculating the water saturation corresponding to pores with different diameters in the target clay sample under the second water content according to the following formula:
Figure BDA0003179817560000121
wherein S isw-in(D) Representing the water saturation of the target clay sample at a pore diameter D,
Figure BDA0003179817560000122
the differential pore volume is expressed in terms of,
Figure BDA0003179817560000123
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content,
Figure BDA0003179817560000124
representing the differential pore volume on the pore size distribution curve of the target clay sample at the second water content.
And the first generation submodule is used for generating a water saturation scatter diagram of the pores in the target clay sample under the second water content according to the water saturation corresponding to the pores with different diameters.
In some embodiments, the second calculation module 40 includes a second calculation submodule and a second generation submodule.
The second calculation submodule is configured to calculate a plurality of differential pore volumes of inorganic matter in the shale sample at the first water cut according to the following equation:
Figure BDA0003179817560000131
wherein S isw-in(D) Representing the water saturation of the target clay sample at a pore diameter D,
Figure BDA0003179817560000132
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content,
Figure BDA0003179817560000133
representing a differential pore volume on a pore size distribution curve of the shale sample at a first water cut,
Figure BDA0003179817560000134
indicating differential porosity on a pore size distribution curve of the shale sample at a second water cutVolume.
A second generation submodule configured to generate a pore size distribution curve for an inorganic substance in the shale sample at the first water cut based on the plurality of differential pore volumes.
In some embodiments, third computing module 50 includes a third computing sub-module and a third generating sub-module.
The third calculation submodule is configured to calculate a plurality of differential pore volumes of organic matter in the shale sample at the first water content according to the following formula:
Figure BDA0003179817560000135
wherein the content of the first and second substances,
Figure BDA0003179817560000136
representing a differential pore volume on a pore size distribution curve of the shale sample at a first water cut,
Figure BDA0003179817560000137
representing a differential pore volume on a pore size distribution curve of the shale sample at a second water cut,
Figure BDA0003179817560000138
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content;
and the third generation submodule is used for generating a pore size distribution curve of organic matters in the shale sample under the first water content according to the plurality of differential pore volumes.
In some embodiments, the shale sample and/or the target clay sample at the second water cut is obtained by: placing the shale sample and/or the target clay sample at the first water content in a sealed environment with a saturated salt solution, and measuring the weight of the sample once every preset time period until the weight difference of the samples is less than a preset percentage when the samples are weighed twice. The saturated salt solution may be K2SO4、CuSO4Or Na2HPO4
In some embodiments, the first acquisition module 10 includes an acquisition submodule and a fourth calculation submodule.
The acquisition submodule is used for acquiring a liquid nitrogen adsorption curve obtained by putting a target clay sample with the first water content into a specific surface area and porosity analyzer to carry out a liquid nitrogen adsorption experiment. And the fourth calculation submodule is used for calculating to obtain a pore size distribution curve of the target clay sample under the first water content by utilizing a BJH method according to the liquid nitrogen adsorption curve.
The description and effects related to the device for determining the distribution of organic matters and inorganic matters in shale can be understood by referring to the method for determining the distribution of organic matters and inorganic matters in shale, and are not described in detail.
An embodiment of the present invention further provides an electronic device, as shown in fig. 8, the electronic device may include a processor 81 and a memory 82, where the processor 81 and the memory 82 may be connected by a bus or in another manner, and fig. 8 takes the connection by the bus as an example.
Processor 81 may be a Central Processing Unit (CPU). The Processor 81 may also be other general purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components, or combinations thereof.
The memory 82, which is a non-transitory computer readable storage medium, may be used to store non-transitory software programs, non-transitory computer executable programs, and modules, such as program instructions/modules corresponding to the method of determining shale organic and inorganic mass distribution in embodiments of the present invention (e.g., the first acquisition module 10, the first calculation module 20, the second acquisition module 30, the second calculation module 40, and the third calculation module 50 shown in fig. 7). The processor 81 executes the non-transitory software programs, instructions and modules stored in the memory 82 to execute various functional applications of the processor and data processing, namely, to implement the method for determining shale organic matter and inorganic matter distribution in the above method embodiment.
The memory 82 may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created by the processor 81, and the like. Further, the memory 82 may include high speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory 82 may optionally include memory located remotely from the processor 81, which may be connected to the processor 81 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The one or more modules are stored in the memory 82 and, when executed by the processor 81, perform a method of determining shale organic and inorganic mass distribution as in the embodiment shown in fig. 1.
The details of the electronic device may be understood with reference to the corresponding related description and effects in the embodiment of fig. 1, and are not described herein again.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD), a Solid State Drive (SSD), or the like; the storage medium may also comprise a combination of memories of the kind described above.
In the 50 s of the 20 th century, improvements in a technology could clearly be distinguished between improvements in hardware (e.g., improvements in circuit structures such as diodes, transistors, switches, etc.) and improvements in software (improvements in process flow). However, as technology advances, many of today's process flow improvements have been seen as direct improvements in hardware circuit architecture. Designers almost always obtain the corresponding hardware circuit structure by programming an improved method flow into the hardware circuit. Thus, it cannot be said that an improvement in the process flow cannot be realized by hardware physical modules. For example, a Programmable Logic Device (PLD), such as a Field Programmable Gate Array (FPGA), is an integrated circuit whose Logic functions are determined by programming the Device by a user. A digital system is "integrated" on a PLD by the designer's own programming without requiring the chip manufacturer to design and fabricate a dedicated integrated circuit chip 2. Furthermore, nowadays, instead of manually making an Integrated Circuit chip, such Programming is often implemented by "logic compiler" software, which is similar to a software compiler used in program development and writing, but the original code before compiling is also written by a specific Programming Language, which is called Hardware Description Language (HDL), and HDL is not only one but many, such as abel (advanced Boolean Expression Language), ahdl (alternate Language Description Language), traffic, pl (core unified Programming Language), HDCal, JHDL (Java Hardware Description Language), langue, Lola, HDL, laspam, hardbyscript Description Language (vhr Description Language), and the like, which are currently used by Hardware compiler-software (Hardware Description Language-software). It will also be apparent to those skilled in the art that hardware circuitry that implements the logical method flows can be readily obtained by merely slightly programming the method flows into an integrated circuit using the hardware description languages described above.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments can be referred to each other, and each embodiment focuses on the differences from the other embodiments.
The systems, devices, modules or units described in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present application may be essentially or partially implemented in the form of software products, which may be stored in a storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and include instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods of some parts of the embodiments of the present application.
The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet-type devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
Although the present application has been described in terms of embodiments, those of ordinary skill in the art will recognize that there are numerous variations and permutations of the present application without departing from the spirit of the application, and it is intended that the appended claims encompass such variations and permutations without departing from the spirit of the application.

Claims (10)

1. A method for determining organic and inorganic substance distribution of shale, which is characterized by comprising the following steps:
acquiring a pore size distribution curve of a target clay sample under a first water content and a pore size distribution curve under a second water content;
calculating a water saturation scatter diagram of pores in the target clay sample at a second water content according to the pore size distribution curve of the target clay sample at a first water content and the pore size distribution curve of the target clay sample at a second water content;
acquiring a pore size distribution curve of the shale sample under a first water content and a pore size distribution curve under a second water content;
calculating a pore size distribution curve of inorganic substances in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content, the pore size distribution curve of the target clay sample under the second water content, and a water saturation scatter diagram of pores in the target clay sample;
and calculating to obtain a pore size distribution curve of organic matters in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content and the pore size distribution curve of inorganic matters in the shale sample under the first water content.
2. The method of claim 1, wherein the second moisture content is greater than the first moisture content.
3. The method of claim 1, wherein calculating a water saturation scatter plot of pores in the target clay sample at a second water cut from the pore size distribution curve of the target clay sample at a first water cut and the pore size distribution curve of the pore size distribution curve at the second water cut comprises:
calculating the water saturation corresponding to pores with different diameters in the target clay sample at the second water content according to the following formula:
Figure FDA0003179817550000011
wherein S isw-in(D) Representing the water saturation of the target clay sample at a pore diameter D,
Figure FDA0003179817550000012
the differential pore volume is expressed in terms of,
Figure FDA0003179817550000013
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content,
Figure FDA0003179817550000014
representing a differential pore volume on a pore size distribution curve of the target clay sample at a second water content;
and generating a water saturation scatter diagram of the pores in the target clay sample under the second water content according to the water saturation corresponding to the pores with different diameters.
4. The method of claim 1, wherein calculating a pore size distribution curve of inorganic matter in the shale sample at the first water cut based on the pore size distribution curve of the shale sample at the first water cut and the pore size distribution curve at the second water cut, the pore size distribution curve of the target clay sample at the first water cut, and a water saturation scatter plot of pores in the target clay sample comprises:
calculating a plurality of differential pore volumes of inorganic matter in the shale sample at the first water content according to the following equation:
Figure FDA0003179817550000021
wherein S isw-in(D) Representing the water saturation of the target clay sample at a pore diameter D,
Figure FDA0003179817550000022
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content,
Figure FDA0003179817550000023
representing a differential pore volume on a pore size distribution curve of the shale sample at a first water cut,
Figure FDA0003179817550000024
representing a differential pore volume on a pore size distribution curve of the shale sample at a second water content;
generating a pore size distribution curve for an inorganic material in the shale sample at the first water cut based on the plurality of differential pore volumes.
5. The method of claim 1, wherein calculating the pore size distribution curve of the organic matter in the shale sample at the first water cut from the pore size distribution curve of the shale sample at the first water cut and the pore size distribution curve of the inorganic matter in the shale sample at the first water cut comprises:
calculating a plurality of differential pore volumes of organic matter in the shale sample at a first water cut according to the following equation:
Figure FDA0003179817550000025
wherein the content of the first and second substances,
Figure FDA0003179817550000026
representing a differential pore volume on a pore size distribution curve of the shale sample at a first water cut,
Figure FDA0003179817550000027
representing a differential pore volume on a pore size distribution curve of the shale sample at a second water cut,
Figure FDA0003179817550000028
representing a differential pore volume on a pore size distribution curve of the target clay sample at a first water content;
and generating a pore size distribution curve of organic matter in the shale sample at the first water content according to the plurality of differential pore volumes.
6. The method of claim 1, wherein the shale sample and/or the target clay sample at a second water cut is obtained by:
placing the shale sample and/or the target clay sample at the first water content in a sealed environment with a saturated salt solution, and measuring the weight of the sample once every preset time period until the weight difference of the samples is less than a preset percentage when the samples are weighed twice.
7. The method of claim 1, wherein obtaining the pore size distribution curve of the target clay sample at the first water cut comprises:
obtaining a liquid nitrogen adsorption curve obtained by putting a target clay sample with a first water content into a specific surface area and porosity analyzer to carry out a liquid nitrogen adsorption experiment;
and according to the liquid nitrogen adsorption curve, calculating by using a BJH (BJH) method to obtain a pore size distribution curve of the target clay sample under the first water content.
8. An apparatus for determining organic and inorganic material distribution of shale, comprising:
the first acquisition module is used for acquiring a pore size distribution curve of the target clay sample under a first water content and a pore size distribution curve under a second water content;
the first calculation module is used for calculating a water saturation scatter diagram of pores in the target clay sample under the second water content according to the pore size distribution curve of the target clay sample under the first water content and the pore size distribution curve of the pore size distribution curve under the second water content;
the second acquisition module is used for acquiring a pore size distribution curve of the shale sample under the first water content and a pore size distribution curve under the second water content;
the second calculation module is used for calculating the inorganic substance pore size distribution curve in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content, the pore size distribution curve of the target clay sample under the second water content, and the water saturation scatter diagram of pores in the target clay sample;
and the third calculation module is used for calculating the pore size distribution curve of the organic matters in the shale sample under the first water content according to the pore size distribution curve of the shale sample under the first water content and the pore size distribution curve of the inorganic matters in the shale sample under the first water content.
9. An electronic device, comprising:
a memory and a processor communicatively coupled to each other, the memory having stored therein computer instructions, the processor implementing the steps of the method of any one of claims 1 to 7 by executing the computer instructions.
10. A computer storage medium storing computer program instructions which, when executed, implement the steps of the method of any one of claims 1 to 7.
CN202110843857.XA 2021-07-26 2021-07-26 Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment Active CN113484220B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110843857.XA CN113484220B (en) 2021-07-26 2021-07-26 Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110843857.XA CN113484220B (en) 2021-07-26 2021-07-26 Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment

Publications (2)

Publication Number Publication Date
CN113484220A true CN113484220A (en) 2021-10-08
CN113484220B CN113484220B (en) 2022-07-22

Family

ID=77942663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110843857.XA Active CN113484220B (en) 2021-07-26 2021-07-26 Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment

Country Status (1)

Country Link
CN (1) CN113484220B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5349528A (en) * 1990-03-12 1994-09-20 Halliburton Logging Services, Inc. Method apparatus for determination of porosity lithological composition
US5992228A (en) * 1996-10-23 1999-11-30 Dunham; Lanny L. Method for determining resistivity derived porosity and porosity derived resistivity
CN106442268A (en) * 2016-10-31 2017-02-22 中国科学技术大学 Method for detecting pore size distribution of shale mesopores
CN107194104A (en) * 2017-06-07 2017-09-22 成都理工大学 A kind of method that use log calculates inanimate matter porosity, total porosity and component in shale reservoir
US20200173902A1 (en) * 2018-03-27 2020-06-04 Min Wang Evaluation method for hydrogen-bearing components, porosity and pore size distribution of organic-rich shale
CN111337410A (en) * 2020-04-14 2020-06-26 中国石油大学(北京) Method, device, equipment and system for determining shale pore size distribution
CN112082918A (en) * 2020-08-04 2020-12-15 中国石油大学(北京) Method, device and equipment for determining porosity
CN112414917A (en) * 2020-11-03 2021-02-26 西安石油大学 Method for dividing and characterizing organic pores and inorganic pores of shale oil reservoir
CN112858136A (en) * 2021-01-29 2021-05-28 中国地质大学(武汉) Quantitative evaluation method for shale organic matter pore structure

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5349528A (en) * 1990-03-12 1994-09-20 Halliburton Logging Services, Inc. Method apparatus for determination of porosity lithological composition
US5992228A (en) * 1996-10-23 1999-11-30 Dunham; Lanny L. Method for determining resistivity derived porosity and porosity derived resistivity
CN106442268A (en) * 2016-10-31 2017-02-22 中国科学技术大学 Method for detecting pore size distribution of shale mesopores
CN107194104A (en) * 2017-06-07 2017-09-22 成都理工大学 A kind of method that use log calculates inanimate matter porosity, total porosity and component in shale reservoir
US20200173902A1 (en) * 2018-03-27 2020-06-04 Min Wang Evaluation method for hydrogen-bearing components, porosity and pore size distribution of organic-rich shale
CN111337410A (en) * 2020-04-14 2020-06-26 中国石油大学(北京) Method, device, equipment and system for determining shale pore size distribution
CN112082918A (en) * 2020-08-04 2020-12-15 中国石油大学(北京) Method, device and equipment for determining porosity
CN112414917A (en) * 2020-11-03 2021-02-26 西安石油大学 Method for dividing and characterizing organic pores and inorganic pores of shale oil reservoir
CN112858136A (en) * 2021-01-29 2021-05-28 中国地质大学(武汉) Quantitative evaluation method for shale organic matter pore structure

Also Published As

Publication number Publication date
CN113484220B (en) 2022-07-22

Similar Documents

Publication Publication Date Title
CN111337410B (en) Method, device, equipment and system for determining shale pore size distribution
Lauritzen et al. A standard test case suite for two-dimensional linear transport on the sphere
Sircar Comments on practical use of Langmuir gas adsorption isotherm model
Too et al. Review of soil water retention characteristic (SWRC) models between saturation and oven dryness
Lapshin et al. On competitive gas adsorption and absorption phenomena in thin films of ionic liquids
CN108008444B (en) A kind of method and device of determining petrophysical model
CN107818379A (en) Absorption tolerance determines method and apparatus
CN109612900B (en) Reservoir rock permeability prediction method and device and storage medium
CN111507537B (en) Method and device for predicting reserves and water invasion of water-flooding gas reservoir
CN113705868B (en) Shale gas state quantitative prediction method, system, equipment and storage medium
CN113484220B (en) Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment
CN110032744A (en) A kind of unconventional oil and gas Reserve Estimation Method and system
Xu et al. Pore scale study of gas sorption hysteresis in shale nanopores using lattice Boltzmann method
Zheng et al. An empirical soil water retention model based on probability laws for pore‐size distribution
Santos et al. A dataset of 3D structural and simulated transport properties of complex porous media
Da Silva et al. Parallel openmp and openacc mixing layer simulation
Ould-Kaddour et al. Determination of the friction coefficient of a Brownian particle by molecular-dynamics simulation
Koniorczyk et al. The influence of salt on hygrothermal behaviour of building materials
Benzi et al. Building cavities in a fluid of spherical or rod‐like particles: A contribution to the solvation free energy in isotropic and anisotropic polarizable continuum model
Barlow et al. WTAQ Version 2: A Computer Program for Analysis of Aquifer Tests in Confined and Water-table Aquifers with Alternative Representations of Drainage from the Unsaturated Zone $ nElektronische Ressource
Ben-Naim Solvation of large molecules: some exact results on the dependence on volume and surface area of the solute
Tsakiroglou et al. Pore‐structure analysis by using nitrogen sorption and mercury intrusion data
CN111999331B (en) Capillary pressure determining method and device based on nuclear magnetic test
Sidelnikov et al. Modeling and study of sorption and separation of acids in solution
Lowe et al. Modelling multi-phase halogen chemistry in the remote marine boundary layer: investigation of the influence of aerosol size resolution on predicted gas-and condensed-phase chemistry

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant