CN116296946B - Shale stratum development degree characterization method and device based on fractal-fluctuation theory - Google Patents
Shale stratum development degree characterization method and device based on fractal-fluctuation theory Download PDFInfo
- Publication number
- CN116296946B CN116296946B CN202310203793.6A CN202310203793A CN116296946B CN 116296946 B CN116296946 B CN 116296946B CN 202310203793 A CN202310203793 A CN 202310203793A CN 116296946 B CN116296946 B CN 116296946B
- Authority
- CN
- China
- Prior art keywords
- scratch
- depth
- distance
- development
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011161 development Methods 0.000 title claims abstract description 101
- 238000012512 characterization method Methods 0.000 title claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000008859 change Effects 0.000 claims description 8
- 239000011435 rock Substances 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 4
- 230000035515 penetration Effects 0.000 claims description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 22
- 229910052683 pyrite Inorganic materials 0.000 description 13
- 239000011028 pyrite Substances 0.000 description 13
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 13
- 239000007789 gas Substances 0.000 description 7
- 239000004973 liquid crystal related substance Substances 0.000 description 6
- 239000003079 shale oil Substances 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 210000000434 stratum corneum Anatomy 0.000 description 4
- 241000258937 Hemiptera Species 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 101100063432 Caenorhabditis elegans dim-1 gene Proteins 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- -1 argon ion Chemical class 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- HBMJWWWQQXIZIP-UHFFFAOYSA-N silicon carbide Chemical compound [Si+]#[C-] HBMJWWWQQXIZIP-UHFFFAOYSA-N 0.000 description 1
- 229910010271 silicon carbide Inorganic materials 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/40—Investigating hardness or rebound hardness
- G01N3/42—Investigating hardness or rebound hardness by performing impressions under a steady load by indentors, e.g. sphere, pyramid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
- G01N3/068—Special adaptations of indicating or recording means with optical indicating or recording means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/40—Investigating hardness or rebound hardness
- G01N3/42—Investigating hardness or rebound hardness by performing impressions under a steady load by indentors, e.g. sphere, pyramid
- G01N3/44—Investigating hardness or rebound hardness by performing impressions under a steady load by indentors, e.g. sphere, pyramid the indentors being put under a minor load and a subsequent major load, i.e. Rockwell system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/40—Investigating hardness or rebound hardness
- G01N3/42—Investigating hardness or rebound hardness by performing impressions under a steady load by indentors, e.g. sphere, pyramid
- G01N3/46—Investigating hardness or rebound hardness by performing impressions under a steady load by indentors, e.g. sphere, pyramid the indentors performing a scratching movement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/48—Analysis of texture based on statistical description of texture using fractals
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
The invention discloses a shale tattoo development degree characterization method and device based on a fractal-fluctuation theory. The method comprises the following steps: acquiring shale scratch test data; drawing a scribing distance-depth curve and a scribing depth-distance curve respectively by using scratch test data; dividing the scribing distance-depth curve by adopting different scales and establishing scratch complexity; dividing the drawn depth-distance curve by adopting different threshold values and establishing a scratch fluctuation step length; and calculating a layer development index according to the fluctuation step length and the scratch complexity, and defining a layer development classification index to represent the layer development degree. The invention solves the problem that the scratch test in the prior art lacks a shale grain development degree characterization method. In addition, the method provided by the invention has the advantages of simple required test equipment, clear test principle and easiness in implementation.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration and material analysis, in particular to a shale tattoo development degree characterization method and device based on fractal-fluctuation theory.
Background
Shale oil gas is an important component of natural gas resources in China, and a large amount of shale oil gas is generated in Sichuan basin, erdos basin, songliao basin and Songyang basin of China. With the continuous deep exploration and development of shale oil and gas, students find that shale oil and gas reservoirs have strong heterogeneity mainly caused by layering stack formed by the change of paleo-deposition environment and the like. The shale stratum can reflect microstructure characteristics and reservoir performance of the shale reservoir, and can reveal the expansion rule and fracturing effect of the pressure fracture network in the shale reservoir reconstruction process. The layers widely develop in the shale layer rich in organic matters, and the quality of the shale oil and gas reservoir layer generated from the self-storage is seriously influenced, so that shale layer research becomes a focus in shale oil and gas exploration and development basic research work.
At present, the stratum characterization method mostly adopts single triaxial rock mechanics experiments with different angles, which are all scatter measurement, and mineral resource exploitation gradually develops to the deep part, so that the geological structure environment is complex, and the coring operation becomes very difficult. And the sample is in a centimeter scale and cannot be characterized for millimeter-scale layers. The continuous section mechanical property can be obtained by adopting a scratch means, the prepared sample size is small, the maximum information about the heterogeneous behavior of the rock can be obtained from a small amount of core analysis, and the advantage is obvious. However, the scratch means can only obtain the parameters such as hardness, fracture toughness and the like, and a shale grain development degree characterization method is lacked.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides the following technical scheme.
The invention provides a shale stratum development degree characterization method based on a fractal-fluctuation theory, which comprises the following steps:
acquiring shale scratch test data;
drawing a scribing distance-depth curve and a scribing depth-distance curve respectively by using scratch test data;
dividing the scribing distance-depth curve by adopting different scales and establishing scratch complexity;
dividing the drawn depth-distance curve by adopting different threshold values and establishing a scratch fluctuation step length;
and calculating a layer development index according to the fluctuation step length and the scratch complexity, and defining a layer development classification index to represent the layer development degree.
Preferably, the shale scratch test data are obtained through a scratch test, wherein the scratch test comprises three stages of pre-scanning, engraving scanning and post-scanning; in the pre-scanning stage, scanning the surface of the sample with constant contact force, and recording the change of the scribing depth along with the scribing distance; in the stage of scribing scanning, linearly loading a scribing sample, and recording the change of loading force, friction force and scribing depth along with the scribing distance; in the post-scanning stage, the scratch morphology is scanned with constant contact force, and the change of the scratch depth along with the scratch distance is recorded.
Preferably, said dividing said scribe distance-depth curve and establishing a scribe complexity using different scales comprises:
dividing the vertical coordinate of the dividing distance-depth curve into depths by using different scales, and generating a plurality of sections corresponding to each scale;
within each section, the corresponding number of boxes is calculated according to the following formula：
If->Is an integer;
if->Is a non-integer;
in the method, in the process of the invention,for the number of boxes->For the scale->Is an integer part of the maximum value of the signal data within the interval,an integer part of a minimum value of the signal data in the section;
adding the box numbers of all intervals corresponding to each scale to obtain the total box number;
establishing a double logarithmic relation curve between each scale and the corresponding total box number;
calculating the slope of a straight line segment in a double-logarithmic relation curve corresponding to each scale as the complexity of a corresponding scribing distance-depth curve;
the mean of the complexities is calculated as the scratch complexity of the rock.
Preferably, the equally dividing the scribing depth-distance curve by using different threshold values and establishing a scratch fluctuation step length includes:
performing equivalence division on the ordinate scratch distances of the scratch depth-distance curve by using different threshold values to generate domains with different numbers;
and calculating the number of the resident points of the scribing depth-distance curve contained in each domain, and defining the corresponding threshold value as the scratch fluctuation step length when the duty ratio of the number of the resident points of 1 reaches the maximum value.
Preferably, the method further comprises, before performing the equivalence division on the ordinate-type distance of the depth-distance curve by using different threshold values:
and amplifying the length of the scribing distance on the ordinate of the scribing depth-distance curve according to a preset proportion, and generating a corrected scribing depth-distance curve while keeping the value of the ordinate unchanged.
Preferably, the calculating the layer development index according to the fluctuation step length and the scratch complexity comprises:
，
wherein, the liquid crystal display device comprises a liquid crystal display device,is the index of development of the schlieren>Takes the value of 100 for the correction coefficient>For fluctuation step size +.>Is scratch complexity.
Preferably, the defining the index of the grade of the development of the tattoos to characterize the degree of development of the tattoos includes:
when (when)The degree of development is defined as the first level of development; when->The degree of development is defined as the secondary level of development; when->The degree of development is defined as the tertiary level of development.
The second aspect of the invention provides a shale grain development degree characterization device based on fractal-fluctuation theory, which comprises:
the data acquisition module is used for acquiring shale scratch test data;
the curve drawing module is used for drawing a drawing distance-depth curve and a drawing depth-distance curve respectively by using scratch test data;
the scratch complexity establishing module is used for dividing the scribing distance-depth curve by adopting different scales and establishing scratch complexity;
the fluctuation step length establishing module is used for equally dividing the dividing depth-distance curve by adopting different threshold values and establishing a scratch fluctuation step length;
the layer development degree characterization module is used for calculating a layer development index according to the fluctuation step length and the scratch complexity, and defining a layer development grading index to characterize the layer development degree.
A third aspect of the invention provides a memory storing a plurality of instructions for implementing the method as described in the first aspect.
A fourth aspect of the invention provides an electronic device comprising a processor and a memory coupled to the processor, the memory storing a plurality of instructions loadable and executable by the processor to enable the processor to perform the method of the first aspect.
The beneficial effects of the invention are as follows: the invention provides a shale tattoo development degree characterization method and device based on a fractal-fluctuation theory. The method is researched by means of shale scratch test, data processing, image analysis and the like, and scratch fluctuation step length is established based on fluctuation theoryThe method is used for quantitatively expressing the fluctuation frequency of the scratch curve of the lamellar shale; establishing a stratum corneum shale development complexity index based on fractal theory>The method is used for quantitatively expressing the complexity of the scratch curve of the lamellar shale; finally combining fractal-fluctuation theory and introducing the stratum corneum development index +.>Definition of the division of the development of the tattooThe grade index is used for comprehensively representing the development degree of the tattoos. The invention solves the problem that the scratch test in the prior art lacks a shale grain development degree characterization method. In addition, the method provided by the invention has the advantages of simple required test equipment, clear test principle and easiness in implementation.
Drawings
FIG. 1 is a flow chart of a shale tattoo development degree characterization method based on fractal-fluctuation theory;
FIG. 2 is a schematic illustration of a scratch test;
figures 3a-3d are graphs of scribe distance versus depth, respectively, for different sample scratches;
FIGS. 4a-4d are plot of the scoring distance versus depth for pyrite divided by scales 50, 100, 150 and 200, respectively;
FIG. 5 is a schematic view of a displayA double logarithmic relationship curve;
FIG. 6a is a scribe depth-distance curve, and FIG. 6b is a modified scribe depth-distance curve;
FIGS. 7a-7d are plot depth-to-distance curves of pyrite divided by threshold values 10, 15, 20, 30, respectively;
FIG. 8 is a diagram illustrating the standing-point duty ratios corresponding to different threshold values;
FIG. 9 is a schematic diagram of a functional module of a shale tattoo development degree characterization device based on fractal-wave theory;
in fig. 2, the meaning of each symbol is:
1. scratch test sample, 2, scratch indenter, FN, vertical load, FT, horizontal load.
Detailed Description
In order to better understand the above technical solutions, the following detailed description will be given with reference to the accompanying drawings and specific embodiments.
Aiming at the problems existing in the prior art, the invention introduces the fractal complexity of the fractal theory and the fluctuation step length of the fluctuation theory based on the fractal theory and the fluctuation theory, quantitatively describes the characteristics of the scratch curve,and introducing a layer development indexAnd defining a grade index of the development of the tattoos to comprehensively characterize the development degree of the tattoos.
Example 1
As shown in fig. 1, the embodiment of the invention provides a shale grain development degree characterization method based on fractal-fluctuation theory, which comprises the following steps: s101, shale scratch test data are obtained; s102, drawing a scribing distance-depth curve and a scribing depth-distance curve respectively by using scratch test data; s103, dividing the scribing distance-depth curve by adopting different scales and establishing scratch complexity; s104, equally dividing the drawn depth-distance curve by adopting different threshold values and establishing a scratch fluctuation step length; s105, calculating a layer development index according to the fluctuation step length and the scratch complexity, and defining a layer development grading index to represent the layer development degree.
Wherein, prior to executing step S101, the shale scratch test may be performed first to obtain relevant data and store the relevant data for subsequent use.
In the actual operation process, sample preparation can be firstly carried out before scratch test is carried out, and the specific preparation method can be as follows: in situ shale sampling, the sample is initially polished parallel to its bedding surface with silicon carbide paper, and then the sample surface is further polished with an argon ion suspension to ensure that the shale sample surface is smooth. After the sample preparation is completed, a scratch test can be performed, and specifically, the following method can be adopted: scratch tests were performed in shale grain page planes. In the test, a rocweil C diamond probe with a half apex angle of 60 degrees was used, with a tip radius r=200 mm. For each test, the load was constant at 50 mN while the scratch speed remained constant, at 30 um/s on a 0.5mm scratch path. A schematic of the scratch test can be shown in fig. 2. Scratch testing is generally divided into three stages: pre-scanning, engraving scanning and post-scanning stages. Wherein, in the pre-scan stage: the sample surface was scanned with a constant contact force of 20uN and the variation of the sample surface relief (i.e., the scribe depth) with scribe distance was recorded. Etching and scanning: the scribing sample is linearly loaded, and the change of loading force, friction force and scribing depth along with the scribing distance is recorded. Post-scanning stage: the scratch morphology was scanned with a constant contact force of 20uN and the variation of the surface relief (i.e. the scratch depth) of the sample with the scratch distance was recorded.
Step S102 is executed, each scratch test generates a group of data, wherein the group of data comprises data of three stages of pre-scanning, engraving scanning and post-scanning, and effective data comes from the engraving scanning stage. Processing the generated data, wherein the drawing depth is set as an abscissa, the drawing distance is set as an ordinate, and a drawing depth-distance curve is generated; the scribing distance is set as the abscissa, and the scribing depth is set as the ordinate, and a scribing distance-depth curve is generated. Different samples produced different scratch penetration depth-distance curves and penetration distance-depth curves. As an example, the scribing distance-depth curve of pyrite sample scratches, sandy sample scratches, carbon dust sample scratches, and scale sample scratches are shown in fig. 3a-3d, respectively. Fig. 6a is a plot of depth of scratch versus distance for a pyrite sample.
Step S103 is performed to divide the scribe distance-depth curves and establish the scratch complexity using different scales, which can be implemented specifically as follows: dividing the ordinate of the dividing distance-depth curve into a plurality of intervals corresponding to each scale. In this step, different scales are selected based on fractal theoryI.e. grid width->Dividing the drawn distance-depth plane image to form a plurality of equidistant square grids, wherein the length of the ordinate interval of each square grid is +.>。
Calculating the corresponding number of boxes in each section. If the integer part of the maximum value of the signal data in the square is divided by the scale +.>When the result of (a) is an integer, subtracting the integer part of the maximum value and the integer part of the minimum value of the signal data in the square, and comparing the result of the subtracting with the scale +.>As the number of boxes; if the integer part of the maximum value of the signal data in the square is divided by the scale +.>When the result of (2) is not an integer, subtracting the integer part of the maximum value and the integer part of the minimum value of the signal data in the square, and comparing the result of the subtracting with the scale +.>The ratio of (2) plus one is used as the number of boxes.
The method can be specifically expressed by the following formula:
if->Is an integer;
if->Is a non-integer;
in the method, in the process of the invention,for the number of boxes->For the scale->Is an integer part of the maximum value of the signal data within the interval,is the integer part of the minimum value of the signal data within the interval.
Adding the box numbers of all intervals corresponding to each scale to obtain the total box number; establishing a double logarithmic relation curve between each scale and the corresponding total box number; calculating the slope of a straight line segment in a double-logarithmic relation curve corresponding to each scale as the complexity of a corresponding scribing distance-depth curve; the mean of the complexities is calculated as the scratch complexity of the rock.
In one embodiment, taking pyrite-based layered shale as an example, different scales are selected based on fractal theoryI.e. grid width->(50, 100, 150, 200) dividing the distance-depth-binned planar image into a plurality of equally spaced bins as shown in fig. 4a-4 d. The total length of the ordinate-drawn depth can be denoted as M if +.>Is an integer, the planar image of the scratch distance-depth is divided into +.>The intervals of the equidistant grids are as follows:
；
if it isNot an integer, then get->The integer part of (2) is denoted as M 1 The sections divided on the ordinate are:
；
the last space grid;
as an example such as when m=5.5,when it is the case, it can be divided into two equally spaced intervals of 2 and one non-equally spaced interval.
In any interval on the ordinate axis, e.g. inIn the interval, the integer part of the maximum value of the signal data is marked as +.>And the integer part of the minimum is denoted +.>. Will->And->Calculate and scale after subtraction>Ratio of (1), if->Divide by the scale->Is an integer, the ratio is taken as the box number +.>If->Divide by the scale->If not, the result of adding one to the ratio is taken as the box number +.>。
The method can be specifically expressed by the following formula:
when (when)When the number is an integer: />
When (when)When not an integer: />。
For example, when the maximum value of signal data in the interval is 5.5, the minimum value is 1.1, the scaleWhen (I)>，The number of boxes is +.>。
Number of boxes of all sections corresponding to the scaleAdding to obtain the total box number of the scale-corresponding scratch-in distance-depth plane image +.>The method can be specifically expressed by the following formula:
superscript in the summation equation above if M is not an integerIs an integer part of (c).
Taking different scalesRepeating the above steps to obtain corresponding +.>Finally get->A double logarithmic relationship is shown in fig. 5. />The shape of the log-log relationship, associated with different shale materials, contains a large amount of information about the rock. If a straight line (or an approximate straight line) exists in the curve, the scale range of the straight line is +.>In this regard, the slope of this line is considered to be the complexity D of the waveform image:
FIG. 5 shows a distance-depth curve of a scratch according to four scratches (corresponding to data sets text1, text2, text3, and text4, respectively)The line is drawn according to the data set corresponding to the scratch) corresponding to the scribing distance-depth curve、Data-rendered +.>Taking straight line segments or approximate straight line segments in the curves respectively, calculating the slope of the straight line segments or approximate straight line segments to obtain the complexity D of the curve (signal waveform) of the distance-depth, and finally calculating the mean value ∈of D>Scratch complexity as the rock>The expression can be expressed as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,for the complexity of the scribing distance-depth curve (signal waveform) corresponding to the ith scratch, i=1, 2,3 … … n, n being the number of scratches.
Step S104 is executed, and the scribing depth-distance curve is divided by using different threshold values and the scratch fluctuation step length is established, which can be implemented in the following manner: performing equivalence division on the ordinate scratch distances of the scratch depth-distance curve by using different threshold values to generate domains with different numbers; and calculating the number of the resident points of the scribing depth-distance curve contained in each domain, and defining the corresponding threshold value as the scratch fluctuation step length when the duty ratio of the number of the resident points of 1 reaches the maximum value.
Further, before the performing the equivalence division on the ordinate of the drawn depth-distance curve by using different threshold values, the method may further include: and amplifying the length of the scribing distance on the ordinate of the scribing depth-distance curve according to a preset proportion, and generating a corrected scribing depth-distance curve while keeping the value of the ordinate unchanged.
In one embodiment, taking pyrite-based layered shale as an example, the scribing depth is set as an abscissa x, the scribing distance is defined as an ordinate y based on the wave theory, and the scribing depth-distance curve is drawn. The ordinate length is scaled up in proportion but the ordinate value remains unchanged, resulting in a modified scribe depth-distance curve as shown in fig. 6 b. The ordinate scribe distance of the modified scribe depth-distance curve is equally divided by different threshold values U (u=10, 15, 20, 30) in which y/U fields will be generated as shown in fig. 7a-7 d. And calculating the number S of standing points of the waveform image contained in each domain, and when the duty ratio of the number S of the standing points is 1 and reaches the maximum value, obtaining the corresponding domain value as the fluctuation step lambda of the scratch. Taking pyrite lamellar shale as an example, when the threshold value u=10, the duty ratio of the number of standing points at this time is 1 is maximum, as shown in fig. 8, and therefore, the fluctuation step λ=10 of pyrite lamellar shale.
Step 105 is executed to calculate the growth index of the layer according to the fluctuation step length and the scratch complexity, and specifically, the following formula can be adopted for calculation:
wherein, the liquid crystal display device comprises a liquid crystal display device,is the index of development of the schlieren>Takes the value of 100 for the correction coefficient>For fluctuation step size +.>Is scratch complexity.
Then, the index of the grade of the development of the tattoo can be defined according to the index of the development of the tattoo to represent the degree of the development of the tattoo, which can comprise: when (when)The degree of development is defined as the first level of development; when->The degree of development is defined as the secondary level of development; when->The degree of development is defined as the tertiary level of development.
In one embodiment, the variation of the surface relief (i.e., the scratch depth) of the sample with the scratch distance as the scratch test is performed, the scratch depth-distance curve and the scratch distance-depth curve are deeply analyzed based on the fractal theory and the wave theory, and the internal cause of the variation of the curve is considered to be due to the strong heterogeneity of the lamellar shale. The fractal complexity of the fractal theory and the fluctuation step length of the fluctuation theory are introduced, and the scratch curve characteristic is quantitatively described, so that the development degree of the tattoo shale is evaluated to a certain extent.
Analysis according to fluctuation theory can be carried out to reduce the development degree of the schlieren when the fluctuation step length is larger, and pyrite is adoptedScale insect->Sandy lamellar shale->For example, by the above method operation, the following results can be obtained:
the development degree of the three shale layers is sequentially from high to low, namely a pyrite layer, a sandy layer and a scale insect layer.
Analysis according to the fractal theory can be carried out when the complexity is higher, and the development degree of the tattoo layer is higher, and by taking pyrite Dim 1, scale insect Dim 2 and sandy tattoo shale Dim 3 as examples, the following results can be obtained through the operation of the method:
the development degree of the three shale layers is sequentially from high to low, namely a pyrite layer, a sandy layer and a scale insect layer.
Finally, combining fractal theory and fluctuation theory, and introducing a stratum corneum development indexAnd according to the index of the growth of the lamina->And defining a grade index of the development of the tattoos to comprehensively characterize the development degree of the tattoos. I.e. the product of the fluctuation step and the complexity is defined as the stratum corneum development index +.>：
Where α is a correction factor, and is typically 100. When (when)The degree of development is defined as the first level of development when +.>The degree of development is defined as the level of secondary development when +.>The degree of development is defined as the tertiary level of development. The greater the BI value, the greater the degree of grain development. For example, pyrite lamellar shale bi=17, the degree of development is tertiary level of developmentThe method comprises the steps of carrying out a first treatment on the surface of the The scale shale bi=6, the degree of development is the first level of development.
Example two
As shown in fig. 9, an embodiment of the present invention provides a shale tattoo development degree characterization device based on fractal-fluctuation theory, including: the data acquisition module 201 is used for acquiring shale scratch test data; the curve drawing module 202 is configured to draw a scribing distance-depth curve and a scribing depth-distance curve respectively using scratch test data; a scratch complexity establishing module 203, configured to divide the scribing distance-depth curves with different scales and establish scratch complexity; the fluctuation step length establishing module 204 is configured to equally divide the scribing depth-distance curve by using different threshold values and establish a scratch fluctuation step length; the layer development degree characterization module 205 is configured to calculate a layer development index according to the fluctuation step size and the scratch complexity, and define a layer development classification index to characterize the layer development degree.
The apparatus may be implemented by the method described in the first embodiment, which is not described herein.
The embodiment of the invention also provides a memory, which stores a plurality of instructions for implementing the method according to the embodiment one.
The embodiment of the invention also provides an electronic device, which comprises a processor and a memory connected with the processor, wherein the memory stores a plurality of instructions, and the instructions can be loaded and executed by the processor so that the processor can execute the method in the embodiment.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. It is therefore intended that the following claims be interpreted as including the preferred embodiments and all such alterations and modifications as fall within the scope of the invention. It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (8)
1. The shale stratum development degree characterization method based on the fractal-fluctuation theory is characterized by comprising the following steps of:
acquiring shale scratch test data;
drawing a scribing distance-depth curve and a scribing depth-distance curve respectively by using scratch test data;
dividing the scribing distance-depth curve by adopting different scales and establishing scratch complexity;
dividing the drawn depth-distance curve by adopting different threshold values and establishing a scratch fluctuation step length;
calculating a layer development index according to the fluctuation step length and the scratch complexity, and defining a layer development classification index to represent the layer development degree;
wherein said dividing said scribe-in distance-depth curve and creating a scribe complexity using different scales comprises:
dividing the ordinate of the dividing distance-depth curve into depths by using different scales, and generating a plurality of intervals corresponding to each scale; within each section, the corresponding number of boxes N is calculated according to the following formula δ (m)：
If->Is an integer;
if->Is a non-integer;
wherein N is δ (m) is the number of cassettes,for the scale->Is an integer part of the maximum value of the signal data within the interval,an integer part of a minimum value of the signal data in the section; adding the box numbers of all intervals corresponding to each scale to obtain the total box number; establishing a double logarithmic relation curve between each scale and the corresponding total box number; calculating the slope of a straight line segment in a double-logarithmic relation curve corresponding to each scale as the complexity of a corresponding scribing distance-depth curve; calculating the average value of the complexity as the scratch complexity of the rock;
the calculating the growth index of the tattoo according to the fluctuation step length and the scratch complexity comprises the following steps:
BI=，
wherein BI is a lamina development index, alpha is a correction coefficient value of 100,for fluctuation step size +.>Is scratch complexity.
2. The shale tattoo development degree characterization method based on the fractal-fluctuation theory as recited in claim 1, wherein the shale scratch test data is obtained through a scratch test, and the scratch test comprises three stages of pre-scanning, engraving scanning and post-scanning; in the pre-scanning stage, scanning the surface of the sample with constant contact force, and recording the change of the scribing depth along with the scribing distance; in the stage of scribing scanning, linearly loading a scribing sample, and recording the change of loading force, friction force and scribing depth along with the scribing distance; in the post-scanning stage, the scratch morphology is scanned with constant contact force, and the change of the scratch depth along with the scratch distance is recorded.
3. The shale tattoo development level characterization method based on fractal-wave theory as recited in claim 1, wherein the equally dividing the scoring depth-distance curve with different threshold values and establishing a scratch fluctuation step size comprises:
performing equivalence division on the ordinate scratch distances of the scratch depth-distance curve by using different threshold values to generate domains with different numbers;
and calculating the number of the resident points of the scribing depth-distance curve contained in each domain, and defining the corresponding threshold value as the scratch fluctuation step length when the duty ratio of the number of the resident points of 1 reaches the maximum value.
4. The shale tattoo development characterization method based on fractal-wave theory as recited in claim 3, further comprising, prior to the performing the equivalence division on the ordinate-type penetration distance of the penetration depth-distance curve with different threshold values:
and amplifying the length of the scribing distance on the ordinate of the scribing depth-distance curve according to a preset proportion, and generating a corrected scribing depth-distance curve while keeping the value of the ordinate unchanged.
5. The method for characterizing the development degree of shale tattoos based on fractal-wave theory as recited in claim 1, wherein defining the development grading index of tattoos to characterize the development degree of tattoos comprises:
when 0< BI is less than or equal to 10, the degree of development is defined as the primary level of development; when 10< BI is less than or equal to 20, the degree of development is defined as the secondary level of development; when 20< BI, the degree of development is defined as tertiary development level.
6. A shale tattoo development level characterization apparatus based on fractal-wave theory for implementing the method as recited in any one of claims 1-5, comprising:
the data acquisition module is used for acquiring shale scratch test data;
the curve drawing module is used for drawing a drawing distance-depth curve and a drawing depth-distance curve respectively by using scratch test data;
the scratch complexity establishing module is used for dividing the scribing distance-depth curve by adopting different scales and establishing scratch complexity;
the fluctuation step length establishing module is used for equally dividing the dividing depth-distance curve by adopting different threshold values and establishing a scratch fluctuation step length;
the layer development degree characterization module is used for calculating a layer development index according to the fluctuation step length and the scratch complexity, and defining a layer development grading index to characterize the layer development degree.
7. A memory, characterized in that a plurality of instructions for implementing the method according to any of claims 1-5 are stored.
8. An electronic device comprising a processor and a memory coupled to the processor, the memory storing a plurality of instructions that are loadable and executable by the processor to enable the processor to perform the method of any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310203793.6A CN116296946B (en) | 2023-03-06 | 2023-03-06 | Shale stratum development degree characterization method and device based on fractal-fluctuation theory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310203793.6A CN116296946B (en) | 2023-03-06 | 2023-03-06 | Shale stratum development degree characterization method and device based on fractal-fluctuation theory |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116296946A CN116296946A (en) | 2023-06-23 |
| CN116296946B true CN116296946B (en) | 2023-08-11 |
Family
ID=86791876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310203793.6A Active CN116296946B (en) | 2023-03-06 | 2023-03-06 | Shale stratum development degree characterization method and device based on fractal-fluctuation theory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116296946B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651966A (en) * | 2016-01-18 | 2016-06-08 | 山东科技大学 | Shale oil and gas high-quality reservoir stratum evaluation method and parameter determination method |
| CN107356489A (en) * | 2017-07-04 | 2017-11-17 | 湖北工业大学 | A kind of test method based on nano indentation test measure granite macromechanical property |
| WO2019223264A1 (en) * | 2018-05-25 | 2019-11-28 | 西南石油大学 | Crushing system for large natural gas hydrate rock sample |
| CN112033866A (en) * | 2020-08-20 | 2020-12-04 | 中国地质调查局油气资源调查中心 | Shale classification method and application thereof and shale lithofacies distribution construction system |
| CN114018697A (en) * | 2021-10-14 | 2022-02-08 | 中国地质大学（武汉） | Method for evaluating compressibility of bedding shale reservoir based on micro-mechanics experiment |
| CN114048640A (en) * | 2022-01-13 | 2022-02-15 | 中国矿业大学(北京) | Method for evaluating shale microcosmic fracture development degree based on nanoindentation experiment |
| CN115078210A (en) * | 2022-05-05 | 2022-09-20 | 中国科学院地质与地球物理研究所 | Shale pore structure testing method |
| CN115170945A (en) * | 2021-03-19 | 2022-10-11 | 中国石油化工股份有限公司 | Quantitative identification method for haircut development degree in continental facies fine-grained sedimentary rock |
| CN115165681A (en) * | 2022-07-05 | 2022-10-11 | 成都理工大学 | Shale reservoir particle structure directional analysis method, system, equipment and terminal |
| CN115356223A (en) * | 2022-10-20 | 2022-11-18 | 中国矿业大学(北京) | Device and method for measuring shale brittleness index continuous section based on high-temperature and high-pressure scratches |
-
2023
- 2023-03-06 CN CN202310203793.6A patent/CN116296946B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651966A (en) * | 2016-01-18 | 2016-06-08 | 山东科技大学 | Shale oil and gas high-quality reservoir stratum evaluation method and parameter determination method |
| CN107356489A (en) * | 2017-07-04 | 2017-11-17 | 湖北工业大学 | A kind of test method based on nano indentation test measure granite macromechanical property |
| WO2019223264A1 (en) * | 2018-05-25 | 2019-11-28 | 西南石油大学 | Crushing system for large natural gas hydrate rock sample |
| CN112033866A (en) * | 2020-08-20 | 2020-12-04 | 中国地质调查局油气资源调查中心 | Shale classification method and application thereof and shale lithofacies distribution construction system |
| CN115170945A (en) * | 2021-03-19 | 2022-10-11 | 中国石油化工股份有限公司 | Quantitative identification method for haircut development degree in continental facies fine-grained sedimentary rock |
| CN114018697A (en) * | 2021-10-14 | 2022-02-08 | 中国地质大学（武汉） | Method for evaluating compressibility of bedding shale reservoir based on micro-mechanics experiment |
| CN114048640A (en) * | 2022-01-13 | 2022-02-15 | 中国矿业大学(北京) | Method for evaluating shale microcosmic fracture development degree based on nanoindentation experiment |
| CN115078210A (en) * | 2022-05-05 | 2022-09-20 | 中国科学院地质与地球物理研究所 | Shale pore structure testing method |
| CN115165681A (en) * | 2022-07-05 | 2022-10-11 | 成都理工大学 | Shale reservoir particle structure directional analysis method, system, equipment and terminal |
| CN115356223A (en) * | 2022-10-20 | 2022-11-18 | 中国矿业大学(北京) | Device and method for measuring shale brittleness index continuous section based on high-temperature and high-pressure scratches |
Non-Patent Citations (1)
| Title |
|---|
| 岩相角度的五峰组页岩小尺度沉积构造系统分类与解析;张妍煜;《非常规油气》;第9卷(第2期);第26-33页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116296946A (en) | 2023-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zargari et al. | Organic maturity, elastic properties, and textural characteristics of self resourcing reservoirs | |
| NO20121394A1 (en) | Procedure for obtaining consistent and integrated physical properties of porous media | |
| SU834514A1 (en) | Method of earth packing quality control | |
| CN106526128B (en) | Method for acquiring strength parameters of fractured rock mass | |
| CN113295561B (en) | Method for evaluating microstructure based on shale nanoindentation curve and electronic equipment | |
| CN103454200A (en) | Method and device for determining pore parameters of rock cores by using computed tomography (CT) technology | |
| CN114460119A (en) | Mechanical parameter acquisition method and device, electronic equipment and computer storage medium | |
| CN116296946B (en) | Shale stratum development degree characterization method and device based on fractal-fluctuation theory | |
| CN112145165A (en) | Dynamic and static permeability conversion method for micro-crack-pore type reservoir | |
| Goldfarb et al. | Segmentationless digital rock physics using different effective medium theories | |
| Su et al. | A comprehensive methodology of evaluation of the fracability of a shale gas play | |
| Bol | The influence of pore pressure gradients in soil classification during piezocone penetration test | |
| CN112711068B (en) | Method and device for predicting effective reservoir of oil gas in sandstone | |
| Bocangel et al. | Acoustic mapping and characterization of organic matter in shales | |
| Moyer et al. | Property changes of oil shale during artificial maturation: The Irati Formation from the Paraná Basin, Brazil | |
| CN114814951A (en) | Method, device, equipment and storage medium for acquiring slowness of fast and slow transverse waves | |
| Aşcı et al. | Correlation of physical and mechanical properties with ultrasonic pulse velocities of sandstones in Çenedağ, Kocaeli-Turkey | |
| CN110258496B (en) | Method for evaluating maximum dry density of calcareous soil based on light dynamic penetration test | |
| CN112346130A (en) | Organic-rich rock transverse wave velocity prediction method, storage medium and system | |
| Horng et al. | Sample Disturbance Effects on Undrained Shear Strengths—Study from Takuhoku Site, Sapporo— | |
| Obaji et al. | Prediction of undrained cohesion obtained from triaxial tests using natural moisture content | |
| Zhao et al. | Stratigraphic sequence and deposition-affected compressibility of fine-grained sediments in the ancient Yellow River Delta during the late Pleistocene and Holocene | |
| Ecemis et al. | Effect of drainage conditions on CPT resistance of silty sand: physical model and field tests | |
| Mukerji et al. | Image Processing of Acoustic Microscopy Data to Estimate Textural Scales and Anistropy in Shales | |
| CN112630019B (en) | Shale brittleness index determination method and device and readable storage medium |
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 |