CN110069844A - A kind of thin sight numerical model generation method considering rock texture feature and mineral composition - Google Patents
A kind of thin sight numerical model generation method considering rock texture feature and mineral composition Download PDFInfo
- Publication number
- CN110069844A CN110069844A CN201910298582.9A CN201910298582A CN110069844A CN 110069844 A CN110069844 A CN 110069844A CN 201910298582 A CN201910298582 A CN 201910298582A CN 110069844 A CN110069844 A CN 110069844A
- Authority
- CN
- China
- Prior art keywords
- rock
- mineral
- particle
- grain
- numerical model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a kind of thin sight numerical model generation methods for considering rock texture feature and mineral composition.It can be considered that the microscopical structure feature of rock and mineral composition, establish the thin sight numerical model of reflection rock texture feature and mineral composition, can more really carry out rock mechanics numerical simulation.The present invention is the following steps are included: obtain the information such as the size, form, mineral type of mineral grain to rock sample progress rock texture signature analysis；The mineral grain of irregular shape is reduced to the ball of a size reduction, makes particle that permutatation occur using partial size plavini, obtains a more closely knit structure；Then the information such as particle space positions and dimensions are extracted, the Spacial domain decomposition based on Voronoi diagram is carried out to aggregates body；Obtained Voronoi polycrystalline structure is finally subjected to FEM meshing, and insertion generates the finite element numerical model that can reflect rock texture and lithofacies structure without thickness boundary element inside mineral grain boundary and particle.
Description
Technical field
The invention belongs to ROCK MECHANICS RESEARCH fields, and in particular to a kind of to consider the thin of rock texture feature and mineral composition
See numerical model generation method.
Background technique
Rock is the basic component of the earth's crust, and a large amount of exposures constitute mountains and rivers valley, be Human dried bloodstains in earth's surface
Underlying carrier and environment.Rock is the solid-state aggregate by one or more of mineral compositions, with stable topography, and mineral
Have certain chemical component and crystal structure.Rock texture refers to the crystallization degree of substance of composition rock, mineral grain
The feature that size, the relative size of mineral grain, the shape of mineral grain and the correlation between them are showed.
Structure can be divided into eucrystalline texture, hypocrystalline texture and glassy struc three categories by the crystallization degree of rock；By in rock
The absolute size of mineral grain can separate coarse grain, middle grain, particulate, the other structure of particulate class；By the relative size of mineral grain
Granulitic texture, inequigranular texture and plaque-like, porphyritic-like texture can be divided again；By the euhedral degree of rock Minerals, can also separate
Idiomorphism structure, hypidiomorphic texture and anhedral texture.
During the Research on Mechanical Properties of rock, size, shape, interlocking degree and the contact type of mineral grain all have
There is different influences.Look down upon greatly from mineral grain, in magmatic rock, metamorphic rock and sedimentary rock, granulitic texture is generally than non-equal grains
Structural strength is high；In granulitic texture, acinose texture is higher than coarse grained texture intensity.In porphyritic texture, particulate matrix compares glass
Matrix strength is high；The acid plutonic rock intensity that coarse grain has phenocryst is minimum；Particulate crystallite and without vitreous mafic effusive rock intensity
Highest.From structure conjoint, chemical rock in most of magmatic rock, metamorphic rock and sedimentary rock is tightly combined between crystal grain, by force
Degree is higher, but in chemical rock, and with soluble crystalline connection, intensity is higher, but water-resistance is poor.Consolidated clay rock some
It is recrystallization connection, intensity is more very different than other solid rocks.
Many scholars both domestic and external have carried out largely the macroscopic deformation and strength characteristics of rock and the relationship of rock texture
Research.
2013, Coggan et al. had studied the mechanical property of southwest England portion granite, it is believed that the increasing of kaolinization
Add is the main reason for causing this area's granite intensity to significantly reduce；Sousa et al. is also to Portugal different regions simultaneously
Granite has carried out comparative analysis, it is believed that quartzy textural characteristics have vital shadow to the mechanical behavior of rock in rock
It rings.
2014, Sajid and Arif et al. investigated the granite of the Pakistani northwestward, they have found rock interior
Lithofacies variation will will lead to granitello strength reduction, while the water absorption rate for also resulting in rock increases.
Advanced observation technologies some in recent years, such as sound emission, scanning electron microscope, x-ray tomography technology, number
Word image correlation technique etc. is introduced in rock mechanics experiment, and the quasi-brittle materials such as study of rocks are under load action, micro-crack
Nucleation, the physical process that extends and converge.Since the mineralogical composition of rock interior, crystal constitute sufficiently complex, mechanical property
It cannot still can be disclosed well by laboratory test.For the deep influence for understanding rock's microstructure to its mechanical property,
Through there is scholar to start to account for the numerical simulation of rock's microstructure, corresponding micromodel is established to reflect the structure of rock
Characteristic.
2014, Ghazvinian et al. proposed a kind of three-dimensional random Voronoi grain model for simulating brittleness rock
The crack propagation of stone inlays acquisition crystal model by Voronoi, and simulates intercrystalline crackle with crystal boundaries.
2017, Peng et al. proposed a kind of grain model and simulates the microcraking behavior of Wu Jizhi horse granite,
The microstructure of crystalline rock is simulated by construction polygon crystal grain.
But the above method can not reflect rock's microstructure completely, be only capable of considering the mineralogical composition and mineral of rock
Content cannot consider the irregular shape and arrangement architecture of rock forming mineral particle, it is therefore desirable to establish and consider rock texture feature
With the thin sight numerical model of mineral composition.
Summary of the invention
The present invention be in order to solve it is above-mentioned insufficient and design, and it is an object of the present invention to provide it is a kind of can reflect rock texture feature and
The thin sight numerical model generation method of mineral composition.Based on rock texture signature analysis as a result, such as rock forming mineral granule content, mine
The spatial position of object size distribution, crystallite dimension, grain shape and crystal grain, more really reflection rock's microstructure is macro to its
See the influence of mechanical property.
The present invention to achieve the goals above, uses following scheme.
A kind of thin sight numerical model generation method considering rock texture feature and mineral composition, which is characterized in that including
Following steps:
Step 1, by rock texture signature analysis, mineral grain distribution and particle geometric information are obtained:
Step 2, discrete element simulation；
It is converged being collected by the mineral grain distribution for obtaining rock texture signature analysis in step 1 and particle geometric information
Always, the mineral grain of irregular shape is reduced to the ball of a size reduction, using in discrete-time epidemic model rock sample
The process that grain diameter gradually expands, until granular system and true rock sample architectural characteristic, that is, grain diameter having the same
When whether reaching the distribution of mineral grain size, stop grain diameter expansion and discrete element numerical simulation, completes rock sample modeling；
Step 3, numerical simulation is generated；
Particle space positions and dimensions information in rock sample model is extracted, aggregates body is carried out based on Voronoi
The Spacial domain decomposition of figure, the mineral grains of the corresponding irregular shape of each Voronoi cellular；According to mines different in rock
The volume fraction of species type distributes mineral type to the mineral grains that each Voronoi cellular represents；
Step 4, finite element model is generated；
Obtained Voronoi polycrystalline structure is subjected to FEM meshing, and in mineral grain boundary and particle
Portion's insertion generates the finite element numerical model that can reflect rock texture and lithofacies structure without thickness boundary element, completes thin sight number
It is worth model to generate.
Preferably, in step 1, rock texture signature analysis method particularly includes:
Rock is sliced first, rock sample slice map is obtained by microexamination, then to rock sample slice map into
Row Digital Image Processing is partitioned into mineral grain using watershed algorithm to the image of binaryzation, identifies the distribution of mineral grain
With particle geometric information.
Preferably, the mineral grain distributed intelligence includes the mineral species of mineral grain, different minerals in step 1
Content, particle size, particle diameter distribution, the position of form center of particle, circularity and the area of type, particle geometric information includes mineral
The size of particle, solid degree, ratio of semi-minor axis length and circularity.
Preferably, the circularity is defined by following formula:
In formula one, SgrainIt is the mineral grain area obtained by rock texture signature analysis, PgrainIt is mineral grain
Perimeter, π are pi.
Preferably, in step 2, the process that is gradually expanded using discrete-time epidemic model rock sample endoparticle partial size
In, particle starts in given boundary according to identical rate volume expansion, in expansion process, the limitation of test sample boundary and
The constraint of condition is not invaded between particle mutually, permutatation occurs during the growth process for particle, tends to a closely knit and isotropism knot
Structure, the architectural characteristic of the above-mentioned aggregates body using discrete-time epidemic model rock sample can using radial distribution function and
Bond-orientational order parameter measurement judges whether grain diameter reaches the distribution of mineral grain size by the two parameters, if being unsatisfactory for
It is required that then continuing partial size expansion；Stop discrete element numerical simulation when demanded.
Preferably, the radial distribution function are as follows:
In formula two, dn (r) show apart from origin particle at r, and that can be found within the scope of width dr
Grain number, N are total mineral grain number in rock sample, and V is the volume of rock sample modeling, and the physical significance of g (r) is and origin
The particle random trajectory model of particle unit volume at r, to obtain the structural information of granular system, π is pi.
Preferably, the bond-orientational order parametric function are as follows:
In formula three,Spheric harmonic function, wherein l indicate symmetry S order parameter, m take-l≤m≤
L,The as direction vector of particle contact pair, andWithPolar angle and azimuth for direction vector, π are circumference
Rate, QlFor bond-orientational order parameter value.
Preferably, carrying out the Spacial domain decomposition based on Voronoi diagram in step 3 to aggregates body and referring to one
The technology for planting subregion in three-dimensional space, for D ∈ R3Region in, R refers to length dimension, R3It indicates in three-dimensional space
There are several seed points in given D, the coordinate of seed point is obtained by partial size plavini in D, for above-mentioned seed point i, is met
{Si(xi) (for i={ 1 ..., N }), each seed point can be assigned a Voronoi unit Ci, wherein Ci, definition
As shown in formula four:
In formula four, d (P, Sj) indicate Euclidean distance.
Preferably, boundary element of the insertion without thickness refers to: rock sample is being divided into finite element in step 4
After grid, the boundary element using conode link is created between the solid element of grid, boundary element itself without thickness,
The geometrical characteristic for not influencing rock sample, after the stress state of boundary element meets fracture criterion, boundary element failure and from
It is deleted in model, it can be with the Micro-fracture of simulation rock.
Advantageous effect of the invention is:
(1) the present invention provides it is a kind of consider rock texture feature and mineral composition thin sight numerical model generation method,
This method not only allows for the mineralogical composition and anisotropy of rock interior, additionally it is possible to consider rock forming mineral granule content, mineral
The spatial position of size distribution, crystallite dimension, grain shape and crystal grain, the microcosmic knot inside simulation rock that can be more accurate
Structure.
(2) numerical model of FEM calculation can be generated in method provided by the present invention, can be had using large scale business
It limits the analogue simulation that first software for calculation carries out mechanical behavior, such as uniaxial compression test, triaxial compression test, is uniaxially stretched examination
It tests, 3 cripping tests etc..
(3) method computational efficiency provided by the present invention is higher, and versatility is good.
Detailed description of the invention
Fig. 1 is the schematic diagram of numerical model generation method involved in the embodiment of the present invention；
Fig. 2 is that the part of rock texture signature analysis involved in the embodiment of the present invention shows that method is intended to；
In figure, 11 indicate that rock sample is sliced, 12 expression slice microexaminations, 13 expression mineral grain Boundary Recognitions, 14
Indicate particle geometric information extraction, 15 indicate that particle geometric information calculates；
Fig. 3 is to have obtained the signal of the mineral grain size regularity of distribution by rock texture signature analysis in the embodiment of the present invention
Figure；
Fig. 4 is the process gradually expanded in the embodiment of the present invention using discrete-time epidemic model rock sample endoparticle partial size
Schematic diagram；Wherein, 24 indicate that ball launches initial stage schematic diagram, 25 indicate that ball expands mid-term schematic diagram, and 26 indicate ball expansion
It completes, reaches dense distribution schematic diagram；
Fig. 5 is the radial distribution function schematic diagram for counting ball stacking system in the embodiment of the present invention in expansion process；
Fig. 6 is the bond-orientational order parameter schematic diagram for counting ball stacking system in the embodiment of the present invention in expansion process；
Fig. 7 is the Spacial domain decomposition method schematic diagram of Voronoi diagram involved in the embodiment of the present invention；
Wherein, 31 obtained ball spatial position and size distribution schematic diagram are indicated, 32 indicate to carry out based on Voronoi diagram
Spacial domain decomposition, generate numerical model schematic diagram；
Fig. 8 is generation finite element model method schematic diagram involved in the embodiment of the present invention；
Wherein, the numerical model schematic diagram generated in the 32 expression embodiment of the present invention, 41 indicate numerical model carrying out net
Lattice divide to obtain grid model schematic diagram, and 42 indicate the boundary element schematic diagram in the boundary of mineral grain insertion zero thickness；
Fig. 9 is algorithm flow chart in the embodiment of the present invention.
Specific embodiment
Referring to the drawings to a kind of thin sight numerical value for considering rock texture feature and mineral composition according to the present invention
Model generating method is elaborated.The part not elaborated in following embodiment belongs to the prior art.
<embodiment>
As shown in Figure 1, the thin sight numerical model according to the present invention for considering rock texture feature and mineral composition generates
Method can be divided into 4 processes: rock texture signature analysis 1, partial size plavini 2, the Spacial domain decomposition 3 based on Voronoi diagram
And generate finite element model 4.
As illustrated in fig. 1 and 2, rock texture signature analysis 1 is specific can include: rock sample slice 11, slice microexamination
12, mineral grain Boundary Recognition 13, particle geometric information extraction 14 and particle geometric information calculate 15.In the present embodiment, rock
Stone sample slice 11 takes diameter 5cm × height 10cm granite to be sliced the amplification for obtaining, being sliced under an optical microscope
Picture carries out Boundary Recognition using open source software ImageJ and obtains boundary profile, and it includes: mineral that particle geometric information, which calculates 15,
The size (as shown in formula five) of grain, the solid degree (shown in formula six) of mineral grain, the ratio of semi-minor axis length of mineral grain are (such as public
Shown in formula seven) and mineral grain circularity (as shown in expression formula eight).
In formula five, PgrainIndicate that the perimeter of mineral grain, π indicate that pi, R are the size for indicating mineral grain.
Formula six, SgrainIndicate the area of mineral grain, SconvexIt indicates to surround mineral with the smallest convex polygon of area
When grain, the area of convex polygon, Grain Solidity indicates the solid degree of mineral grain.
In formula seven, LMajorWhen indicating with ellipse fitting mineral grain boundary, the long axial length of the ellipse, LMinorAs
The short axle of the ellipse is long, and Aspect ratio indicates the ratio of semi-minor axis length of mineral grain.
In formula eight, SgrainIndicate the area of mineral grain, LMajorWhen indicating with ellipse fitting mineral grain boundary,
The long axial length of the ellipse, π indicate that pi, Grain Roundness indicate the circularity of mineral grain.
As shown in Fig. 1, Fig. 3, Fig. 4, Fig. 5 and Fig. 6, the distribution of mineral grain size has been obtained by rock texture signature analysis
After rule 21, ball 23 is launched according to the size regularity of distribution in cylindrical boundary 22 using open source software LIGGGHTS,
Random distribution rule is presented at dispensing initial stage 24 in ball, and subsequent ball starts with the expansion of identical rate, and expansion mid-term 25 is divided
The stacking that cloth rule shows as ball gradually starts closely knit, is finally reached a kind of more dense distribution 26, unites in expansion process
The radial distribution function 27 and bond-orientational order parameter 28 of ball stacking system are counted, whether the two parameters can measure whole system
The size distribution of true rock sample chats composition granule is met.In this example, observes and represent warp in radial distribution function 27
Cross the radial distribution function curve (blue curve) after partial size expansion peak value whether with the radial direction that represents true rock style
The first peak value of distribution function curve (black curve) coincide, and first peak value in radial distribution function indicates entire granular system
Order degree, first peak value coincide the aggregates body and true rock chats composition granule collection for indicating that numerical method generates
It matches in fit queueing discipline degree；The key represented after partial size expansion in bond-orientational order parameter 28 is observed simultaneously
Orientational order parameter number (black scatterplot) whether close to HCP (body-centered cubic arrangement architecture) and FCC (face-centered cubic arrangement architecture) friendship
Point.When above-mentioned parameter is met the requirements, partial size, which expands this process, to be terminated.
As illustrated in figures 1 and 7, it after obtained ball spatial position and size distribution 31 being extracted, carries out based on Voronoi
The Spacial domain decomposition of (Thiessen polygon) figure generates numerical model 32.
As shown in figs. 1 and 8, subsequent progress mechanical property analogue simulation for convenience, this method is by the numerical model of generation
32 carry out grid dividing, the grid model 41 that can be used for finite element software (such as ANSYS software) calculating are divided into, in order to guarantee
Wherein each mineral grain 45 can be separated into about 100 tetrahedron elements to computational accuracy, and on the boundary of mineral grain
It is inserted into the boundary element 42 of zero thickness, while in order to simulate the cracking situation inside mineral grain, also can inside mineral grain
It is inserted into boundary element 44.
It is sliced as shown in figure 9, this method passes through microexamination rock sample first, rock sample is sliced and carries out digitized map
As processing, mineral grain is partitioned into using watershed algorithm to the image of binaryzation, and then obtain the size of mineral grain, shape
The information such as state, mineral type；The mineral grain of irregular shape is reduced to the ball of a size reduction, using discrete element method
The process that simulation rock sample endoparticle partial size gradually expands does not invade the pact of condition mutually between the limitation of test sample boundary and particle
Permutatation occurs during the growth process for beam, particle, tends to a more closely knit and isotropic structure；Whether judge grain diameter
Reach the distribution of mineral grain size, continues partial size expansion if being unsatisfactory for requiring；Stop discrete element number when demanded
Value simulation, extracts the information such as particle space positions and dimensions, carries out the area of space based on Voronoi diagram to aggregates body and draws
Point, the mineral grains of the corresponding irregular shape of each Voronoi cellular；According to the volume of different minerals type in rock
Score distributes mineral type to the mineral grains that each Voronoi cellular represents；By obtained Voronoi polycrystalline structure into
Row FEM meshing, and without thickness boundary element, generation can reflect rock for insertion inside mineral grain boundary and particle
The finite element numerical model of structure and lithofacies structure.
Above embodiments are only the illustrations done to technical solution of the present invention.A kind of polycrystalline rock according to the present invention
The generation method of the microcosmic numerical model of stone is not merely defined in described process in the embodiment above, but is wanted with right
It asks subject to limited range.Any modification or benefit that those skilled in the art of the invention are made on the basis of the embodiment
It fills or equivalence replacement, all in scope of the present invention.
Claims (9)
1. it is a kind of consider rock texture feature and mineral composition thin sight numerical model generation method, which is characterized in that including with
Lower step:
Step 1, by rock texture signature analysis, mineral grain distribution and particle geometric information are obtained:
Step 2, discrete element simulation；
It will be summarized by collecting the mineral grain distribution that rock texture signature analysis obtains with particle geometric information in step 1, it will
The mineral grain of irregular shape is reduced to the ball of a size reduction, using discrete-time epidemic model rock sample endoparticle grain
The process that diameter gradually expands, until whether granular system and true rock sample architectural characteristic, that is, grain diameter having the same reach
To when the distribution of mineral grain size, stops grain diameter expansion and discrete element numerical simulation, complete rock sample modeling；
Step 3, numerical simulation is generated；
Particle space positions and dimensions information in rock sample model is extracted, aggregates body is carried out based on Voronoi diagram
Spacial domain decomposition, the mineral grains of the corresponding irregular shape of each Voronoi cellular；According to different minerals class in rock
The volume fraction of type distributes mineral type to the mineral grains that each Voronoi cellular represents；
Step 4, finite element model is generated；
Obtained Voronoi polycrystalline structure is subjected to FEM meshing, and is inserted inside mineral grain boundary and particle
Enter no thickness boundary element, generate the finite element numerical model that can reflect rock texture and lithofacies structure, completes thin sight Numerical-Mode
Type generates.
2. thin sight numerical model generation method according to claim 1, it is characterised in that: in step 1, rock texture feature
Analysis method particularly includes:
Rock is sliced first, rock sample slice map is obtained by microexamination, then rock sample slice map is counted
Word image procossing is partitioned into mineral grain using watershed algorithm to the image of binaryzation, identify mineral grain distribution and
Grain geological information.
3. thin sight numerical model generation method according to claim 2, it is characterised in that: in step 1, the mineral grain
Distributed intelligence includes the mineral species of mineral grain, the content of different minerals type, particle size, particle diameter distribution, the shape of particle
Heart position, circularity and area, particle geometric information include size, solid degree, ratio of semi-minor axis length and the circularity of mineral grain.
4. thin sight numerical model generation method according to claim 3, it is characterised in that: the circularity passes through following formula
Definition:
In formula one, SgrainIt is the mineral grain area obtained by rock texture signature analysis, PgrainIt is mineral grain perimeter,
π is pi.
5. thin sight numerical model generation method according to claim 1, it is characterised in that: in step 2, using discrete unit
During method simulation rock sample endoparticle partial size gradually expands, particle starts in given boundary according to identical rate
Volume expansion does not invade the constraint of condition mutually, particle is in growth course in expansion process between the limitation of test sample boundary and particle
Middle generation permutatation tends to a closely knit and isotropic structure, the above-mentioned particle using discrete-time epidemic model rock sample
The architectural characteristic of aggregate can be measured using radial distribution function and bond-orientational order parameter, judge particle by the two parameters
Whether partial size reaches the distribution of mineral grain size, continues partial size expansion if being unsatisfactory for requiring；Stop when demanded
Discrete element numerical simulation.
6. thin sight numerical model generation method according to claim 5, it is characterised in that: the radial distribution function are as follows:
In formula two, dn (r) show apart from origin particle at r, and the granule number that can be found within the scope of width dr,
N is that total mineral grain number, V are the volume of rock sample modeling in rock sample, and the physical significance of g (r) is and origin particle phase
Particle random trajectory model away from unit volume at r, to obtain the structural information of granular system, π is pi.
7. thin sight numerical model generation method according to claim 5, it is characterised in that: the bond-orientational order parametric function
Are as follows:
In formula three,It is spheric harmonic function, wherein l indicates that the S order parameter of symmetry, m take-l≤m≤l,I.e.
For particle contact pair direction vector, andWithPolar angle and azimuth for direction vector, π are pi, QlFor
Bond-orientational order parameter value.
8. thin sight numerical model generation method according to claim 1, it is characterised in that: in step 3, to aggregates body
The technology that the Spacial domain decomposition based on Voronoi diagram refers to a kind of subregion in three-dimensional space is carried out, for D ∈ R3Area
In domain, R refers to length dimension, R3Indicate there are several seed points in the D given in three-dimensional space, the coordinate of seed point in D
It is obtained by partial size plavini, for above-mentioned seed point i, meets { Si(xi) (for i={ 1 ..., N }), each seed point
A Voronoi unit C can be assignedi, wherein Ci, it defines as shown in formula four:
In formula four, d (P, Sj) indicate Euclidean distance.
9. thin sight numerical model generation method according to claim 1, it is characterised in that: in step 4, be inserted into without thickness
Boundary element refers to: after rock sample is divided into finite element grid, creating between the solid element of grid using altogether
The boundary element of node link, boundary element itself do not influence the geometrical characteristic of rock sample, when boundary element without thickness
After stress state meets fracture criterion, boundary element fails and deletes from model, can be with the Micro-fracture of simulation rock.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910298582.9A CN110069844B (en) | 2019-04-15 | 2019-04-15 | Mesoscopic model generation method considering rock structure characteristics and mineral composition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910298582.9A CN110069844B (en) | 2019-04-15 | 2019-04-15 | Mesoscopic model generation method considering rock structure characteristics and mineral composition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110069844A true CN110069844A (en) | 2019-07-30 |
| CN110069844B CN110069844B (en) | 2020-10-30 |
Family
ID=67367579
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910298582.9A Active CN110069844B (en) | 2019-04-15 | 2019-04-15 | Mesoscopic model generation method considering rock structure characteristics and mineral composition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110069844B (en) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110765572A (en) * | 2019-09-12 | 2020-02-07 | 中国科学院武汉岩土力学研究所 | Continuous discontinuous numerical simulation method for single triaxial test of almond-shaped basalt |
| CN110990911A (en) * | 2019-10-29 | 2020-04-10 | 武汉大学 | Rock crystal structure rapid modeling method based on particle cluster substitution |
| CN111159927A (en) * | 2019-11-24 | 2020-05-15 | 浙江大学 | Numerical modeling method for three-dimensional irregular-shaped particle throwing based on voxel matrix |
| CN111307669A (en) * | 2020-02-29 | 2020-06-19 | 上海健康医学院 | Method for measuring particle local structure in sparse two-phase flow |
| CN111353237A (en) * | 2020-03-23 | 2020-06-30 | 河海大学 | Anisotropic rock modeling method based on oriented development of mineral grains |
| CN111428359A (en) * | 2020-03-23 | 2020-07-17 | 河海大学 | Anisotropic rock modeling method considering grain occlusion |
| CN111537699A (en) * | 2020-05-29 | 2020-08-14 | 煤炭科学技术研究院有限公司 | Method for distinguishing hard rock stratum of coal mine through multiple indexes |
| CN111583358A (en) * | 2020-03-21 | 2020-08-25 | 长沙迪迈数码科技股份有限公司 | Method and device for automatically delineating mine rock boundary of surface mine and storage medium |
| CN112150430A (en) * | 2020-09-21 | 2020-12-29 | 中国矿业大学(北京) | Numerical analysis method utilizing rock microscopic structure digital image |
| CN112257326A (en) * | 2020-10-27 | 2021-01-22 | 中国矿业大学(北京) | Numerical simulation method and system for stress state of surrounding rock |
| CN112287563A (en) * | 2020-11-18 | 2021-01-29 | 华能澜沧江水电股份有限公司 | Method and system for detecting microstructure of crystalline phase precursor of monodisperse particle system |
| CN112326524A (en) * | 2020-10-22 | 2021-02-05 | 中国石油大学(华东) | Rock pore permeability measurement method based on CT scanning image |
| CN112444534A (en) * | 2020-12-11 | 2021-03-05 | 中国地质大学（北京） | Method for rapidly measuring and calculating fire passing temperature of stone cultural relics |
| CN112507581A (en) * | 2020-11-18 | 2021-03-16 | 华能澜沧江水电股份有限公司 | Rock mesoscopic numerical model generation method and system capable of avoiding malformation computing unit |
| CN113866389A (en) * | 2021-09-27 | 2021-12-31 | 中国电建集团成都勘测设计研究院有限公司 | Method for generating mineral crystal model containing bedding rock |
| CN115078362A (en) * | 2022-07-01 | 2022-09-20 | 西南石油大学 | Drillability characterization method based on rock mesoscopic structure |
| CN115618700A (en) * | 2022-09-26 | 2023-01-17 | 长安大学 | Mineral identification technology-based rock modeling method and system for containing physical surface |
| CN115937294A (en) * | 2022-12-30 | 2023-04-07 | 重庆大学 | Method for predicting height of goaf collapse zone after ground fracturing of coal mine |
| CN116011302A (en) * | 2023-02-16 | 2023-04-25 | 广西大学 | Continuous-discontinuous numerical model building method considering real microscopic characteristics of rock |
| CN116467901A (en) * | 2023-03-02 | 2023-07-21 | 中国地质大学（北京） | Discrete element modeling method and device for broken stone |
| CN118278195A (en) * | 2024-04-07 | 2024-07-02 | 长江大学 | Soil shearing and slope simulation method for block stone polygonal and random block stone structure |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2686714A1 (en) * | 2007-06-01 | 2008-12-11 | Exxonmobil Upstream Research Company | Generation of constrained voronoi grid in a plane |
| CN104036538A (en) * | 2014-06-25 | 2014-09-10 | 清华大学 | Method and system for reconstruction and analysis of earth-rock aggregate three-dimensional microstructure |
| CN106294984A (en) * | 2016-08-08 | 2017-01-04 | 中国科学院武汉岩土力学研究所 | A kind of method determining Rock Damage threshold based on micro-mechanical model rate of increase |
| CN108170959A (en) * | 2017-12-28 | 2018-06-15 | 天地科技股份有限公司 | Mechanical response of the rock mass numerical analysis method and device based on discrete element |
| CN108629147A (en) * | 2018-05-02 | 2018-10-09 | 西安电子科技大学 | A kind of polycrystal Geometric Modeling Method |
-
2019
- 2019-04-15 CN CN201910298582.9A patent/CN110069844B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2686714A1 (en) * | 2007-06-01 | 2008-12-11 | Exxonmobil Upstream Research Company | Generation of constrained voronoi grid in a plane |
| CN104036538A (en) * | 2014-06-25 | 2014-09-10 | 清华大学 | Method and system for reconstruction and analysis of earth-rock aggregate three-dimensional microstructure |
| CN106294984A (en) * | 2016-08-08 | 2017-01-04 | 中国科学院武汉岩土力学研究所 | A kind of method determining Rock Damage threshold based on micro-mechanical model rate of increase |
| CN108170959A (en) * | 2017-12-28 | 2018-06-15 | 天地科技股份有限公司 | Mechanical response of the rock mass numerical analysis method and device based on discrete element |
| CN108629147A (en) * | 2018-05-02 | 2018-10-09 | 西安电子科技大学 | A kind of polycrystal Geometric Modeling Method |
Non-Patent Citations (6)
| Title |
|---|
| XIANG JI;WEI ZHOU;YUAN CHEN;GANG MA: "Generation of the Polycrystalline Rock Microstructure by a Novel Voronoi Grain-Based Model With Particle Growth", 《52ND U.S. ROCK MECHANICS/GEOMECHANICS SYMPOSIUM》 * |
| ZHOU W , LIU J , MA G , ET AL.: "Macroscopic and microscopic behaviors of granular materials under proportional strain path: a DEM study", 《INTERNATIONAL JOURNAL FOR NUMERICAL & ANALYTICAL METHODS IN GEOMECHANICS》 * |
| 徐琨，周伟，马刚，常晓林等: "基于离散元法的颗粒破碎模拟研究进展", 《岩土工程学报》 * |
| 朱亮，陈明等: "轮廓爆破下柱状节理岩体开裂过程的数值模拟", 《爆炸与冲击》 * |
| 马刚，周伟，常晓林，周创兵: "考虑颗粒破碎的堆石体三维随机多面体细观数值模拟", 《岩石力学与工程学报》 * |
| 马刚，周伟，常晓林等: "堆石料缩尺效应的细观机制研究", 《岩石力学与工程学报》 * |
Cited By (31)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110765572A (en) * | 2019-09-12 | 2020-02-07 | 中国科学院武汉岩土力学研究所 | Continuous discontinuous numerical simulation method for single triaxial test of almond-shaped basalt |
| CN110765572B (en) * | 2019-09-12 | 2021-04-09 | 中国科学院武汉岩土力学研究所 | Continuous discontinuous numerical simulation method for single triaxial test of almond-shaped basalt |
| CN110990911A (en) * | 2019-10-29 | 2020-04-10 | 武汉大学 | Rock crystal structure rapid modeling method based on particle cluster substitution |
| CN111159927A (en) * | 2019-11-24 | 2020-05-15 | 浙江大学 | Numerical modeling method for three-dimensional irregular-shaped particle throwing based on voxel matrix |
| CN111307669A (en) * | 2020-02-29 | 2020-06-19 | 上海健康医学院 | Method for measuring particle local structure in sparse two-phase flow |
| CN111307669B (en) * | 2020-02-29 | 2022-07-05 | 上海健康医学院 | Method for measuring particle local structure in sparse two-phase flow |
| CN111583358A (en) * | 2020-03-21 | 2020-08-25 | 长沙迪迈数码科技股份有限公司 | Method and device for automatically delineating mine rock boundary of surface mine and storage medium |
| CN111583358B (en) * | 2020-03-21 | 2024-02-20 | 长沙迪迈科技股份有限公司 | Automatic surface mine rock boundary line encircling method, device and storage medium |
| CN111353237A (en) * | 2020-03-23 | 2020-06-30 | 河海大学 | Anisotropic rock modeling method based on oriented development of mineral grains |
| CN111428359A (en) * | 2020-03-23 | 2020-07-17 | 河海大学 | Anisotropic rock modeling method considering grain occlusion |
| CN111428359B (en) * | 2020-03-23 | 2022-03-08 | 河海大学 | Anisotropic rock modeling method considering grain occlusion |
| CN111353237B (en) * | 2020-03-23 | 2022-03-08 | 河海大学 | Anisotropic rock modeling method based on oriented development of mineral grains |
| CN111537699A (en) * | 2020-05-29 | 2020-08-14 | 煤炭科学技术研究院有限公司 | Method for distinguishing hard rock stratum of coal mine through multiple indexes |
| CN112150430A (en) * | 2020-09-21 | 2020-12-29 | 中国矿业大学(北京) | Numerical analysis method utilizing rock microscopic structure digital image |
| CN112150430B (en) * | 2020-09-21 | 2023-10-10 | 中国矿业大学(北京) | Numerical analysis method for digital image by using rock microstructure |
| CN112326524A (en) * | 2020-10-22 | 2021-02-05 | 中国石油大学(华东) | Rock pore permeability measurement method based on CT scanning image |
| CN112257326A (en) * | 2020-10-27 | 2021-01-22 | 中国矿业大学(北京) | Numerical simulation method and system for stress state of surrounding rock |
| CN112257326B (en) * | 2020-10-27 | 2023-05-30 | 中国矿业大学(北京) | Numerical simulation method and system for stress state of surrounding rock |
| CN112287563B (en) * | 2020-11-18 | 2022-08-16 | 华能澜沧江水电股份有限公司 | Method and system for detecting microstructure of crystalline phase precursor of monodisperse particle system |
| CN112507581A (en) * | 2020-11-18 | 2021-03-16 | 华能澜沧江水电股份有限公司 | Rock mesoscopic numerical model generation method and system capable of avoiding malformation computing unit |
| CN112287563A (en) * | 2020-11-18 | 2021-01-29 | 华能澜沧江水电股份有限公司 | Method and system for detecting microstructure of crystalline phase precursor of monodisperse particle system |
| CN112444534A (en) * | 2020-12-11 | 2021-03-05 | 中国地质大学（北京） | Method for rapidly measuring and calculating fire passing temperature of stone cultural relics |
| CN113866389A (en) * | 2021-09-27 | 2021-12-31 | 中国电建集团成都勘测设计研究院有限公司 | Method for generating mineral crystal model containing bedding rock |
| CN115078362A (en) * | 2022-07-01 | 2022-09-20 | 西南石油大学 | Drillability characterization method based on rock mesoscopic structure |
| CN115618700A (en) * | 2022-09-26 | 2023-01-17 | 长安大学 | Mineral identification technology-based rock modeling method and system for containing physical surface |
| CN115937294A (en) * | 2022-12-30 | 2023-04-07 | 重庆大学 | Method for predicting height of goaf collapse zone after ground fracturing of coal mine |
| CN115937294B (en) * | 2022-12-30 | 2023-10-27 | 重庆大学 | Prediction method for goaf caving zone height after coal mine ground fracturing |
| CN116011302A (en) * | 2023-02-16 | 2023-04-25 | 广西大学 | Continuous-discontinuous numerical model building method considering real microscopic characteristics of rock |
| CN116011302B (en) * | 2023-02-16 | 2024-06-25 | 广西大学 | Continuous-discontinuous numerical model building method considering real microscopic characteristics of rock |
| CN116467901A (en) * | 2023-03-02 | 2023-07-21 | 中国地质大学（北京） | Discrete element modeling method and device for broken stone |
| CN118278195A (en) * | 2024-04-07 | 2024-07-02 | 长江大学 | Soil shearing and slope simulation method for block stone polygonal and random block stone structure |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110069844B (en) | 2020-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110069844A (en) | A kind of thin sight numerical model generation method considering rock texture feature and mineral composition | |
| Kong et al. | Microstructure characteristics and fractal analysis of 3D-printed sandstone using micro-CT and SEM-EDS | |
| CN105487121B (en) | Multiple dimensioned digital cores method is built based on CT scan image and Electrical imaging image co-registration | |
| Steiner et al. | The developmental cycles of early Cambrian Olivooidae fam. nov.(? Cycloneuralia) from the Yangtze Platform (China) | |
| CN107292919B (en) | A kind of multiple dimensioned settlement place matching process of space-oriented subdivision | |
| Latypov et al. | Application of chord length distributions and principal component analysis for quantification and representation of diverse polycrystalline microstructures | |
| CN107844757A (en) | A kind of method using river width in digital elevation model extraction basin | |
| CN112507581B (en) | Rock mesoscopic numerical model generation method and system capable of avoiding malformation computing unit | |
| CN109444973B (en) | Gravity forward modeling acceleration method under spherical coordinate system | |
| CN115018986A (en) | Strong-structure active region rock microstructure information interpretation and three-dimensional modeling method | |
| Rollett et al. | Modeling polycrystalline microstructures in 3D | |
| CN113223631B (en) | Crystal structure analysis method and system | |
| CN113743027B (en) | Method and device for drawing wind resource map based on CFD technology | |
| CN107240103A (en) | Boundary processing method in a kind of digital volume related algorithm split based on image | |
| CN108875260B (en) | Thermal spraying coating performance parameter calculation method | |
| CN111291497A (en) | Improved digital core modeling method | |
| Groeber et al. | Developing a robust 3-D characterization-representation framework for modeling polycrystalline materials | |
| CN114528729B (en) | Method for predicting gas reservoir yield of buried hill cracks based on multi-scale coupling | |
| CN114544286A (en) | Particle random generation method based on particle size distribution | |
| CN113514899A (en) | Self-adaptive subdivision inversion method for gravity gradient | |
| Zeng et al. | When did the large-scale extensional tectonics begin in North China Craton? | |
| CN108827987B (en) | Method for evaluating distribution uniformity of spatial positions of bubbles in gas-containing soil sample | |
| CN111415346A (en) | Particle model damage characterization method based on gridding | |
| Helminiak et al. | Uniaxial wave propagation through copper mountain sandstone | |
| Wang et al. | Peridynamic Simulation of Fracturing in Heterogeneous Rock Based on X-Ray Diffraction and Scanning Electron Microscope Tests |
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 |