CN106228558B - A method of layer being carried out into formula breakdown diagnosis to rock - Google Patents

A method of layer being carried out into formula breakdown diagnosis to rock Download PDF

Info

Publication number
CN106228558B
CN106228558B CN201610597318.1A CN201610597318A CN106228558B CN 106228558 B CN106228558 B CN 106228558B CN 201610597318 A CN201610597318 A CN 201610597318A CN 106228558 B CN106228558 B CN 106228558B
Authority
CN
China
Prior art keywords
rock
group
sample
ring
carried out
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.)
Active
Application number
CN201610597318.1A
Other languages
Chinese (zh)
Other versions
CN106228558A (en
Inventor
刘杰
雷岚
颜媛媛
梁佳丽
兰钰
支永艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ANHUI PROVINCE WANJIANG MECHANICAL AND ELECTRICAL EQUIPMENT MANUFACTURING Co.,Ltd.
Original Assignee
China Three Gorges University CTGU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN201610597318.1A priority Critical patent/CN106228558B/en
Publication of CN106228558A publication Critical patent/CN106228558A/en
Application granted granted Critical
Publication of CN106228558B publication Critical patent/CN106228558B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0004—Industrial image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30108—Industrial image inspection
    • G06T2207/30132—Masonry; Concrete

Abstract

A method of layer being carried out into formula breakdown diagnosis to rock, this approach includes the following steps：Step 1：Region division is carried out to the cross section；Step 2：Calculate the CT average values CT of various kinds ring under original stateTave(N)；Step 3：The freeze-thaw cycle experiment that the different numbers of n groups are carried out to the rock sample, calculates separately the CT average values CT of various kinds ring after each group freeze-thaw cycle experimentTave (N) n-th group, step 4：The data that step 2 and step 3 obtain are arranged, breakdown diagnosis is carried out to rock, judges the concentration position that freeze-thaw cycle experiment is damaged caused by rock sample.It is provided by the invention it is a kind of to rock carry out layer into formula breakdown diagnosis method, it can solve existing method and be unable to get the CT average values for covering all the points, as a result inaccurate problem, the extent of damage of accurate reflection rock, the extent of damage in infinitely small damage range can be obtained, it is easy to operate, do not have to frequent operation.

Description

A method of layer being carried out into formula breakdown diagnosis to rock
Technical field
It is especially a kind of that layer being carried out into formula damage to rock the present invention relates to a kind of method carrying out breakdown diagnosis to rock sample The method of analysis.
Background technology
With going deep into for scientific research, the method for carrying out breakdown diagnosis to rock sample now is to carry out CT to impaired rock sample mostly Scanning analyzes the extent of damage of rock interior by the size of the bright-dark degree of the imaging after CT scan and CT values.And existing rank Section by CT scan analyze rock method it is main there are three：
1, comparison of light and shade with the naked eye is carried out to CT images；
The CT scan figure of a certain level of rock, figure is brighter, indicates that rock compactness is better.A small number of stains at middle part indicate There are localized voids, although can intuitively observe that the hole of rock interior compares using the method, substantially obtains the damage of rock Hinder rule, the analysis of a quantization can not be but carried out to it, can not accurately determine damage range.
2, the foundation that limited circle takes the CT average values in face to analyze as a whole is chosen in level
Whether there is or not several circumference from the center of circle to periphery, it was therefore concluded that accuracy can be taken with the circle of selection face quantity difference And it is different, but all the points can not be covered always, so being unable to get an accurate conclusion.
3, criterion of the CT average values of rounding scanning level as scanning level compactness.But this method can not be anti- The damage comparison inside rock sample is reflected, using the CT values as basis for estimation, can not reflect the inhomogeneities of scanning level.
Invention content
Layer is carried out into the method for formula breakdown diagnosis to rock technical problem to be solved by the invention is to provide a kind of, it can be with It solves existing method and is unable to get the CT average values for covering all the points, as a result inaccurate problem is accurate to reflect the impaired of rock The extent of damage in infinitely small damage range can be obtained in degree, easy to operate, does not have to frequent operation.
In order to solve the above technical problems, the technical solution adopted in the present invention is：It is a kind of that layer being carried out into formula damage to rock The method of analysis, this approach includes the following steps：
Step 1：The a certain cross section on rock sample is chosen, region division is carried out to the cross section：From the cross section Central point sets out, and marks off the sample ring of N number of equal semidiameters；
Step 2：Calculate the CT average values CT of various kinds ring under original stateTave(N)：Respectively to the circle where each sample ring Initial CT scan is carried out, the CT values for obtaining the circle where each sample ring are CT1 is initial、CT2 is initial、CT3 is initial…CTN is initial, measure each The area value of circle where sample ring is S1、S2、S3…SN,
The CT average values CT of various kinds ring under original stateTave (N) is initial=(CTN-CTN-1)/(SN-SN-1)；
Step 3：It is real to calculate separately each group Frozen-thawed cycled for the freeze-thaw cycle experiment that the different numbers of n groups are carried out to the rock sample Test the CT average values CT of rear various kinds ringTave (N) n-th group, calculation formula is
CTThe 1st group of Tave (N)=(CTThe 1st group of N-CT(N-1) the 1st group)/(SN-SN-1),
CTThe 2nd group of Tave (N)=(CTThe 1st group of N-CT(N-1) the 2nd group)/(SN-SN-1),
…
CTTave (N) n-th group=(CTN n-th groups-CT(N-1) n-th group)/(SN-SN-1),
In formula, CTN n-th groupsFor the CT values of the circle where n-th sample ring after n-th group freeze-thaw cycle experiment；
Step 4：The data that step 2 and step 3 obtain are arranged, breakdown diagnosis is carried out to rock, judge Frozen-thawed cycled reality Test the concentration position damaged caused by rock sample.
In step 1, rock sample is the cylindrical rock sample that radius is 2.5 centimetres, from the central point of the cross section, is divided Go out the sample ring of 5 equal semidiameters, semidiameter is 0.5 centimetre.
The method further includes step 5：
Step 5：Region is carried out again the concentration position that the freeze-thaw cycle experiment that step 4 obtains is damaged caused by rock sample It divides, from the central point of the cross section, the concentration position of damage is marked off into the M sample ring for waiting semidiameters；According still further to Above-mentioned steps 2-4 carries out breakdown diagnosis to rock, judges the accurately concentration that freeze-thaw cycle experiment is damaged caused by rock sample Position.
It is provided by the invention it is a kind of to rock carry out layer into the method for formula breakdown diagnosis, have the beneficial effect that：
1, accurately reflect the extent of damage of rock.All the points can not be covered when previous method CT scan, only with limited The CT average values reflection of point is so the CT average values put, can not accurately reflect the extent of damage of rock.And this method passes through division The mode of sample ring including covering all the points, is no longer limited to take the mean CT-number in face to damage rock stratum with limited circle Wound analysis, then utilizes mathematical formulae, obtains the mean CT-number of various kinds ring, (is declined by the changing value of CT values in comparative sample ring Value is bigger, and the extent of damage is bigger) it can accurately determine the extent of damage of rock sample various kinds loop section.
2, ring can be divided again to sample ring, the extent of damage in infinitely small damage range can be obtained.Primarily determine various kinds ring The extent of damage after, ring can be divided again to arbitrary sample ring, obtain several thin rings, by similarly calculating, obtained in each thin ring CT average values, then by comparing the CT average values in each thin ring, obtain each thin ring extent of damage, further reduce rock by Damage range.
3, easy to operate, it is not required to frequent operation.Commonsense method obtains a diameter of 0.5 centimetre of cylinder rock sample damage range The minimum time needed for CT values is 10 minutes, and is only needed using this method 2 minutes, and the operating time is greatly saved.
Description of the drawings
The invention will be further described with reference to the accompanying drawings and examples：
Fig. 1 is 5 sample ring schematic diagrames for waiting semidiameters that step 1 marks off in the embodiment of the present invention one；
Fig. 2 is to be followed according to every sample ring mean CT-number-freeze thawing that the data that step 4 sorts out are made in the embodiment of the present invention one Ring experiment number curve graph；
Fig. 3 is to be carried out again the concentration position that freeze-thaw cycle experiment is damaged caused by rock sample in the embodiment of the present invention two The sample ring schematic diagram of region division；
Fig. 4 is to be followed according to every sample ring mean CT-number-freeze thawing that the data that step 4 sorts out are made in the embodiment of the present invention two Ring experiment number curve graph；
Fig. 5 is that the sample ring of the n equal semidiameters marked off to the rock sample that cross section is rectangle in the embodiment of the present invention three shows It is intended to.
Specific implementation mode
Embodiment one
Step 1：Rock sample is the cylindrical rock sample that radius is 2.5 centimetres, a certain cross section on rock sample is chosen, to the cross Section carries out region division：From the central point of the cross section, the sample ring of 5 equal semidiameters is marked off, is labeled as first Sample ring, the second sample ring, third sample ring, the 4th sample ring, the 5th sample ring, semidiameter are 0.5 centimetre, as shown in Figure 1；
Step 2：Calculate the CT average values CT of various kinds ring under original stateTave(N)：Respectively to the circle where each sample ring Initial CT scan is carried out, the CT values for obtaining the circle where each sample ring are CT1 is initial、CT2 is initial、CT3 is initial…CTN is initial, measure each The area value of circle where sample ring is S1、S2、S3…SN,
The CT average values CT of various kinds ring under original stateTave (N) is initial=(CTN-CTN-1)/(SN-SN-1)
N=5；
Step 3：It (is respectively 10 freeze-thaw cycle experiments, 30 Frozen-thawed cycled realities to carry out 3 groups of difference numbers to the rock sample Test and 50 freeze-thaw cycle experiments) freeze-thaw cycle experiment, the CT for calculating separately various kinds ring after each group freeze-thaw cycle experiment is average Value CTTave (N) n-th group, calculation formula is
CTThe 1st group of Tave (N)=(CTThe 1st group of N-CT(N-1) the 1st group)/(SN-SN-1),
CTThe 2nd group of Tave (N)=(CTThe 1st group of N-CT(N-1) the 2nd group)/(SN-SN-1),
…
CTTave (N) n-th group=(CTN n-th groups-CT(N-1) n-th group)/(SN-SN-1),
In formula, CTN n-th groupsFor the CT values of the circle where n-th sample ring after n-th group freeze-thaw cycle experiment, N=5, n=3, step The data of 2 and 3 gained are as shown in table 1：
Table 1
Step 4：The data that step 2 and step 3 obtain are arranged, as shown in table 2
Table 2
According to the data of table 2 to every sample ring mean CT-number-freeze-thaw cycle experiment frequency curve figure for making, as shown in Fig. 2,
From figure 2 it can be seen that with the increase of number of freezing and thawing, the CT values amplitude of variation of the 5th sample ring is most apparent, gap It is more maximum than variation, so Frozen-thawed cycled is concentrated mainly on the 5th sample ring, i.e. layer of the radius at 2.0-2.5 centimetres to the damage of rock Face range.
Embodiment two
The method further includes step 5：
Step 5：The 5th sample ring mainly concentrated to the damage of rock to the Frozen-thawed cycled that step 4 obtains carries out region again It divides, from the central point of the cross section, the concentration position of damage is marked off to the sample ring of semidiameters such as M (M=5), Respectively 5-1 samples ring, 5-2 samples ring, 5-3 samples ring, 5-4 samples ring and 5-5 sample rings, as shown in figure 3, semidiameter is 0.1 centimetre；According still further to the step 2-4 of above-described embodiment one, the data sorted out are as shown in table 3：
Table 3
According to the data of table 3 to every sample ring mean CT-number-freeze-thaw cycle experiment frequency curve figure for making, as shown in figure 4, Figure 4, it is seen that the CT values amplitude of variation of 5-5 sample rings is most apparent, void ratio variation is maximum, so Frozen-thawed cycled is to rock The damage of stone is concentrated mainly on 5-5 sample rings, i.e. level range of the radius at 2.4-2.5 centimetres.
This method carries out region division with reference to certain mathematical principle to rock sample, divides the injuring rule one of rule and rock sample It causes, the degree of injury of rock stratum each section is determined using the changing value of each region CT values, show that rock sample is damaged most serious region, then The region is divided again, Query refinement, damage range can be accurate to arbitrary small range.
Embodiment three
It can be also rectangle by sample ring design, the sample ring of semidiameters is waited such as to n that the rock sample that cross section is rectangle marks off Shown in Fig. 5, analyzed according to the step of above-described embodiment one and two.

Claims (2)

1. it is a kind of to rock carry out layer into formula breakdown diagnosis method, it is characterised in that this approach includes the following steps：
Step 1：The a certain cross section on rock sample is chosen, region division is carried out to the cross section：Center from the cross section Point sets out, and marks off N number of sample ring for waiting semidiameters；
Step 2：Calculate the CT average values CT of various kinds ring under original stateTave(N)：The circle where each sample ring is carried out just respectively Beginning CT scan, the CT values for obtaining the circle where each sample ring are CT1 is initial、CT2 is initial、CT3 is initial…CTN is initial, measure each sample ring institute Circle area value be S1、S2、S3…SN,
The CT average values CT of various kinds ring under original stateTave (N) is initial=(CTN is initial-CTN-1 is initial)/(SN-SN-1)；
Step 3：The freeze-thaw cycle experiment that the different numbers of n groups are carried out to the rock sample, after calculating separately each group freeze-thaw cycle experiment The CT average values CT of various kinds ringTave (N) n-th group, calculation formula is
CTThe 1st group of Tave (N)=(CTThe 1st group of N-CT(N-1) the 1st group)/(SN-SN-1),
CTThe 2nd group of Tave (N)=(CTThe 1st group of N-CT(N-1) the 2nd group)/(SN-SN-1),
…
CTTave (N) n-th group=(CTN n-th groups-CT(N-1) n-th group)/(SN-SN-1),
In formula, CTN n-th groupsFor the CT values of the circle where n-th sample ring after n-th group freeze-thaw cycle experiment；
Step 4：The data that step 2 and step 3 obtain are arranged, breakdown diagnosis is carried out to rock, judges freeze-thaw cycle experiment pair The concentration position damaged caused by rock sample；
It further include step 5：
Step 5：Region is carried out again the concentration position that the freeze-thaw cycle experiment that step 4 obtains is damaged caused by rock sample to draw Divide, from the central point of the cross section, the concentration position of damage is marked off to the sample ring of M equal semidiameters；According still further to upper Step 2 is stated to 4, breakdown diagnosis is carried out to rock, judges the accurately concentration that freeze-thaw cycle experiment is damaged caused by rock sample Position.
2. it is according to claim 1 it is a kind of to rock carry out layer into formula breakdown diagnosis method, it is characterised in that：Step 1 In, rock sample is the cylindrical rock sample that radius is 2.5 centimetres, from the central point of the cross section, marks off 5 equal semidiameters Sample ring, semidiameter be 0.5 centimetre.
CN201610597318.1A 2016-07-26 2016-07-26 A method of layer being carried out into formula breakdown diagnosis to rock Active CN106228558B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610597318.1A CN106228558B (en) 2016-07-26 2016-07-26 A method of layer being carried out into formula breakdown diagnosis to rock

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610597318.1A CN106228558B (en) 2016-07-26 2016-07-26 A method of layer being carried out into formula breakdown diagnosis to rock

Publications (2)

Publication Number Publication Date
CN106228558A CN106228558A (en) 2016-12-14
CN106228558B true CN106228558B (en) 2018-11-09

Family

ID=57532919

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610597318.1A Active CN106228558B (en) 2016-07-26 2016-07-26 A method of layer being carried out into formula breakdown diagnosis to rock

Country Status (1)

Country Link
CN (1) CN106228558B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106680129B (en) * 2017-02-24 2019-08-02 新华水力发电有限公司 Rock sample recycles the simulation of dry and wet freezing-thawing and damage deterioration test device
CN107462492B (en) * 2017-07-17 2019-07-12 三峡大学 The anchor rod anchored research method of large dosage swelling agent is carried out for CT scan technology
CN109211949B (en) * 2018-09-30 2021-06-29 三峡大学 Test method for measuring compactness of anchoring body and surrounding soil layer of expansive cement paste based on CT scanning
CN110297003A (en) * 2019-07-01 2019-10-01 三峡大学 A kind of CT scan post-processing approach being layered based on axis, diameter ring-type

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103512816A (en) * 2013-10-15 2014-01-15 三峡大学 Real-time measuring method for rock sample damage in circulation loading process
CN104330474A (en) * 2014-10-22 2015-02-04 中国地质大学（武汉） Calibration method for rock damage fracture system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9903826B2 (en) * 2014-09-17 2018-02-27 Saudi Arabian Oil Company Multi-objective core-flood test system for oil recovery evaluation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103512816A (en) * 2013-10-15 2014-01-15 三峡大学 Real-time measuring method for rock sample damage in circulation loading process
CN104330474A (en) * 2014-10-22 2015-02-04 中国地质大学（武汉） Calibration method for rock damage fracture system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
冻融循环下混凝土力学性能试验及损伤演化研究;邢凯;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20160115;第四章，第2.3节 *
混凝土单轴动力压缩CT试验过程定量化分区研究;方建银等;《振动与冲击》;20160531;第2-3节 *

Also Published As

Publication number Publication date
CN106228558A (en) 2016-12-14

Similar Documents

Publication Publication Date Title
CN106228558B (en) A method of layer being carried out into formula breakdown diagnosis to rock
Schurink et al. Diffusion-weighted imaging in rectal cancer: current applications and future perspectives
CN106874956A (en) The construction method of image classification convolutional neural networks structure
CN104020135B (en) Calibration model modeling method based near infrared spectrum
CN105844228A (en) Remote sensing image cloud detection method based on convolution nerve network
CN101061951A (en) Method and apparatus for classifying tissue using image data
CN1311247C (en) Multi-measuring NMR analysing based on maximum entropy
CN108446616B (en) Road extraction method based on full convolution neural network ensemble learning
CN105528638A (en) Method for grey correlation analysis method to determine number of hidden layer characteristic graphs of convolutional neural network
CN108229550A (en) A kind of cloud atlas sorting technique that network of forests network is cascaded based on more granularities
CN109671063B (en) Image quality evaluation method based on importance among deep network features
CN105513066B (en) It is a kind of that the generic object detection method merged with super-pixel is chosen based on seed point
CN102177435A (en) Histopathologic diagnosis support system, histopathologic diagnosis support program, and histopathologic diagnosis support method
CN108122221A (en) The dividing method and device of diffusion-weighted imaging image midbrain ischemic area
CN108629772A (en) Image processing method and device, computer equipment and computer storage media
CN109472798A (en) Live pig fat content detection model training method and live pig fat content detection method
CN103399852A (en) Multi-channel spectrum clustering method based on local density estimation and neighbor relation spreading
CN105004745B (en) A method of viscosity of crude is predicted by nuclear magnetic resoance spectrum
CN103985110A (en) Methods for positioning and displaying artifact or intensity inhomogeneity in digital image
CN104126188A (en) Key word detection device, control method and control program for same, and display apparatus
Makowski et al. De novo radiomics approach using image augmentation and features from T1 mapping to predict Gleason scores in prostate cancer
CN105046265B (en) Iris image intestines ring region area detecting method based on texture difference
CN108665546A (en) A kind of multiple spot geological statistics three-dimensional modeling method of combination deep learning
CN106023188A (en) Breast tumor feature selection method based on Relief algorithm
CN104637060A (en) Image partition method based on neighbor-hood PCA (Principal Component Analysis)-Laplace

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201209

Address after: Unit 2414-2416, main building, no.371, Wushan Road, Tianhe District, Guangzhou City, Guangdong Province

Patentee after: GUANGDONG GAOHANG INTELLECTUAL PROPERTY OPERATION Co.,Ltd.

Address before: 443002 No. 8, University Road, Yichang, Hubei

Patentee before: CHINA THREE GORGES University

Effective date of registration: 20201209

Address after: No.005, Wanshui Road, Qianshan Economic Development Zone, Anqing City, Anhui Province

Patentee after: ANHUI PROVINCE WANJIANG MECHANICAL AND ELECTRICAL EQUIPMENT MANUFACTURING Co.,Ltd.

Address before: Unit 2414-2416, main building, no.371, Wushan Road, Tianhe District, Guangzhou City, Guangdong Province

Patentee before: GUANGDONG GAOHANG INTELLECTUAL PROPERTY OPERATION Co.,Ltd.

TR01 Transfer of patent right