CN117370718B - Shale brittleness index acquisition method, shale brittleness index acquisition device, shale brittleness index acquisition terminal and shale brittleness index storage medium - Google Patents
Shale brittleness index acquisition method, shale brittleness index acquisition device, shale brittleness index acquisition terminal and shale brittleness index storage medium Download PDFInfo
- Publication number
- CN117370718B CN117370718B CN202311355279.0A CN202311355279A CN117370718B CN 117370718 B CN117370718 B CN 117370718B CN 202311355279 A CN202311355279 A CN 202311355279A CN 117370718 B CN117370718 B CN 117370718B
- Authority
- CN
- China
- Prior art keywords
- wave velocity
- representing
- shale sample
- shale
- acoustic parameters
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 239000000463 material Substances 0.000 claims abstract description 88
- 239000000126 substance Substances 0.000 claims abstract description 32
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 48
- 239000011707 mineral Substances 0.000 claims description 48
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 45
- 239000010453 quartz Substances 0.000 claims description 43
- 229910052683 pyrite Inorganic materials 0.000 claims description 42
- 239000011028 pyrite Substances 0.000 claims description 42
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 42
- 229910021532 Calcite Inorganic materials 0.000 claims description 41
- 239000002734 clay mineral Substances 0.000 claims description 41
- 229910000514 dolomite Inorganic materials 0.000 claims description 41
- 239000010459 dolomite Substances 0.000 claims description 41
- 239000010433 feldspar Substances 0.000 claims description 41
- 239000005416 organic matter Substances 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 20
- 239000011435 rock Substances 0.000 claims description 14
- 238000010606 normalization Methods 0.000 claims description 13
- 239000000203 mixture Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 22
- 239000000523 sample Substances 0.000 description 198
- 239000007789 gas Substances 0.000 description 12
- 239000003245 coal Substances 0.000 description 11
- 238000011156 evaluation Methods 0.000 description 10
- 238000011161 development Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000002441 X-ray diffraction Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000007822 coupling agent Substances 0.000 description 2
- 235000021185 dessert Nutrition 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 239000010445 mica Substances 0.000 description 2
- 229910052618 mica group Inorganic materials 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 239000003079 shale oil Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 229910052925 anhydrite Inorganic materials 0.000 description 1
- 229910052586 apatite Inorganic materials 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- VSIIXMUUUJUKCM-UHFFFAOYSA-D pentacalcium;fluoride;triphosphate Chemical compound [F-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O VSIIXMUUUJUKCM-UHFFFAOYSA-D 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- 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 provides a shale brittleness index acquisition method, a shale brittleness index acquisition device, a shale brittleness index acquisition terminal and a shale brittleness index storage medium. The method comprises the following steps: acquiring acoustic parameters corresponding to a plurality of typical shale samples in a target area layer system and characteristic parameters of different substance components in each typical shale sample; determining acoustic parameters of different material components of the representative shale samples in the target zone layer according to the acoustic parameters of the representative shale samples and the percentage of the different material components in each representative shale sample; and calculating the acoustic parameters of the shale sample to be tested in the target zone layer system based on the acoustic parameters of different material components of the typical shale sample in the target zone layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested. According to the method, acoustic parameters of different material components can be rapidly calculated, the shale brittleness index is obtained based on the elastic wave theory, and the calculation efficiency of the shale brittleness index is greatly improved.
Description
Technical Field
The invention relates to the technical field of shale exploration, in particular to a shale brittleness index acquisition method, device, terminal and storage medium.
Background
Shale is a fine-grained sedimentary rock with main body composed of grains with granularity smaller than 0.0625mm and developed in page, and the rock components mainly comprise clay, quartz, feldspar, carbonate rock, organic matters and the like. Shale formed under different sedimentary conditions has different types and contents of material compositions, so that different types of shale are developed. The shale forms various sedimentary environments, and the organic shale mainly develops in the canopy and large Liu Xiepo of the sea, the deep lake-semi-deep lake of the land, the bay between river diversion, swamps, the front delta of the sea-land transition phase, the mud lawn, the lagoon and the like. For a conventional oil and gas reservoir, shale can play a role in sealing physical properties of a top plate and a bottom plate by the compactness of the shale, and the shale rich in organic matters can play a role in sealing the concentration of hydrocarbons when the shale contains higher hydrocarbon content.
Shale brittleness index is a key index for evaluating the compressibility of shale reservoirs, and higher brittleness index of shale is the manifestation of good compressibility of the shale. The evaluation of shale brittleness index is not only applied to the development engineering evaluation of shale oil gas, but also is involved in the development evaluation of coal bed gas. In the development practice of the coal bed gas for constructing the development of the coal, the development strategy of the fracturing coal bed is often avoided, the exploitation of the coal bed gas is generally realized through fracturing the top and bottom plates, and under the condition, the calculation result of the brittleness index is needed for evaluating the shale compressibility of the top and bottom plates of the coal bed. In the practice of co-exploring and co-mining of three gases (coal bed gas, shale gas and tight sandstone gas) in a coal-based stratum, the brittleness and compressibility of shale are one of key factors for realizing simultaneous transformation of shale, coal bed and sandstone reservoirs, and are also important evaluation contents for evaluating the feasibility of co-mining of the coal-based stratum. Therefore, the accurate acquisition of the shale brittleness index has extremely important scientific value and practical significance for evaluating the feasibility of unconventional oil and gas development.
Currently, the most commonly used methods for obtaining shale brittleness index mainly comprise two types: (1) According to the shale brittleness index calculation method based on the brittle minerals, the shale brittleness index is calculated through analysis of the content of the brittle minerals such as quartz, feldspar, carbonate rock and pyrite in the shale; (2) According to the method, the shale mechanical parameter elastic modulus and the poisson ratio are mainly adopted, and the brittleness index of shale based on the rock mechanical parameter is obtained by normalizing the shale mechanical parameter elastic modulus and the poisson ratio and calculating the average value of the normalized quantity parameter. In the whole, the existing shale brittleness index calculation method generally needs to calculate shale brittleness index based on a large amount of experimental test results or engineering practice data, so that the method has strong dependence on experience of experimenters, needs longer time and is low in efficiency.
Disclosure of Invention
The embodiment of the invention provides a method, a device, a terminal and a storage medium for acquiring a shale brittleness index, which are used for solving the problem of low calculation efficiency of the shale brittleness index.
In a first aspect, an embodiment of the present invention provides a method for obtaining a shale brittleness index, including:
Acquiring acoustic parameters corresponding to a plurality of typical shale samples in a target area layer system and characteristic parameters of different substance components in each typical shale sample; the characteristic parameters comprise percentage content;
determining acoustic parameters of different material components of the representative shale sample in the target zone layer according to the acoustic parameters of a plurality of representative shale samples and the percentage content of the different material components in each representative shale sample based on the principle that the acoustic parameters of the shale sample are controlled by the acoustic parameters of the different material components in the shale sample;
and calculating the acoustic parameters of the shale sample to be tested in the target zone layer system based on the acoustic parameters of different material components of the typical shale sample in the target zone layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested.
In a second aspect, an embodiment of the present invention provides an apparatus for obtaining a shale brittleness index, including:
the data acquisition module is used for acquiring acoustic parameters corresponding to a plurality of typical shale samples in the target area layer system and characteristic parameters of different material components in each typical shale sample; the characteristic parameters comprise percentage content;
The component parameter acquisition module is used for determining the acoustic parameters of the different material components of the typical shale samples in the target area layer according to the acoustic parameters of the plurality of typical shale samples and the percentage content of the different material components in each typical shale sample based on the principle that the acoustic parameters of the shale samples are controlled by the acoustic parameters of the different material components in the shale samples;
the brittleness index acquisition module is used for calculating the acoustic parameters of the shale sample to be tested in the target area layer system based on the acoustic parameters of different substance components of the typical shale sample in the target area layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested.
In a third aspect, an embodiment of the present invention provides a terminal, including a memory, a processor, and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method according to any one of the possible implementations of the first aspect above when the computer program is executed.
In a fourth aspect, embodiments of the present invention provide a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method as described in any one of the possible implementations of the first aspect above.
The embodiment of the invention provides a method, a device, a terminal and a storage medium for acquiring shale brittleness index, wherein the method firstly acquires acoustic parameters corresponding to a plurality of typical shale samples in a target area layer system and characteristic parameters of different substance components in each typical shale sample; then determining the acoustic parameters of the different material components of the typical shale samples in the target zone layer according to the acoustic parameters of a plurality of typical shale samples and the percentage content of the different material components in each typical shale sample based on the principle that the acoustic parameters of the shale samples are controlled by the acoustic parameters of the different material components in the shale samples; and finally, calculating the acoustic parameters of the shale sample to be tested in the target area layer system based on the acoustic parameters of different material components of the typical shale sample in the target area layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested. According to the method, the acoustic parameters of different material components can be rapidly calculated on the basis of obtaining the acoustic parameters of the shale and the percentage contents of the material components, and the brittleness index of the shale to be measured is calculated on the basis of the elastic wave theory, so that the calculation efficiency of the shale brittleness index is greatly improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of an implementation of a method for obtaining a shale brittleness index provided by an embodiment of the invention;
FIG. 2 is a schematic structural diagram of an apparatus for obtaining a shale brittleness index according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a terminal according to an embodiment of the present invention.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth such as the particular system architecture, techniques, etc., in order to provide a thorough understanding of the embodiments of the present invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the following description will be made by way of specific embodiments with reference to the accompanying drawings.
Referring to fig. 1, a flowchart of an implementation of the method for obtaining the shale brittleness index according to the embodiment of the invention is shown, and the details are as follows:
s101: acquiring acoustic parameters corresponding to a plurality of typical shale samples in a target area layer system and characteristic parameters of different substance components in each typical shale sample; the characteristic parameters include percentage content.
In this embodiment, during the stage of collecting the typical shale samples, the collected typical shale samples are required to be free of weathering, and are preferably newly drilled well samples, the number of which is denoted as M. Meanwhile, the deterioration degree, the sedimentation environment, the construction condition and other geological background conditions of a plurality of typical shale samples are similar, and the shale samples which are in the same area and the same layer and are not influenced by special geological events or phenomena (such as magma activities, formation water and the like) are better, so that reasonable condition constraint is provided for the application of a subsequent method.
After the typical shale samples are collected, XRD (X-ray diffraction) is firstly carried out on all the collected typical shale samples, and component analysis is carried out on the typical shale samples to obtain a substance component; and then performing TOC (Total Organic Carbon, total organic carbon content) analysis on the typical shale sample to obtain the percentage content of each substance component in the typical shale sample, and finally performing density analysis on the typical shale sample to determine the density of each substance component. It should be noted that the percentages of the material components mentioned in the present application are mass percentages.
In this embodiment, the material composition includes quartz, feldspar, clay minerals, calcite, dolomite, pyrite, accessory minerals, organic matter, and the like. The percentage of material components of different shale samples in the same target zone are different. The species of material composition of shale samples in different target zones also vary. The present example is explained below taking the general composition of matter present in shale (including quartz, feldspar, clay minerals, calcite, dolomite, pyrite, by-minerals and organic matter) as an example.
Specifically, the mineral components in different shale samples have larger differences, for example, dolomite or calcite in some lake-phase shale is relatively developed, but the quartz content in the shallow lake-phase is relatively high, and all brittle minerals are obviously unreasonable in a summary; in addition, the rock density and acoustic parameters of certain minerals can change under different conditions, and particularly, the acoustic parameters of soluble minerals (such as calcite, dolomite, quartz and the like) can be different according to the difference of the development degree of the corrosion holes; secondly, organic matters serving as special objects in the shale rich in organic matters can influence the density and acoustic parameter characteristics of the shale to a great extent, and the density or acoustic parameter is also commonly used in practice for predicting the organic matter content in the shale. Thus, considering the differences in acoustic parameters of different minerals, the embodiment finely considers the contribution of various different mineral components to the acoustic parameters of the shale to more accurately reflect the overall acoustic characteristics of the shale.
S102: based on the principle that the acoustic parameters of the shale samples are controlled by the acoustic parameters of different material components in the shale samples, the acoustic parameters of the different material components of the typical shale samples in the target zone layer are determined according to the acoustic parameters of a plurality of typical shale samples and the percentage content of the different material components in each typical shale sample.
In the present embodiment, the acoustic parameters include a shear wave velocity and a longitudinal wave velocity. In order to obtain the transverse wave velocity and the longitudinal wave velocity of different material components in shale, the percentage content of the material components of a selected typical shale sample is firstly combed. The shale contains mainly quartz, feldspar, clay minerals, calcite, dolomite, pyrite, organic matters, other auxiliary minerals and the like, and the percentages of the components in different shale samples are different.
It will be appreciated that the sum of the percentages of the individual material components in a typical shale sample is equal to 100%. The method comprises the following steps:
and is also provided with
w Q-i +w F-i +w C-i +w Ca-i +w D-i +w P-i +w AM-i +w OM-i ＝100％
Wherein w is Q-i Representing the percentage of quartz in an ith typical shale sample; w (w) F-i Representing the percentage of feldspar in the ith typical shale sample; w (w) C-i Representing the percentage of clay minerals in an ith typical shale sample; w (w) Ca-i Representing the percentage of calcite in the ith typical shale sample; w (w) D-i Representing the percentage of dolomite in an ith typical shale sample; w (w) P-i Representing the percentage of pyrite in the ith typical shale sample; w (w) AM-i Representing the percentage of secondary minerals in an ith typical shale sample; w (w) OM-i Representing the percentage of organic matter in an ith typical shale sample.
Different material components in shaleThe wave velocities of the shale are different, so that the wave velocities of different components are needed to be analyzed for better understanding and analyzing the wave velocity and brittleness characteristics of the shale, and a longitudinal wave velocity matrix I of different material components in the shale vp Can be expressed as:
wherein v is p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; v p-D Representing the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is shown.
Transverse wave velocity matrix I of different material components in shale vs Can be expressed as:
wherein v is s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; v s-D Represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
The wave velocity of each material component of shale is closely related to the density thereof, and the wave velocity corresponding to the higher the density of the component is also higher. Generally, the densities of quartz, feldspar, clay minerals, calcite, dolomite, pyrite and organic matters are respectively 2.65, 2.55-2.75, 2.00-2.9, 2.71, 2.8-2.9, 4.9-5.2 and 1.3g/cm < 3 >, and the pyrite is found to be a heavy mineral, the density of the pyrite is far higher than that of other components, and the density of the organic matters is far lower than that of other components, and the densities of quartz, feldspar, clay minerals, mica, calcite and dolomite are relatively close. Other accessory minerals (such as apatite, gypsum and anhydrite) have a lower density than pyrite and higher density than organic matter. Notably, clay minerals have a severely reduced density after water and shale with a high clay mineral content has a poor brittleness. Thus, pyrite has a significantly higher wave velocity than other components, while organic matter has a much lower wave velocity than other components, quartz, feldspar, clay minerals, calcite, dolomite, and a relatively close wave velocity.
In this embodiment, it may be determined based on the physical properties of the shale that the sum of the products of the acoustic parameters and percentages corresponding to the different material components in the same typical shale sample approximates the acoustic parameters of the typical shale sample.
Based on the principle, the embodiment collects at least 30 typical shale samples from the rock stratum in the same area, and selects 80% of the typical shale samples to construct an equation set to obtain the longitudinal wave velocity and the transverse wave velocity of the typical shale samples. The selected typical shale sample simultaneously needs to meet the condition that the length is not less than 5cm, and no artificial damage phenomenon exists. And (3) obtaining a columnar sample with the diameter of 2.5cm and the length of not less than 5cm by adopting a linear cutting technology, drying a typical shale sample at the temperature of 90 ℃ for 12 hours before an experiment, and accurately measuring the length of the sample by using a vernier caliper and marking the length as L (in cm).
Specifically, the selected representative shale samples are tested for longitudinal and transverse wave velocities using the rock acoustic parameter testing system for all representative shale samples therein. To meet the requirements of accurate testing and eliminate the prominence of shale mineral particle scattering, the selected ultrasonic wave wavelength should be less than one tenth of the sample length and greater than 0.5mm. Before testing, a layer of coupling agent is coated on the probes, then the probes are closed and pressed tightly, and the propagation time t of sound waves between the probes can be obtained by adjusting the waveform 0 (in s). The transmitting probe is placed at both ends of the sample and a coupling agent is applied between the sample and the probe. When the ultrasonic transmitting probe transmits sound wave information and passes through the sample, the ultrasonic receiving probe receives the sound wave information, and the time scale is adjusted to the initial motion of the waveform, so that the longitudinal wave measurement time t can be obtained p (in s). By adjusting the bandwidth and attenuation of the acoustic wave deviceThe button is subtracted, the transverse wave is identified on the display, the time scale is adjusted to the point of the jump of the transverse wave, and the measuring time t of the transverse wave can be read s (in s). To ensure accuracy and reliability in acquiring the wave velocity data, each representative shale sample is tested at least three times to take its average as a measurement.
The calculation formulas of the longitudinal wave velocity and the transverse wave velocity of a typical sample are as follows:
wherein v is p Representing longitudinal wave velocity, v s Representing the transverse wave velocity, wherein the unit is m/s; l is the length of the shale sample in cm; t is t p And t s The unit is us respectively for longitudinal wave measurement time and transverse wave measurement time; t is t 0 The unit is us for travel time between probes.
Specifically, the specific implementation flow of S102 includes:
s201: for each representative shale sample, constructing an equation with the acoustic parameters of the individual material components of the representative shale sample as unknowns and the acoustic parameters of the representative shale sample and the percentages of the individual material components of the representative shale sample as known quantities;
S202: and solving equations corresponding to the plurality of typical shale samples to obtain acoustic parameters of each material component of the typical shale samples in the target zone strata.
In one possible embodiment, the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the equation corresponding to the longitudinal wave velocity is as follows:
wherein v is p-i Representing the longitudinal wave velocity of an ith typical shale sample; w (w) Q-i Representing the percentage of quartz in an ith typical shale sample; w (w) F-i Representing the percentage of feldspar in the ith typical shale sample; w (w) C-i Represents the ithThe percentage of clay minerals in a typical shale sample; w (w) Ca-i Representing the percentage of calcite in the ith typical shale sample; w (w) D-i Representing the percentage of dolomite in an ith typical shale sample; w (w) P-i Representing the percentage of pyrite in the ith typical shale sample; w (w) AM-i Representing the percentage of secondary minerals in an ith typical shale sample; w (w) OM-i Representing the percentage of organic matter in an ith typical shale sample; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; v p-D Representing the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
the equation corresponding to the transverse wave velocity is as follows:
wherein v is s-i Representing the shear wave velocity of an ith typical shale sample; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; v s-D Represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
In this embodiment, the wave velocity of the shale is controlled by the wave velocity of each material component, so an equation set with the longitudinal wave velocity of different material components as an unknown quantity and the acoustic parameters of a typical shale sample and the percentage content of different material components as a known quantity can be constructed, and then the wave velocities of different material components of the typical shale sample can be obtained by solving the equation set.
In one possible embodiment, the characteristic parameter comprises density; the specific implementation flow of S202 includes:
combining equations corresponding to a plurality of typical shale samples to form an equation set; taking the sum of the percentages of all the substance components of the same typical shale sample as 100%, taking the acoustic parameters of all the substance components as non-negative numbers, and taking the acoustic parameters of the substance components with relatively higher density in any two substance components as constraint conditions, wherein the acoustic parameters of the substance components with relatively lower density are larger than those of the substance components with relatively lower density;
And solving the equation set by adopting a least square method to obtain acoustic parameters of different substance components in the typical shale sample.
In this embodiment, the equation set of the longitudinal wave velocity solved by the least square method is:
the system of equations for transverse wave velocity is:
wherein,transposed matrix representing longitudinal wave velocities of different material components +.>A transposed matrix representing transverse wave velocities of different material components; n is the number of typical shale samples that construct the equation set.
The constraints of the equation set include:
1. the sum of the percentages of all the material components of the same typical shale sample is 100%, namely:
w Q-i +w F-i +w C-i +w Ca-i +w D-i +w P-i +w AM-i +w OM-i ＝100％；
2. all of the material componentsThe acoustic parameters are all non-negative numbers, namely: i vs ≥0，I ps ≥0；
3. The acoustic parameters of the relatively denser material component of any two material components are greater than the acoustic parameters of the relatively less dense material component.
After acoustic parameters of different material components in a typical shale sample are calculated, in order to verify the accuracy of the acoustic parameters, in the embodiment, acoustic parameters of the remaining 20% of the typical shale sample are calculated by adopting the calculated acoustic parameters of each material component, then the acoustic parameters calculated by the corresponding typical shale sample are verified by adopting the acoustic parameters measured by experiments of the remaining 20% of the typical shale sample, and the calculation error of the acoustic parameters is determined. If the calculation error of the acoustic parameters exceeds 15%, the number of the typical shale samples is increased, the number of equations in the equation set is increased, the equation set is solved again by adopting the specific method of S202, and the process is circulated until the error between the calculated acoustic parameters of the typical shale samples and the actual acoustic parameters of the corresponding typical shale samples is less than 15%.
From the above embodiments, it can be seen that, according to the embodiment of the present invention, based on the acoustic parameters actually measured by the typical shale sample, an equation set of longitudinal wave velocities and transverse wave velocities of acoustic parameters of different material components is established, and in the process of solving the equation set, actual situations such as non-negativity of wave velocities, approximate relation between wave velocities of different material components, and sum of mineral component contents being 100% are considered as constraint conditions. The acoustic parameters of each shale sample can be more accurately calculated. According to the processing mode, the fact that the acoustic parameters of shale samples under different geological backgrounds are different is considered, so that the shale acoustic characteristics in actual conditions can be reflected better based on the acoustic parameter solving of shale samples at the same horizon in the same region, and the processing mode is more specific and adaptive.
S103: and calculating the acoustic parameters of the shale sample to be tested in the target zone layer system based on the acoustic parameters of different material components of the typical shale sample in the target zone layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested.
In one possible implementation manner, the implementation procedure for calculating the acoustic parameters of the shale sample to be measured in S103 includes:
Based on the formulaCalculating to obtain the longitudinal wave velocity of the shale sample to be measured;
wherein v is p-x Representing the longitudinal wave velocity of a shale sample x to be measured; w (w) Q-x Representing the percentage content of quartz in the shale sample x to be measured; w (w) F-x Representing the percentage content of feldspar in the shale sample x to be tested; w (w) C-x Representing the percentage content of clay minerals in the shale sample x to be measured; w (w) Ca-x Representing the percentage content of calcite in the shale sample x to be tested; w (w) D-x Representing the percentage content of dolomite in the shale sample x to be measured; w (w) P-x Representing the percentage content of pyrite in the shale sample x to be tested; w (w) AM-x Representing the percentage content of accessory minerals in the shale sample x to be measured; w (w) OM-x Representing the percentage content of organic matters in the shale sample x to be measured; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; v p-D Representing the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
based on the formulaCalculating to obtain the transverse wave velocity of the shale sample to be measured;
wherein v is s-x Representing the longitudinal wave velocity of a shale sample x to be measured; w (w) Q-x Representing the percentage content of quartz in the shale sample x to be measured; w (w) F-x Representing the percentage content of feldspar in the shale sample x to be tested; w (w) C-x Representing the percentage content of clay minerals in the shale sample x to be measured; w (w) Ca-x Representing the percentage content of calcite in the shale sample x to be tested; w (w) D-x Representation of dolomite in shale sample x to be testedIs a percentage of (1); w (w) P-x Representing the percentage content of pyrite in the shale sample x to be tested; w (w) AM-x Representing the percentage content of accessory minerals in the shale sample x to be measured; w (w) OM-x Representing the percentage content of organic matters in the shale sample x to be measured; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; v s-D Represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
In one possible embodiment, the specific implementation procedure for calculating the brittleness index includes:
s301: calculating Young modulus and Poisson ratio of the shale sample to be tested according to the acoustic parameters of the shale sample to be tested;
s302: and calculating the brittleness index of the shale sample to be tested based on the Young modulus and the Poisson ratio of the shale sample to be tested.
Specifically, the mutual calculation of the acoustic parameters and the rock mechanical parameters of the rock can be realized by the elastic wave theory. Longitudinal wave propagation velocity v of shale sample x to be measured p-x And transverse wave propagation velocity v s-x Can be obtained from the mechanical parameters of rock modulus of elasticity E x And Poisson ratio u x Expressed as:
wherein v is p-x Represents the longitudinal wave velocity of the shale sample x to be measured, and the unit is m/s, v s-x The transverse wave velocity of the shale sample x to be measured is m/s; e (E) x The elastic modulus of the shale sample x to be measured is expressed in GPa; mu (mu) x The poisson ratio of the shale sample x to be measured is represented, and the shale sample x is dimensionless; ρ x The rock density of the shale sample x to be measured is expressed in kg/m 3 ；
The calculating the brittleness index of the shale sample to be measured based on the Young modulus and the Poisson ratio of the shale sample to be measured comprises:
carrying out normalization treatment on Young modulus and Poisson ratio of the shale sample to be measured;
based on the formulaCalculating the brittleness index of the shale sample to be measured;
wherein BI x Indicating the brittleness index, E, of the shale sample x to be tested Brit-x Represents the Young modulus and mu of the shale sample x to be measured after normalization Brit-x And (5) representing the poisson ratio of the shale sample x to be tested after normalization.
In one possible embodiment, the normalization process includes:
based on the formulaCarrying out normalization treatment on Young modulus and Poisson ratio of the shale sample to be measured;
wherein E is Brit-x Represents the Young modulus and mu of the shale sample x to be measured after normalization Brit-x Represents the Poisson's ratio of the shale sample x to be measured after normalization, E x The elastic modulus of the shale sample x to be measured is represented; mu (mu) x The poisson ratio of the shale sample x to be measured is represented; e (E) min Represents the minimum value of Young's modulus in all shale samples to be tested, E max Represents the Young's modulus maximum value, mu, in all shale samples to be tested min Represents the minimum value of poisson ratio, mu, in all shale samples to be tested max The poisson's ratio maximum in all shale samples tested is shown.
According to the embodiment of the invention, the acoustic parameters of the shale sample are firstly determined through different material components (comprising quartz, feldspar, clay minerals, mica, calcite, dolomite, pyrite, accessory minerals and organic matters), then the corresponding shale mechanical parameters are calculated according to the acoustic parameters, and finally the integral brittleness index of the shale sample is calculated through the calculated shale mechanical parameters, so that the dependence of the whole calculation process on human experience can be reduced, the evaluation efficiency of the brittleness index is greatly improved, and the input cost is reduced; the method provided by the embodiment can fully consider the action of different substance components on the shale brittleness, avoid the limitation of a single evaluation parameter, and integrally improve the scientificity, rationality and credibility of the shale brittleness calculation result.
Furthermore, the technical scheme provided by the invention not only can be applied to calculation of the brittleness index of the shale sample in the core scale, but also has strong applicability to evaluation of the shale brittleness of the whole well section, and further can provide basis and direction for optimization and development layout of the favorable interval of shale oil gas development.
Furthermore, the method is not only suitable for calculating the brittleness index of shale, but also has strong heuristics and comparability to the brittleness index calculation method of reservoirs such as tight sandstone, coal seam and the like, and has wider extension application range and field.
Furthermore, the method can perform large-scale and efficient evaluation on the shale brittleness index of the same shale layer in the same evaluation area or in the areas with comparable geological conditions on the basis of obtaining a certain amount of typical shale sample acoustic parameters, and simultaneously can perform organic shale rich in organic matters, gas content or oil content distribution prediction, dessert section or dessert area evaluation and the like by utilizing the acoustic parameters obtained through calculation.
It should be understood that the sequence number of each step in the foregoing embodiment does not mean that the execution sequence of each process should be determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present invention.
The following are device embodiments of the invention, for details not described in detail therein, reference may be made to the corresponding method embodiments described above.
Fig. 2 is a schematic structural diagram of an apparatus for obtaining a shale brittleness index according to an embodiment of the present invention, and for convenience of explanation, only a portion related to the embodiment of the present invention is shown, which is described in detail below:
as shown in fig. 2, the shale brittleness index obtaining apparatus 100 includes:
the data acquisition module 110 is configured to acquire acoustic parameters corresponding to a plurality of typical shale samples in the target zone layer, and characteristic parameters of different material components in each typical shale sample; the material components comprise quartz, feldspar, clay minerals, calcite, dolomite, pyrite, accessory minerals and organic matters; the characteristic parameters comprise percentage content;
a component parameter obtaining module 120, configured to determine acoustic parameters of different material components of the typical shale samples in the target zone layer according to the acoustic parameters of the plurality of typical shale samples and the percentage content of the different material components of each typical shale sample based on a principle that the acoustic parameters of the shale samples are controlled by the acoustic parameters of the different material components of the shale samples;
The brittleness index obtaining module 130 is configured to calculate acoustic parameters of a shale sample to be tested in the target zone layer system based on acoustic parameters of different material components of a typical shale sample in the target zone layer system, and determine the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested.
In one possible implementation, the component parameter acquisition module 120 includes:
an equation construction unit for constructing, for each representative shale sample, an equation having an acoustic parameter of each of the material components of the representative shale sample as an unknown quantity and having an acoustic parameter of the representative shale sample and a percentage content of each of the material components of the representative shale sample as a known quantity;
and the acoustic parameter calculation unit is used for solving the acoustic parameters of each material component of the typical shale samples in the target area layer system by combining equations corresponding to the plurality of typical shale samples.
In one possible embodiment, the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the equation corresponding to the longitudinal wave velocity is as follows:
wherein v is p-i Representing the longitudinal wave velocity of an ith typical shale sample; w (w) Q-i Representing the percentage of quartz in an ith typical shale sample; w (w) F-i Representing the percentage of feldspar in the ith typical shale sample; w (w) C-i Representing the percentage of clay minerals in an ith typical shale sample; w (w) Ca-i Representing the percentage of calcite in the ith typical shale sample; w (w) D-i Representing the percentage of dolomite in an ith typical shale sample; w (w) P-i Representing the percentage of pyrite in the ith typical shale sample; w (w) AM-i Representing the percentage of secondary minerals in an ith typical shale sample; w (w) OM-i Representing the percentage of organic matter in an ith typical shale sample; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; v p-D Representing the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
the equation corresponding to the transverse wave velocity is as follows:
wherein v is s-i Representing the shear wave velocity of an ith typical shale sample; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; v s-D Represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
In one possible embodiment, the characteristic parameter comprises density;
the acoustic parameter calculation unit includes:
combining equations corresponding to a plurality of typical shale samples to form an equation set; taking the sum of the percentages of all the substance components of the same typical shale sample as 100%, taking the acoustic parameters of all the substance components as non-negative numbers, and taking the acoustic parameters of the substance components with relatively higher density in any two substance components as constraint conditions, wherein the acoustic parameters of the substance components with relatively lower density are larger than those of the substance components with relatively lower density;
and solving the equation set by adopting a least square method to obtain acoustic parameters of different substance components in the typical shale sample.
In one possible embodiment, the brittleness index obtaining module includes:
the acoustic parameter calculation unit of the sample to be measured is used for calculating the acoustic parameter of the sample to be measured based on the formulaCalculating to obtain the longitudinal wave velocity of the shale sample to be measured;
wherein v is p-x Representing the longitudinal wave velocity of a shale sample x to be measured; w (w) Q-x Representing the percentage content of quartz in the shale sample x to be measured; w (w) F-x Representing the percentage content of feldspar in the shale sample x to be tested; w (w) C-x Representing the percentage content of clay minerals in the shale sample x to be measured; w (w) Ca-x Representing the percentage content of calcite in the shale sample x to be tested; w (w) D-x Representing the percentage content of dolomite in the shale sample x to be measured; w (w) P-x Representing the percentage content of pyrite in the shale sample x to be tested; w (w) AM-x Representing the percentage content of accessory minerals in the shale sample x to be measured; w (w) OM-x Representing the percentage content of organic matters in the shale sample x to be measured; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; v p-D Representing the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
based on the formulaCalculating to obtain the transverse wave velocity of the shale sample to be measured;
wherein v is s-x Representing the longitudinal wave velocity of a shale sample x to be measured; w (w) Q-x Representing the percentage content of quartz in the shale sample x to be measured; w (w) F-x Representing the percentage content of feldspar in the shale sample x to be tested; w (w) C-x Representing the percentage content of clay minerals in the shale sample x to be measured; w (w) Ca-x Representing the percentage content of calcite in the shale sample x to be tested; w (w) D-x Representing the percentage content of dolomite in the shale sample x to be measured; w (w) P-x Representing the percentage content of pyrite in the shale sample x to be tested; w (w) AM-x Representing the percentage content of accessory minerals in the shale sample x to be measured; w (w) OM-x Representing the percentage content of organic matters in the shale sample x to be measured; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; v s-D Represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
In one possible embodiment, the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the brittleness index obtaining module further comprises a brittleness index calculating unit for:
the mechanical parameter calculation subunit is used for calculating Young modulus and Poisson ratio of the shale sample to be measured according to the acoustic parameters of the shale sample to be measured;
and the brittleness index calculating subunit is used for calculating the brittleness index of the shale sample to be measured based on the Young modulus and the Poisson ratio of the shale sample to be measured.
In one possible embodiment, the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the mechanical parameter calculation subunit includes:
based on the formulaCalculating Young modulus and Poisson ratio of the shale sample to be tested;
wherein v is p-x Representing the longitudinal wave velocity, v, of a shale sample x to be measured s-x The transverse wave velocity of the shale sample x to be measured; e (E) x The elastic modulus of the shale sample x to be measured is represented; mu (mu) x The poisson ratio of the shale sample x to be measured is represented; ρ x Representing the rock density of the shale sample x to be measured;
the brittleness index calculation subunit includes:
carrying out normalization treatment on Young modulus and Poisson ratio of the shale sample to be measured;
based on the formulaCalculating the brittleness index of the shale sample to be measured;
wherein BI x Indicating the brittleness index, E, of the shale sample x to be tested Brit-x Represents the Young modulus and mu of the shale sample x to be measured after normalization Brit-x And (5) representing the poisson ratio of the shale sample x to be tested after normalization.
According to the shale brittleness index acquisition device provided by the embodiment of the invention, characteristic parameters of substance components corresponding to a plurality of shale samples in the same regional layer are acquired firstly; the material components comprise quartz, feldspar, clay minerals, calcite, dolomite, pyrite, accessory minerals and organic matters; then determining the acoustic parameters of the different material components of the typical shale samples in the target zone layer according to the acoustic parameters of the plurality of typical shale samples and the percentage of the different material components in each typical shale sample; and calculating the acoustic parameters of the shale sample to be tested in the target zone layer system based on the acoustic parameters of different material components of the typical shale sample in the target zone layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested. According to the device, on the basis of obtaining typical shale acoustic parameters and the percentage contents of the material components thereof, the acoustic parameters of different material components can be rapidly calculated, and the shale brittleness index to be measured is calculated based on the elastic wave theory, so that the calculation efficiency of the shale brittleness index is greatly improved.
Fig. 3 is a schematic diagram of a terminal according to an embodiment of the present invention. As shown in fig. 3, the terminal 3 of this embodiment includes: a processor 30 and a memory 31. The memory 31 is used for storing a computer program 32, and the processor 30 is used for calling and running the computer program 32 stored in the memory 31 to execute the steps in the above-mentioned method embodiment for obtaining the shale brittleness index, such as the steps S101 to S103 shown in fig. 1. Alternatively, the processor 30 is configured to invoke and run the computer program 32 stored in the memory 31 to implement the functions of the modules/units in the above-described device embodiments, such as the functions of the modules 110 to 130 shown in fig. 2.
Illustratively, the computer program 32 may be partitioned into one or more modules/units that are stored in the memory 31 and executed by the processor 30 to complete the present invention. The one or more modules/units may be a series of computer program instruction segments capable of performing specific functions for describing the execution of the computer program 32 in the terminal 3.
The terminal 3 may be a computing device such as a desktop computer, a notebook computer, a palm computer, a cloud server, etc. The terminal 3 may include, but is not limited to, a processor 30, a memory 31. It will be appreciated by those skilled in the art that fig. 3 is merely an example of the terminal 3 and does not constitute a limitation of the terminal 3, and may include more or less components than illustrated, or may combine certain components, or different components, e.g., the terminal may further include an input-output device, a network access device, a bus, etc.
The processor 30 may be a central processing unit (Central Processing Unit, CPU), other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field-programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 31 may be an internal storage unit of the terminal 3, such as a hard disk or a memory of the terminal 3. The memory 31 may be an external storage device of the terminal 3, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) or the like, which are provided on the terminal 3. Further, the memory 31 may also include both an internal storage unit and an external storage device of the terminal 3. The memory 31 is used for storing the computer program as well as other programs and data required by the terminal. The memory 31 may also be used for temporarily storing data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and in part, not described or illustrated in any particular embodiment, reference is made to the related descriptions of other embodiments.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus/terminal and method may be implemented in other manners. For example, the apparatus/terminal embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via interfaces, devices or units, which may be in electrical, mechanical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated modules/units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present invention may be implemented in whole or in part by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, where the computer program, when executed by a processor, may implement the steps of the method embodiments for obtaining the shale brittleness index. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. It should be noted that the computer readable medium may include content that is subject to appropriate increases and decreases as required by jurisdictions in which such content is subject to legislation and patent practice, such as in certain jurisdictions in which such content is not included as electrical carrier signals and telecommunication signals.
The above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention.
Claims (9)
1. A method for obtaining a shale brittleness index, comprising:
acquiring acoustic parameters corresponding to a plurality of typical shale samples in a target area layer system and characteristic parameters of different substance components in each typical shale sample; the characteristic parameters comprise percentage content;
determining acoustic parameters of different material components of the representative shale sample in the target zone layer according to the acoustic parameters of a plurality of representative shale samples and the percentage content of the different material components in each representative shale sample based on the principle that the acoustic parameters of the shale sample are controlled by the acoustic parameters of the different material components in the shale sample;
Calculating acoustic parameters of a shale sample to be tested in the target zone layer system based on acoustic parameters of different material components of a typical shale sample in the target zone layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested;
wherein the material components comprise quartz, feldspar, clay minerals, calcite, dolomite, pyrite, accessory minerals and organic matters; the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the equation corresponding to the longitudinal wave velocity is as follows:
wherein v is p-i Representing the longitudinal wave velocity of an ith typical shale sample; w (w) Q-i Representing the percentage of quartz in an ith typical shale sample; w (w) F-i Representing the percentage of feldspar in the ith typical shale sample; w (w) C-i Representing the percentage of clay minerals in an ith typical shale sample; w (w) Ca-i Representing the percentage of calcite in the ith typical shale sample; w (w) D-i Representing the percentage of dolomite in an ith typical shale sample; w (w) P-i Representing the percentage of pyrite in the ith typical shale sample; w (w) AM-i Representing the percentage of secondary minerals in an ith typical shale sample; w (w) OM-i Representing the percentage of organic matter in an ith typical shale sample; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; vp-D represents the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
the equation corresponding to the transverse wave velocity is as follows:
wherein v is s-i Representing the shear wave velocity of an ith typical shale sample; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; vs-D represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
2. The method of claim 1, wherein determining the acoustic parameters of the different material components of the representative shale sample in the target zone layer based on the principle that the acoustic parameters of the shale sample are controlled by the acoustic parameters of the different material components of the representative shale sample according to the percentage of the acoustic parameters of the representative shale samples to the different material components of the respective representative shale samples comprises:
For each representative shale sample, constructing an equation with the acoustic parameters of the individual material components of the representative shale sample as unknowns and the acoustic parameters of the representative shale sample and the percentages of the individual material components of the representative shale sample as known quantities;
and solving equations corresponding to the plurality of typical shale samples to obtain acoustic parameters of each material component of the typical shale samples in the target zone strata.
3. The method of claim 2, wherein the characteristic parameter comprises density;
the equation corresponding to the plurality of typical shale samples is combined, and the acoustic parameters of each material component of the typical shale samples in the target area layer system are obtained by solving the equation, wherein the equation comprises the following components:
combining equations corresponding to a plurality of typical shale samples to form an equation set; taking the sum of the percentages of all the substance components of the same typical shale sample as 100%, taking the acoustic parameters of all the substance components as non-negative numbers, and taking the acoustic parameters of the substance components with relatively higher density in any two substance components as constraint conditions, wherein the acoustic parameters of the substance components with relatively lower density are larger than those of the substance components with relatively lower density;
and solving the equation set by adopting a least square method to obtain acoustic parameters of different substance components in the typical shale sample.
4. The method of claim 1, wherein the composition of matter comprises quartz, feldspar, clay minerals, calcite, dolomite, pyrite, accessory minerals, and organic matter; the calculating the acoustic parameters of the shale sample to be measured in the target zone layer system based on the acoustic parameters of different material components of the typical shale sample in the target zone layer system comprises the following steps:
based on the formulaCalculating to obtain the longitudinal wave velocity of the shale sample to be measured;
wherein v is p-x Representing the longitudinal wave velocity of a shale sample x to be measured; w (w) Q-x Representing the percentage content of quartz in the shale sample x to be measured; w (w) F-x Representing the percentage content of feldspar in the shale sample x to be tested; w (w) C-x Representing the percentage content of clay minerals in the shale sample x to be measured; w (w) Ca-x Representing the percentage content of calcite in the shale sample x to be tested; wD-x represents the percentage content of dolomite in the shale sample x to be measured; w (w) P-x Representing the percentage content of pyrite in the shale sample x to be tested; w (w) AM-x Representing the percentage content of accessory minerals in the shale sample x to be measured; w (w) OM-x Representing the percentage content of organic matters in the shale sample x to be measured; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; vp-D represents the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
based on the formulaCalculating to obtain the transverse wave velocity of the shale sample to be measured;
wherein v is s-x Representing the longitudinal wave velocity of a shale sample x to be measured; w (w) Q-x Representing the percentage content of quartz in the shale sample x to be measured; w (w) F-x Representing the percentage content of feldspar in the shale sample x to be tested; w (w) C-x Representing the percentage content of clay minerals in the shale sample x to be measured; w (w) Ca-x Representing the percentage content of calcite in the shale sample x to be tested; wD-x represents the percentage content of dolomite in the shale sample x to be measured; w (w) P-x Representing the percentage content of pyrite in the shale sample x to be tested; w (w) AM-x Representing the percentage content of accessory minerals in the shale sample x to be measured; w (w) OM-x Representing the percentage content of organic matters in the shale sample x to be measured; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; vs-D represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
5. The method of claim 1, wherein the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested comprises:
calculating Young modulus and Poisson ratio of the shale sample to be tested according to the acoustic parameters of the shale sample to be tested;
and calculating the brittleness index of the shale sample to be tested based on the Young modulus and the Poisson ratio of the shale sample to be tested.
6. The method of claim 1, wherein the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
calculating young modulus and poisson ratio of the shale sample to be measured according to the acoustic parameters of the shale sample to be measured, wherein the calculating comprises the following steps:
based on the formulaCalculating Young modulus and Poisson ratio of the shale sample to be tested;
wherein v is p-x Representing the longitudinal wave velocity, v, of a shale sample x to be measured s-x The transverse wave velocity of the shale sample x to be measured; e (E) x The elastic modulus of the shale sample x to be measured is represented; mu (mu) x The poisson ratio of the shale sample x to be measured is represented; ρ x Representing the rock density of the shale sample x to be measured;
The calculating the brittleness index of the shale sample to be measured based on the Young modulus and the Poisson ratio of the shale sample to be measured comprises:
carrying out normalization treatment on Young modulus and Poisson ratio of the shale sample to be measured;
based on the formulaCalculating the brittleness index of the shale sample to be measured;
wherein BI x Indicating the brittleness index, E, of the shale sample x to be tested Brit-x Represents the Young modulus and mu of the shale sample x to be measured after normalization Brit-x Represents the x normalization of shale samples to be measuredPoisson's ratio after conversion.
7. An acquisition device of shale brittleness index, characterized by comprising:
the data acquisition module is used for acquiring acoustic parameters corresponding to a plurality of typical shale samples in the target area layer system and characteristic parameters of different material components in each typical shale sample; the characteristic parameters comprise percentage content;
the component parameter acquisition module is used for determining the acoustic parameters of the different material components of the typical shale samples in the target area layer according to the acoustic parameters of the plurality of typical shale samples and the percentage content of the different material components in each typical shale sample based on the principle that the acoustic parameters of the shale samples are controlled by the acoustic parameters of the different material components in the shale samples;
The brittleness index acquisition module is used for calculating acoustic parameters of the shale sample to be tested in the target area layer system based on acoustic parameters of different material components of the typical shale sample in the target area layer system, and determining the brittleness index of the shale sample to be tested based on the acoustic parameters of the shale sample to be tested;
wherein the material components comprise quartz, feldspar, clay minerals, calcite, dolomite, pyrite, accessory minerals and organic matters; the acoustic parameters include longitudinal wave velocity and transverse wave velocity;
the equation corresponding to the longitudinal wave velocity is as follows:
wherein v is p-i Representing the longitudinal wave velocity of an ith typical shale sample; w (w) Q-i Representing the percentage of quartz in an ith typical shale sample; w (w) F-i Representing the percentage of feldspar in the ith typical shale sample; w (w) C-i Representing the percentage of clay minerals in an ith typical shale sample; w (w) Ca-i Representing the percentage of calcite in the ith typical shale sample; w (w) D-i Representing the percentage of dolomite in an ith typical shale sample; w (w) P-i Representing the percentage of pyrite in the ith typical shale sample; w (w) AM-i Representing the percentage of secondary minerals in an ith typical shale sample; w (w) OM-i Representing the percentage of organic matter in an ith typical shale sample; v p-Q Representing the longitudinal wave velocity of quartz; v p-F Representing the longitudinal wave velocity of feldspar; v p-C Representing the longitudinal wave velocity of clay minerals; v p-Ca Representing the longitudinal wave velocity of calcite; vp-D represents the longitudinal wave velocity of dolomite; v p-P Representing the longitudinal wave velocity of pyrite; v p-AM The longitudinal wave velocity of the accessory mineral; v p-OM The longitudinal wave velocity of the organic matter is represented;
the equation corresponding to the transverse wave velocity is as follows:
wherein v is s-i Representing the shear wave velocity of an ith typical shale sample; v s-Q Representing the transverse wave velocity of quartz; v s-F Representing the transverse wave velocity of feldspar; v s-C Representing the transverse wave velocity of clay minerals; v s-Ca Represents the transverse wave velocity of calcite; vs-D represents the transverse wave velocity of dolomite; v s-P Represents the transverse wave velocity of pyrite; v s-AM The transverse wave velocity of the accessory minerals; v s-OM Represents the transverse wave velocity of the organic matter.
8. A terminal comprising a processor and a memory, the memory for storing a computer program, the processor for invoking and running the computer program stored in the memory, performing the shale brittleness index acquisition method according to any of claims 1-6.
9. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the method according to any of the preceding claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311355279.0A CN117370718B (en) | 2023-10-18 | 2023-10-18 | Shale brittleness index acquisition method, shale brittleness index acquisition device, shale brittleness index acquisition terminal and shale brittleness index storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311355279.0A CN117370718B (en) | 2023-10-18 | 2023-10-18 | Shale brittleness index acquisition method, shale brittleness index acquisition device, shale brittleness index acquisition terminal and shale brittleness index storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117370718A CN117370718A (en) | 2024-01-09 |
| CN117370718B true CN117370718B (en) | 2024-03-29 |
Family
ID=89400835
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311355279.0A Active CN117370718B (en) | 2023-10-18 | 2023-10-18 | Shale brittleness index acquisition method, shale brittleness index acquisition device, shale brittleness index acquisition terminal and shale brittleness index storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117370718B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111898065A (en) * | 2020-06-18 | 2020-11-06 | 长江大学 | Shale brittleness index calculation method and device and readable storage medium |
| US11282183B1 (en) * | 2021-04-29 | 2022-03-22 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Rock brittleness analysis method and system based on mineral content and distribution and device |
| CN116774309A (en) * | 2022-03-09 | 2023-09-19 | 中国石油化工股份有限公司 | Shale oil reservoir brittleness evaluation method based on bedding structure |
-
2023
- 2023-10-18 CN CN202311355279.0A patent/CN117370718B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111898065A (en) * | 2020-06-18 | 2020-11-06 | 长江大学 | Shale brittleness index calculation method and device and readable storage medium |
| US11282183B1 (en) * | 2021-04-29 | 2022-03-22 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Rock brittleness analysis method and system based on mineral content and distribution and device |
| CN116774309A (en) * | 2022-03-09 | 2023-09-19 | 中国石油化工股份有限公司 | Shale oil reservoir brittleness evaluation method based on bedding structure |
Non-Patent Citations (2)
| Title |
|---|
| 矿物弹性参数约束的脆性指数数字岩石物理分析;韩金鑫等;《第五届油气地球物理学术年会论文集》;20230419;全文 * |
| 页岩可压裂性声学模型及应用;李帆等;《应用声学》;20200131;第第39卷卷(第第1期期);第46-52页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117370718A (en) | 2024-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Adams et al. | Climate controls on erosion in tectonically active landscapes | |
| Aboutaleb et al. | Using non-destructive tests for estimating uniaxial compressive strength and static Young’s modulus of carbonate rocks via some modeling techniques | |
| Scarponi et al. | Sequence stratigraphy and the resolution of the fossil record | |
| Croke et al. | Use of the IHACRES rainfall-runoff model in arid and semi arid regions | |
| Fan et al. | Effects of hydromechanical loading history and antecedent soil mechanical damage on shallow landslide triggering | |
| Morey et al. | Assessment of numerical simulations of deep circulation and variability in the Gulf of Mexico using recent observations | |
| Longman et al. | Runoff events and related rainfall variability in the Southern Carpathians during the last 2000 years | |
| Ma et al. | A probabilistic method for evaluating wedge stability based on blind data theory | |
| Peña et al. | Evaluation of idealized large-eddy simulations performed with the Weather Research and Forecasting Model using turbulence measurements from a 250-m meteorological mast | |
| De Haas et al. | Fan‐surface evidence for debris‐flow avulsion controls and probabilities, Saline Valley, California | |
| Cooper et al. | Wind waves in sea ice of the western Arctic and a global coupled wave-ice model | |
| Ebeling et al. | Long‐term nitrate trajectories vary by season in Western European catchments | |
| CN117370718B (en) | Shale brittleness index acquisition method, shale brittleness index acquisition device, shale brittleness index acquisition terminal and shale brittleness index storage medium | |
| Koshurnikov et al. | The influence of salinity on the acoustic and electrical properties of frozen soils | |
| Sajinkumar et al. | Picturing beach erosion and deposition trends using PSInSAR: an example from the non-barred southern west coast of India | |
| CN105372631A (en) | Polarizing radar inversion method based on genetic-partial least square algorithm, and application of polarizing radar inversion method | |
| Nowell et al. | High energy benthic boundary layer experiment: HEBBLE | |
| Wilson | Comparative taphonomy and paleoecological reconstruction of two microvertebrate accumulations from the Late Cretaceous Hell Creek Formation (Maastrichtian), eastern Montana | |
| Martin et al. | Quantifying downstream, vertical and lateral variation in fluvial deposits: implications from the Huesca Distributive Fluvial System | |
| CN115308205B (en) | Fine-grained mixed rock dividing method, system, medium and terminal | |
| Lashkaripour et al. | Correlations between index parameters and mechanical properties of shales | |
| Zhou et al. | Parametrical Model for Estimating VS 30 from Shallow Borehole Profiles Using a Database for China | |
| Persson et al. | Estimating soil solution electrical conductivity from time domain reflectometry measurements using neural networks | |
| Chen et al. | Mutation analysis of annual sediment discharge at Wu Long station in Wu Jiang River Basin from 1960 to 2016 | |
| CN114548520A (en) | Shale oil fracturing potential prediction method, system, equipment and terminal |
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 |