CN111810133A - Stratum brittleness evaluation method - Google Patents
Stratum brittleness evaluation method Download PDFInfo
- Publication number
- CN111810133A CN111810133A CN201910285962.9A CN201910285962A CN111810133A CN 111810133 A CN111810133 A CN 111810133A CN 201910285962 A CN201910285962 A CN 201910285962A CN 111810133 A CN111810133 A CN 111810133A
- Authority
- CN
- China
- Prior art keywords
- content
- depth point
- brittleness
- silicon
- quartz
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 67
- 229910052710 silicon Inorganic materials 0.000 claims abstract description 66
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims abstract description 63
- 239000010703 silicon Substances 0.000 claims abstract description 63
- 238000000034 method Methods 0.000 claims abstract description 54
- 239000010453 quartz Substances 0.000 claims abstract description 31
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims abstract description 31
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 20
- 239000011707 mineral Substances 0.000 claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 19
- 239000004927 clay Substances 0.000 claims description 18
- 239000011575 calcium Substances 0.000 claims description 13
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 12
- BHPQYMZQTOCNFJ-UHFFFAOYSA-N Calcium cation Chemical compound [Ca+2] BHPQYMZQTOCNFJ-UHFFFAOYSA-N 0.000 claims description 8
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 8
- 239000011593 sulfur Substances 0.000 claims description 8
- 229910052717 sulfur Inorganic materials 0.000 claims description 8
- 229910052791 calcium Inorganic materials 0.000 claims description 7
- 229910052782 aluminium Inorganic materials 0.000 claims description 6
- 229910052742 iron Inorganic materials 0.000 claims description 6
- 238000005553 drilling Methods 0.000 abstract description 4
- 239000011777 magnesium Substances 0.000 description 11
- 229910052749 magnesium Inorganic materials 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 229910021532 Calcite Inorganic materials 0.000 description 4
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 4
- 229910000514 dolomite Inorganic materials 0.000 description 4
- 239000010459 dolomite Substances 0.000 description 4
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 3
- 238000002441 X-ray diffraction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052655 plagioclase feldspar Inorganic materials 0.000 description 2
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 229910052900 illite Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical compound [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011028 pyrite Substances 0.000 description 1
- 229910052683 pyrite Inorganic materials 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Physics & Mathematics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Geology (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
A formation brittleness evaluation method comprising: the method comprises the steps of firstly, obtaining the content of each element in a preset element set at each depth point in a well section to be analyzed, and determining the content of silicon element in quartz at each depth point according to the content of each element; secondly, respectively determining the content of brittle characteristic elements at each depth point in the well section to be analyzed according to the content of silicon elements in quartz at each depth point; and step three, determining the brittleness state of each depth point according to the brittleness characteristic element content. The method can accurately and continuously evaluate the brittleness of the stratum in real time through element logging information on a drilling site, and further provides reliable technical support for evaluating the drillability and the compressibility of the stratum. Meanwhile, the method only needs to calculate the content of the characteristic elements capable of accurately representing the brittle minerals in the implementation process, and has the advantages of simple operation, accurate evaluation and low cost, so that the method has a wider application range.
Description
Technical Field
The invention relates to the technical field of petroleum engineering, in particular to a formation brittleness evaluation method.
Background
In the prior art, there are twenty-thirty models for evaluating formation brittleness and brittleness index. In the field of shale gas, three common brittleness evaluation methods are available.
The first method is to use dipole sonic logging information to calculate brittleness index according to rock mechanics parameters. However, there are three disadvantages to this approach: firstly, only a few straight pilot hole wells can carry out measurement, and the horizontal section is not measured; secondly, the acoustic data is greatly influenced by borehole factors such as borehole diameter and the like; and thirdly, the difference between the rock mechanical parameters calculated by the dipole sound waves and the rock mechanical parameters actually measured in a laboratory is larger.
The second method is to use mineral data to determine the brittleness index or total amount of brittle minerals. The disadvantages of this method mainly include two points: firstly, minerals are generally measured through X-ray diffraction (XRD) whole rock analysis, and are relative content, the results given by instruments with different resolutions and different software are very different, the precision of an XRD instrument used on site is insufficient, and the sampling distance of an instrument used in a laboratory is insufficient; secondly, the components used by the molecule are controversial when calculating the brittleness index.
The third method uses the element data to calculate the brittleness index, and uses Si (silicon), Ca (calcium) and Mg (magnesium) to represent quartz, calcite and dolomite respectively. The disadvantages of this method mainly include two points: firstly, silicon and magnesium are also present in other minerals (in particular clay minerals) and do not accurately represent brittle minerals; and pyrite is not considered.
Fig. 1 to 3 show the distribution of Si (silicon), Ca (calcium) and Mg (magnesium) in different minerals. As can be seen from fig. 1 to 3, Si (silicon) is distributed in quartz, plagioclase feldspar and illite, Mg (magnesium) is distributed in calcite and chlorite, and Ca (calcium) is distributed in calcite, plagioclase feldspar and dolomite. So Si does not correlate well with quartz and Mg does with dolomite, i.e. Si does not represent quartz well, Mg does not represent dolomite well, nor Ca only calcite.
Disclosure of Invention
In order to solve the above problems, the present invention provides a formation brittleness evaluation method, including:
the method comprises the steps of firstly, obtaining the content of each element in a preset element set at each depth point in a well section to be analyzed, and determining the content of silicon element in quartz at each depth point according to the content of each element;
secondly, respectively determining the content of brittle characteristic elements at each depth point in the well section to be analyzed according to the content of silicon elements in quartz at each depth point;
and step three, determining the brittleness state of each depth point according to the brittleness characteristic element content.
According to an embodiment of the invention, the preset element set comprises silicon element, aluminum element and iron element, and in the step one, the content of the silicon element in the clay is determined according to the content of the silicon element, the content of the aluminum element and the content of the iron element in the clay, and the content of the silicon element in the quartz is determined according to the content of the silicon element in the clay.
According to one embodiment of the present invention, in said step one,
respectively calculating the total content of the three elements of the silicon element, the aluminum element and the iron element at each depth point;
and calculating the ratio of the total content of the silicon element and the three elements at each depth point, and selecting the minimum value from the ratio corresponding to each depth point to obtain the content of the silicon element in the clay.
According to one embodiment of the invention, in the second step, the content of the silicon element in the quartz is determined according to the following expression:
wherein, k'jRepresents the content of silicon element, Si, in the quartz at the jth depth pointjRepresents the content of silicon element, Al, at the j-th depth pointjDenotes the content of aluminum element, Fe, at the j-th depth pointjRepresents the content of iron element at the j depth point, and m represents the content of silicon element in the clay.
According to one embodiment of the invention, in the second step, according to the obtained contents of the calcium element and the sulfur element at each depth point, the content of the brittleness characteristic element at each depth point is respectively determined by combining the content of the silicon element in the quartz at each depth point.
According to one embodiment of the invention, in the second step, the content of brittle characteristic elements at each depth point is determined according to the following expression:
BEj＝k′j*100+Caj+Sj
wherein BEjDenotes the content, k ', of brittle characteristic element at the j-th depth point'jDenotes the content of silicon element, Ca, in the quartz at the j-th depth pointjAnd SjRespectively representing the contents of calcium element and sulfur element at the j depth point.
According to an embodiment of the invention, in the third step, a brittleness evaluation parameter is determined according to the content of the brittleness characteristic element at each depth point and a preset deviation threshold, and a brittleness state at each depth point is determined according to the value of the brittleness evaluation parameter.
According to one embodiment of the invention, in the third step, the brittleness evaluation parameter is determined by calculating the sum of the content of the brittleness characteristic element at each depth point and a preset deviation threshold value, and the brittleness state at each depth point is determined according to the brittleness evaluation parameter by utilizing a preset BM brittleness pricing standard.
According to one embodiment of the invention, the preset deviation threshold value characterizes the deviation between the total amount of brittle minerals and the content of brittle characteristic elements, which is determined by means of data fitting.
According to an embodiment of the invention, in the third step, the content of the brittle characteristic element at each depth point is matched with a preset BE brittleness evaluation standard, and the brittleness state at each depth point is determined according to the matching result.
The stratum brittleness evaluation method provided by the invention adopts the element logging technology which is relatively popular on site and relatively high in precision, utilizes the characteristic element relational expression which can accurately represent the total amount of the brittle minerals to evaluate the stratum brittleness, and can effectively overcome the defects of the existing brittleness evaluation method in accuracy, continuity and timeliness.
The method can accurately and continuously evaluate the brittleness of the stratum in real time through element logging information on a drilling site, and further provides reliable technical support for evaluating the drillability and the compressibility of the stratum. Meanwhile, the method only needs to calculate the content of the characteristic elements capable of accurately representing the brittle minerals in the implementation process, and has the advantages of simple operation, accurate evaluation and low cost, so that the method has a wider application range.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required in the description of the embodiments or the prior art:
FIGS. 1 to 3 are schematic diagrams of the distribution of Si, Ca and Mg in different minerals;
FIG. 4 is a schematic flow chart illustrating an implementation of a formation brittleness evaluation method according to one embodiment of the present invention;
FIG. 5 is a schematic diagram of an implementation flow for determining the elemental silicon content of a clay according to one embodiment of the present invention;
FIGS. 6 and 7 are graphs comparing the brittleness evaluation of the X-well 4012.5-4096.5m well segments BM to BE according to one embodiment of the present invention.
Detailed Description
The following detailed description of the embodiments of the present invention will be provided with reference to the drawings and examples, so that how to apply the technical means to solve the technical problems and achieve the technical effects can be fully understood and implemented. It should be noted that, as long as there is no conflict, the embodiments and the features of the embodiments of the present invention may be combined with each other, and the technical solutions formed are within the scope of the present invention.
In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced without some of these specific details or with other methods described herein.
Additionally, the steps illustrated in the flow charts of the figures may be performed in a computer system such as a set of computer-executable instructions and, although a logical order is illustrated in the flow charts, in some cases, the steps illustrated or described may be performed in an order different than here.
Aiming at the problems in the prior art, the invention provides a novel stratum brittleness evaluation method. The method is used for evaluating the formation brittleness by calculating the content of characteristic elements capable of accurately representing brittle minerals.
Fig. 4 shows a schematic implementation flow chart of the formation brittleness evaluation method provided by the embodiment.
As shown in fig. 4, in the method for evaluating formation brittleness provided in this embodiment, first, in step S401, the content of each element in the preset element set at each depth point in the well section to be analyzed is obtained, and in step S402, the content of silicon element in quartz at each depth point is determined according to the content of each element obtained in step S401.
In this embodiment, the method preferably obtains the content of each element in the preset element set at each depth point according to the element logging data in step S401. The elements contained in the preset element set preferably include silicon, aluminum and iron. In step S402, the method preferably determines the content of silicon in the clay according to the contents of silicon, aluminum and iron, and then determines the content of silicon in the quartz according to the content of silicon in the clay.
Specifically, as shown in fig. 5, in the present embodiment, when determining the content of silicon in the clay, the method preferably calculates the total content of three elements, namely, silicon, aluminum and iron, at each depth point in step S501, calculates the ratio of silicon to the total content of the three elements at each depth point in step S502, and selects the minimum value from the ratios corresponding to each depth point in step S503, where the minimum value is the content of silicon in the clay.
For example, in this embodiment, for any depth point j in the well section to be analyzed, the ratio of the corresponding silicon element to the total content of the three elements may be represented as Sij/(Sij+Alj+Fej). Wherein, SijRepresents the content of silicon element, Al, at the j-th depth pointjDenotes the content of aluminum element, Fe, at the j-th depth pointjIndicating the content of iron element at the j-th depth point. Assuming that the well section to be analyzed contains N depth points, N ratios can be correspondingly obtained, and a value m with the minimum value is selected from the N ratios, wherein the m can more accurately reflect the content of the silicon element in the clay of the well section to be analyzed.
In this embodiment, after the content m of the silicon element in the clay in the well section to be analyzed is obtained, the method can determine the content of the silicon element in the quartz at each depth point according to the content m of the silicon element in the clay. For example, in the present embodiment, the method may determine the content of silicon element in quartz at each depth point according to the following expression:
wherein, k'jThe content of silicon element in the quartz at the j-th depth point is shown.
Of course, in other embodiments of the present invention, according to actual needs, the method may further determine the content of the silicon element in the quartz at each depth point according to the obtained content of each element in step S402.
As shown in fig. 4 again, in this embodiment, after obtaining the content of the silicon element in the quartz at each depth point, the method preferably determines the content of the brittle characteristic element at each depth point in the well section to be analyzed according to the content of the silicon element in the quartz at each depth point in step S403.
Specifically, in this embodiment, in step S403, the method preferably determines the content of the brittle characteristic element at each depth point according to the obtained content of the calcium element and the sulfur element at each depth point, in combination with the content of the silicon element in the quartz at each depth point.
For example, the method may specifically determine the brittle characteristic element content at each depth point according to the following expression:
BEj＝k′j*100+Caj+Sj(2)
wherein BEjDenotes the content, k ', of brittle characteristic element at the j-th depth point'jDenotes the content of silicon element, Ca, in the quartz at the j-th depth pointjAnd SjRespectively representing the contents of calcium element and sulfur element at the j depth point. Content Ca of calcium element at jth depth pointjAnd the content S of sulfur elementjThe units of (c) are the same (all%).
Of course, in other embodiments of the present invention, the method may also use other reasonable manners to determine the content of the brittle characteristic element at each depth point according to the content of the calcium element and the sulfur element at each depth point and the content of the silicon element in the quartz in step S403.
In this embodiment, after obtaining the content of the brittle characteristic element at each depth point in the well section to be analyzed, the method may determine the brittle state at each depth point according to the content of the brittle characteristic element at each depth point in step S404.
Since the brittleness evaluation is mostly performed by using the evaluation standard of the total amount of brittle minerals in the prior art, in this embodiment, in step S404, the method preferably obtains the brittleness evaluation parameter corresponding to each depth point by calculating the sum of the content of the brittle characteristic element at each depth point and the preset deviation threshold, and then determines the brittleness state at each depth point according to the brittleness evaluation parameter by using the preset BM brittleness evaluation standard.
Wherein the preset deviation threshold value is characterized by the deviation between the total amount of the brittle minerals and the content of the brittle characteristic elements, and the deviation can be determined in advance by means of data fitting based on the drilled data.
For example, in this embodiment, with a preset BM brittleness evaluation criterion, for any depth point, the method may determine in step S404 whether the brittleness evaluation parameter at the depth point is smaller than the first BM brittleness evaluation threshold in the preset BM brittleness evaluation criterion. If the brittleness evaluation parameter at the depth point is smaller than the first BM brittleness evaluation threshold value, the method can judge that the brittleness state at the depth point belongs to the first brittleness state.
If the brittleness evaluation parameter at the depth point is greater than or equal to the first BM brittleness evaluation threshold, the method further determines whether the brittleness evaluation parameter is less than a second BM brittleness evaluation threshold. If the brittleness evaluation parameter at the depth point is greater than or equal to the first BM brittleness evaluation threshold value but less than the second BM brittleness evaluation threshold value, the method determines that the brittleness state at the depth point belongs to a second type of brittleness state.
And if the brittleness evaluation parameter at the depth point is greater than or equal to the second BM brittleness evaluation threshold value, the method judges that the brittleness state at the depth point belongs to a third brittleness state.
Wherein the first BM brittleness evaluation threshold is less than a second BM brittleness evaluation threshold, and the first, second, and third brittle states are indicative of low, medium, and high brittleness, respectively.
It should be noted that, in different embodiments of the present invention, specific values of the first BM brittleness evaluation threshold and the second BM brittleness evaluation threshold may be configured to be different reasonable values according to a preset BM brittleness evaluation criterion, and the present invention does not limit the specific values of the first BM brittleness evaluation threshold and the second BM brittleness evaluation threshold.
Of course, in other embodiments of the present invention, according to actual needs, the method may also establish a set of independent BE evaluation brittleness evaluation criteria in advance. In step S404, the method may determine the brittleness state at each depth point according to the matching result by matching the content of the brittle characteristic element at each depth point with a preset BE brittleness evaluation criterion. For example, the preset BE evaluation criterion may also include a first BE brittleness evaluation threshold value and a second BE brittleness evaluation threshold value, and a difference between the first BE brittleness evaluation threshold value and the first BM brittleness evaluation threshold value and a difference between the second BE brittleness evaluation threshold value and the second BM brittleness evaluation threshold value may BE preset deviation threshold values.
FIG. 6 shows a graph comparing the brittleness evaluation of the X-well 4012.5-4096.5m interval BM with BE. Wherein, BIMeasuringIs brittleness index, BI, calculated by dipole sonic loggingRecording deviceThe brittleness index calculated by (Si + Ca + Mg) recorded on the well at present shows that the two brittleness index curves have larger trend difference with the laboratory measured brittleness mineral total amount (BM) curve, and the brittleness characteristic element content (BE) determined by the method provided by the embodiment is consistent with the trend of the brittleness mineral total amount (BM).
Taking the well sections of the X wells 4012.5-4096.5m as an example, XRD whole rock analysis of 87 cores was performed on the core drilling of the X wells 4012.5-4096.5m, and the total amount of the brittle minerals counted is shown in BM curve in the attached figure 7. The BE curve in fig. 7 is a brittle characteristic element curve obtained by the method provided in this embodiment, and as can BE seen from fig. 7, the trend consistency or parallelism between the BM curve and the BE curve is high. If the criterion of low, medium and high brittleness is divided by using not less than 50 percent, not less than 50 percent to not less than 60 percent and not less than 60 percent as BM, and the criterion of low, medium and high brittleness is divided by using not less than 15 percent, not less than 15 percent to not less than 25 percent and not less than 25 percent as BE, or the criterion of low, medium and high brittleness is divided by using not less than 50 percent, not less than 50 percent to not less than 60 percent and not less than 60 percent as BE +35, the evaluation results of the two are very consistent.
From the description, the stratum brittleness evaluation method provided by the invention adopts the element logging technology which is relatively popular on site and relatively high in precision, utilizes the characteristic element relational expression capable of accurately representing the total amount of the brittle minerals to evaluate the stratum brittleness, and can effectively overcome the defects of the existing brittleness evaluation method in accuracy, continuity and timeliness.
The method can accurately and continuously evaluate the brittleness of the stratum in real time through element logging information on a drilling site, and further provides reliable technical support for evaluating the drillability and the compressibility of the stratum. Meanwhile, the method only needs to calculate the content of the characteristic elements capable of accurately representing the brittle minerals in the implementation process, and has the advantages of simple operation, accurate evaluation and low cost, so that the method has a wider application range.
It is to be understood that the disclosed embodiments of the invention are not limited to the particular structures or process steps disclosed herein, but extend to equivalents thereof as would be understood by those skilled in the relevant art. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" or "an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
While the above examples are illustrative of the principles of the present invention in one or more applications, it will be apparent to those of ordinary skill in the art that various changes in form, usage and details of implementation can be made without departing from the principles and concepts of the invention. Accordingly, the invention is defined by the appended claims.
Claims (10)
1. A formation brittleness evaluation method, comprising:
the method comprises the steps of firstly, obtaining the content of each element in a preset element set at each depth point in a well section to be analyzed, and determining the content of silicon element in quartz at each depth point according to the content of each element;
secondly, respectively determining the content of brittle characteristic elements at each depth point in the well section to be analyzed according to the content of silicon elements in quartz at each depth point;
and step three, determining the brittleness state of each depth point according to the brittleness characteristic element content.
2. The method according to claim 1, wherein the preset element set comprises silicon element, aluminum element and iron element, and in the step one, the content of the silicon element in the clay is determined according to the content of the silicon element, the content of the aluminum element and the content of the iron element in the clay, and the content of the silicon element in the quartz is determined according to the content of the silicon element in the clay.
3. The method of claim 2, wherein, in step one,
respectively calculating the total content of the three elements of the silicon element, the aluminum element and the iron element at each depth point;
and calculating the ratio of the total content of the silicon element and the three elements at each depth point, and selecting the minimum value from the ratio corresponding to each depth point to obtain the content of the silicon element in the clay.
4. The method according to claim 2 or 3, wherein in the second step, the content of silicon element in the quartz is determined according to the following expression:
wherein, k'jRepresents the content of silicon element, Si, in the quartz at the jth depth pointjRepresents the content of silicon element, Al, at the j-th depth pointjDenotes the content of aluminum element, Fe, at the j-th depth pointjRepresents the content of iron element at the j depth point, and m represents the content of silicon element in the clay.
5. The method according to any one of claims 1 to 4, wherein in the second step, according to the obtained contents of the calcium element and the sulfur element at each depth point, the content of the brittle characteristic element at each depth point is respectively determined by combining the contents of the silicon element in the quartz at each depth point.
6. The method of claim 5, wherein in step two, the brittle characteristic element content at each depth point is determined according to the expression:
BEj＝k′j*100+Caj+Sj
wherein BEjDenotes the content, k ', of brittle characteristic element at the j-th depth point'jDenotes the content of silicon element, Ca, in the quartz at the j-th depth pointjAnd SjRespectively representing the contents of calcium element and sulfur element at the j depth point.
7. The method according to any one of claims 1 to 6, wherein in the third step, a brittleness evaluation parameter is determined according to the content of the brittleness characteristic element at each depth point and a preset deviation threshold value, and the brittleness state at each depth point is determined according to the value of the brittleness evaluation parameter.
8. The method according to claim 7, wherein in step three, the brittleness evaluation parameter is determined by calculating the sum of the content of the brittle characteristic element at each depth point and a preset deviation threshold, and the brittleness state at each depth point is determined from the brittleness evaluation parameter using a preset BM brittleness pricing criterion.
9. The method according to claim 7 or 8, wherein the preset deviation threshold value is indicative of a deviation between the total amount of brittle minerals and the content of brittle characteristic elements, which is determined by means of data fitting.
10. The method according to any one of claims 1 to 6, wherein in the third step, the content of the brittle characteristic element at each depth point is matched with a preset BE brittleness evaluation standard, and the brittleness state at each depth point is determined according to the matching result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910285962.9A CN111810133B (en) | 2019-04-10 | 2019-04-10 | Stratum brittleness evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910285962.9A CN111810133B (en) | 2019-04-10 | 2019-04-10 | Stratum brittleness evaluation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111810133A true CN111810133A (en) | 2020-10-23 |
| CN111810133B CN111810133B (en) | 2023-05-09 |
Family
ID=72843758
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910285962.9A Active CN111810133B (en) | 2019-04-10 | 2019-04-10 | Stratum brittleness evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111810133B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5786595A (en) * | 1996-03-29 | 1998-07-28 | Schlumberger Technology Corporation | Method for estimating lithological fractions using nuclear spectroscopy measurements |
| US20130182819A1 (en) * | 2012-01-13 | 2013-07-18 | Ingrain, Inc. | Method Of Determining Reservoir Properties And Quality With Multiple Energy X-Ray Imaging |
| CN103982178A (en) * | 2014-04-16 | 2014-08-13 | 孙赞东 | Mineral content-based brittleness evaluation method for shale gas reservoir |
| CN104564037A (en) * | 2013-10-10 | 2015-04-29 | 中国石油天然气股份有限公司 | Shale gas reservoir brittle mineral content logging calculation method |
| CN104615802A (en) * | 2014-12-09 | 2015-05-13 | 中国石油集团川庆钻探工程有限公司 | Method for calculating content of shale brittle minerals by utilizing rock debris element analysis data |
| CN105114047A (en) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | Method utilizing logging information for determining optimal perforation layer for SRV fracturing |
| CN105986817A (en) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | Method for recognizing engineering sweet spots in shale stratum |
| CN106600436A (en) * | 2016-11-28 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司 | Method for Calculating Mineral Component Content and Porosity of Shale Gas Formation |
| CN107620593A (en) * | 2017-08-08 | 2018-01-23 | 中国科学院地质与地球物理研究所 | Scanning system and method are considered to be worth doing with rock drilling in a kind of oil gas field scene |
| CN107766662A (en) * | 2017-10-26 | 2018-03-06 | 中国石油化工股份有限公司 | A kind of horizontal well test sectional evaluation method of shale gas |
| US20180334897A1 (en) * | 2015-12-31 | 2018-11-22 | Landmark Graphics Corporation | Drilling control based on brittleness index correlation |
-
2019
- 2019-04-10 CN CN201910285962.9A patent/CN111810133B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5786595A (en) * | 1996-03-29 | 1998-07-28 | Schlumberger Technology Corporation | Method for estimating lithological fractions using nuclear spectroscopy measurements |
| US20130182819A1 (en) * | 2012-01-13 | 2013-07-18 | Ingrain, Inc. | Method Of Determining Reservoir Properties And Quality With Multiple Energy X-Ray Imaging |
| CN104564037A (en) * | 2013-10-10 | 2015-04-29 | 中国石油天然气股份有限公司 | Shale gas reservoir brittle mineral content logging calculation method |
| CN103982178A (en) * | 2014-04-16 | 2014-08-13 | 孙赞东 | Mineral content-based brittleness evaluation method for shale gas reservoir |
| CN104615802A (en) * | 2014-12-09 | 2015-05-13 | 中国石油集团川庆钻探工程有限公司 | Method for calculating content of shale brittle minerals by utilizing rock debris element analysis data |
| CN105986817A (en) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | Method for recognizing engineering sweet spots in shale stratum |
| CN105114047A (en) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | Method utilizing logging information for determining optimal perforation layer for SRV fracturing |
| US20180334897A1 (en) * | 2015-12-31 | 2018-11-22 | Landmark Graphics Corporation | Drilling control based on brittleness index correlation |
| CN106600436A (en) * | 2016-11-28 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司 | Method for Calculating Mineral Component Content and Porosity of Shale Gas Formation |
| CN107620593A (en) * | 2017-08-08 | 2018-01-23 | 中国科学院地质与地球物理研究所 | Scanning system and method are considered to be worth doing with rock drilling in a kind of oil gas field scene |
| CN107766662A (en) * | 2017-10-26 | 2018-03-06 | 中国石油化工股份有限公司 | A kind of horizontal well test sectional evaluation method of shale gas |
Non-Patent Citations (6)
| Title |
|---|
| TAEYOUN KIM等: "Petrophysical approach for S-wave velocity prediction based on brittleness index and total organic carbon of shale gas reservoir: A case study from Horn River Basin, Canada", 《JOURNAL OF APPLIED GEOPHYSICS》 * |
| 张新华;邹筱春;赵红艳;李芳;秦黎明;: "利用X荧光元素录井资料评价页岩脆性的新方法", 石油钻探技术 * |
| 李芳;: "利用录井技术计算页岩脆性含量的方法", 江汉石油职工大学学报 * |
| 牛强等: "X射线元素录井技术在胜利油区泥页岩脆性评价中的应用", 《油气地质与采收率》 * |
| 王淑芳等: "四川盆地富有机质页岩硅质生物成因及对页岩气开发的意义", 《北京大学学报（自然科学版）》 * |
| 赵建华等: "四川盆地五峰组-龙马溪组页岩岩相类型与沉积环境", 《石油学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111810133B (en) | 2023-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102175832B (en) | Method for determining optimal saturation calculation model of typical reservoir | |
| CN104989392B (en) | Lithology identification method | |
| CN104239703A (en) | Quantitative analogical evaluation method for multiple parameters of shale gas reservoir | |
| CN108194077B (en) | Gas logging total hydrocarbon correction method | |
| CN109655909B (en) | Method and system for calculating shale full-angle anisotropy brittleness index | |
| CN110219644B (en) | Method for determining spatial distribution of reservoir fracability index values | |
| CN112145165B (en) | Microcrack-pore type reservoir dynamic and static permeability conversion method | |
| CN106054254B (en) | A kind of log curve correction method and device | |
| CN103306672B (en) | A kind of method predicting the abrasiveness of the different drilling direction of shale formation | |
| CN105350959A (en) | Method for determining gas saturation of shale gas reservoir through well-logging lithologic density | |
| CN107545079A (en) | A kind of wellbore trace offsets monitoring method | |
| CN111810133B (en) | Stratum brittleness evaluation method | |
| CN107103165A (en) | A kind of method of CALCULATING OILFIELD reservoir steep dip shearing crack line density | |
| CN114060015B (en) | Method and device for evaluating gas content of compact sandstone | |
| CN108647417B (en) | Simple method for determining gas saturation of shale gas reservoir | |
| CN109632546B (en) | Rapid testing system and method for quartz content of rock of Tunnel Boring Machine (TBM) | |
| CN108627878A (en) | The crack identification method and system of tight sand formation | |
| CN112147713B (en) | Shale total organic carbon content segmented prediction method | |
| CN108240999A (en) | Method for identifying lithology of clastic rock | |
| CN113775329A (en) | Method and system for acquiring deep carbonate rock mechanical parameters | |
| Amanipoor | Providing a subsurface reservoir quality maps in oil fields by geostatistical methods | |
| CN115479540B (en) | Karst water bursting disaster magnitude assessment method and system based on laser scanning | |
| CN118167280A (en) | Method and device for evaluating EUR of tight sandstone gas reservoir sidetracking horizontal well | |
| CN115983694A (en) | Method and system for determining surrounding rock level of deep-buried extra-long tunnel | |
| CN105484732A (en) | Method for processing well depth in drilling geosteering construction process of horizontal well |
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 |