CN110057846A - Rock core internal tiny crack recognition methods based on digital picture, system, device - Google Patents
Rock core internal tiny crack recognition methods based on digital picture, system, device Download PDFInfo
- Publication number
- CN110057846A CN110057846A CN201910290206.5A CN201910290206A CN110057846A CN 110057846 A CN110057846 A CN 110057846A CN 201910290206 A CN201910290206 A CN 201910290206A CN 110057846 A CN110057846 A CN 110057846A
- Authority
- CN
- China
- Prior art keywords
- point
- slice
- rock
- obtains
- sub
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/751—Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Artificial Intelligence (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Immunology (AREA)
- Radiology & Medical Imaging (AREA)
- Pulmonology (AREA)
- Chemical & Material Sciences (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
Abstract
The invention belongs to rock behavio(u)r analysis fields, and in particular to a kind of rock core internal tiny crack recognition methods based on digital picture, system, device, it is intended in order to solve the problems, such as that rock core internal tiny crack is unrecognized.The method of the present invention includes: to obtain the first rock sample, obtains initial slice sequence by CT scan, and obtain each initial slice in initial slice sequence with reference to figure point；Crushing test is carried out to the first rock sample based on three axis pressurized equipments, each process Slice Sequence corresponding with initial slice sequence under different rate of loading is obtained by CT scan；Based on the matched method of local area image, the corresponding points with reference to figure point in each process slice are obtained, obtain the first displacement field of each process slice；Based on each process slice the first displacement field, obtain sub-pixel positioning after matching figure point, obtain second displacement field, and extract rock deformation destroy after micro-crack.The present invention effectively improves the accuracy of identification of crackle, is conducive to the extension and evolution of simulation rock crackle.
Description
Technical field
The invention belongs to rock behavio(u)r analysis fields, and in particular to a kind of rock core internal tiny crack knowledge based on digital picture
Other method, system, device.
Background technique
The extension and distribution of crackle and the morphological feature of crackle are indispensable one of study of rocks destruction characteristic
Point.At present detection rock crackle forming distribution there are two ways to: one is based on rock face crack observation, another kind be by
CT scan obtains the underbead crack distribution of rock.The internal modification of the available rock of second method destroys, but due to setting
Standby hardware limitation can only identify that width is higher than the crackle of device resolution size.Currently, the available rock of high energy acclerator CT
Internal crackle.Energy be 6MeV high energy acclerator while load while three axis pressure process of dynamic scan in, for diameter 100mm's
The crackle minimum widith that rock sample can identify is 1mm, and in hydraulic fracturing test, effective crack width is largely less than
1mm.The information such as the crack density obtained accordingly are typically only capable to the crack information of reflection large scale, for being lower than device resolution
The crack distribution of size is difficult to obtain.Inventor it has been investigated that, to the unconventional oil and gas such as shale gas exploitation design, geological disaster
Prevention and treatment etc., effectively simulate and observe the pressure break micro-crack of rock to be formed and developed be necessary.
Summary of the invention
In order to solve the above problem in the prior art, in order to solve the problems, such as that rock core internal tiny crack is unrecognized,
First aspect of the present invention it is proposed a kind of rock core internal tiny crack recognition methods based on digital picture, comprising:
Step S10 obtains the first rock sample, obtains initial slice sequence by CT scan, and obtain the initial slice sequence
Each initial slice with reference to figure point in column；
Step S20 carries out crushing test to the first rock sample based on three axis pressurized equipments, obtains different loads by CT scan
Each process Slice Sequence corresponding with the initial slice sequence under intensity；
Step S30 is based on the matched method of local area image, obtains pair with reference to figure point in each process slice
Ying Dian obtains the first displacement field of each process slice；
Step S40, based on the first displacement field of each process slice, matching figure point after obtaining sub-pixel positioning obtains the
Two displacement fields, and extract the micro-crack after rock deformation destroys.
In some preferred embodiments, step S30 " is based on the matched method of local area image, obtains the ginseng
Figure point is examined in the corresponding points in each process slice " in, with reference to the acquisition methods of corresponding points of the figure point P in process slice (N, M)
Are as follows:
To include the sub-district of the point based on point P building；The size of the sub-district is pre-set dimension；
With the size building sliding window of the sub-district in process slice (N, M) on piece traversal, one group of image to be matched is obtained, point
The related coefficient of each image to be matched Yu the sub-district is not calculated；
Based on related coefficient calculated, acquisition and the highest image to be matched of sub-district matching degree are schemed as matching
Picture；
Based on point P in the positional relationship of the sub-district, the point D of same location relationship in matching image is calculated, is existed as point P
Process is sliced the Corresponding matching figure point in (N, M)；
Wherein, N is the serial number of rate of loading, and M is the serial number of process slice in process Slice Sequence.
In some preferred embodiments, positional relationship of the point P in the sub-district are as follows: point P is the center of the sub-district
Point.
In some preferred embodiments, the related coefficient of each image to be matched and the sub-district, calculating side
Method are as follows:
Wherein, C is related coefficient, and f (x, y) is the grey value profile g (x with reference to the corresponding sub-district of figure point*,y*) cut for process
The grey value profile of image to be matched in piece.
In some preferred embodiments, the displacement field of the process slice, acquisition methods are as follows:
All Corresponding matching figure points with reference to figure point in process slice (N, M) are obtained, and calculate each position with reference to figure point
Shifting amount obtains the displacement field of process slice (N, M).
In some preferred embodiments, right during step S40 " extracting the micro-crack after rock deformation destroys "
The matching figure point of process slice carries out sub-pixel positioning, method are as follows:
Wherein,For the matching figure point after sub-pixel positioning；(x0,y0) it is with reference to figure point；(u, ν) indicates (x0,y0)
The displacement of point；Δ x and Δ y are the corresponding Displacement of whole pixel displacement respectively；gxFor the First-order Gradient of X-direction；gyFor the side y
To First-order Gradient；F is the pixel of original image；G is the pixel of image after pressure break.
In some preferred embodiments, " different rate of loading ", acquisition methods in step S20 are as follows:
Crushing test is carried out to multiple second rock samples respectively, obtains the peak strength of each second rock sample respectively；
Mean value computation is carried out to the peak strength of each second rock sample, obtains peak strength mean value；
Based on preset rate of loading ratio, the corresponding rate of loading of each ratio is obtained for triggering CT scan；
Wherein, second rock sample and first rock sample are derived from same rock, and the stratification drilled through is identical.
In some preferred embodiments, this method further includes step S50:
Based on each process slice the first displacement field or second displacement field, obtain each slice of rock horizontal displacement field and
Vertical displacement field judges the type of crackle.
The second aspect of the present invention proposes a kind of rock core internal tiny crack identifying system based on digital picture, comprising:
Initial slice obtains module, process slice obtains module, displacement field obtains module, micro-crack identification module；
The initial slice obtains module, is configured to obtain the first rock sample, obtains initial slice sequence by CT scan, and
Obtain each initial slice in the initial slice sequence with reference to figure point；
The process slice obtains module, is configured to three axis pressurized equipments and carries out crushing test to the first rock sample, leads to
It crosses CT scan and obtains each process Slice Sequence corresponding with the initial slice sequence under different rate of loading；
The displacement field obtains module, is configured to the matched method of local area image, obtains described with reference to figure point
Corresponding points in each process slice obtain the first displacement field of each process slice；
The micro-crack identification module is configured to the first displacement field of each process slice, after obtaining sub-pixel positioning
Matching figure point, obtain second displacement field, and extract rock deformation destroy after micro-crack.
The third aspect of the present invention proposes a kind of processing unit, including processor, storage device；Processor, suitable for holding
Each program of row；Storage device is suitable for storing a plurality of program；Described program is suitable for being loaded by processor and being executed above-mentioned to realize
The rock core internal tiny crack recognition methods based on digital picture.
The present invention has effectively improved the accuracy of identification of crackle, is conducive to the extension and evolution of simulation rock crackle.
Detailed description of the invention
By reading a detailed description of non-restrictive embodiments in the light of the attached drawings below, the application's is other
Feature, objects and advantages will become more apparent upon:
Fig. 1 is the process signal of the rock core internal tiny crack recognition methods based on digital picture of an embodiment of the present invention
Figure；
Fig. 2 is the local displacement field amplification exemplary diagram that process is sliced in an embodiment of the present invention.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing to the embodiment of the present invention
In technical solution be clearly and completely described, it is clear that described embodiments are some of the embodiments of the present invention, without
It is whole embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not before making creative work
Every other embodiment obtained is put, shall fall within the protection scope of the present invention.
The application is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining related invention, rather than the restriction to the invention.It also should be noted that in order to
Convenient for description, part relevant to related invention is illustrated only in attached drawing.
It should be noted that in the absence of conflict, the features in the embodiments and the embodiments of the present application can phase
Mutually combination.
A kind of rock core internal tiny crack recognition methods based on digital picture of the invention is based on digital picture correlation skill
Art obtains the distribution of rock interior micro-crack, displacement field, the method for determining crackle property.To improve the accuracy of identification of crackle, have
Conducive to the extension and evolution of simulation rock crackle.As shown in Figure 1, the method for the present invention includes:
Step S10 obtains the first rock sample, obtains initial slice sequence by CT scan, and obtain the initial slice sequence
Each initial slice with reference to figure point in column；
Step S20 carries out crushing test to the first rock sample based on three axis pressurized equipments, obtains different loads by CT scan
Each process Slice Sequence corresponding with the initial slice sequence under intensity；
Step S30 is based on the matched method of local area image, obtains pair with reference to figure point in each process slice
Ying Dian obtains the first displacement field of each process slice；
Step S40, based on the first displacement field of each process slice, matching figure point after obtaining sub-pixel positioning obtains the
Two displacement fields, and extract the micro-crack after rock deformation destroys.
Can also include in some embodiments further includes step S50: the first displacement field or the based on each process slice
Two displacement fields obtain horizontal displacement field and the vertical displacement field of each slice of rock, judge the type of crackle.
In order to be more clearly illustrated to the rock core internal tiny crack recognition methods the present invention is based on digital picture, below
Expansion detailed description is carried out to each step in a kind of embodiment of our inventive method in conjunction with attached drawing.
The embodiment of the present invention is based on same rock, and the identical stratification drilled through obtains multiple sizes one of rock production
The rock sample of cause, and selection is one of as the first rock sample progress sweep test in load, remaining is multiple as the second rock sample
Carry out peak strength test.Rock sample is the cylindrical body of diameter 100mm, high 200mm in the present embodiment.And used device is
The high energy acclerator CT multi- scenarios method rock mechanics testing system of Institute of Geophysics, Academia Sinica.Herein to rock sample ruler
Purpose that is very little, being described using testing equipment is only to keep technical solution of the present invention apparent, be should not be understood as to the technology of the present invention side
The restriction of case.
Step S10 obtains the first rock sample, obtains initial slice sequence by CT scan, and obtain the initial slice sequence
Each initial slice with reference to figure point in column.
In the present embodiment, the first rock sample can be put into three axis pressurized equipments of high energy acclerator CT turntable, opened
Initial slice sequence is obtained by CT scan before beginning crushing test, and obtains each initial slice in the initial slice sequence
With reference to figure point.
Step S20 carries out crushing test to the first rock sample based on three axis pressurized equipments, obtains different loads by CT scan
Each process Slice Sequence corresponding with the initial slice sequence under intensity.
To the first rock sample in the three axis pressurized equipments for being put into high energy acclerator CT turntable, carries out the dynamic in load and sweep
Test is retouched, CT image reconstruction is carried out after scanning, scan slice of the rock core in different distortion failure stage is obtained, identifies not same order
The underbead crack of section.
Used different rate of loading, acquisition methods are as follows: crushing test is carried out to multiple second rock samples respectively, respectively
Obtain the peak strength of each second rock sample；Mean value computation is carried out to the peak strength of each second rock sample, obtains peak strength mean value；
Based on preset rate of loading ratio, the corresponding rate of loading of each ratio is obtained for triggering CT scan.
In one example, peak strength can be obtained by the test to multiple second rock samples, and calculates its average value
σca, rate of loading is obtained according to 20%, 40%, 60%, 80%, then 20% σca, 40% σca, 60% σca, 80% σcaIt is corresponding
Rate of loading σciIt (i=1,2,3,4) is the acquired rate of loading in crushing test progress CT scan.
Then, the first rock sample is subjected to the dynamic scan in load, reaches the load of above-mentioned calculating in rate of loading respectively
Intensity σciWhen carry out CT dynamic scan, until sample destroy.After the load of first rock sample destroys, it is strong to obtain corresponding peak value
Spend σc, can recalculate to the corresponding true stress ratio of each rate of loading of the rock sample, i.e. σci/σc。
In the present embodiment, obtained under each load scanning intensity rock sample process slice (can be with linear array scanning, can also be with
Planar array scanning, the linear array scanning of example of the present invention), slice thickness 2mm is corresponding to obtain 100 CT slices.Sample after destruction
Underbead crack using VG Studio software extract crack, can recognize measurement crackle minimum widith be 1mm.Process is cut
Chip resolution is 1024*1024.
Step S30 is based on the matched method of local area image, obtains pair with reference to figure point in each process slice
Ying Dian obtains the first displacement field of each process slice.
Rock has situations such as deformation displacement in pressure process, therefore the slice obtained needs to carry out matching primitives: false
If axis direction maximum distortion is a%, height of specimen d, then m layers of slice computer capacity of this programme are [m-d*a%n/um, m+
D*a%/num], num is total number of plies of slice.Similarity measures are carried out with deformed slice to before deformation respectively, obtain water
First displacement field of square tangential section, the micro-crack after extracting rock deformation and failure.For being vertically sliced, radial deformation and sample
Diameter similarly calculates, and can obtain vertical cross section crackle.
In the present embodiment, be based on the matched method of local area image, obtain it is described with reference to figure point each process be sliced in
Corresponding points method, essence be exactly with reference to figure point P process be sliced (N, M) in corresponding points acquisition methods (N be load
The serial number of intensity, M are the serial number of process slice in process Slice Sequence), this method are as follows:
Step S31, to include the sub-district of the point based on point P building；The size of the sub-district is pre-set dimension.
In the present embodiment, point P is the central point of the sub-district, that is, the point centered on point P, according to pre-set dimension structure
Build sub-district.Pre-set dimension can be the pixel of 50*50.More information can be provided in this way in sub-district, so that unique match
It comes true.
Step S32 obtains one group and waits for the size of sub-district building sliding window in process slice (N, M) on piece traversal
With image, the related coefficient of each image to be matched Yu the sub-district is calculated separately.
Sub-district traverses on process slice, constantly calculates the related coefficient of the sub-district, the related coefficient meter that this patent uses
Shown in calculation method such as formula (1):
Wherein, C is related coefficient, and f (x, y) is the grey value profile g (x with reference to the corresponding sub-district of figure point*,y*) cut for process
The grey value profile of image to be matched in piece.
Step S33 is based on related coefficient calculated, obtains and the highest image to be matched of sub-district matching degree, work
For matching image.
Above-mentioned related coefficient indicates that two modules are perfectly correlated when C=1, when C=0, two modules
It is uncorrelated.The maximum image to be matched of related coefficient can be chosen in the present embodiment as matching image.
Step S34 calculates the point D of same location relationship in matching image based on point P in the positional relationship of the sub-district,
As Corresponding matching figure point of the point P in process slice (N, M).
Step S31-S34 is repeated, after completing traversal with reference to figure point, all matching figure points are obtained, based on all with reference to figure
Corresponding matching figure point of the point in process slice (N, M), and each displacement with reference to figure point is calculated, obtain process slice (N, M)
The first displacement field (example for being illustrated in figure 2 a local displacement field), so as to obtain the expansion of rock interior microcrack
Exhibition.
Since the first displacement field is the vector representation being displaced in process slice with reference to figure point, by each in the first displacement field
Displacement vector, can be apparent show displacement vector it is inconsistent refer to figure point, using the region of these figure points composition as
The crack to be extracted.
Step S40, based on the first displacement field of each process slice, matching figure point after obtaining sub-pixel positioning obtains the
Two displacement fields, and extract the micro-crack after rock deformation destroys.
The minimum unit of digital picture is pixel, in order to improve the accuracy of identification of crackle, using sub-pix matching algorithm, is calculated
The accuracy of identification of method and the accuracy of identification of crackle are positively correlated.It is right during " extracting the micro-crack after rock deformation destroys "
The matching figure point of process slice carries out sub-pixel positioning, shown in method such as formula (2), (3):
Wherein,For the matching figure point after sub-pixel positioning；(x0,y0) it is with reference to figure point；(u, ν) indicates (x0,y0)
The displacement of point；Δ x and Δ y are the corresponding Displacement of whole pixel displacement respectively.
This paper patent calculates the gray scale and gray scale derivative of sub-pixel location using bicubic spline interpolation, and here is sub- picture
Shown in the calculation method such as formula (4) of element displacement.
Wherein, gxFor the First-order Gradient of X-direction；gyFor the First-order Gradient in the direction y；F is the pixel of original image；G is pressure break
The pixel of image afterwards.
It is similar with the extracting method of crackle in the first displacement field, it, can be very bright by each displacement vector in second displacement field
It is aobvious show displacement vector it is inconsistent with reference to figure point, the region that these figure points are constituted is as the micro-crack to be extracted.
The present invention only gives the calculating of a slice, and to analyze rock seam net form state can repeat, and obtains multiple
Crack/micro-crack of scan period multiple slices is to analyze rock seam net form state.
Step S50, the first displacement field or second displacement field based on each process slice, obtains the level of each slice of rock
Displacement field and vertical displacement field, judge the type of crackle.
The micro-crack that the crack obtained based on step S30/step S40 is obtained, to crack/micro-crack in corresponding displacement field
The displacement vector in region is judged:
It is tensile crack/micro-crack if having relative displacement with reference to figure point perpendicular to crack/micro-crack direction；
It is shearing crack/micro-crack if being parallel to crack/micro-crack direction to have relative displacement with reference to figure point；
If above-mentioned two situations all include, to draw shear crack/micro-crack.
The rock core internal tiny crack identifying system based on digital picture of an embodiment of the present invention, including initial slice obtain
Modulus block, process slice obtain module, displacement field obtains module, micro-crack identification module；
The initial slice obtains module, is configured to obtain the first rock sample, obtains initial slice sequence by CT scan, and
Obtain each initial slice in the initial slice sequence with reference to figure point；
The process slice obtains module, is configured to three axis pressurized equipments and carries out crushing test to the first rock sample, leads to
It crosses CT scan and obtains each process Slice Sequence corresponding with the initial slice sequence under different rate of loading；
The displacement field obtains module, is configured to the matched method of local area image, obtains described with reference to figure point
Corresponding points in each process slice obtain the first displacement field of each process slice；
The micro-crack identification module is configured to the first displacement field of each process slice, after obtaining sub-pixel positioning
Matching figure point, obtain second displacement field, and extract rock deformation destroy after micro-crack.
Person of ordinary skill in the field can be understood that, for convenience and simplicity of description, foregoing description
The specific work process of system and related explanation, can refer to corresponding processes in the foregoing method embodiment, details are not described herein.
It should be noted that the rock core internal tiny crack identifying system provided by the above embodiment based on digital picture, only
The example of the division of the above functional modules, in practical applications, it can according to need and by above-mentioned function distribution
Completed by different functional modules, i.e., by the embodiment of the present invention module or step decompose or combine again, for example, on
The module for stating embodiment can be merged into a module, multiple submodule can also be further split into, to complete above description
All or part of function.For module involved in the embodiment of the present invention, the title of step, it is only for distinguish each
Module or step, are not intended as inappropriate limitation of the present invention.
The storage device of an embodiment of the present invention, wherein being stored with a plurality of program, described program is suitable for being added by processor
It carries and executes to realize the above-mentioned rock core internal tiny crack recognition methods based on digital picture.
The processing unit of an embodiment of the present invention, including processor, storage device；Processor is adapted for carrying out each journey
Sequence；Storage device is suitable for storing a plurality of program；Described program is suitable for being loaded by processor and being executed above-mentioned based on number to realize
The rock core internal tiny crack recognition methods of word image.
Person of ordinary skill in the field can be understood that, for convenience and simplicity of description, foregoing description
The specific work process and related explanation of storage device, processing unit, can refer to corresponding processes in the foregoing method embodiment,
Details are not described herein.
Those skilled in the art should be able to recognize that, mould described in conjunction with the examples disclosed in the embodiments of the present disclosure
Block, method and step, can be realized with electronic hardware, computer software, or a combination of the two, software module, method and step pair
The program answered can be placed in random access memory (RAM), memory, read-only memory (ROM), electrically programmable ROM, electric erasable and can compile
Any other form of storage well known in journey ROM, register, hard disk, moveable magnetic disc, CD-ROM or technical field is situated between
In matter.In order to clearly demonstrate the interchangeability of electronic hardware and software, in the above description according to function generally
Describe each exemplary composition and step.These functions are executed actually with electronic hardware or software mode, depend on technology
The specific application and design constraint of scheme.Those skilled in the art can carry out using distinct methods each specific application
Realize described function, but such implementation should not be considered as beyond the scope of the present invention.
Term " first ", " second " etc. are to be used to distinguish similar objects, rather than be used to describe or indicate specific suitable
Sequence or precedence.
Term " includes " or any other like term are intended to cover non-exclusive inclusion, so that including a system
Process, method, article or equipment/device of column element not only includes those elements, but also including being not explicitly listed
Other elements, or further include the intrinsic element of these process, method, article or equipment/devices.
So far, it has been combined preferred embodiment shown in the drawings and describes technical solution of the present invention, still, this field
Technical staff is it is easily understood that protection scope of the present invention is expressly not limited to these specific embodiments.Without departing from this
Under the premise of the principle of invention, those skilled in the art can make equivalent change or replacement to the relevant technologies feature, these
Technical solution after change or replacement will fall within the scope of protection of the present invention.
Claims (10)
1. a kind of rock core internal tiny crack recognition methods based on digital picture characterized by comprising
Step S10 obtains the first rock sample, obtains initial slice sequence by CT scan, and obtain in the initial slice sequence
Each initial slice with reference to figure point；
Step S20 carries out crushing test to the first rock sample based on three axis pressurized equipments, obtains different rate of loading by CT scan
Lower each process Slice Sequence corresponding with the initial slice sequence；
Step S30 is based on the matched method of local area image, obtains the correspondence with reference to figure point in each process slice
Point obtains the first displacement field of each process slice；
Step S40, based on each process slice the first displacement field, obtain sub-pixel positioning after matching figure point, obtain second
Field is moved, and extracts the micro-crack after rock deformation destroys.
2. the rock core internal tiny crack recognition methods according to claim 1 based on digital picture, which is characterized in that step
In S30 " being based on the matched method of local area image, obtain the corresponding points with reference to figure point in each process is sliced ", reference
Scheme the acquisition methods of corresponding points of the point P in process slice (N, M) are as follows:
To include the sub-district of the point based on point P building；The size of the sub-district is pre-set dimension；
With the size building sliding window of the sub-district in process slice (N, M) on piece traversal, one group of image to be matched is obtained, is counted respectively
Calculate the related coefficient of each image to be matched Yu the sub-district；
Based on related coefficient calculated, acquisition and the highest image to be matched of sub-district matching degree, as matching image；
Based on point P in the positional relationship of the sub-district, the point D of same location relationship in matching image is calculated, as point P in process
The Corresponding matching figure point being sliced in (N, M)；
Wherein, N is the serial number of rate of loading, and M is the serial number of process slice in process Slice Sequence.
3. the rock core internal tiny crack recognition methods according to claim 2 based on digital picture, which is characterized in that point P
In the positional relationship of the sub-district are as follows: point P is the central point of the sub-district.
4. the rock core internal tiny crack recognition methods according to claim 2 based on digital picture, which is characterized in that described
The related coefficient of each image to be matched and the sub-district, calculation method are as follows:
Wherein, C is related coefficient, and f (x, y) is the grey value profile g (x with reference to the corresponding sub-district of figure point*,y*) it is in process slice
The grey value profile of image to be matched.
5. the rock core internal tiny crack recognition methods according to claim 2 based on digital picture, which is characterized in that described
The displacement field of process slice, acquisition methods are as follows:
All Corresponding matching figure points with reference to figure point in process slice (N, M) are obtained, and calculate each displacement with reference to figure point,
Acquisition process is sliced the displacement field of (N, M).
6. the rock core internal tiny crack recognition methods according to claim 2 based on digital picture, which is characterized in that step
During S40 " extracting the micro-crack after rock deformation destroys ", sub-pixel positioning is carried out to the matching figure point of process slice,
Method are as follows:
Wherein,For the matching figure point after sub-pixel positioning；(x0,y0) it is with reference to figure point；(u, ν) indicates (x0,y0) point
Displacement；Δ x and Δ y are the corresponding Displacement of whole pixel displacement respectively；gxFor the First-order Gradient of X-direction；gyFor the direction y
First-order Gradient；F is the pixel of original image；G is the pixel of image after pressure break.
7. the rock core internal tiny crack recognition methods according to claim 1-6 based on digital picture, feature
It is, " different rate of loading ", acquisition methods in step S20 are as follows:
Crushing test is carried out to multiple second rock samples respectively, obtains the peak strength of each second rock sample respectively；
Mean value computation is carried out to the peak strength of each second rock sample, obtains peak strength mean value；
Based on preset rate of loading ratio, the corresponding rate of loading of each ratio is obtained for triggering CT scan；
Wherein, second rock sample and first rock sample are derived from same rock, and the stratification drilled through is identical.
8. the rock core internal tiny crack recognition methods according to claim 1-6 based on digital picture, feature
It is, this method further includes step S50:
The first displacement field or second displacement field based on each process slice obtain the horizontal displacement field of each slice of rock and vertical
Displacement field judges the type of crackle.
9. a kind of rock core internal tiny crack identifying system based on digital picture characterized by comprising initial slice obtains mould
Block, process slice obtain module, displacement field obtains module, micro-crack identification module；
The initial slice obtains module, is configured to obtain the first rock sample, obtains initial slice sequence by CT scan, and obtain
Each initial slice with reference to figure point in the initial slice sequence；
The process slice obtains module, is configured to three axis pressurized equipments and carries out crushing test to the first rock sample, passes through CT
Scanning obtains each process Slice Sequence corresponding with the initial slice sequence under different rate of loading；
The displacement field obtains module, is configured to the matched method of local area image, obtain it is described with reference to figure point each
Corresponding points in process slice obtain the first displacement field of each process slice；
The micro-crack identification module is configured to the first displacement field of each process slice, after obtaining sub-pixel positioning
Figure point obtains second displacement field, and extracts the micro-crack after rock deformation destroys.
10. a kind of processing unit, including processor, storage device；Processor is adapted for carrying out each program；Storage device is suitable for
Store a plurality of program；It is characterized in that, described program is suitable for being loaded by processor and being executed to realize any one of claim 1-8
The rock core internal tiny crack recognition methods based on digital picture.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910290206.5A CN110057846A (en) | 2019-04-11 | 2019-04-11 | Rock core internal tiny crack recognition methods based on digital picture, system, device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910290206.5A CN110057846A (en) | 2019-04-11 | 2019-04-11 | Rock core internal tiny crack recognition methods based on digital picture, system, device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110057846A true CN110057846A (en) | 2019-07-26 |
Family
ID=67318782
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910290206.5A Pending CN110057846A (en) | 2019-04-11 | 2019-04-11 | Rock core internal tiny crack recognition methods based on digital picture, system, device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110057846A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110686971A (en) * | 2019-10-14 | 2020-01-14 | 中国科学院武汉岩土力学研究所 | Method for accurately judging cracking stress and cracking type of fractured rock mass |
| CN111754500A (en) * | 2020-06-30 | 2020-10-09 | 中国科学院地质与地球物理研究所 | Rock fracturing network topological structure depicting system |
| CN111754471A (en) * | 2020-06-12 | 2020-10-09 | 中国科学院地质与地球物理研究所 | Rock fracture detection method under in-situ real-time pressurization of soil-rock mixture |
| CN111767631A (en) * | 2019-08-02 | 2020-10-13 | 中国石油大学(华东) | Method and system for simulating rock crack propagation based on multiphase digital rock core |
| CN113538473A (en) * | 2021-07-08 | 2021-10-22 | 南京航空航天大学 | Random grid and special-shaped subarea division method for cracks in digital image correlation |
| CN114266776A (en) * | 2022-03-03 | 2022-04-01 | 南京航空航天大学 | Digital image correlation method applying composite crack displacement field function |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104330023A (en) * | 2014-10-15 | 2015-02-04 | 浙江大学 | Acquisition system and identification method of concrete surface initial crack information |
-
2019
- 2019-04-11 CN CN201910290206.5A patent/CN110057846A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104330023A (en) * | 2014-10-15 | 2015-02-04 | 浙江大学 | Acquisition system and identification method of concrete surface initial crack information |
Non-Patent Citations (5)
| Title |
|---|
| 何向 等: "基于CT 扫描技术和数字图像相关法的岩石变形场分析", 《北京力学会第20届学士年会论文集》 * |
| 潘兵 等: "数字图像相关中的亚像素位移定位算法进展", 《力学进展》 * |
| 王开福 等: "《散斑计量》", 31 December 2012, 北京理工大学出版社 * |
| 袁则循: "基于CT图像的混凝土损伤演化及数值模拟研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 韩依颖: "基于DIC方法的混凝土压缩试件裂纹识别", 《土木建筑与环境工程》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111767631A (en) * | 2019-08-02 | 2020-10-13 | 中国石油大学(华东) | Method and system for simulating rock crack propagation based on multiphase digital rock core |
| CN111767631B (en) * | 2019-08-02 | 2022-10-18 | 中国石油大学(华东) | Method and system for simulating rock crack propagation based on multiphase digital rock core |
| CN110686971A (en) * | 2019-10-14 | 2020-01-14 | 中国科学院武汉岩土力学研究所 | Method for accurately judging cracking stress and cracking type of fractured rock mass |
| CN111754471A (en) * | 2020-06-12 | 2020-10-09 | 中国科学院地质与地球物理研究所 | Rock fracture detection method under in-situ real-time pressurization of soil-rock mixture |
| CN111754500A (en) * | 2020-06-30 | 2020-10-09 | 中国科学院地质与地球物理研究所 | Rock fracturing network topological structure depicting system |
| CN111754500B (en) * | 2020-06-30 | 2023-06-27 | 中国科学院地质与地球物理研究所 | Rock fracture network topology structure describing system |
| CN113538473A (en) * | 2021-07-08 | 2021-10-22 | 南京航空航天大学 | Random grid and special-shaped subarea division method for cracks in digital image correlation |
| CN114266776A (en) * | 2022-03-03 | 2022-04-01 | 南京航空航天大学 | Digital image correlation method applying composite crack displacement field function |
| CN114266776B (en) * | 2022-03-03 | 2022-06-10 | 南京航空航天大学 | Digital image correlation method applying composite crack displacement field function |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110057846A (en) | Rock core internal tiny crack recognition methods based on digital picture, system, device | |
| CN110415167B (en) | Rough surface crack generation method and test system based on digital image technology | |
| Benn | Clast morphology | |
| Loeffler et al. | Detection and segmentation of mechanical damage in concrete with X-Ray microtomography | |
| Panahi et al. | A 4D synchrotron X-ray-tomography study of the formation of hydrocarbon-migration pathways in heated organic-rich shale | |
| CN106023158B (en) | The fresh water pipless pearl pearly layer defect identification method of SD-OCT images | |
| US20030174877A1 (en) | Method for inspecting patterns | |
| JP2016508295A (en) | Defect detection on wafers using template image matching | |
| KR101787161B1 (en) | Design-based inspection using repeating structures | |
| KR20170124984A (en) | Method for processing data of ground penetrating radar | |
| Adhikari et al. | Tracking of defects in reinforced concrete bridges using digital images | |
| Rivera et al. | Automated detection and measurement of cracks in reinforced concrete components | |
| Chen et al. | Digital image based approach for three-dimensional mechanical analysis of heterogeneous rocks | |
| CN111507972A (en) | Tunnel surface defect detection method combining convolutional neural network and support vector machine | |
| Lindquist | Quantitative analysis of three-dimensional X-ray tomographic images | |
| CN111624146A (en) | Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir | |
| Peinke et al. | Experimental study of cone penetration in snow using x-ray tomography | |
| CN114821079A (en) | Rock crack information identification method and system based on variational self-encoder | |
| CN111462056A (en) | Workpiece surface defect detection method, device, equipment and storage medium | |
| Liu et al. | Automatic fracture segmentation and detection from image logging using mask R-CNN | |
| Elbischger et al. | Automatic analysis of collagen fiber orientation in the outermost layer of human arteries | |
| CN115147340A (en) | Rock crack deformation information acquisition method, device, equipment and storage medium | |
| Wang | An edge based segmentation algorithm for rock fracture tracing | |
| CN115100402A (en) | Small sample segmentation model, small target detection method, system, medium, and apparatus | |
| Diogenes et al. | Porous media microstructure reconstruction using pixel-based and object-based simulated annealing: comparison with other reconstruction methods |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190726 |