CN104535475A - Carbonate rock microstructure determination method and device - Google Patents
Carbonate rock microstructure determination method and device Download PDFInfo
- Publication number
- CN104535475A CN104535475A CN201510009254.4A CN201510009254A CN104535475A CN 104535475 A CN104535475 A CN 104535475A CN 201510009254 A CN201510009254 A CN 201510009254A CN 104535475 A CN104535475 A CN 104535475A
- Authority
- CN
- China
- Prior art keywords
- hole
- image
- carbonatite
- segment
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title abstract description 21
- 239000011435 rock Substances 0.000 title abstract description 21
- 239000002245 particle Substances 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 17
- 239000010430 carbonatite Substances 0.000 claims description 103
- 238000010586 diagram Methods 0.000 claims description 97
- 239000011148 porous material Substances 0.000 claims description 43
- 238000006243 chemical reaction Methods 0.000 claims description 34
- 238000009826 distribution Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000029777 axis specification Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 abstract description 3
- 150000004649 carbonic acid derivatives Chemical class 0.000 abstract 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 19
- 229910052753 mercury Inorganic materials 0.000 description 19
- 238000005266 casting Methods 0.000 description 10
- 239000002734 clay mineral Substances 0.000 description 5
- 241000272525 Anas platyrhynchos Species 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 230000035699 permeability Effects 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000006378 damage Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004090 dissolution Methods 0.000 description 2
- 239000003822 epoxy resin Substances 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 238000002493 microarray Methods 0.000 description 2
- 229920000647 polyepoxide Polymers 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 230000007096 poisonous effect Effects 0.000 description 1
- 229920003217 poly(methylsilsesquioxane) Polymers 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000037452 priming Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- 235000002639 sodium chloride Nutrition 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The embodiment of the invention provides a carbonate rock microstructure determination method and device. The method includes the steps of obtaining an electronic computer X-ray CT image of a carbonate rock sample and conducting binarization processing on the CT image, wherein the pixel number of holes of the carbonate rock sample in the image obtained after binarization processing is different from the pixel number of particles of the carbonate rock sample; calculating a total area of all holes according to the pixel number of the holes, and calculating the area of the CT image of the carbonate rock sample according to the pixel number of the holes and the pixel number of the particles; determining a specific value of the total area of the holes to the area of the CT image of the carbonate rock sample to be the porosity. The scheme is based on the CT image of the carbonate rock sample, the area of the holes and the area of the CT image of the carbonate rock sample can be calculated in a quantitative mode, the porosity of the carbonate rock sample can be further calculated accurately and reliability, and quantitative and reliable analysis on a microstructure of a carbonate rock reservoir is better achieved.
Description
Technical field
The present invention relates to carbonate reservoir Study on Microstructure technical field, particularly a kind of defining method of carbonatite micromechanism and device.
Background technology
At present, to the method for the hole of carbonate reservoir, hole, the research of seam microstructure characteristic be familiar with less, the means such as casting body flake, scanning electron microscope and pressure mercury that main employing is conventional, but, all there are some defects in these means, makes to carry out reliable analysis to the hole of carbonate reservoir, hole, seam micromechanism quantitatively.
Such as, based on the microarrays techniques in the reservoir study of casting body flake and scanning electron-microscopy means.Pore cast thin slice is dried through washing oil by underground oil water rock, to be put into colored organic glass or epoxy resin in body of casting instrument by priming by vacuum in the hole crack of rock, again after rock particles is removed in High Temperature High Pressure solidification, acidifying, only retain pore framework, be placed on basis of microscopic observation to after the process of pore framework plated surface mould, directly can observe the microscopic appearance information such as the Areal porosity in core wafer, hole, venturi, pore throat coordination number and detrital component.The method is simple, with low cost, is the common method in current RESERVOIR PORE STRUCTURE research.But there is following defect in this analytical approach:
1) utilize casting body flake image to ask for reservoir porosity, need professional technique personage to go interpretation, large to the technical merit experience dependence of operating personnel, the thickness of casting body flake, the difference of flatness directly can affect the quality of qualification result; And in real work, because of the difference that operating personnel, camera parameter are arranged etc., may there is difference in analysis result;
2) recognition image represents two dimension local hole, representative poor.Because casting body flake yardstick is too little, the ken after amplification is very little, and the nonuniformity of rock is very strong, therefore, is difficult to make evaluation to the micromechanism of carbonatite based on the very little ken;
3) casting body flake is in the colored organic glass of hole crack grouting of rock or the process of epoxy resin, destruction can be produced to the clay mineral etc. of original hole and inside thereof, produce artificial induction's hole and crack, make the pore framework formed not reflect real carbonate porosity architectural feature.
Such as, based on the microarrays techniques in the reservoir study of the means of mercury intrusion method.Conventional mercury intrusion method, based on bundle of capillary tubes pack, supposes that porous medium is made up of the capillary bundle that diameter is different.The mercury of non-wetted phase is injected rarefied rock sample space when going, the capillary resistance of blowhole system to mercury must be overcome.Mercury is as non-wetted phase, and in the injection process of mercury, first mercury enter in larger pore constriction; Along with injection pressure constantly increases, namely mercury constantly enter less hole.The capillary pressure of rock sample is all relevant with pore diameter distribution, pore volume distribution, pore throat connected relation, factor of porosity, permeability and fluid saturation etc.Therefore, utilize hollow billet curve can analyze throat distribution and the permeability contribution rate distribution characteristics of different permeability rock core, inquire into permeability, the average parameter of pore structure such as throat radius, sorting coefficient.But also there is following shortcoming in it:
1) press in mercury experimentation, the mercury of use is poisonous, can work the mischief to human body and environment;
2), when pressing mercury experiment, can produce destruction to the clay mineral etc. of original hole and inside thereof in the injection process of mercury, make not reflect real rock pore structure feature, and injection pressure is higher, this impact is larger.
Summary of the invention
Embodiments provide a kind of defining method of carbonatite micromechanism, to realize carrying out quantitatively, reliably analyzing to the micromechanism of carbonate reservoir.The method comprises: the CT CT image obtaining carbonatite sample, and binary conversion treatment is carried out to CT image, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample is different from the pixel count of the particle of carbonatite sample; Calculate the total area of all holes according to the pixel count of hole, calculate the CT image area of described carbonatite sample according to the pixel count of hole and the pixel count of particle; By the ratio of the CT image area of the total area of hole and carbonatite sample, be defined as factor of porosity.
In one embodiment, also comprise: will be profile diagram through the image Hole of binary conversion treatment and the graphics process of hole, wherein, profile diagram is the figure closed that the outline line of the figure of hole or hole is formed, and the pixel count through the image Hole of binary conversion treatment is identical with the pixel count of hole; According to the distribution mode of distribution mode determination hole and pore gap on carbonatite sample of profile diagram.
In one embodiment, also comprise: for each profile diagram, profile diagram is divided into multiple segment, in each segment, set up the ellipse tangent with figure block edge, short-axis direction oval in multiple segment is identical; The length of each segment minor axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, obtain and be the minor axis length of the ellipse tangent with this profile diagram edge, be defined as the minor axis length of hole or hole.
In one embodiment, also comprise: the length of each segment major axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, obtain and be the long axis length of the ellipse tangent with this profile diagram edge; According to the long axis length of profile diagram and the magnitude relationship of minor axis length, determine the areal shape of hole or hole.
In one embodiment, also comprise: according to the profile diagram of hole and pore gap, obtain the skeleton diagram of hole and pore gap, the figure of the expression hole that skeleton diagram is made up of line segment or the topological structure of hole; According to the coordinate of skeleton diagram on the CT image of carbonatite sample, determine hole or the position of hole on carbonatite sample.
The embodiment of the present invention additionally provides a kind of determining device of carbonatite micromechanism, to realize carrying out quantitatively, reliably analyzing to the micromechanism of carbonate reservoir.This device comprises: image processing module, for obtaining the CT CT image of carbonatite sample, and binary conversion treatment is carried out to CT image, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample is different from the pixel count of the particle of carbonatite sample; Areal calculation module, for calculating the total area of all holes according to the pixel count of hole, calculates the CT image area of described carbonatite sample according to the pixel count of hole and the pixel count of particle; Factor of porosity determination module, for the ratio of the CT image area by the total area of hole and carbonatite sample, is defined as factor of porosity.
In one embodiment, described device also comprises: profile diagram processing module, for will be profile diagram through the image Hole of binary conversion treatment and the graphics process of hole, wherein, profile diagram is the figure closed that the outline line of the figure of hole or hole is formed, and the pixel count through the image Hole of binary conversion treatment is identical with the pixel count of hole; Distribution mode determination module, for the distribution mode of distribution mode determination hole and pore gap on carbonatite sample according to profile diagram.
In one embodiment, also comprise: segment divides module, for for each profile diagram, profile diagram is divided into multiple segment, in each segment, set up the ellipse tangent with figure block edge, short-axis direction oval in multiple segment is identical; Minor axis determination module, the length of each segment minor axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, that obtain and be the minor axis length of the ellipse tangent with this profile diagram edge, be defined as the minor axis length of hole or hole.
In one embodiment, also comprising: major axis determination module, for the length of each segment major axis being multiplied with this segment area, obtaining the product of each segment, by the product addition of multiple segment, obtain and be the long axis length of the ellipse tangent with this profile diagram edge; Areal shape determination module, for according to the long axis length of profile diagram and the magnitude relationship of minor axis length, determines the areal shape of hole or hole.
In one embodiment, also comprise: skeleton diagram processing module, for the profile diagram according to hole and pore gap, obtain the skeleton diagram of hole and pore gap, wherein, the figure of the expression hole that skeleton diagram is made up of line segment or the topological structure of hole; Position determination module, for according to the coordinate of skeleton diagram on the CT image of carbonatite sample, determines hole or the position of hole on carbonatite sample.
In embodiments of the present invention, by carrying out CT scan to carbonatite sample, obtain the CT image of carbonatite sample, and binary conversion treatment is carried out to CT image, pixel count based on the image mesoporosity through binary conversion treatment is different from the pixel count of particle, calculate the total area of hole and the CT image area of carbonatite sample, and then can based on the CT image area determining hole porosity of the total area of hole and carbonatite sample.Determining hole porosity is carried out owing to adopting the CT image of carbonatite sample, as compared to means such as casting body flake, scanning electron microscope and pressure mercury in prior art, can avoid damaging the clay mineral etc. of original hole and inside thereof, to reflect real carbonate porosity architectural feature; Avoid the dependence to operating personnel's technical merit, avoid using the objectionable impuritiess such as mercury in operating process; Also can obtain the CT image of carbonatite sample entirety, improve the representativeness of CT image; Simultaneously, based on the CT image of carbonatite sample, the area of hole and the CT image area of carbonatite sample can be calculated quantitatively, and then accurately, reliably can calculate the factor of porosity of carbonatite sample, be conducive to carrying out quantitatively, reliably analyzing to the micromechanism of carbonate reservoir.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the process flow diagram of the defining method of a kind of carbonatite micromechanism that the embodiment of the present invention provides;
Fig. 2 is the schematic diagram of the profile diagram of a kind of hole and pore gap that the embodiment of the present invention provides;
Fig. 3 is the classification schematic diagram of a kind of hole and pore gap that the embodiment of the present invention provides;
Fig. 4 is the schematic diagram of the skeleton diagram of a kind of hole and pore gap that the embodiment of the present invention provides;
Fig. 5 is the profile diagram of a kind of hole and pore gap that the embodiment of the present invention provides and the schematic diagram of skeleton diagram overlap;
Fig. 6 is the structured flowchart of the determining device of a kind of carbonatite micromechanism that the embodiment of the present invention provides.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with embodiment and accompanying drawing, the present invention is described in further details.At this, exemplary embodiment of the present invention and illustrating for explaining the present invention, but not as a limitation of the invention.
In embodiments of the present invention, provide a kind of defining method of carbonatite micromechanism, as shown in Figure 1, the method comprises:
Step 101: the CT CT image obtaining carbonatite sample, and binary conversion treatment is carried out to image, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample is different from the pixel count of the particle of carbonatite sample;
Step 102: the total area calculating all holes according to the pixel count of hole, calculates the CT image area of described carbonatite sample according to the pixel count of hole and the pixel count of particle;
Step 103: by the ratio of the CT image area of the total area of hole and carbonatite sample, be defined as factor of porosity.
Flow process is as shown in Figure 1 known, in embodiments of the present invention, by carrying out CT scan to carbonatite sample, obtain the CT image of carbonatite sample, and binary conversion treatment is carried out to CT image, pixel count based on the image mesoporosity through binary conversion treatment is different from the pixel count of particle, calculates the total area of hole and the CT image area of carbonatite sample, and then can based on the CT image area determining hole porosity of the total area of hole and carbonatite sample.Determining hole porosity is carried out owing to adopting the CT image of carbonatite sample, as compared to means such as casting body flake, scanning electron microscope and pressure mercury in prior art, can avoid damaging the clay mineral etc. of original hole and inside thereof, to reflect real rock salt pore structure characteristic; Avoid the dependence to operating personnel's technical merit, avoid using the objectionable impuritiess such as mercury in operating process; Also can obtain the CT image of carbonatite sample entirety, improve the representativeness of CT image; Simultaneously, based on the CT image of carbonatite sample, the area of hole and the CT image area of carbonatite sample can be calculated quantitatively, and then accurately, reliably can calculate the factor of porosity of carbonatite sample, be conducive to carrying out quantitatively, reliably analyzing to the micromechanism of carbonate reservoir.
During concrete enforcement, when binary conversion treatment is carried out to CT image, the graythresh method calculated threshold in image processing techniques can be adopted, then the binary conversion treatment that im2bw carries out image is called, gray level image by 256 brightness degrees is chosen (such as by suitable threshold values, it can be 60 that threshold value is chosen) and obtain the binary image that still can reflect integral image and local feature, filtering and Edge contrast can also be carried out to image.
Concrete, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample and the particle of carbonatite sample are (in the image of binary conversion treatment, all particles are shown as the form of matrix) pixel count different, in binary image, the gray-scale value of image only has 2 kinds, 255 and 0, can look be (such as a setting by particle, white), namely gray-scale value is 0, and so the gray-scale value of hole is exactly 255, like this in binary image, white portion is the matrix of particle composition, and black region is hole.When calculating the CT image area of the total area of hole and carbonatite sample, can calculate according to pixel count, such as, the total area of hole is calculated according to the pixel count of hole, the area of particle is calculated according to the pixel count of particle, the total area of hole is added with the area of particle, the area sum obtained can be the CT image area of carbonatite sample, achieve the CT image area of the total area and the carbonatite sample quantitatively, reliably calculating hole, and then factor of porosity can be calculated exactly.
During concrete enforcement, in order to the distribution situation of hole and pore gap can be determined, in the present embodiment, said method also comprises: will be profile diagram (as shown in Figure 2) through the image Hole of binary conversion treatment and the graphics process of hole, wherein, profile diagram is the figure closed that the outline line of the figure of hole or hole is formed, and the pixel count through the image Hole of binary conversion treatment is identical with the pixel count of hole; According to the distribution mode of distribution mode determination hole and pore gap on carbonatite sample of profile diagram.Concrete, bwboundaries method in image processing techniques can be adopted to obtain profile diagram, irregular figure as shown in Figure 2 is the profile diagram of hole or hole, specifically according to the diameter of figure, sorting technique with reference to Fig. 3 determines hole or hole, such as, the diameter of figure be less than 20 millimeters be hole, the diameter of figure be greater than 20 millimeters be hole.
During concrete enforcement, in order to determine the size of the hole and pore gap of different scale quantitatively, in the present embodiment, for each profile diagram, profile diagram is divided into multiple segment, in each segment, sets up the ellipse tangent with figure block edge, short-axis direction oval in multiple segment is identical; The length of each segment minor axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, obtain and be the minor axis length of the ellipse tangent with this profile diagram edge, be defined as the minor axis length of hole or hole.Concrete, after determining the minor axis length of hole or hole, namely can the size of determining hole hole or hole, and then can come hole and pore gap classification according to the diameter length in Fig. 3, such as, diameter be less than 0.025 millimeter for aperture, diameter be less than 0.5 millimeter be greater than 0.025 millimeter for middle dissolution pore, the like, hole is divided into aperture, middle dissolution pore and macropore, hole is divided into duck eye, Zhong Dong and great Dong.After hole or Pore classification, can also count the frequency of occurrences of certain hole or hole on carbonatite sample, such as, the frequency of statistics aperture, can calculate according to following formula:
wherein η represents the frequency size of aperture, and n represents the number of aperture, and N represents the total number of hole; Can also add up the frequency of certain hole according to this formula, such as, when adding up the frequency of duck eye, η represents the frequency size of duck eye, and n represents the number of duck eye, and N represents the total number of hole.
Concrete, in order to determine the form of hole and pore gap quantitatively, in the present embodiment, said method also comprises: the length of each segment major axis be multiplied with this segment area, obtain the product of each segment, by the product addition of multiple segment, obtain and be the long axis length of the ellipse tangent with this profile diagram edge; According to the long axis length of profile diagram and the magnitude relationship of minor axis length, determine the areal shape of hole or hole.Such as, when the long axis length of profile diagram is greater than minor axis length, then the figure of profile diagram is oval, and the form of the hole that profile diagram is corresponding or hole is oval; When the long axis length of profile diagram equals minor axis length, then the figure of profile diagram is circular, and the form of the hole that profile diagram is corresponding or hole is circular.
During concrete enforcement, in order to determine the particular location of hole and pore gap on carbonatite sample, in the present embodiment, according to the profile diagram of hole and pore gap, obtain the skeleton diagram (as shown in Figure 4,5) of hole and pore gap, according to the coordinate of skeleton diagram on the CT image of carbonatite sample, determine hole or the position of hole on carbonatite sample.Concrete, when scanning carbonatite sample and obtaining CT image, set up the corresponding relation of position on carbonatite sample and CT image coordinate, according to the coordinate of skeleton diagram on the CT image of carbonatite sample, determine hole or the position of hole on carbonatite sample.
The defining method of above-mentioned carbonatite micromechanism can solve the quantitative description problem of micro-nano hole to different scale rock core and venturi well, extensively can use the exploratory development of the unconventionaloil pool such as shale gas, coal-seam gas, fine and close oil gas.
Based on same inventive concept, additionally provide a kind of determining device of carbonatite micromechanism in the embodiment of the present invention, as described in the following examples.The principle of dealing with problems due to the determining device of carbonatite micromechanism is similar to the defining method of carbonatite micromechanism, therefore the enforcement of the determining device of carbonatite micromechanism see the enforcement of the defining method of carbonatite micromechanism, can repeat part and repeats no more.Following used, term " unit " or " module " can realize the software of predetermined function and/or the combination of hardware.Although the device described by following examples preferably realizes with software, hardware, or the realization of the combination of software and hardware also may and conceived.
Fig. 6 is a kind of structured flowchart of the determining device of the carbonatite micromechanism of the embodiment of the present invention, as shown in Figure 6, comprising: image processing module 601, areal calculation module 602 and factor of porosity determination module 603, be described this structure below.
Image processing module 601, for obtaining the CT CT image of carbonatite sample, and binary conversion treatment is carried out to CT image, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample is different from the pixel count of the particle of carbonatite sample;
Areal calculation module 602, is connected with image processing module 601, for calculating the total area of all holes according to the pixel count of hole, calculates the CT image area of described carbonatite sample according to the pixel count of hole and the pixel count of particle;
Factor of porosity determination module 603, is connected with areal calculation module 602, for the ratio of the CT image area by the total area of hole and carbonatite sample, is defined as factor of porosity.
In one embodiment, described device also comprises: profile diagram processing module, for will be profile diagram through the image Hole of binary conversion treatment and the graphics process of hole, wherein, profile diagram is the figure closed that the outline line of the figure of hole or hole is formed, and the pixel count through the image Hole of binary conversion treatment is identical with the pixel count of hole; Distribution mode determination module, is connected with profile diagram processing module, for the distribution mode of distribution mode determination hole and pore gap on carbonatite sample according to profile diagram.
In one embodiment, also comprise: segment divides module, for for each profile diagram, profile diagram is divided into multiple segment, in each segment, set up the ellipse tangent with figure block edge, short-axis direction oval in multiple segment is identical; Minor axis determination module, model calling is divided with segment, the length of each segment minor axis is multiplied with this segment area, obtain the product of each segment, by the product addition of multiple segment, that obtain and be the minor axis length of the ellipse tangent with this profile diagram edge, be defined as the minor axis length of hole or hole.
In one embodiment, also comprising: major axis determination module, for the length of each segment major axis being multiplied with this segment area, obtaining the product of each segment, by the product addition of multiple segment, obtain and be the long axis length of the ellipse tangent with this profile diagram edge; Areal shape determination module, for according to the long axis length of profile diagram and the magnitude relationship of minor axis length, determines the areal shape of hole or hole.
In one embodiment, also comprise: skeleton diagram processing module, for the profile diagram according to hole and pore gap, obtain the skeleton diagram of hole and pore gap, wherein, the figure of the expression hole that skeleton diagram is made up of line segment or the topological structure of hole; Position determination module, is connected with skeleton diagram processing module, for according to the coordinate of skeleton diagram on the CT image of carbonatite sample, determines hole or the position of hole on carbonatite sample.
In embodiments of the present invention, by carrying out CT scan to carbonatite sample, obtain the CT image of carbonatite sample, and binary conversion treatment is carried out to CT image, pixel count based on the image mesoporosity through binary conversion treatment is different from the pixel count of the matrix that particle forms, calculate the total area of hole and the CT image area of carbonatite sample, and then can based on the CT image area determining hole porosity of the total area of hole and carbonatite sample.Determining hole porosity is carried out owing to adopting the CT image of carbonatite sample, as compared to means such as casting body flake, scanning electron microscope and pressure mercury in prior art, can avoid damaging the clay mineral etc. of original hole and inside thereof, to reflect real carbonate porosity architectural feature; Avoid the dependence to operating personnel's technical merit, avoid using the objectionable impuritiess such as mercury in operating process; Also can obtain the CT image of carbonatite sample entirety, improve the representativeness of CT image; Simultaneously, based on the CT image of carbonatite sample, the area of hole and the CT image area of carbonatite sample can be calculated quantitatively, and then accurately, reliably can calculate the factor of porosity of carbonatite sample, be conducive to analyzing reliably the micromechanism of carbonate reservoir.
Obviously, those skilled in the art should be understood that, each module of the above-mentioned embodiment of the present invention or each step can realize with general calculation element, they can concentrate on single calculation element, or be distributed on network that multiple calculation element forms, alternatively, they can realize with the executable program code of calculation element, thus, they can be stored and be performed by calculation element in the storage device, and in some cases, step shown or described by can performing with the order be different from herein, or they are made into each integrated circuit modules respectively, or the multiple module in them or step are made into single integrated circuit module to realize.Like this, the embodiment of the present invention is not restricted to any specific hardware and software combination.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. a defining method for carbonatite micromechanism, is characterized in that, comprising:
Obtain the CT CT image of carbonatite sample, and carry out binary conversion treatment to CT image, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample is different from the pixel count of the particle of carbonatite sample;
Calculate the total area of all holes according to the pixel count of hole, calculate the CT image area of described carbonatite sample according to the pixel count of hole and the pixel count of particle;
By the ratio of the CT image area of the total area of hole and carbonatite sample, be defined as factor of porosity.
2. the method for claim 1, is characterized in that, also comprises:
To be profile diagram through the image Hole of binary conversion treatment and the graphics process of hole, wherein, profile diagram be the figure closed that the outline line of the figure of hole or hole is formed, and the pixel count through the image Hole of binary conversion treatment is identical with the pixel count of hole;
According to the distribution mode of distribution mode determination hole and pore gap on carbonatite sample of profile diagram.
3. method as claimed in claim 2, is characterized in that, also comprise:
For each profile diagram, profile diagram is divided into multiple segment, in each segment, sets up the ellipse tangent with figure block edge, short-axis direction oval in multiple segment is identical;
The length of each segment minor axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, obtain and be the minor axis length of the ellipse tangent with this profile diagram edge, be defined as the minor axis length of hole or hole.
4. method as claimed in claim 3, is characterized in that, also comprise:
The length of each segment major axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, obtain and be the long axis length of the ellipse tangent with this profile diagram edge;
According to the long axis length of profile diagram and the magnitude relationship of minor axis length, determine the areal shape of hole or hole.
5. the method according to any one of claim 2 to 4, is characterized in that, also comprises:
According to the profile diagram of hole and pore gap, obtain the skeleton diagram of hole and pore gap, the figure of the expression hole that skeleton diagram is made up of line segment or the topological structure of hole;
According to the coordinate of skeleton diagram on the CT image of carbonatite sample, determine hole or the position of hole on carbonatite sample.
6. a determining device for carbonatite micromechanism, is characterized in that, comprising:
Image processing module, for obtaining the CT CT image of carbonatite sample, and binary conversion treatment is carried out to CT image, in the image of binary conversion treatment, the pixel count of the hole of carbonatite sample is different from the pixel count of the particle of carbonatite sample;
Areal calculation module, for calculating the total area of all holes according to the pixel count of hole, calculates the CT image area of described carbonatite sample according to the pixel count of hole and the pixel count of particle;
Factor of porosity determination module, for the ratio of the CT image area by the total area of hole and carbonatite sample, is defined as factor of porosity.
7. device as claimed in claim 6, it is characterized in that, described device also comprises:
Profile diagram processing module, for will be profile diagram through the image Hole of binary conversion treatment and the graphics process of hole, wherein, profile diagram is the figure closed that the outline line of the figure of hole or hole is formed, and the pixel count through the image Hole of binary conversion treatment is identical with the pixel count of hole;
Distribution mode determination module, for the distribution mode of distribution mode determination hole and pore gap on carbonatite sample according to profile diagram.
8. device as claimed in claim 7, is characterized in that, also comprise:
Segment divides module, for for each profile diagram, profile diagram is divided into multiple segment, in each segment, set up the ellipse tangent with figure block edge, and short-axis direction oval in multiple segment is identical;
Minor axis determination module, the length of each segment minor axis is multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, that obtain and be the minor axis length of the ellipse tangent with this profile diagram edge, be defined as the minor axis length of hole or hole.
9. device as claimed in claim 8, is characterized in that, also comprise:
Major axis determination module, for the length of each segment major axis being multiplied with this segment area, obtains the product of each segment, by the product addition of multiple segment, obtain and be the long axis length of the ellipse tangent with this profile diagram edge;
Areal shape determination module, for according to the long axis length of profile diagram and the magnitude relationship of minor axis length, determines the areal shape of hole or hole.
10. the device according to any one of claim 7 to 9, is characterized in that, also comprises:
Skeleton diagram processing module, for the profile diagram according to hole and pore gap, obtains the skeleton diagram of hole and pore gap, wherein, and the figure of the expression hole that skeleton diagram is made up of line segment or the topological structure of hole;
Position determination module, for according to the coordinate of skeleton diagram on the CT image of carbonatite sample, determines hole or the position of hole on carbonatite sample.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510009254.4A CN104535475B (en) | 2015-01-08 | 2015-01-08 | Method and device for determining microstructure of carbonate rock |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510009254.4A CN104535475B (en) | 2015-01-08 | 2015-01-08 | Method and device for determining microstructure of carbonate rock |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104535475A true CN104535475A (en) | 2015-04-22 |
| CN104535475B CN104535475B (en) | 2018-02-06 |
Family
ID=52851049
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510009254.4A Active CN104535475B (en) | 2015-01-08 | 2015-01-08 | Method and device for determining microstructure of carbonate rock |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104535475B (en) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105021511A (en) * | 2015-07-28 | 2015-11-04 | 昆明理工大学 | Perforated tipping paper air permeability detection method applied to tobacco industry |
| CN105092616A (en) * | 2015-09-07 | 2015-11-25 | 中国兵器科学研究院宁波分院 | Method for measuring minor detail feature sizes in industrial CT (computerized tomography) detection |
| CN105300868A (en) * | 2015-11-24 | 2016-02-03 | 昆明理工大学 | Detection method for air permeability of perforated tipping paper used in tobacco industry |
| CN105510204A (en) * | 2015-12-03 | 2016-04-20 | 大连理工大学 | Permeability predication method based on CT (computed tomography) images |
| CN106204463A (en) * | 2015-05-05 | 2016-12-07 | 中国石油化工股份有限公司 | A kind of processing method of Core Scanning Image |
| CN106383078A (en) * | 2016-09-20 | 2017-02-08 | 中国石油天然气股份有限公司 | Method and device for determining water drive efficiency of rock |
| CN106442271A (en) * | 2016-11-18 | 2017-02-22 | 杨康 | Rock core permeability simulating method and rock core permeability simulating device |
| CN107133630A (en) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of method that carbonate porosity type is judged based on scan image |
| CN108061697A (en) * | 2017-12-06 | 2018-05-22 | 石家庄铁道大学 | Soil body three-dimensional pore space rate computational methods |
| CN108734751A (en) * | 2018-05-08 | 2018-11-02 | 中山大学 | A kind of porous media microscopical structure model generating method based on random fractal |
| CN108961332A (en) * | 2017-05-23 | 2018-12-07 | 中国石油化工股份有限公司 | Sandstone image representative region area size determines method and device |
| CN108956416A (en) * | 2018-06-06 | 2018-12-07 | 中国地质大学（北京） | A method of based on Matlab analysis tight sandstone reservoir hole characterization |
| CN109374624A (en) * | 2018-12-12 | 2019-02-22 | 山东大学 | Porous pavement choke detecting method and system |
| CN109447944A (en) * | 2018-09-21 | 2019-03-08 | 中国石油天然气股份有限公司 | Lithofacies identification method and system for carbonate rock |
| CN109490165A (en) * | 2018-11-02 | 2019-03-19 | 中国石油化工股份有限公司 | The method for characterizing the non-fabric selectivity reservoir space of carbonate rock |
| CN111257188A (en) * | 2018-11-30 | 2020-06-09 | 中国石油天然气股份有限公司 | Multi-scale characterization method and device for carbonate rock pore system |
| CN111738978A (en) * | 2020-03-27 | 2020-10-02 | 中国石油化工股份有限公司 | Evaluation method and device for connectivity of pore throats of reservoir, electronic equipment and storage medium |
| CN112084884A (en) * | 2020-08-19 | 2020-12-15 | 厦门理工学院 | Scanning electron microscope image pore identification method, terminal device and storage medium |
| CN113447419A (en) * | 2021-06-28 | 2021-09-28 | 哈尔滨工业大学 | Dividing system for porous medium communicated pore structure unit |
| CN113670791A (en) * | 2021-08-03 | 2021-11-19 | 中国地质大学（北京） | Quantitative analysis method for pore space and particle surface area of reservoir |
| CN117593299A (en) * | 2024-01-18 | 2024-02-23 | 北京大学 | Method, device, equipment and medium for evaluating space effectiveness of lamellar shale reservoir |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101403683B (en) * | 2008-11-17 | 2010-12-01 | 长安大学 | Method for analyzing porous asphalt mixture gap structure by using CT technology |
| CN101556703B (en) * | 2009-05-16 | 2011-01-26 | 中国石油大学(华东) | Method for establishing network model based on serial section image |
| JP5588713B2 (en) * | 2010-03-31 | 2014-09-10 | イビデン株式会社 | Simulation apparatus and computer program therefor |
| CN102175832B (en) * | 2011-01-10 | 2013-10-16 | 中国石油天然气股份有限公司 | Method for determining optimal saturation calculation model of typical reservoir |
| CN102222359B (en) * | 2011-05-24 | 2012-09-26 | 中国石油天然气股份有限公司 | Method for reconstructing three-dimensional pore structure of rock core |
| CN103698803B (en) * | 2012-09-27 | 2017-02-08 | 中国石油天然气股份有限公司 | Rock pore structure characterization method and device |
| CN103575751A (en) * | 2013-09-10 | 2014-02-12 | 济南大学 | Characterization method of pore structure of porous concrete based on CT (Computerized Tomography) technology |
-
2015
- 2015-01-08 CN CN201510009254.4A patent/CN104535475B/en active Active
Cited By (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106204463A (en) * | 2015-05-05 | 2016-12-07 | 中国石油化工股份有限公司 | A kind of processing method of Core Scanning Image |
| CN105021511A (en) * | 2015-07-28 | 2015-11-04 | 昆明理工大学 | Perforated tipping paper air permeability detection method applied to tobacco industry |
| CN105092616B (en) * | 2015-09-07 | 2017-09-12 | 中国兵器科学研究院宁波分院 | Industry CT detects medium and small minutia dimension measurement method |
| CN105092616A (en) * | 2015-09-07 | 2015-11-25 | 中国兵器科学研究院宁波分院 | Method for measuring minor detail feature sizes in industrial CT (computerized tomography) detection |
| CN105300868A (en) * | 2015-11-24 | 2016-02-03 | 昆明理工大学 | Detection method for air permeability of perforated tipping paper used in tobacco industry |
| CN105510204A (en) * | 2015-12-03 | 2016-04-20 | 大连理工大学 | Permeability predication method based on CT (computed tomography) images |
| CN107133630A (en) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of method that carbonate porosity type is judged based on scan image |
| CN107133630B (en) * | 2016-02-29 | 2020-04-03 | 中国石油化工股份有限公司 | Method for judging carbonate rock pore type based on scanned image |
| CN106383078A (en) * | 2016-09-20 | 2017-02-08 | 中国石油天然气股份有限公司 | Method and device for determining water drive efficiency of rock |
| CN106383078B (en) * | 2016-09-20 | 2019-05-07 | 中国石油天然气股份有限公司 | Method and device for determining water drive efficiency of rock |
| CN106442271A (en) * | 2016-11-18 | 2017-02-22 | 杨康 | Rock core permeability simulating method and rock core permeability simulating device |
| CN108961332A (en) * | 2017-05-23 | 2018-12-07 | 中国石油化工股份有限公司 | Sandstone image representative region area size determines method and device |
| CN108961332B (en) * | 2017-05-23 | 2021-06-15 | 中国石油化工股份有限公司 | Method and device for determining area size of representative region of sandstone image |
| CN108061697A (en) * | 2017-12-06 | 2018-05-22 | 石家庄铁道大学 | Soil body three-dimensional pore space rate computational methods |
| CN108061697B (en) * | 2017-12-06 | 2020-09-25 | 石家庄铁道大学 | Method for calculating three-dimensional porosity of soil body |
| CN108734751B (en) * | 2018-05-08 | 2022-02-11 | 中山大学 | Porous medium mesoscopic structure model generation method based on random fractal |
| CN108734751A (en) * | 2018-05-08 | 2018-11-02 | 中山大学 | A kind of porous media microscopical structure model generating method based on random fractal |
| CN108956416A (en) * | 2018-06-06 | 2018-12-07 | 中国地质大学（北京） | A method of based on Matlab analysis tight sandstone reservoir hole characterization |
| CN109447944A (en) * | 2018-09-21 | 2019-03-08 | 中国石油天然气股份有限公司 | Lithofacies identification method and system for carbonate rock |
| CN109447944B (en) * | 2018-09-21 | 2020-08-11 | 中国石油天然气股份有限公司 | Lithofacies identification method and system for carbonate rock |
| CN109490165A (en) * | 2018-11-02 | 2019-03-19 | 中国石油化工股份有限公司 | The method for characterizing the non-fabric selectivity reservoir space of carbonate rock |
| CN111257188A (en) * | 2018-11-30 | 2020-06-09 | 中国石油天然气股份有限公司 | Multi-scale characterization method and device for carbonate rock pore system |
| CN109374624A (en) * | 2018-12-12 | 2019-02-22 | 山东大学 | Porous pavement choke detecting method and system |
| CN111738978A (en) * | 2020-03-27 | 2020-10-02 | 中国石油化工股份有限公司 | Evaluation method and device for connectivity of pore throats of reservoir, electronic equipment and storage medium |
| CN112084884A (en) * | 2020-08-19 | 2020-12-15 | 厦门理工学院 | Scanning electron microscope image pore identification method, terminal device and storage medium |
| CN113447419A (en) * | 2021-06-28 | 2021-09-28 | 哈尔滨工业大学 | Dividing system for porous medium communicated pore structure unit |
| CN113447419B (en) * | 2021-06-28 | 2022-10-25 | 哈尔滨工业大学 | Dividing system for porous medium communicated pore structure unit |
| CN113670791A (en) * | 2021-08-03 | 2021-11-19 | 中国地质大学（北京） | Quantitative analysis method for pore space and particle surface area of reservoir |
| CN117593299A (en) * | 2024-01-18 | 2024-02-23 | 北京大学 | Method, device, equipment and medium for evaluating space effectiveness of lamellar shale reservoir |
| CN117593299B (en) * | 2024-01-18 | 2024-04-19 | 北京大学 | Method, device, equipment and medium for evaluating space effectiveness of lamellar shale reservoir |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104535475B (en) | 2018-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104535475A (en) | Carbonate rock microstructure determination method and device | |
| CN104751473A (en) | Device and method for determining multiscale porous characteristics of carbonatite | |
| CN107817199B (en) | A kind of construction method of tight sand multi-scale porosity model and application | |
| US11506650B2 (en) | Method for automatic quantitative statistical distribution characterization of dendrite structures in a full view field of metal materials | |
| CN105487121B (en) | Multiple dimensioned digital cores method is built based on CT scan image and Electrical imaging image co-registration | |
| CN107507173B (en) | No-reference definition evaluation method and system for full-slice image | |
| CN104732227B (en) | A kind of Location Method of Vehicle License Plate based on definition and luminance evaluation | |
| CN113227755A (en) | Assay accuracy improvement | |
| CN109285222A (en) | The building of organic shale high-resolution digital rock core and analysis method | |
| CN103984939B (en) | A kind of sample visible component sorting technique and system | |
| CN110021030A (en) | A kind of segmentation threshold of material of rock and soil digital picture determines method | |
| Liu et al. | A new method for threshold determination of gray image | |
| CN106446952A (en) | Method and apparatus for recognizing score image | |
| CN114926463A (en) | Production quality detection method suitable for chip circuit board | |
| CN108734714B (en) | Method for analyzing carbonate rock structure based on Matlab | |
| CN112150430B (en) | Numerical analysis method for digital image by using rock microstructure | |
| CN103400139A (en) | Method for identifying concrete crack characteristic information | |
| JP2023511869A (en) | Method and system for automatic identification and grading of low multiplication acid etching defects by machine vision | |
| KR101245923B1 (en) | Method of analyzing cell structures and their components | |
| Xin et al. | A fine extraction algorithm for image-based surface cracks in underwater dams | |
| CN103196404A (en) | Electronic probe rock slice sample positioning device | |
| CN106093035B (en) | A kind of microspur video image identification method that the soil body develops | |
| CN103093231A (en) | Automatic identification and quantitative analysis method of second phase in alloy | |
| CN104732555A (en) | Image edge detection method based on Sobel operator | |
| CN113191330A (en) | Method for identifying area growth pores by fusing secondary electrons and back scattering electron images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20151013 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Applicant after: PetroChina Company Limited Address before: 102249 Beijing city Changping District Road No. 18 Applicant before: China University of Petroleum (Beijing) |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |