CN107045580A - A kind of shale mechanics parameter quick calculation method based on digital cores - Google Patents
A kind of shale mechanics parameter quick calculation method based on digital cores Download PDFInfo
- Publication number
- CN107045580A CN107045580A CN201710288425.0A CN201710288425A CN107045580A CN 107045580 A CN107045580 A CN 107045580A CN 201710288425 A CN201710288425 A CN 201710288425A CN 107045580 A CN107045580 A CN 107045580A
- Authority
- CN
- China
- Prior art keywords
- mrow
- shale
- msub
- formula
- rock core
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 29
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 84
- 239000011707 mineral Substances 0.000 claims abstract description 84
- 239000011435 rock Substances 0.000 claims abstract description 80
- 238000005381 potential energy Methods 0.000 claims abstract description 13
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 35
- 238000010008 shearing Methods 0.000 claims description 16
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 10
- 239000010453 quartz Substances 0.000 claims description 10
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 10
- 229910021532 Calcite Inorganic materials 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 6
- 229910000514 dolomite Inorganic materials 0.000 claims description 6
- 239000010459 dolomite Substances 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 239000004927 clay Substances 0.000 claims description 4
- 239000011028 pyrite Substances 0.000 claims description 4
- 229910052683 pyrite Inorganic materials 0.000 claims description 4
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 4
- 238000002083 X-ray spectrum Methods 0.000 claims description 3
- DLHONNLASJQAHX-UHFFFAOYSA-N aluminum;potassium;oxygen(2-);silicon(4+) Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[Al+3].[Si+4].[Si+4].[Si+4].[K+] DLHONNLASJQAHX-UHFFFAOYSA-N 0.000 claims description 3
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 claims description 3
- 229910052786 argon Inorganic materials 0.000 claims description 3
- -1 argon ion Chemical class 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000002939 conjugate gradient method Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000001000 micrograph Methods 0.000 claims description 3
- 229910052655 plagioclase feldspar Inorganic materials 0.000 claims description 3
- 238000005498 polishing Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 229910052570 clay Inorganic materials 0.000 claims description 2
- 239000000203 mixture Substances 0.000 abstract description 3
- 239000000470 constituent Substances 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000004876 x-ray fluorescence Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910001748 carbonate mineral Inorganic materials 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000724 energy-dispersive X-ray spectrum Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The present invention relates to a kind of shale mechanics parameter quick calculation method based on digital cores, including:(1) shale minerals component 3-dimensional digital rock core is built;A, scanning shale structure, obtain shale minerals structure chart, calculate brittle mineral percentage composition;B, structure shale minerals component 3-dimensional digital rock core;(2) it is based on shale minerals component 3-dimensional digital rock core FEM parallel computation;C, Region Decomposition, discretization are carried out successively to 3-dimensional digital rock core, obtain some voxels;D, one macro-strain of application so that the 3-dimensional digital rock core elastic potential energy is minimum, asks for being tried to achieve the rock mechanics parameters of overall shale as mean stress caused by the macro-strain.The present invention considers influence of the shale minerals component to shale mechanics parameter, and result of calculation is more accurate.Region Decomposition is carried out to digital cores, finite element modelling is carried out to 3-dimensional digital rock core based on parallel computation, digital cores information is not lost, rock mechanics parameters are accurately calculated, the calculating time is saved.
Description
Technical field
The present invention relates to a kind of shale mechanics parameter quick calculation method based on digital cores, belong to shale oil-gas reservoir and survey
Visit development technique field.
Background technology
The acquisition of shale mechanics parameter plays an important roll for horizontal well drilling and hydraulic fracture mining, and research shows bullet
Property modulus has certain correlation with shale fragility, therefore Accurate Determining shale modulus of elasticity also has for look for oil and gas dessert
Play an important role.Actual mining site generally determines shale mechanics parameter using Rock Mechanics Test, high yet with shale fragility, right
Shale carries out Rock Mechanics Test difficulty greatly, and has the shortcomings that the costly, time is long.
With the development of CT and scanning electron-microscopy, the research to shale microstructure achieves great progress.Pass through
Binary image segmentation method, is converted into digital cores image to characterize the binary map of interstitial space and grain structure, the method pair
In traditional hypertonic sandrock effect preferably, but for fine and close shale, mineral constituent information is complicated, and bianry image can not be accurate
Ground describes shale mechanical characteristics, and being directly used in measure shale elastic parameter has larger error.Therefore, shale minerals component is set up
3-dimensional digital rock core is significant for calculating its rock mechanics parameters.
The content of the invention
In view of the shortcomings of the prior art, quickly calculated the invention provides a kind of shale mechanics parameter based on digital cores
Method;
Shale minerals constitutional diagram can be set up according to x-ray fluorescence analysis (or EDS energy spectrum analysis), digital rock is then based on
Heart construction method, can build the 3-dimensional digital rock core for considering mineral constituent.It is each when carrying out finite element modelling to digital cores
Individual pixel is expressed as a three linear finite elements, in order to calculate the mechanics parameter of shale exactly, it is necessary to take what structure was obtained
The larger overlay area of 3-dimensional digital rock core, but this will cause the finite element free degree rapidly to increase, it is therefore necessary to use area
Domain decomposition method is based on parallel computation and carries out finite element modelling to 3-dimensional digital rock core, finally asks for overall shale rock mechanics ginseng
Number.
Not only calculating speed is fast for the method for the invention, but also spends few.
Term is explained
1st, shale mechanics parameter, refers to the rock mechanics parameters of shale, mainly including Young's modulus, bulk modulus, shearing
The parameters such as modulus, Poisson's ratio, tensile strength and compression strength.
2nd, brittle mineral, refers to be broken or mineralogical composition few or without plastic deformation is shown before destroying, shale is crisp
Property mineral mainly include quartz and carbonate mineral.
3rd, voxel, is volume element (volumepixel) abbreviation, the solid comprising voxel can be rendered by solid or
The polygon contour surface that person extracts given threshold value profile is showed.
4th, shale volume modulus, refers to the ratio of shale mean stress and body strain.
5th, shale modulus of shearing, refers to shale shear stress and the ratio of shear strain.
6th, shale Young's modulus, refers to the ratio between Rock Under Uniaxial Compression stress and uniaxial strain.
7th, shale Poisson's ratio, refers to rock when unidirectional tension or in compression, transverse normal strain is absolute with axial direction normal strain
The ratio of value.
The technical scheme is that:
A kind of shale mechanics parameter quick calculation method based on digital cores, including:
(1) shale minerals component 3-dimensional digital rock core is built
A, scanning shale structure, analysis obtains shale minerals structure chart, and calculates the volumn concentration of brittle mineral;
B, based on shale minerals structure chart, shale minerals component 3-dimensional digital rock is built using digital cores construction method
The heart;Each voxel of shale minerals component 3-dimensional digital rock core is made up of shale minerals component;
(2) FEM parallel computation based on shale minerals component 3-dimensional digital rock core
C, Region Decomposition, discretization are carried out successively to shale minerals component 3-dimensional digital rock core, obtain several three-dimensional images
Element, each voxel is considered as three linear finites of 8 nodes;
D, to step C processing after shale minerals component 3-dimensional digital rock core apply a macro-strain so that the three-dimensional
Digital cores elastic potential energy is minimum, is asked for based on parallel computation using quick conjugate gradient method as average caused by the macro-strain
Stress, tries to achieve the rock mechanics parameters of overall shale, and rock mechanics parameters include shale volume modulus, shale modulus of shearing, page
Rock Young's modulus and shale Poisson's ratio.Basic data is provided for shale pressure break.
According to currently preferred, the step A, shale structure is scanned by CT or FIB-SEM.
According to currently preferred, the step A, shale structure is scanned, analysis obtains shale minerals structure chart, including:
A, shale samples surface is polished successively, argon ion polishing;For reducing sample surfaces roughness, improve into
As quality;
B, shale samples are scanned using CT or FIB-SEM, carry out electronic imaging;
C, collect the X-ray spectrum that produces in scanning shale samples formation Electron Microscope images, and by itself and mineral
Matter characteristic spectrum is compared, and obtains the corresponding mineral matter of each pixel, builds shale minerals structure chart;
D, the volumn concentration for calculating brittle mineral.By counting brittle mineral (quartz and carbonate) voxel number,
Divided by total voxel number, produce the volumn concentration of brittle mineral.
According to currently preferred, the shale minerals component includes quartz, dolomite, calcite, potassium feldspar, plagioclase
Stone, pyrite and clay, the brittle mineral include quartz, dolomite and calcite.
According to currently preferred, the step B, the digital cores construction method includes process simulation method, Ma Er can
The random reconstruction method of husband.
Process simulation method sets up digital rock by the geology diagenetic process (including deposition, compacting and diagenesis) of simulation rock
The heart, not only allows for the particle size distribution of rock, but also by some other rock physicses got by thin section analysis
Matter is incorporated, and can set up anisotropic digital cores;The digital cores that the random reconstruction method of markov is set up have good
The pore communication and modeling speed got well quickly, therefore mainly carry out the structure of digital cores from this two methods.
According to currently preferred, the step C, region point is carried out successively to shale minerals component 3-dimensional digital rock core
Solution, discretization, including:
E, shale minerals component 3-dimensional digital rock core divided equally according to size, obtain many sub-regions, each sub-district
Domain interface meets displacement continuously, shown in the weak form such as formula (I) per sub-regions:
∫Γ(u (l)-v (l)) f (l) dl=0 (I)
In formula (I), l is coordinate variable, and u (l) and v (l) represent the displacement in adjacent subarea domain, and f (l) is finite element shape letter
Number;
F, the every sub-regions obtained to step e carry out discretization, according to pixel cell that digital cores are directly discrete, from
Each voxel obtained after dispersion regards three linear finites of 8 nodes as, and the unit is square, by 8
Individual node composition, finite element shape function represents that X, Y, Z are three directions of digital cores coordinate system respectively using linear function,
Each three linear finite is occupied by a kind of mineral matter, and the modulus of elasticity and Poisson's ratio of various mineral matters are by existing experimental data
Provide.
According to currently preferred, the step D, including step are as follows:
G, 3-dimensional digital rock core elastic potential energy EnExpression formula such as formula (II) shown in:
In formula (II), p, q, rs be digital cores coordinate system x, y, tri- directions of z, d is integral sign, εpqAnd εrsPoint
Not Wei pq directions and rs directions dependent variable, CpqrsFor modulus of elasticity tensor;
Due to symmetry, strain tensor ε includes 6 independent variable (εxx、εyy、εzz、εxy、εxz、εyz), εxx、εyy、εzz、
εxy、εxz、εyzRefer respectively to the dependent variable of all directions, CpqrsIt is expressed as Cαβ, α, β is x, y, tri- directions of z, such as formula (III), formula
(IV) shown in:
In formula (III), formula (IV), upRefer to the displacement in p directions, xpRefer to the coordinate in p directions, uqRefer to the displacement in q directions,
xqRefer to the coordinate in q directions;
CpqrsIt is expressed as Cαβ, εpqIt is expressed as εα, εrsIt is expressed as εβ, the elastic potential energy E of 3-dimensional digital rock corenExpression formula formula
(II) it is converted into formula (V):
Make 3-dimensional digital rock core elastic potential energy minimum, it is necessary to meet formula (VI):
In formula (VI), umIt is displacement, m refers to all nodes and all directions;
Formula (VI) is solved using Fast Field method, finite element is carried out to all subregion based on MPI concurrent techniques
Calculate, calculate the stress and strain for obtaining all subregion, the mean stress peace of overall digital rock core is tried to achieve by arithmetic average
Strain;
K, based on 3-dimensional digital rock core mean strain stress, calculating obtains shale volume modulus, shale modulus of shearing, shale
Young's modulus and shale Poisson's ratio:
Shown in the calculation formula of shale volume modulus such as formula (VII):
In formula (VII), σ0For mean stress, σxx、σyy、σzzThe respectively direct stress in x, y, z direction, K is shale volume mould
Amount, εxx、εyy、εzzThe respectively normal strain in x, y, z direction;
Shown in the calculation formula of shale modulus of shearing such as formula (VIII):
σij=2 μ εij (Ⅷ)
In formula (VIII), σijRefer to i, the direct stress in j directions, i, j is respectively x, y direction, and μ is shale modulus of shearing, εijFor
Shearing stress;
Shown in the calculation formula of shale Young's modulus such as formula (Ⅸ):
In formula (Ⅸ), E is shale Young's modulus;
Shown in the calculation formula of shale Poisson's ratio such as formula (Ⅹ):
In formula (Ⅹ), ν is shale Poisson's ratio.
Beneficial effects of the present invention are:
1. the present invention can build shale minerals structure chart exactly using EDS/X ray fluorescence analysis, so as to set up page
Rock mineral constituent 3-dimensional digital rock core, the present invention considers influence of the shale minerals component to shale mechanics parameter, therefore calculates knot
Fruit is more accurate.
2. by carrying out Region Decomposition to digital cores, finite element mould is carried out to 3-dimensional digital rock core based on parallel computation
Intend, due to taking the entirety or large area of 3-dimensional digital rock core, digital cores information is not lost, therefore, it is possible to accurately calculate its rock
Stone mechanics parameter, in addition parallel computation can save many calculating times.
3rd, the present invention has the advantages that speed is fast, expense is low, and the exploration and development for shale gas is significant.
Brief description of the drawings
Fig. 1 is a kind of shale mechanics parameter quick calculation method flow chart based on digital cores of the present invention;
Fig. 2 is embodiment shale core CT scan image schematic diagram;
Fig. 3 is embodiment 3-dimensional digital rock core Region Decomposition schematic diagram;
Fig. 4 is embodiment pixel cell node serial number sequential schematic.
Embodiment
The present invention is further qualified with reference to Figure of description and embodiment, but not limited to this.
Embodiment
A kind of shale mechanics parameter quick calculation method based on digital cores, as shown in figure 1, including:
(1) shale minerals component 3-dimensional digital rock core is built
A, shale structure, obtained wherein piece image scanned as shown in Fig. 2 analysis obtains page by CT or FIB-SEM
Rock mineral structure figure, and calculate the volumn concentration of brittle mineral;Shale minerals component include quartz, dolomite, calcite,
Potassium feldspar, plagioclase, pyrite and clay, brittle mineral include quartz, dolomite and calcite.
B, based on shale minerals structure chart, shale minerals group is built using process simulation method or the random reconstruction method of markov
Divide 3-dimensional digital rock core;Each voxel of shale minerals component 3-dimensional digital rock core is made up of shale minerals component;Process simulation method is led to
Cross the geology diagenetic process (including deposition, compacting and diagenesis) of simulation rock to set up digital cores, not only allow for rock
Particle size distribution, but also some other petrophysical property got by thin section analysis is incorporated, it can set up
Anisotropic digital cores;The digital cores that the random reconstruction method of markov is set up have good pore communication and modeling
Speed quickly, therefore mainly carries out the structure of digital cores from this two methods.
(2) FEM parallel computation based on shale minerals component 3-dimensional digital rock core
C, Region Decomposition, discretization are carried out successively to shale minerals component 3-dimensional digital rock core, obtain several three-dimensional images
Element, each voxel is considered as three linear finites of 8 nodes;
D, to step C processing after shale minerals component 3-dimensional digital rock core apply a macro-strain so that the three-dimensional
Digital cores elastic potential energy is minimum, is asked for based on parallel computation using quick conjugate gradient method as average caused by the macro-strain
Stress, tries to achieve the rock mechanics parameters of overall shale, and rock mechanics parameters include shale volume modulus, shale modulus of shearing, page
Rock Young's modulus and shale Poisson's ratio.Basic data is provided for shale pressure break.
Step A, scans shale structure, and analysis obtains shale minerals structure chart, including:
A, shale samples surface is polished successively, argon ion polishing;For reducing sample surfaces roughness, improve into
As quality;
B, shale samples are scanned using CT or FIB-SEM, carry out electronic imaging;
C, collect the EDS/X ray spectrums that produce in scanning shale samples formation Electron Microscope images, and by its with
Mineral matter characteristic spectrum is compared, and obtains the corresponding mineral matter of each pixel, builds shale minerals structure chart;
D, the volumn concentration for calculating brittle mineral.By counting brittle mineral (quartz and carbonate) voxel number,
Divided by total voxel number, produce the volumn concentration of brittle mineral.
Step C, Region Decomposition, discretization are carried out to shale minerals component 3-dimensional digital rock core successively, including:
E, shale minerals component 3-dimensional digital rock core divided equally according to size, many sub-regions are obtained, such as Fig. 3 institutes
Show, each sub-region interface meets displacement continuously, shown in the weak form such as formula (I) per sub-regions:
∫Γ(u (l)-v (l)) f (l) dl=0 (I)
In formula (I), l is coordinate variable, and u (l) and v (l) represent the displacement in adjacent subarea domain, and f (l) is finite element shape letter
Number;
F, the every sub-regions obtained to step e carry out discretization, according to pixel cell that digital cores are directly discrete, from
Each voxel obtained after dispersion regards three linear finites of 8 nodes as, and pixel cell node serial number is suitable
Sequence figure is as shown in figure 4, for any pixel unit of digital cores, according to being from bottom to top counterclockwise numbered, the list
Member is square, is made up of 8 nodes, finite element shape function is represented using linear function, X, Y, Z are digital cores coordinate respectively
Three directions of system, each three linear finite occupies by a kind of mineral matter, the modulus of elasticity and Poisson's ratio of various mineral matters by
Existing experimental data is provided.
Step D, including step are as follows:
G, 3-dimensional digital rock core elastic potential energy EnExpression formula such as formula (II) shown in:
In formula (II), p, q, rs be digital cores coordinate system x, y, tri- directions of z, d is integral sign, εpqAnd εrsPoint
Not Wei pq directions and rs directions dependent variable, CpqrsFor modulus of elasticity tensor;
Due to symmetry, strain tensor ε includes 6 independent variable (εxx、εyy、εzz、εxy、εxz、εyz), εxx、εyy、εzz、
εxy、εxz、εyzRefer respectively to the dependent variable of all directions, CpqrsIt is expressed as Cαβ, α, β is x, y, tri- directions of z, such as formula (III), formula
(IV) shown in:
In formula (III), formula (IV), upRefer to the displacement in p directions, xpRefer to the coordinate in p directions, uqRefer to the displacement in q directions,
xqRefer to the coordinate in q directions;
CpqrsIt is expressed as Cαβ, εpqIt is expressed as εα, εrsIt is expressed as εβ, the elastic potential energy E of 3-dimensional digital rock corenExpression formula formula
(II) it is converted into formula (V):
Make 3-dimensional digital rock core elastic potential energy minimum, it is necessary to meet formula (VI):
In formula (VI), umIt is displacement, m refers to all nodes and all directions;
Formula (VI) is solved using Fast Field method, finite element is carried out to all subregion based on MPI concurrent techniques
Calculate, calculate the stress and strain for obtaining all subregion, the mean stress peace of overall digital rock core is tried to achieve by arithmetic average
Strain;
K, based on 3-dimensional digital rock core mean strain stress, calculating obtains shale volume modulus, shale modulus of shearing, shale
Young's modulus and shale Poisson's ratio:
Shown in the calculation formula of shale volume modulus such as formula (VII):
In formula (VII), σ0For mean stress, σxx、σyy、σzzThe respectively direct stress in x, y, z direction, K is shale volume mould
Amount, εxx、εyy、εzzThe respectively normal strain in x, y, z direction;
Shown in the calculation formula of shale modulus of shearing such as formula (VIII):
σij=2 μ εij(Ⅷ)
In formula (VIII), σijRefer to i, the direct stress in j directions, i, j is respectively x, y direction, and μ is shale modulus of shearing, εijFor
Shearing stress;
Shown in the calculation formula of shale Young's modulus such as formula (Ⅸ):
In formula (Ⅸ), E is shale Young's modulus;
Shown in the calculation formula of shale Poisson's ratio such as formula (Ⅹ):
In formula (Ⅹ), ν is shale Poisson's ratio.
Be respectively adopted the method for the present embodiment with Experimental Method in Laboratory in the prior art, do not consider the meter of shale minerals component
Calculation method, the calculating of rock mechanics parameters is carried out to certain shale core, and the shale component information of shale core as shown in table 1, is counted
Calculate result as shown in table 2.
Table 1
Sequence number | Component | Percent by volume | E(GPa) | v |
1 | Hole (oil/gas) | 4 | 0 | 0.5 |
2 | Clay | 33 | 27 | 0.35 |
3 | Calcite | 20 | 55 | 0.3 |
4 | Quartz | 40 | 47 | 0.12 |
5 | Pyrite | 3 | 85 | 0.17 |
Table 2
Can be obtained by table 1, table 2, the present embodiment methods described computational solution precision is high, the error with test value within 5%,
Without considering that the error calculated of shale minerals component is larger.Do not consider that the computational methods of shale minerals component refer to numeral
Rock core is regarded as to be only made up of hole and skeleton, and the rock mechanics parameters of skeleton are provided by the arithmetic mean of instantaneous value of mineral constituent.
In addition, compared for the calculating time of the present embodiment method and conventional method, sub-district for different pixels number of unit
Domain number and parallel computation check figure are all 8, as shown in table 3, and research finds using Region Decomposition and uses parallel computation can
Save many calculating times.Table 3
Pixel cell | Conventional method/h | The present embodiment method/h |
60*60*60 | 0.35 | 0.06 |
80*80*80 | 4.4 | 0.72 |
100*100*100 | 8.22 | 1.55 |
Claims (7)
1. a kind of shale mechanics parameter quick calculation method based on digital cores, it is characterised in that including:
(1) shale minerals component 3-dimensional digital rock core is built
A, scanning shale structure, analysis obtains shale minerals structure chart, and calculates the volumn concentration of brittle mineral;
B, based on shale minerals structure chart, shale minerals component 3-dimensional digital rock core is built using digital cores construction method;
(2) FEM parallel computation based on shale minerals component 3-dimensional digital rock core
C, Region Decomposition, discretization are carried out successively to shale minerals component 3-dimensional digital rock core, obtain several voxels, often
Individual voxel is considered as three linear finites of 8 nodes;
D, to step C processing after shale minerals component 3-dimensional digital rock core apply a macro-strain so that the 3-dimensional digital
Rock core elastic potential energy is minimum, is asked for based on parallel computation using quick conjugate gradient method as averagely should caused by the macro-strain
Power, tries to achieve the rock mechanics parameters of overall shale, and rock mechanics parameters include shale volume modulus, shale modulus of shearing, shale
Young's modulus and shale Poisson's ratio.
2. a kind of shale mechanics parameter quick calculation method based on digital cores according to claim 1, its feature exists
In the step A scans shale structure by CT or FIB-SEM.
3. a kind of shale mechanics parameter quick calculation method based on digital cores according to claim 1, its feature exists
In, the step A scans shale structure, and analysis obtains shale minerals structure chart, including:
A, shale samples surface is polished successively, argon ion polishing;
B, shale samples are scanned using CT or FIB-SEM, carry out electronic imaging;
The X-ray spectrum that c, collection are produced in scanning shale samples formation Electron Microscope images, and itself and mineral matter is special
Levy spectrum to be compared, obtain the corresponding mineral matter of each pixel, build shale minerals structure chart;
D, the volumn concentration for calculating brittle mineral.
4. a kind of shale mechanics parameter quick calculation method based on digital cores according to claim 1, its feature exists
In the shale minerals component includes quartz, dolomite, calcite, potassium feldspar, plagioclase, pyrite and clay, the fragility
Mineral include quartz, dolomite and calcite.
5. a kind of shale mechanics parameter quick calculation method based on digital cores according to claim 1, its feature exists
In the step B, the digital cores construction method includes process simulation method, the random reconstruction method of markov.
6. a kind of shale mechanics parameter quick calculation method based on digital cores according to claim 1, its feature exists
In, the step C carries out Region Decomposition, discretization successively to shale minerals component 3-dimensional digital rock core, including:
E, shale minerals component 3-dimensional digital rock core divided equally according to size, many sub-regions are obtained, per sub-regions circle
Face meets displacement continuously, shown in the weak form such as formula (I) per sub-regions:
∫Γ(u (l)-v (l)) f (l) dl=0 (I)
In formula (I), l is coordinate variable, and u (l) and v (l) represent the displacement in adjacent subarea domain, and f (l) is finite element shape function;
Each voxel that f, the every sub-regions obtained to step e obtained after discretization, discretization is regarded as
Three linear finites of one 8 node, each three linear finite is occupied by a kind of mineral matter, the modulus of elasticity of various mineral matters
Provided with Poisson's ratio by existing experimental data.
7. according to a kind of any described shale mechanics parameter quick calculation methods based on digital cores of claim 1-6, its
It is characterised by, the step D, including step are as follows:
G, 3-dimensional digital rock core elastic potential energy EnExpression formula such as formula (II) shown in:
<mrow>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&Integral;</mo>
<msup>
<mi>d</mi>
<mn>3</mn>
</msup>
<msub>
<mi>l&epsiv;</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msub>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
In formula (II), p, q, rs be digital cores coordinate system x, y, tri- directions of z, d is integral sign, εpqAnd εrsRespectively pq
Direction and the dependent variable in rs directions, CpqrsFor modulus of elasticity tensor;
Due to symmetry, strain tensor ε includes 6 independent variable (εxx、εyy、εzz、εxy、εxz、εyz), εxx、εyy、εzz、εxy、εxz、
εyzRefer respectively to the dependent variable of all directions, CpqrsIt is expressed as Cαβ, α, β is x, y, tri- directions of z, such as formula (III), formula (IV) institute
Show:
<mrow>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>p</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>I</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>q</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>q</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>V</mi>
<mo>)</mo>
</mrow>
</mrow>
In formula (III), formula (IV), upRefer to the displacement in p directions, xpRefer to the coordinate in p directions, uqRefer to the displacement in q directions, xqIt is
Refer to the coordinate in q directions;
CpqrsIt is expressed as Cαβ, εpqIt is expressed as εα, εrsIt is expressed as εβ, the elastic potential energy E of 3-dimensional digital rock corenExpression formula formula (II)
It is converted into formula (V):
<mrow>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&Integral;</mo>
<msup>
<mi>d</mi>
<mn>3</mn>
</msup>
<msub>
<mi>r&epsiv;</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>&alpha;</mi>
<mi>&beta;</mi>
</mrow>
</msub>
<msub>
<mi>&epsiv;</mi>
<mi>&beta;</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>)</mo>
</mrow>
</mrow>
Make 3-dimensional digital rock core elastic potential energy minimum, it is necessary to meet formula (VI):
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
In formula (VI), umIt is displacement, m refers to all nodes and all directions;
Formula (VI) is solved using Fast Field method, FEM calculation is carried out to all subregion based on MPI concurrent techniques,
The stress and strain for obtaining all subregion is calculated, the mean stress peace for trying to achieve overall digital rock core by arithmetic average all should
Become;
K, based on 3-dimensional digital rock core mean strain stress, calculating obtains shale volume modulus, shale modulus of shearing, shale Young
Modulus and shale Poisson's ratio:
Shown in the calculation formula of shale volume modulus such as formula (VII):
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mi>I</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
In formula (VII), σ0For mean stress, σxx、σyy、σzzThe respectively direct stress in x, y, z direction, K is shale volume modulus,
εxx、εyy、εzzThe respectively normal strain in x, y, z direction;
Shown in the calculation formula of shale modulus of shearing such as formula (VIII):
σij=2 μ εij (Ⅷ)
In formula (VIII), σijRefer to i, the direct stress in j directions, i, j is respectively x, y direction, and μ is shale modulus of shearing, εijShould to cut
Power;
Shown in the calculation formula of shale Young's modulus such as formula (Ⅸ):
<mrow>
<mi>E</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>9</mn>
<mi>K</mi>
<mi>&mu;</mi>
</mrow>
<mrow>
<mn>3</mn>
<mi>K</mi>
<mo>+</mo>
<mi>&mu;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
In formula (Ⅸ), E is shale Young's modulus;
Shown in the calculation formula of shale Poisson's ratio such as formula (Ⅹ):
<mrow>
<mi>&nu;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>K</mi>
<mo>-</mo>
<mn>2</mn>
<mi>&mu;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mi>K</mi>
<mo>+</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
In formula (Ⅹ), ν is shale Poisson's ratio.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710288425.0A CN107045580B (en) | 2017-04-27 | 2017-04-27 | Shale mechanical parameter rapid calculation method based on digital rock core |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710288425.0A CN107045580B (en) | 2017-04-27 | 2017-04-27 | Shale mechanical parameter rapid calculation method based on digital rock core |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107045580A true CN107045580A (en) | 2017-08-15 |
CN107045580B CN107045580B (en) | 2020-04-21 |
Family
ID=59546577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710288425.0A Active CN107045580B (en) | 2017-04-27 | 2017-04-27 | Shale mechanical parameter rapid calculation method based on digital rock core |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107045580B (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107908913A (en) * | 2017-12-22 | 2018-04-13 | 中国海洋大学 | Geodynamic digital-to-analogue algorithm based on parallel computer |
CN108376422A (en) * | 2018-01-05 | 2018-08-07 | 中山大学 | A kind of continuous across the scale measurement characterizing method of rock-soil material |
CN109472112A (en) * | 2018-12-06 | 2019-03-15 | 长江大学 | A kind of shale digital cores modeling method |
CN110473597A (en) * | 2019-09-05 | 2019-11-19 | 中国石油大学(北京) | Conglomerate mechanical property evaluation analysis method and system |
CN111680440A (en) * | 2020-06-03 | 2020-09-18 | 中国地质大学(北京) | Method for analyzing acoustic characteristics of digital rock cores at different temperatures |
CN111767650A (en) * | 2020-06-28 | 2020-10-13 | 中国石油大学(华东) | Method and device for estimating equivalent elastic parameters of digital rock core |
CN112131714A (en) * | 2020-08-28 | 2020-12-25 | 中国石油天然气集团有限公司 | Shale matrix modulus calculation method and system |
CN112392472A (en) * | 2019-08-19 | 2021-02-23 | 中国石油天然气股份有限公司 | Method and device for determining integrated development mode of shale and adjacent oil layer |
CN112485278A (en) * | 2020-11-24 | 2021-03-12 | 洛阳理工学院 | Method for researching microstructure and properties of rock based on molecular layer |
CN113533104A (en) * | 2020-04-17 | 2021-10-22 | 中国石油化工股份有限公司 | Method for acquiring elastic parameters before and after shale water-rock action |
CN116844651A (en) * | 2023-03-05 | 2023-10-03 | 东北石油大学 | Method for solving mineral microscopic elastic parameters based on molecular dynamics simulation |
CN116912458A (en) * | 2023-09-05 | 2023-10-20 | 煤炭科学研究总院有限公司 | Finite element model generation method and device for coal rock and electronic equipment |
CN118013784A (en) * | 2024-01-22 | 2024-05-10 | 中国石油大学(华东) | Digital rock core reconstruction and compressibility evaluation method based on compact sandstone |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020276A (en) * | 2014-06-24 | 2014-09-03 | 中国石油大学(北京) | Determination method for mechanical parameters of transverse isotropy shale reservoir rocks |
CN104729904A (en) * | 2015-03-31 | 2015-06-24 | 中国石油大学(华东) | Complicated rock core preparation method based on CT scanning and 3D printing |
CN105221141A (en) * | 2014-06-23 | 2016-01-06 | 中国石油化工股份有限公司 | A kind of mud shale brittleness index Forecasting Methodology |
US20160025895A1 (en) * | 2014-07-28 | 2016-01-28 | Schlumberger Technology Corporation | Well Treatment with Digital Core Analysis |
CN106547034A (en) * | 2016-11-09 | 2017-03-29 | 西南石油大学 | A kind of method for calculating compact reservoir rock brittleness index |
CN206132572U (en) * | 2016-10-18 | 2017-04-26 | 中国石油大学(华东) | Hyposmosis rock core time domain developments CT scan imbibition experimental apparatus |
-
2017
- 2017-04-27 CN CN201710288425.0A patent/CN107045580B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105221141A (en) * | 2014-06-23 | 2016-01-06 | 中国石油化工股份有限公司 | A kind of mud shale brittleness index Forecasting Methodology |
CN104020276A (en) * | 2014-06-24 | 2014-09-03 | 中国石油大学(北京) | Determination method for mechanical parameters of transverse isotropy shale reservoir rocks |
US20160025895A1 (en) * | 2014-07-28 | 2016-01-28 | Schlumberger Technology Corporation | Well Treatment with Digital Core Analysis |
CN104729904A (en) * | 2015-03-31 | 2015-06-24 | 中国石油大学(华东) | Complicated rock core preparation method based on CT scanning and 3D printing |
CN206132572U (en) * | 2016-10-18 | 2017-04-26 | 中国石油大学(华东) | Hyposmosis rock core time domain developments CT scan imbibition experimental apparatus |
CN106547034A (en) * | 2016-11-09 | 2017-03-29 | 西南石油大学 | A kind of method for calculating compact reservoir rock brittleness index |
Non-Patent Citations (1)
Title |
---|
刘学锋: "基于数字岩心的岩石声电特性微观数值模拟研究", 《中国博士学位论文全文数据库基础科学辑》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107908913A (en) * | 2017-12-22 | 2018-04-13 | 中国海洋大学 | Geodynamic digital-to-analogue algorithm based on parallel computer |
CN107908913B (en) * | 2017-12-22 | 2020-12-11 | 中国海洋大学 | Earth power digital-analog method based on parallel computer |
CN108376422A (en) * | 2018-01-05 | 2018-08-07 | 中山大学 | A kind of continuous across the scale measurement characterizing method of rock-soil material |
CN109472112A (en) * | 2018-12-06 | 2019-03-15 | 长江大学 | A kind of shale digital cores modeling method |
CN109472112B (en) * | 2018-12-06 | 2023-05-30 | 长江大学 | Shale digital rock core modeling method |
CN112392472B (en) * | 2019-08-19 | 2022-08-02 | 中国石油天然气股份有限公司 | Method and device for determining integrated development mode of shale and adjacent oil layer |
CN112392472A (en) * | 2019-08-19 | 2021-02-23 | 中国石油天然气股份有限公司 | Method and device for determining integrated development mode of shale and adjacent oil layer |
CN110473597A (en) * | 2019-09-05 | 2019-11-19 | 中国石油大学(北京) | Conglomerate mechanical property evaluation analysis method and system |
CN113533104A (en) * | 2020-04-17 | 2021-10-22 | 中国石油化工股份有限公司 | Method for acquiring elastic parameters before and after shale water-rock action |
CN111680440A (en) * | 2020-06-03 | 2020-09-18 | 中国地质大学(北京) | Method for analyzing acoustic characteristics of digital rock cores at different temperatures |
CN111680440B (en) * | 2020-06-03 | 2024-03-19 | 中国地质大学(北京) | Method for analyzing acoustic characteristics of digital rock core at different temperatures |
CN111767650A (en) * | 2020-06-28 | 2020-10-13 | 中国石油大学(华东) | Method and device for estimating equivalent elastic parameters of digital rock core |
CN111767650B (en) * | 2020-06-28 | 2022-05-03 | 中国石油大学(华东) | Method and device for estimating equivalent elastic parameters of digital rock core |
CN112131714A (en) * | 2020-08-28 | 2020-12-25 | 中国石油天然气集团有限公司 | Shale matrix modulus calculation method and system |
CN112485278A (en) * | 2020-11-24 | 2021-03-12 | 洛阳理工学院 | Method for researching microstructure and properties of rock based on molecular layer |
CN112485278B (en) * | 2020-11-24 | 2024-08-27 | 洛阳理工学院 | Research method for rock microstructure and property based on molecular level |
CN116844651A (en) * | 2023-03-05 | 2023-10-03 | 东北石油大学 | Method for solving mineral microscopic elastic parameters based on molecular dynamics simulation |
CN116844651B (en) * | 2023-03-05 | 2024-01-02 | 东北石油大学 | Method for solving mineral microscopic elastic parameters based on molecular dynamics simulation |
CN116912458A (en) * | 2023-09-05 | 2023-10-20 | 煤炭科学研究总院有限公司 | Finite element model generation method and device for coal rock and electronic equipment |
CN118013784A (en) * | 2024-01-22 | 2024-05-10 | 中国石油大学(华东) | Digital rock core reconstruction and compressibility evaluation method based on compact sandstone |
Also Published As
Publication number | Publication date |
---|---|
CN107045580B (en) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107045580A (en) | A kind of shale mechanics parameter quick calculation method based on digital cores | |
Liu et al. | Quantification of the microstructures of Bakken shale reservoirs using multi-fractal and lacunarity analysis | |
Liu et al. | Multi-scale fractal analysis of pores in shale rocks | |
Yu et al. | Numerical modeling of jointed rock under compressive loading using X-ray computerized tomography | |
Volland et al. | Anisotropy quantification: the application of fractal geometry methods on tectonic fracture patterns of a Hercynian fault zone in NW Sardinia | |
CN105651966A (en) | Shale oil and gas high-quality reservoir stratum evaluation method and parameter determination method | |
CN106600553A (en) | DEM super-resolution method based on convolutional neural network | |
NO20121394A1 (en) | Procedure for obtaining consistent and integrated physical properties of porous media | |
CN104885124A (en) | Method for producing a three-dimensional characteristic model of a porous material sample for analysis of permeability characteristics | |
CN111027818A (en) | Shale oil classification evaluation method | |
Zhang et al. | Permeability prediction of carbonate rocks based on digital image analysis and rock typing using random forest algorithm | |
Zhang et al. | Shale lamina thickness study based on micro-scale image processing of thin sections | |
CN113075731B (en) | Deep reservoir continuity wellbore digital modeling method and device | |
Pavičić et al. | Geometric and fractal characterization of pore systems in the Upper Triassic dolomites based on image processing techniques (example from Žumberak Mts, NW Croatia) | |
Wu et al. | Fractal characteristics of low-permeability sandstone reservoirs | |
Kahl et al. | Microfabric and anisotropy of elastic waves in sandstone–An observation using high-resolution X-ray microtomography | |
CN111507988A (en) | Broken rock mass modeling and seepage test method based on digital image processing | |
Kaczmarek et al. | Numerical approach in recognition of selected features of rock structure from hybrid hydrocarbon reservoir samples based on microtomography | |
Tien et al. | Uncertainty in estimation of volumetric block proportion of bimrocks by using scanline method | |
CN114781151A (en) | Continental facies shale multi-cluster fracturing cluster spacing optimization method based on longitudinal fracture crossing | |
CN113588722A (en) | High-temperature rock electrical property analysis method and system based on digital rock physics | |
Pyrak-Nolte et al. | Fractures: finite-size scaling and multifractals | |
Bishwal | Scope of 3-D printing in mining and geology: An overview | |
Jacques et al. | The use of microtomography in structural geology: A new methodology to analyse fault faces | |
Zhou et al. | The shape of moving boundary of fluid flow in sandstone: Video microscopic investigation and stochastic modeling approach |
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 |