CN108804849A - A kind of rock mechanics parameters evaluation method based on structure complexity - Google Patents

A kind of rock mechanics parameters evaluation method based on structure complexity Download PDF

Info

Publication number
CN108804849A
CN108804849A CN201810648437.4A CN201810648437A CN108804849A CN 108804849 A CN108804849 A CN 108804849A CN 201810648437 A CN201810648437 A CN 201810648437A CN 108804849 A CN108804849 A CN 108804849A
Authority
CN
China
Prior art keywords
rock
hugger
analysis
mechanics parameters
sample
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810648437.4A
Other languages
Chinese (zh)
Other versions
CN108804849B (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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201810648437.4A priority Critical patent/CN108804849B/en
Publication of CN108804849A publication Critical patent/CN108804849A/en
Application granted granted Critical
Publication of CN108804849B publication Critical patent/CN108804849B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Abstract

The rock mechanics parameters evaluation method based on structure complexity that the invention discloses a kind of, is related to petrol resources exploitation field.The present invention includes the following steps:Obtain mining site coal rock cleat development information to be evaluated, and measure complex degree of structure, hugger density and its complex degree of structure are subjected to correlation analysis, face cleat quantity and its complex degree of structure are subjected to correlation analysis, obtain the lithostructure complexity evaluation model under two kinds of analysis modes;Respectively by the correlation analysis of the rock mechanics parameters and its complex degree of structure of the coal petrography rock sample under two kinds of analysis modes, obtained rock mechanics parameters evaluation model integrate preferably, the rock mechanics parameters evaluation model based on structure complexity is obtained, is then evaluated.The lithostructure complexity that the present invention is characterized using Fractal has obtained the evaluation method that a large amount of coal rock parameters are obtained by a small amount of standard sample as point of penetration.

Description

A kind of rock mechanics parameters evaluation method based on structure complexity
Technical field
The present invention relates to petrol resources exploitation field, specially a kind of rock mechanics parameters evaluation based on structure complexity Method.
Background technology
Rock mechanics parameters are through the underlying parameter of entire petrol resources exploitation process, and it is general to obtain rock mechanics parameters There are laboratory experiment acquisition and mining site to obtain two kinds of approach.Indoor mechanical test is to obtain one of rock strength most efficient method, But its limitation essentially consist in indoor experiment of machanics cannot simulation underground high-temperature and high-pressure conditions very well, meanwhile, laboratory test A large amount of underground coal samples are generally required, the sampling of big section of underground is very difficult for coal seam reservoirs, to be difficult to a large amount of indoor physics Experiment is to obtain the rock mechanics parameters of coal petrography.Furthermore the coal petrography being crushed for structure is difficult to prepare Standard rock sample, even if sample preparation The success rate that successful coal sample carries out indoor experiment of machanics is also relatively low.A small amount of certified reference coal is only needed therefore, it is necessary to find And the method for obtaining a large amount of rock mechanics parameters.According to existing literature, Fractal is typically used to describe complexity Complexity, in petroleum industry field, Fractal is typically used to evaluate the anisotropism of rock, is typically used to study hollow billet The variation of pressure curve, and then for studying flow unit etc., and application achievements are more.Coal rock cleat is developed, and structure is multiple It is miscellaneous, and complicated structure is affected to the rock mechanics parameters of coal petrography, current existing technology can not accurately be come to this It evaluates.
Invention content
Technical problem to be solved by the invention is to provide a kind of rock mechanics parameters evaluation side based on structure complexity Method, the present invention is directed to the lithostructure complexities characterized using Fractal as point of penetration, obtain through a small amount of standard sample To obtain the evaluation method of a large amount of coal rock parameters.
In order to solve the above technical problems, the technical solution used in the present invention is:A kind of rock based on structure complexity Mechanics parameter evaluation method, which is characterized in that it includes the following steps:
A, to mining site to be evaluated, require to obtain a large amount of coal petrography rock samples according to physics laboratory, and collect include hugger density, The hugger of hugger quantity and hugger angle develops information, and measures the complex degree of structure of coal petrography rock sample;
B, according to the hugger developmental state of coal petrography rock sample in step a, the hugger density of its coal petrography rock sample is answered with its structure Miscellaneous degree carries out correlation analysis, and the hugger quantity of its coal petrography rock sample and its complex degree of structure are carried out correlation analysis;Point It Huo Qu not evaluation model of the coal petrography rock sample based on the structure complexity under hugger density analysis and hugger quantitative analysis two ways;
C, the coal petrography rock sample hugger in mining site to be evaluated is identified and rock is carried out to hugger and base portion respectively Value estimations of mechanical parameters obtains coal rock parameter by numerical simulation mode;It will be based on hugger density analysis and hugger number Structure complexity under amount analysis two ways carries out correlation analysis with its rock mechanics parameters respectively;By two kinds of analysis modes Under correlation analysis integrate preferably, obtain the coal rock parameter evaluation model based on structure complexity, then Coal rock parameter can be evaluated.
Further technical solution is, the hugger development information of the rock sample include face cleat angle, face cleat density, Face cleat quantity, end hugger angle, end hugger density, end hugger quantity.
Further technical solution is that the rock mechanics parameters include uniaxial compressive strength, elasticity modulus and Poisson Than three rock mechanics parameters.
Further technical solution is, the rock mechanics parameters evaluation method based on structure complexity, specific steps For:
A, to mining site to be evaluated, require to obtain a large amount of coal petrography rock samples according to physics laboratory, and collect include hugger density, The hugger of hugger quantity and angle develops information, prepares and screens core sample;
B, according to the hugger developmental state of core sample in step a, with reference to indoor Physical Experiment require to core sample into Row data test, the box counting dimension FD of program calculation two-dimensional mathematical model, then the complex degree of structure of core sample is FD, Jiang Qiyan The hugger density of core sample carries out correlation analysis with its complex degree of structure, and structure changes the coal under hugger density analysis mode Rock structure complexity evaluation model, the following formula of model:
FD=f0(ρ) (1)
In formula:ρ is hugger density, item/5cm;FD is the complex degree of structure of coal sample model, dimensionless;
The hugger density of its core sample and its complex degree of structure are subjected to correlation analysis, structure changes face cleat number Lithostructure complexity evaluation model under analysis mode, model are as follows:
FD=f1(M) (2)
In formula:M is face cleat number, item;FD is the complex degree of structure of coal sample model, dimensionless;
C, the coal petrography rock sample hugger in mining site to be evaluated is identified and rock is carried out to hugger and base portion respectively Value estimations of mechanical parameters;Logarithm simulation result of calculation carries out each computation model can be obtained in processing uniaxial compressive strength, elasticity Three rock mechanics parameters of modulus and Poisson's ratio;By three rock mechanics parameters respectively with based under hugger density analysis mode Structure complexity carries out correlation analysis, show that the rock mechanics parameters based on the coal petrography rock sample under hugger density analysis are tied with it The correlation analysis of structure complexity builds the rock mechanics parameters evaluation based on rock sample complexity under hugger density analysis mode Model is as follows:
UCS=f2(FD) (3)
E=f3(FD) (4)
μ=f4(FD) (5)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, nothing Dimension;FD is rock sample box counting dimension, dimensionless;
Three rock mechanics parameters are subjected to correlation with based on the structure complexity under hugger quantitative analysis mode respectively Analysis show that the rock mechanics based on the coal petrography rock sample under face cleat quantitative analysis participates in the correlation point of its complex degree of structure Analysis, rock mechanics parameters evaluation model of the structure based on rock sample complexity under hugger quantitative analysis mode are as follows:
UCS=f5(FD) (6)
E=f6(FD) (7)
μ=f7(FD) (8)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, nothing Dimension;FD is rock sample box counting dimension, dimensionless.
D, the rock mechanics parameters evaluation model that two kinds of analysis modes obtain in comparison step c;
The evaluation model obtained according to two kinds of analysis modes is close, then by the Data Integration of two kinds of analysis modes together into Row correlation analysis obtains the rock mechanics parameters evaluation model based on structure complexity;
Or, it is larger according to the computation model difference that two kinds of analysis modes obtain, then with practical indoor experiment of machanics data difference The result of calculation of the evaluation model obtained with two kinds of analysis modes is compared, and is preferably gone out and practical indoor experiment of machanics result phase The evaluation model smaller to error, the rock mechanics parameters based on structure complexity as the rock mechanics parameters evaluate mould Type.
Further technical solution is that the core sample is to drill through acquirement, is cylinder and core axis and face Plane normal angular separation is between 0 °~90 ° where hugger, and core grinding two end surfaces, two end faces are perpendicular to core axis.
Further technical solution is that steps are as follows by the box counting dimension FD of program calculation two-dimensional mathematical model:
The hugger developmental state collection of core sample is handled first, difference is denoted as according to the combination of different number pixel " box " of size, note " box " size are εi, wherein i is different subdivision number;It is repeatedly cutd open with various sizes of " box " Sub-model counts " box " number N passed through by hugger under each sizei, it is one group available that a subdivision is often carried out to model Data (εi, Ni), Q group data (ε then can be obtained in Q divisioni, Ni), with logNiFor ordinate, log εiFor abscissa, make logNi-logεiRelational graph and carry out linear fit, then the negative value for being fitted the slope of formula is the box counting dimension FD of corresponding model, i.e., For the complex degree of structure value FD of corresponding model.
Further technical solution is, carries out the correlation analysis under two kinds of analysis modes to integrate preferred method For:
The rock mechanics parameters evaluation model that two kinds of analysis modes of comparison obtain;
The evaluation model obtained according to two kinds of analysis modes is close, then by the Data Integration of two kinds of analysis modes together into Row correlation analysis obtains the rock mechanics parameters evaluation model based on structure complexity;
Or, it is larger according to the computation model difference that two kinds of analysis modes obtain, then with practical indoor experiment of machanics data difference The result of calculation of the evaluation model obtained with two kinds of analysis modes is compared, and is preferably gone out and practical indoor experiment of machanics result phase The evaluation model smaller to error, the rock mechanics parameters based on structure complexity as the rock mechanics parameters evaluate mould Type.
It is using advantageous effect caused by above-mentioned technical proposal:The lithostructure that the present invention is characterized with Fractal Complexity is point of penetration, has obtained the evaluation method that a large amount of coal rock parameters are obtained by a small amount of standard sample.
Description of the drawings
The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
Fig. 1 is to change face cleat quantity mode to build different structure complexity coal petrography model;
Fig. 2 is to change face cleat density manner structure different structure complexity coal petrography model;
Fig. 3 is to change face cleat quantity mode model box counting dimension calculating figure;
Fig. 4 is to change face cleat density manner model box counting dimension calculating figure;
Fig. 5 a- Fig. 5 b are the relationships of face cleat quantity and face cleat density and box counting dimension;
Fig. 6 a- Fig. 6 c are rock mechanics parameters and box counting dimension relationship analysis;
Fig. 7 a- Fig. 7 c are rock mechanics parameters and box counting dimension relationship analysis.
Specific implementation mode
With reference to the attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete Ground describes, it is clear that described embodiment is only a part of the embodiment of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
Many details are elaborated in the following description to facilitate a thorough understanding of the present invention, still the present invention can be with Implemented different from other manner described here using other, those skilled in the art can be without prejudice to intension of the present invention In the case of do similar popularization, therefore the present invention is not limited by following public specific embodiment.
The rock mechanics parameters evaluation method based on structure complexity that the present invention provides a kind of, which is characterized in that it is wrapped Include following steps:
A, to mining site to be evaluated, require to obtain a large amount of coal petrography rock samples according to physics laboratory, and collect include hugger density, The hugger of hugger quantity and hugger angle develops information, and measures the complex degree of structure of coal petrography rock sample;
B, according to the hugger developmental state of coal petrography rock sample in step a, the hugger density of its coal petrography rock sample is answered with its structure Miscellaneous degree carries out correlation analysis, and the hugger quantity of its coal petrography rock sample and its complex degree of structure are carried out correlation analysis;Point It Huo Qu not evaluation model of the coal petrography rock sample based on the structure complexity under hugger density analysis and hugger quantitative analysis two ways;
C, the coal petrography rock sample hugger in mining site to be evaluated is identified and rock is carried out to hugger and base portion respectively Value estimations of mechanical parameters obtains coal rock parameter by numerical simulation mode;It will be based on hugger density analysis and hugger number Structure complexity under amount analysis two ways carries out correlation analysis with its rock mechanics parameters respectively;By two kinds of analysis modes Under correlation analysis integrate preferably, obtain the coal rock parameter evaluation model based on structure complexity, then Coal rock parameter can be evaluated.
Carry out the correlation analysis under two kinds of analysis modes to integrate preferred method be:
The rock mechanics parameters evaluation model that two kinds of analysis modes of comparison obtain, if the evaluation mould that two kinds of analysis modes obtain Type is close, then the Data Integration of two kinds of analysis modes is carried out correlation analysis together, obtains the rock based on structure complexity Stone mechanics parameter evaluation model, if the computation model difference that two kinds of analysis modes obtain is larger, with practical indoor experiment of machanics number It is compared, is preferably gone out and practical indoor experiment of machanics according to the result of calculation of the evaluation model obtained respectively with two kinds of analysis modes As a result the smaller evaluation model of relative error, the rock mechanics parameters based on structure complexity as the rock mechanics parameters are commented Valence model.
Preferably, the hugger development information of the rock sample includes face cleat angle, face cleat density, face cleat quantity, end Hugger angle, end hugger density, end hugger quantity.
Preferably, the rock mechanics parameters include three uniaxial compressive strength, elasticity modulus and Poisson's ratio rock mechanics Parameter.
Preferably, the rock mechanics parameters evaluation method based on structure complexity, which is characterized in that its specific steps are:
A, to mining site to be evaluated, require to obtain a large amount of coal petrography rock samples according to physics laboratory, and collect include hugger density, The hugger of hugger quantity and angle develops information, prepares and screens core sample;
B, according to the hugger developmental state of core sample in step a, with reference to indoor Physical Experiment require to core sample into Row data test, the box counting dimension FD of program calculation two-dimensional mathematical model, then the complex degree of structure of core sample is FD, Jiang Qiyan The hugger density of core sample carries out correlation analysis with its complex degree of structure, and structure changes the coal under hugger density analysis mode Rock structure complexity evaluation model, the following formula of model:
FD=f0(ρ) (1)
In formula:ρ is hugger density, item/5cm;FD is the complex degree of structure of coal sample model, dimensionless;
The hugger density of its core sample and its complex degree of structure are subjected to correlation analysis, structure changes face cleat number Lithostructure complexity evaluation model under analysis mode, model are as follows:
FD=f1(M) (2)
In formula:M is face cleat number, item;FD is the complex degree of structure of coal sample model, dimensionless;
C, the coal petrography rock sample hugger in mining site to be evaluated is identified and rock is carried out to hugger and base portion respectively Value estimations of mechanical parameters;Logarithm simulation result of calculation carries out each computation model can be obtained in processing uniaxial compressive strength, elasticity Three rock mechanics parameters of modulus and Poisson's ratio;By three rock mechanics parameters respectively with based under hugger density analysis mode Structure complexity carries out correlation analysis, show that the rock mechanics parameters based on the coal petrography rock sample under hugger density analysis are tied with it The correlation analysis of structure complexity builds the rock mechanics parameters evaluation based on rock sample complexity under hugger density analysis mode Model is as follows:
UCS=f2(FD) (3)
E=f3(FD) (4)
μ=f4(FD) (5)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, nothing Dimension;FD is rock sample box counting dimension, dimensionless;
Three rock mechanics parameters are subjected to correlation with based on the structure complexity under hugger quantitative analysis mode respectively Analysis show that the rock mechanics based on the coal petrography rock sample under face cleat quantitative analysis participates in the correlation point of its complex degree of structure Analysis, rock mechanics parameters evaluation model of the structure based on rock sample complexity under hugger quantitative analysis mode are as follows:
UCS=f5(FD) (6)
E=f6(FD) (7)
μ=f7(FD) (8)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, nothing Dimension;FD is rock sample box counting dimension, dimensionless.
D, the rock mechanics parameters evaluation model that two kinds of analysis modes obtain in comparison step c;
The evaluation model obtained according to two kinds of analysis modes is close, then by the Data Integration of two kinds of analysis modes together into Row correlation analysis obtains the rock mechanics parameters evaluation model based on structure complexity;
Or, it is larger according to the computation model difference that two kinds of analysis modes obtain, then with practical indoor experiment of machanics data difference The result of calculation of the evaluation model obtained with two kinds of analysis modes is compared, and is preferably gone out and practical indoor experiment of machanics result phase The evaluation model smaller to error, the rock mechanics parameters based on structure complexity as the rock mechanics parameters evaluate mould Type.
Preferably, the core sample is to drill through acquirement, is cylinder and core axis and planar process where face cleat Line angular separation is between 0 °~90 °, core grinding two end surfaces, and two end faces are perpendicular to core axis.
Preferably, steps are as follows by the box counting dimension FD of program calculation two-dimensional mathematical model:
The hugger developmental state collection of core sample is handled first, difference is denoted as according to the combination of different number pixel " box " of size, note " box " size are εi, wherein i is different subdivision number;It is repeatedly cutd open with various sizes of " box " Sub-model counts " box " number N passed through by hugger under each sizei, it is one group available that a subdivision is often carried out to model Data (εi, Ni), Q group data (ε then can be obtained in Q divisioni, Ni), with logNiFor ordinate, log εiFor abscissa, make logNi-logεiRelational graph and carry out linear fit, then the negative value for being fitted the slope of formula is the box counting dimension FD of corresponding model, i.e., For the complex degree of structure value FD of corresponding model.
Embodiment
Using the rock mechanics parameters evaluation method based on structure complexity, which is characterized in that this method includes following step Suddenly:
The first step carries out observation description to mining site coal petrography to be evaluated, and the hugger development information for obtaining coal petrography rock sample (including is cut Reason density and angle) it prepares and screens core sample;
Coal rock cleat is developed, and the core sample drilled through is cylinder and core axis and plane normal direction where face cleat Angle is between 0 °~90 °, core grinding two end surfaces, and the hugger of rock sample is observed and counted in two end faces perpendicular to core axis Developmental state, including face cleat angle, face cleat density, end hugger angle, end hugger density.
Second step, according to practical rock sample size, with reference to practical rock sample hugger developmental state from change face cleat quantity and Change two modes of face cleat density to establish the computation model (as shown in Figure 1 and Figure 2) of different structure complexity.
Third walks, the complex degree of structure of each model of program calculation.
Fig. 1, Fig. 2 institute representation model are read in into memory first, the combination of different number pixel is denoted as various sizes of " box Son ", note " box " size are εi(i is different subdivision number), with various sizes of " box " multiple subdivision model, statistics is every " box " number N passed through by hugger under a sizei, a subdivision is often carried out to model, and one group of data (ε can be obtainedi, Ni), Q times Q group data (ε then can be obtained in divisioni, Ni), with logNiFor ordinate, log εiFor abscissa, logN is madei-logεiRelationship Scheme and carry out linear fit, then the negative value for being fitted the slope of formula is the box counting dimension FD of corresponding model, as corresponds to the structure of model Complexity value FD.
4th step calculates each model using RFPA numerical simulation softwares and carries out rock mechanics parameters.
Fig. 1, Fig. 2 institute representation model are read into RFPA softwares first, hugger is identified and respectively to hugger and matrix Part carries out rock mechanics parameters assignment, logarithm simulation result of calculation carry out processing can be obtained each computation model single shaft it is anti- Three Compressive Strength, elasticity modulus and Poisson's ratio rock mechanics parameters.
5th step builds the rock mechanics parameters evaluation model based on rock sample complexity.
The rock mechanics that third is walked to the box counting dimension for calculating each model of gained and each model of four-step calculation gained is joined Numerical value carries out correlation analysis, can obtain the rock mechanics parameters evaluation model based on rock sample complexity, based on hugger density point Under analysis, rock mechanics parameters evaluation model following formula of the structure based on rock sample complexity is shown below:
UCS=f2(FD) (3)
E=f3(FD) (4)
μ=f4(FD) (5)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, nothing Dimension;FD is rock sample box counting dimension, dimensionless;
Based under face cleat quantitative analysis, rock mechanics parameters evaluation model of the structure based on rock sample complexity is public as follows Formula:
UCS=f5(FD) (6)
E=f6(FD) (7)
μ=f7(FD) (8)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, nothing Dimension;FD is rock sample box counting dimension, dimensionless.
D, the rock mechanics parameters evaluation model that two kinds of analysis modes obtain in comparison step c;
The evaluation model obtained according to two kinds of analysis modes is close, then by the Data Integration of two kinds of analysis modes together into Row correlation analysis obtains the rock mechanics parameters evaluation model based on structure complexity;
Or, it is larger according to the computation model difference that two kinds of analysis modes obtain, then with practical indoor experiment of machanics data difference The result of calculation of the evaluation model obtained with two kinds of analysis modes is compared, and is preferably gone out and practical indoor experiment of machanics result phase The evaluation model smaller to error, the rock mechanics parameters based on structure complexity as the rock mechanics parameters evaluate mould Type.
Calculated examples
A kind of rock mechanics parameters evaluation method based on structure complexity
The first step, coal petrography rock sample used in this example are the Standard rock sample of length 50mm, diameter 25mm, and face cleat angle is about For 80 ° (angles for being defined as bed plane and one end section of rock sample), face cleat density is 9-15 items/5cm, and end hugger density is 7-14 items/5cm.
Second step, according to practical rock sample size, with reference to practical rock sample hugger developmental state from change face cleat quantity and Change two modes of face cleat density to establish the computation model (Fig. 1, Fig. 2) of different structure complexity.
Change under face cleat quantity mode, face cleat quantity variation range is 3-14 items, and face cleat quantity change step is 0.5, amount to 23 computation models (Fig. 1).
Change under face cleat density manner, face cleat variable density ranging from 6-19 items/5cm, face cleat variable density step A length of 1/5cm amounts to 16 computation models (Fig. 2).
Third walks, the complex degree of structure of each model of program calculation.
The box counting dimension of each model of program calculation, as shown in Figure 3, Figure 4, to logNi-logεiRelational graph go forward side by side line Fitting, the then negative value for being fitted the slope of formula are the box counting dimension FD of corresponding model, as correspond to the complex degree of structure value FD of model. The box counting dimension value of each model is as shown in table 1.
Fig. 1, Fig. 2 institute representation model are read in into memory first, the combination of different number pixel is denoted as various sizes of " box Son ", note " box " size are εi(i is different subdivision number, and value is respectively:2 pixels, 4 pixels, 6 pixel ..., 46 pixel, 48 Pixel, 50 pixels), with various sizes of " box " multiple subdivision model, count " box " passed through by hugger under each size Number Ni, a subdivision is often carried out to model, and one group of data (ε can be obtainedi, Ni), 25 groups of data (ε then can be obtained in 25 divisionsi, Ni), with logNiFor ordinate, log εiFor abscissa, logN is madei-logεiRelational graph and carry out linear fit (Fig. 3, figure 4) negative value for, being then fitted the slope of formula is the box counting dimension FD of corresponding model, as corresponds to the complex degree of structure value FD (tables of model 1, table 2).
Table 1 changes face cleat quantity mode model box counting dimension
Pattern number Hugger quantity (item) Box counting dimension
M1 3.00 1.1939
M2 3.50 1.2905
M3 4.00 1.2929
M4 4.50 1.3102
M5 5.00 1.3583
M6 5.50 1.4069
M7 6.00 1.3936
M8 6.50 1.3937
M9 7.00 1.4319
M10 7.50 1.4662
M11 8.00 1.4544
M12 8.50 1.4533
M13 9.00 1.4786
M14 9.50 1.5043
M15 10.00 1.4926
M16 10.50 1.492
M17 11.00 1.5105
M18 11.50 1.5181
M19 12.00 1.5346
M20 12.50 1.5514
M21 13.00 1.5656
M22 13.50 1.579
M23 14.00 1.5881
Table 2 changes face cleat density manner model box counting dimension
Correlation analysis is carried out with the box counting dimension of corresponding model to face cleat quantity, face cleat density, such as Fig. 5 a- Fig. 5 b institutes Show, coal petrography sample box counting dimension with hugger quantity and hugger density with having preferable correlation, with hugger quantity or hugger The increase of density, coal petrography sample box counting dimension gradually increase, and coal petrography sample box counting dimension can be used in characterizing its complex degree of structure, then coal petrography Structure complexity evaluation model is:
4th step calculates each model using RFPA numerical simulation softwares and carries out rock mechanics parameters.
Fig. 1 models are read into RFPA softwares first, hugger is identified and hugger and base portion are carried out respectively Rock mechanics parameters assignment (table 3), logarithm simulation result of calculation are handled, and each calculating can be obtained from stress-strain diagram Uniaxial compressive strength, three rock mechanics parameters of elasticity modulus and Poisson's ratio (table 4, table 5) of model.
3 model matrix calculating parameter of table
Table 4 changes hugger quantity mode drag mechanics parameter result of calculation
Pattern number Hugger quantity (item) Uniaxial compressive strength (MPa) Elasticity modulus (MPa) Poisson's ratio
M1 3.00 19.093 7220.8 0.3462
M2 3.50 18.120 6945.7 0.3468
M3 4.00 17.287 6897.0 0.3482
M4 4.50 16.960 6827.0 0.3562
M5 5.00 16.662 6790.2 0.3552
M6 5.50 17.231 6741.8 0.3538
M7 6.00 15.976 6683.2 0.361
M8 6.50 15.124 6608.4 0.3716
M9 7.00 15.647 6574.3 0.3724
M10 7.50 15.841 6525.0 0.3686
M11 8.00 15.489 6477.4 0.3812
M12 8.50 14.227 6414.2 0.3956
M13 9.00 14.447 6377.6 0.3954
M14 9.50 15.227 6332.5 0.3902
M15 10.00 15.146 6278.3 0.4018
M16 10.50 14.330 6199.8 0.42
M17 11.00 14.506 6178.2 0.4152
M18 11.50 14.335 6121.3 0.4326
M19 12.00 14.544 6071.9 0.4266
M20 12.50 14.310 6036.2 0.4196
M21 13.00 14.221 6006.3 0.4198
M22 13.50 13.984 5952.8 0.4202
M23 14.00 14.029 5945.8 0.4094
Table 5 changes hugger density manner model mechanics parameter result of calculation
5th step builds the rock mechanics parameters evaluation model based on rock sample complexity.
The rock mechanics that third is walked to the box counting dimension for calculating each model of gained and each model of four-step calculation gained is joined Numerical value carries out correlation analysis, can obtain the rock mechanics parameters evaluation model based on rock sample complexity, correlation analysis is as follows Shown in Fig. 6.
Table 7 is shown under two kinds of complexity makes, the relational model system of coal sample box counting dimension and rock mechanics parameters Meter.3 groups of fit correlations that two kinds of different modes of comparison obtain, two uniaxial compressive strength relational model formulas, elasticity modulus relationships Model consistency is higher, and two Poisson's ratio relational models differ greatly, this is because changing under hugger quantity mode, the face of model Hugger does not run through, and causes the box counting dimension value of corresponding model relatively small, and the effect of model " weak structure face " dies down, " micro- There are more radial deformations before the stage of development of crack ", and then Poisson ratio is increased.Therefore, in follow-up study, single shaft is anti- The acquisition modes of Compressive Strength and elasticity modulus are:Face cleat quantity will be changed and change two groups of data of face cleat density and merge progress Obtained by correlation analysis (Fig. 7 a, Fig. 7 b).Example evaluates the Poisson's ratio of coal petrography under two kinds of analysis modes, and by evaluation result and room Interior experimental result is compared (table 6), and as shown in Table 6, the error for changing hugger density analysis mode is relatively small, therefore, pool The acquisition modes of loose ratio are:Data to change hugger density manner carry out obtained by correlation analysis (Fig. 7 c).
The Poisson's ratio of 6 two kinds of analysis modes of table and measured result contrast table
The relational model of coal sample box counting dimension and rock mechanics parameters under 7 two kinds of complexity makes of table
That is, the rock mechanics parameters of different structure complexity shale and the relationship of box counting dimension are respectively:
UCS=-12.3531FD+33.345
E=-4983.7FD+13608
μ=0.0669FD+0.2466
Wherein, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, immeasurable Guiding principle;FD is rock sample box counting dimension, dimensionless.
Above formula is the rock sample rock mechanics parameters evaluation model based on structure complexity.
The rock mechanics parameters of same regional 10 lump coal rock are evaluated using above-mentioned evaluation model, and are pressed with practical single shaft Contracting test data is compared, and data are as shown in table 7.
8 prediction rock mechanics parameter of table and actual measurement rock mechanics parameters contrast table
As seen from the above table, the error of prediction three obtained rock mechanics parameters is all within 9%, therefore, institute of the present invention The coal rock parameter evaluation method of foundation has preferable applicability.

Claims (7)

1. a kind of rock mechanics parameters evaluation method based on structure complexity, which is characterized in that it includes the following steps:
A, it to mining site to be evaluated, requires to obtain a large amount of coal petrography rock samples according to physics laboratory, and it includes hugger density, hugger to collect The hugger of quantity and hugger angle develops information, and measures the complex degree of structure of coal petrography rock sample;
B, according to the hugger developmental state of coal petrography rock sample in step a, by the hugger density of its coal petrography rock sample and its complicated journey Degree carries out correlation analysis, and the hugger quantity of its coal petrography rock sample and its complex degree of structure are carried out correlation analysis;It obtains respectively Take coal petrography rock sample based on the evaluation model of the structure complexity under hugger density analysis and hugger quantitative analysis two ways;
C, the coal petrography rock sample hugger in mining site to be evaluated is identified and rock mechanics is carried out to hugger and base portion respectively Parameter assignment obtains coal rock parameter by numerical simulation mode;It will be based on hugger density analysis and hugger quantity point The structure complexity analysed under two ways carries out correlation analysis with its rock mechanics parameters respectively;It will be under two kinds of analysis modes Correlation analysis integrate preferably, obtains the coal rock parameter evaluation model based on structure complexity, then can be right Coal rock parameter is evaluated.
2. a kind of rock mechanics parameters evaluation method based on structure complexity according to claim 1, which is characterized in that The hugger development information of the rock sample includes face cleat angle, face cleat density, face cleat quantity, end hugger angle, end hugger Density, end hugger quantity.
3. a kind of rock mechanics parameters evaluation method based on structure complexity according to claim 1, which is characterized in that The rock mechanics parameters include three uniaxial compressive strength, elasticity modulus and Poisson's ratio rock mechanics parameters.
4. a kind of rock mechanics parameters evaluation method based on structure complexity according to claim 1 or 2 or 3, special Sign is, the specific steps are:
A, it to mining site to be evaluated, requires to obtain a large amount of coal petrography rock samples according to physics laboratory, and it includes hugger density, hugger to collect The hugger of quantity and angle develops information, prepares and screens core sample;
B, it according to the hugger developmental state of core sample in step a, is required to core sample into line number with reference to indoor Physical Experiment According to test, the box counting dimension FD of program calculation two-dimensional mathematical model, then the complex degree of structure of core sample is FD, its core is tried The hugger density of sample carries out correlation analysis with its complex degree of structure, and structure changes the coal petrography knot under hugger density analysis mode Structure complexity evaluation model, the following formula of model:
FD=f0(ρ) (1)
In formula:ρ is hugger density, item/5cm;FD is the complex degree of structure of coal sample model, dimensionless;
The hugger density of its core sample and its complex degree of structure are subjected to correlation analysis, structure changes the analysis of face cleat number Lithostructure complexity evaluation model under mode, model are as follows:
FD=f1(M) (2)
In formula:M is face cleat number, item;FD is the complex degree of structure of coal sample model, dimensionless;
C, the coal petrography rock sample hugger in mining site to be evaluated is identified and rock mechanics is carried out to hugger and base portion respectively Parameter assignment;Logarithm simulation result of calculation carries out each computation model can be obtained in processing uniaxial compressive strength, elasticity modulus And three rock mechanics parameters of Poisson's ratio;By three rock mechanics parameters respectively with based on the structure under hugger density analysis mode Complexity carries out correlation analysis, show that the rock mechanics parameters based on the coal petrography rock sample under hugger density analysis are answered with its structure The correlation analysis of miscellaneous degree builds the rock mechanics parameters evaluation model based on rock sample complexity under hugger density analysis mode It is as follows:
UCS=f2(FD) (3)
E=f3(FD) (4)
μ=f4(FD) (5)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, immeasurable Guiding principle;FD is rock sample box counting dimension, dimensionless;
Three rock mechanics parameters are subjected to correlation analysis with based on the structure complexity under hugger quantitative analysis mode respectively, Show that the rock mechanics based on the coal petrography rock sample under face cleat quantitative analysis participates in the correlation analysis of its complex degree of structure, structure It is as follows to build the rock mechanics parameters evaluation model based on rock sample complexity under hugger quantitative analysis mode:
UCS=f5(FD) (6)
E=f6(FD) (7)
μ=f7(FD) (8)
In above formula, UCS is rock specimen in uniaxial compression strength, MPa;E is rock sample elasticity modulus, MPa;μ is rock sample Poisson's ratio, immeasurable Guiding principle;FD is rock sample box counting dimension, dimensionless;
D, the rock mechanics parameters evaluation model that two kinds of analysis modes obtain in comparison step c;
The evaluation model obtained according to two kinds of analysis modes is close, then the Data Integration of two kinds of analysis modes is carried out phase together The analysis of closing property, obtains the rock mechanics parameters evaluation model based on structure complexity;
Or, the computation model difference obtained according to two kinds of analysis modes is larger, then use practical indoor experiment of machanics data respectively with two The result of calculation for the evaluation model that kind analysis mode obtains is compared, and is preferably gone out opposite accidentally with practical indoor experiment of machanics result The smaller evaluation model of difference, the rock mechanics parameters evaluation model based on structure complexity as the rock mechanics parameters.
5. a kind of rock mechanics parameters evaluation method based on structure complexity according to claim 3, which is characterized in that The core sample is to drill through acquirements, be cylinder and core axis with face cleat place plane normal angular separation 0 °~ Between 90 °, core grinding two end surfaces, two end faces are perpendicular to core axis.
6. a kind of rock mechanics parameters evaluation method based on structure complexity according to claim 3, which is characterized in that Steps are as follows by the box counting dimension FD of program calculation two-dimensional mathematical model:
The hugger developmental state collection of core sample is handled first, different sizes are denoted as according to the combination of different number pixel " box ", note " box " size be εi, wherein i is different subdivision number;With various sizes of " box " multiple subdivision mould Type counts " box " number N passed through by hugger under each sizei, a subdivision is often carried out to model, and one group of data can be obtained (εi, Ni), Q group data (ε then can be obtained in Q divisioni, Ni), with logNiFor ordinate, log εiFor abscissa, logN is madei- logεiRelational graph and carry out linear fit, then the negative value for being fitted the slope of formula is the box counting dimension FD of corresponding model, as corresponding The complex degree of structure value FD of model.
7. a kind of rock mechanics parameters evaluation method based on structure complexity according to claim 1, which is characterized in that Carry out the correlation analysis under two kinds of analysis modes to integrate preferred method be:
The rock mechanics parameters evaluation model that two kinds of analysis modes of comparison obtain;
The evaluation model obtained according to two kinds of analysis modes is close, then the Data Integration of two kinds of analysis modes is carried out phase together The analysis of closing property, obtains the rock mechanics parameters evaluation model based on structure complexity;
Or, the computation model difference obtained according to two kinds of analysis modes is larger, then use practical indoor experiment of machanics data respectively with two The result of calculation for the evaluation model that kind analysis mode obtains is compared, and is preferably gone out opposite accidentally with practical indoor experiment of machanics result The smaller evaluation model of difference, the rock mechanics parameters evaluation model based on structure complexity as the rock mechanics parameters.
CN201810648437.4A 2018-06-22 2018-06-22 Rock mechanical parameter evaluation method based on structural complexity Active CN108804849B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810648437.4A CN108804849B (en) 2018-06-22 2018-06-22 Rock mechanical parameter evaluation method based on structural complexity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810648437.4A CN108804849B (en) 2018-06-22 2018-06-22 Rock mechanical parameter evaluation method based on structural complexity

Publications (2)

Publication Number Publication Date
CN108804849A true CN108804849A (en) 2018-11-13
CN108804849B CN108804849B (en) 2022-04-26

Family

ID=64084542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810648437.4A Active CN108804849B (en) 2018-06-22 2018-06-22 Rock mechanical parameter evaluation method based on structural complexity

Country Status (1)

Country Link
CN (1) CN108804849B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109596713A (en) * 2018-12-29 2019-04-09 西南石油大学 A kind of coal petrography ultrasonic wave attenuation coefficient evaluation method based on Fractal
CN109933844A (en) * 2019-01-28 2019-06-25 西南石油大学 A kind of rock rupture complexity characterizing method based on fractal dimension
AU2019429806B2 (en) * 2019-04-18 2021-11-04 China University Of Mining And Technology Method for evaluating degree of complexity of regional geological structure
CN113758839A (en) * 2021-07-19 2021-12-07 山东大学 Large-scale rock simulation method and system based on coarse grained bonding model
CN115407045A (en) * 2022-08-02 2022-11-29 西南石油大学 Rock mechanical parameter evaluation model construction method and rock mechanical characteristic evaluation method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014205248A2 (en) * 2013-06-19 2014-12-24 Conocophillips Company Mechanical characterization of core samples
CN104314563A (en) * 2014-10-21 2015-01-28 西安科技大学 Logging quantitative evaluation method of coal bed methane reservoir fracturing capability
CN104453872A (en) * 2014-10-22 2015-03-25 西安科技大学 Sound wave time difference hole enlargement affect correction method based on coal rock industrial components
CN106295042A (en) * 2016-08-17 2017-01-04 西安科技大学 A kind of coal seam top rock stability Quantitative Evaluation with Well Logging method
CN107132334A (en) * 2017-04-28 2017-09-05 中国科学院地质与地球物理研究所 Rock physical and mechanic parameter intelligent integral test system and its method of testing

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014205248A2 (en) * 2013-06-19 2014-12-24 Conocophillips Company Mechanical characterization of core samples
CN104314563A (en) * 2014-10-21 2015-01-28 西安科技大学 Logging quantitative evaluation method of coal bed methane reservoir fracturing capability
CN104453872A (en) * 2014-10-22 2015-03-25 西安科技大学 Sound wave time difference hole enlargement affect correction method based on coal rock industrial components
CN106295042A (en) * 2016-08-17 2017-01-04 西安科技大学 A kind of coal seam top rock stability Quantitative Evaluation with Well Logging method
CN107132334A (en) * 2017-04-28 2017-09-05 中国科学院地质与地球物理研究所 Rock physical and mechanic parameter intelligent integral test system and its method of testing

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
和志浩: "煤岩井壁稳定性的岩石力学基础研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109596713A (en) * 2018-12-29 2019-04-09 西南石油大学 A kind of coal petrography ultrasonic wave attenuation coefficient evaluation method based on Fractal
CN109933844A (en) * 2019-01-28 2019-06-25 西南石油大学 A kind of rock rupture complexity characterizing method based on fractal dimension
AU2019429806B2 (en) * 2019-04-18 2021-11-04 China University Of Mining And Technology Method for evaluating degree of complexity of regional geological structure
CN113758839A (en) * 2021-07-19 2021-12-07 山东大学 Large-scale rock simulation method and system based on coarse grained bonding model
CN115407045A (en) * 2022-08-02 2022-11-29 西南石油大学 Rock mechanical parameter evaluation model construction method and rock mechanical characteristic evaluation method

Also Published As

Publication number Publication date
CN108804849B (en) 2022-04-26

Similar Documents

Publication Publication Date Title
CN108804849A (en) A kind of rock mechanics parameters evaluation method based on structure complexity
CN109594968A (en) Fracture parameters evaluation method and system after a kind of shale gas multistage pressure break horizontal well pressure
Martinez et al. Bio-geochemical reactive transport modeling of microbial induced calcite precipitation to predict the treatment of sand in one-dimensional flow
Lemenkov et al. Using TeX markup language for 3D and 2D geological plotting
Manouchehrian et al. Application of artificial neural networks and multivariate statistics to estimate UCS using textural characteristics
CN105484739B (en) Pore Pressure on Carbonate Rock Formation test method and device
CN104533400B (en) Method for reconstructing logging curve
CN106841001A (en) A kind of tight sand porosity based on reservoir quality Analysis The Main Control Factor, Permeability Prediction method
Ozturk et al. Strength classification of rock material based on textural properties
CN109214026A (en) A kind of shale gas horizontal well initial productivity prediction technique
Elbadry Simplified reliable prediction method for determining the volume change of expansive soils based on simply physical tests
Liu et al. Estimating the joint roughness coefficient of rock joints from translational overlapping statistical parameters
Aligholi et al. Evaluating the relationships between NTNU/SINTEF drillability indices with index properties and petrographic data of hard igneous rocks
Taboada et al. Evaluation of the quality of a granite quarry
Li et al. A novel approach of two-dimensional representation of rock fracture network characterization and connectivity analysis
Yesiloglu-Gultekin et al. A comparison among some non-linear prediction tools on indirect determination of uniaxial compressive strength and modulus of elasticity of basalt
Singer et al. Quarry Analysis at Bodie Hills, Mono County, California: A Case Study
CN105372631A (en) Polarizing radar inversion method based on genetic-partial least square algorithm, and application of polarizing radar inversion method
Geng et al. Application of full-scale experimental cutterhead system to study penetration performance of tunnel boring machines (TBMs)
CN112946780B (en) Method and device for determining running and sliding fracture activity period
CN106383053A (en) Engineering mechanical parameter related brittleness index prediction method
Li et al. A coupling model based on grey relational analysis and stepwise discriminant analysis for wood defect area identification by stress wave
Molina-Gómez et al. Assessment of groundwater level variations using multivariate statistical methods
Wang et al. Reference strain γr in hyperbolic modeling of dynamic shear modulus of soils
CN112523748B (en) Multi-stage and multi-dimensional evaluation method and device for volume fracturing effect of compact oil

Legal Events

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