[go: up one dir, main page]

CN111458767B - A method and system for identifying lithology based on intersection map method - Google Patents

A method and system for identifying lithology based on intersection map method Download PDF

Info

Publication number
CN111458767B
CN111458767B CN202010281356.2A CN202010281356A CN111458767B CN 111458767 B CN111458767 B CN 111458767B CN 202010281356 A CN202010281356 A CN 202010281356A CN 111458767 B CN111458767 B CN 111458767B
Authority
CN
China
Prior art keywords
lithology
logging
amplitude
mudstone
sandstone
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.)
Expired - Fee Related
Application number
CN202010281356.2A
Other languages
Chinese (zh)
Other versions
CN111458767A (en
Inventor
滕娟
夏宇
邓虎成
陈文玲
王园园
解馨慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202010281356.2A priority Critical patent/CN111458767B/en
Publication of CN111458767A publication Critical patent/CN111458767A/en
Application granted granted Critical
Publication of CN111458767B publication Critical patent/CN111458767B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种基于交汇图法识别岩性的方法及系统,其中方法包括,S01测井解释系列的选择,SO2测井数据归一化处理,S03利用两步岩性识别法精细识别岩性,S04通过对不同岩性类型的数据进行统计分析,对其频率进行对比,确定岩性分界线、岩性类型,S05判断识别岩性。本发明通过对研究区岩性进行分析,以一个新的思路,方法和模型快速确定岩性方法,在保留了交汇图法使用简单,将数据图示化,快速直观识别岩性的基础上,解决了只有少量岩心、薄片及录井资料时,交汇图法难以运用和无法达到精细识别岩性的标准的问题,并利用模型对大量测井数据进行岩性识别。

Figure 202010281356

The invention provides a method and system for identifying lithology based on a cross-graph method, wherein the method includes: S01 selection of logging interpretation series, SO2 logging data normalization processing, S03 using two-step lithology identification method to finely identify lithology , S04 by performing statistical analysis on the data of different lithology types and comparing their frequencies to determine the lithologic boundary and lithologic types, and S05 to judge and identify the lithology. By analyzing the lithology of the research area, the present invention uses a new idea, method and model to quickly determine the lithology method, and on the basis of retaining the simple use of the intersection diagram method, graphing the data, and quickly and intuitively identifying the lithology, It solves the problem that the intersection map method is difficult to use and cannot reach the standard of fine identification of lithology when there is only a small amount of core, thin section and logging data, and the model is used to identify the lithology of a large amount of logging data.

Figure 202010281356

Description

Method and system for identifying lithology based on intersection graph method
Technical Field
The invention belongs to the technical field of oil-gas exploration, and particularly relates to a method and a system for identifying lithology based on a cross plot method.
Background
The lithology characteristic identification is the basis for recognizing stratum and reservoir prediction, the lithology identification is the first step of reservoir prediction, the lithology characteristic identification has very important significance for predicting high-yield enriched areas of oil gas, and the lithology characteristic identification is beneficial to exploration and development of the oil gas. The rock core observation is the most direct and effective method for identifying lithology, but the method cannot be widely applied due to high cost and complex operation flow, and the method for explaining lithology by logging information is promoted.
At present, various methods have been proposed at home and abroad for explaining lithology of logging data, such as conventional logging series, imaging logging, element logging and the like. However, due to cost, most oil fields still use conventional well logging series, and the method for explaining lithology by using conventional well logging is widely applied. At present, the intersection graph method is the method which explains lithology application most by conventional well logging information, but the conventional well logging information used by the intersection graph method is the comprehensive reflection of the combined action of all factors of a stratum, and the problem of low lithology depicting precision exists.
At present, the lithology is explained by using an intersection map at home and abroad mainly by a direct intersection map method. The direct intersection graph method is used for identifying lithology, two kinds of lithology-sensitive logging curve data are intersected directly, a standard graph is built, and lithology is identified. And realizing lithology division according to the characteristic that different lithology projects in different areas on a plane after intersection. The intersection graph method has the advantages that data are converted into graphs, and the purpose of more intuitionistic lithology identification is achieved. However, the method has the following disadvantages and drawbacks:
1. the intersection graph method adopts the conventional well logging series curves, has the condition of low recognition degree and cannot reach the standard of fine lithology recognition.
2. The direct intersection graph method requires a lot of time for data processing.
Disclosure of Invention
In order to overcome the problems in the background art, the invention provides a method and a system for identifying lithology based on a cross plot method, which can quickly identify various lithologies, have high discrimination, can effectively identify the lithology and achieve the semi-quantitative standard.
The invention adopts the following technical scheme:
a lithology recognition method based on a junction graph method comprises the following steps: the method comprises the following steps:
s01, selecting a logging interpretation series;
selecting a well with existing drilling information, researching and analyzing the well logging, well logging and lithology of the well, selecting any two of GR, SP, DEN, AC and RILD well logging curves to perform cluster analysis in SPSS, and selecting a well logging curve with high lithology recognition degree;
s02, carrying out normalization processing on logging data;
determining a standard well and a standard layer, wherein the selection principle of the standard well and the standard well is determined according to the selection principle of the standard well and the standard layer in the petroleum logging data standardization;
correcting logging curves of all wells in a research area by adopting a peak value method by taking the selected standard well as a standard, and calculating to obtain a difference value of the logging curves of the standard well and the well to be detected until the logging data standardization processing of the rest wells in the work area is completed;
s03, finely identifying the lithology by using a two-step lithology identification method,
performing Fisher discriminant analysis by SPSS software to divide lithology into three major categories including carbonate rock, sandstone rock and mudstone;
on the basis of completing the lithological major classification judgment, further utilizing a Fisher judgment model to continuously divide the lithological properties in each major classification into subclasses, such as: mudstone, grey mudstone, sandy mudstone, oil shale limestone, dolomite, marl, sandy limestone, sandstone, grey sandstone and argillaceous sandstone;
s04, performing statistical analysis on data of different lithological types, comparing frequencies of the data, and determining a lithological boundary and a rock type;
and S05, inputting a correction value (difference value) obtained by the well logging data standardization processing module according to the divided lithological boundary to finish the identification of lithological judgment.
The further technical proposal is that in S03, the DEN value of carbonate rocks is more than 2.6g/cm3The DEN value amplitude of sandstone and mudstone is less than 2.6g/cm3And the GR amplitude of the sandstone is smaller than 62API, the GR amplitude of the mudstone is larger than 62API, and the DEN amplitude is 2.6g/cm3And GR amplitude 62API as a boundary value for discriminating the three major lithologies of the study area.
The further technical scheme is that in S03, the sandstone includes sandstone, gray sandstone and argillaceous sandstone, the GR amplitude of the gray sandstone is smaller than 44API, the sandstone is 44-56API, the argillaceous sandstone is larger than 56API, and the GR amplitudes of 44API and 56API are used as the boundary values of the sandstone, the gray sandstone and the argillaceous sandstone.
The further technical scheme is that the mudstone comprises mudstone, gray mudstone, sandy mudstone and oil shale, and the DEN amplitude of the gray mudstone is more than 2.4g/cm3The sandy mudstone is 2.28-2.4g/cm3And mudstone is less than 2.28g/cm3The DEN amplitude was set to 2.28g/cm3And 2.4g/cm3As the boundary value of the gray mud rock, the sandy mud rock and the mudstone.
The further technical scheme is that carbonate rocks comprise limestone, dolomite, marl and sandy limestone, and DEN amplitudes of the limestone and the dolomite are both more than 2.6g/cm3When the GR amplitude is larger than 32API, the GR amplitude is dolostone and is smaller than 32API, the limestone amplitude is lower than 2.6g/cm, and the DEN amplitude of impurity limestone is lower than3The DEN amplitude is 2.6g/cm3As the boundary value of carbonate rock and impurity limestone, the GR amplitude 32API is used as the boundary value of limestone and dolomite.
A system for identifying lithology based on a junction graph method comprises the following modules:
and the well logging data curve distinguishing module selects a well with existing well drilling data, carries out research and analysis on well logging, well logging and lithology of the well, selects any two curves from GR, SP, DEN, AC and RILD well logging curves to carry out cluster analysis in SPSS, and selects a well logging curve with high intersection map recognition degree.
And the well logging data standardization processing module is used for determining the standard wells and the standard layers according to the selection principle of the standard wells and the standard layers in the petroleum well logging data standardization, inputting data of the rest wells, correcting by adopting a peak method, calculating to obtain a difference value of the well logging curves of the standard wells and the wells to be measured by utilizing the curves screened out by the well logging data curve discrimination module, and completing the well logging data standardization of the rest wells in the research area.
And the lithology boundary generation module is used for carrying out Fisher discriminant analysis by using SPSS software, counting data of different lithologies, analyzing frequency distribution and comparing to obtain a lithology boundary.
And the lithology judging module is used for inputting the original logging data and the correction value (difference value) obtained by the logging data standardization processing module according to the divided lithology boundary so as to finish the identification of the lithology judgment.
The invention has the beneficial effects that:
1. the method selects key wells and standard layers, carries out well logging data standardization processing on other wells, reduces data errors, provides a basis for fine interpretation of lithology, preferably selects well logging curves sensitive to lithology to be intersected on coordinate axes, projects existing data points on a coordinate system, lays a cushion for dividing lithology boundaries, determines areas with different lithologies based on different falling points of different lithologies in the coordinate system, compiles an intersection diagram, and predicts the lithologies of other wells by using an established chart.
2. The method firstly determines the lithologic boundary by using a small amount of data, then establishes the identification chart, establishes the two-step method for identifying the lithologic property, identifies the large lithologic property firstly, identifies various lithologic properties, has high discrimination, can effectively identify the lithologic property and reaches the semi-quantitative standard.
Compared with the prior direct intersection graph method:
the direct intersection graph method adopts a conventional well logging series curve, has low recognition degree, can only achieve the condition of qualitative recognition, and cannot achieve the standard of fine lithology recognition. The lithology can be directly identified by establishing a lithology identification system, and a large amount of time for repeatedly utilizing plate casting points and then processing data by a direct intersection mapping method is reduced.
Drawings
1(a) -1 (b) are amplitude comparison graphs of logs of a thick mudstone section according to the invention;
FIG. 2(a) -FIG. 2(d) are schematic diagrams of lithology identification plates;
FIG. 3 is a diagram of a lithology identification procedure;
FIG. 4 is a schematic diagram of lithology judgment and identification;
FIG. 5 is a schematic diagram of lithology identification verification;
FIG. 6 is a flowchart of the process of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention are described below clearly and completely, and it is obvious that the described embodiments are some, not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
GR-natural gamma logging, SP-natural potential logging, DEN-density logging, AC-sonic time difference logging, RILD-deep induction logging;
API-natural gamma logging unit.
A system for identifying lithology based on a junction graph method comprises the following modules:
and the well logging data curve distinguishing module selects a well with existing well drilling data, carries out research and analysis on well logging, well logging and lithology of the well, selects any two curves from GR, SP, DEN, AC and RILD well logging curves to carry out cluster analysis in SPSS, and selects a well logging curve with high intersection map recognition degree.
And the well logging data standardization processing module is used for determining the standard wells and the standard layers according to the selection principle of the standard wells and the standard layers in the petroleum well logging data standardization, inputting data of the rest wells, correcting by adopting a peak method, calculating to obtain a difference value of the well logging curves of the standard wells and the wells to be measured by utilizing the curves screened out by the well logging data curve discrimination module, and completing the well logging data standardization of the rest wells in the research area.
And the lithology boundary generation module is used for carrying out Fisher discriminant analysis by using SPSS software, counting data of different lithologies, analyzing frequency distribution and comparing to obtain a lithology boundary.
And the lithology judging module is used for inputting the original logging data and the correction value (difference value) obtained by the logging data standardization processing module according to the divided lithology boundary so as to finish the identification of the lithology judgment.
As shown in fig. 6, a method for identifying lithology based on a junction graph method includes the following steps:
s01, selecting a logging interpretation series;
selecting a well with existing drilling information, researching and analyzing the well logging, well logging and lithology of the well, selecting any two of GR, SP, DEN, AC and RILD well logging curves to perform cluster analysis in SPSS, and selecting a well logging curve with high lithology recognition degree;
s02, carrying out normalization processing on logging data;
determining a standard well and a standard layer, wherein the selection principle of the standard well and the standard well is determined according to the selection principle of the standard well and the standard layer in the petroleum logging data standardization;
correcting logging curves of all wells in a research area by adopting a peak value method by taking the selected standard well as a standard, and calculating to obtain a difference value of the logging curves of the standard well and the well to be detected until the logging data standardization processing of the rest wells in the work area is completed;
s03, finely identifying the lithology by using a two-step lithology identification method,
performing Fisher discriminant analysis by SPSS software to divide lithology into three major categories including carbonate rock, sandstone rock and mudstone;
on the basis of completing the lithological major classification judgment, further utilizing a Fisher judgment model to continuously divide the lithological properties in each major classification into subclasses, such as: mudstone, grey mudstone, sandy mudstone, oil shale limestone, dolomite, marl, sandy limestone, sandstone, grey sandstone and argillaceous sandstone;
s04, performing statistical analysis on data of different lithological types, comparing frequencies of the data, and determining a lithological boundary and a lithological type;
and S05, inputting a correction value (difference value) obtained by the well logging data standardization processing module according to the divided lithological boundary to finish the identification of lithological judgment.
The further technical proposal is that in S03, the DEN value of carbonate rocks is more than 2.6g/cm3The DEN value amplitude of sandstone and mudstone is less than 2.6g/cm3And the GR amplitude of the sandstone is smaller than 62API, the GR amplitude of the mudstone is larger than 62API, and the DEN amplitude is 2.6g/cm3And GR amplitude 62API as a boundary value for discriminating the three major lithologies of the study area.
The further technical scheme is that in S03, the sandstone includes sandstone, gray sandstone and argillaceous sandstone, the GR amplitude of the gray sandstone is smaller than 44API, the sandstone is 44-56API, the argillaceous sandstone is larger than 56API, and the GR amplitudes of 44API and 56API are used as the boundary values of the sandstone, the gray sandstone and the argillaceous sandstone.
The further technical scheme is that the mudstone comprises mudstone, gray mudstone, sandy mudstone and oil shale, and the DEN amplitude of the gray mudstone is more than 2.4g/cm3The sandy mudstone is 2.28-2.4g/cm3And mudstone is less than 2.28g/cm3Amplitude of DEN is set to 2.28g/cm3And 2.4g/cm3As the boundary value of the gray mud rock, the sandy mud rock and the mudstone.
The further technical scheme is that carbonate rocks comprise limestone, dolomite, marl and sandy limestone, and DEN amplitudes of the limestone and the dolomite are both more than 2.6g/cm3When the GR amplitude is larger than 32API, the GR amplitude is dolostone and is smaller than 32API, the limestone amplitude is lower than 2.6g/cm, and the DEN amplitude of impurity limestone is lower than3The DEN amplitude is 2.6g/cm3As the boundary value of carbonate rock and impurity limestone, the GR amplitude 32API is used as the boundary value of limestone and dolomite.
Example 1.
A3-section sample of a Shahechu block in Bohai Bay basin Boxing hollow pure beam region is selected, and the sample has a small amount of core and slice data and a large amount of logging data. The lithology is judged by lacking visual data, the lithology needs to be identified through well logging interpretation, and the lithology is finely identified by adopting the method.
Taking fine lithology identification of samples at 3 sections of a sand river street group in Bohai Bay basin Boxing hollow pure beam region as an example.
(1) The method comprises the steps of selecting a well with existing well logging, well logging and lithology information, and conducting research and analysis on lithology, wherein the lithology of a research area is considered to be mainly mudstone and is inferior to sandstone and carbonate. Aiming at the main lithology of a research area, two curves in GR, SP, DEN, AC and RILD are selected in turn to carry out cluster analysis in SPSS software, and then the GR-AC and DEN-GR intersection graph is found to have high recognition degree, so that the lithology can be effectively distinguished. And selecting three logging curves of GR, AC and DEN to carry out logging data standardization. Here, an F120(F120 well is one of the wells to be selected) well is selected as a standard well and the thickest mudstone is selected as a standard layer according to the standard well and standard layer selection principle in the standardization of the petroleum logging data.
Correcting a well logging curve of an FY1 well to be measured (an FY1 well is a well to be measured) by using a peak value method with an F120 well as a standard well, and respectively reading GR and AC characteristic peak values of the well logging curves of the F120 well and the FY1 well (the FY1 well is one well in the research area) corresponding to a standard layer, as shown in fig. 1(a) -1 (b), calculating to obtain a difference value of the well logging curves of the FY1 well and the F120 well, wherein the GR is +10, the AC is +10, and the DEN is-0.09. The method is used for finally completing the standardization of logging data of the rest wells in the research area.
After discrimination is carried out through SPSS software, Fisher discrimination analysis is carried out to divide lithology into three major classes, including carbonate rock class, sandstone rock class and mudstone class;
on the basis of completing the lithological major classification judgment, further utilizing a Fisher judgment model to continuously divide the lithological properties in each major classification into minor classifications. Mudstone can be divided into mudstone, gray mudstone, sandy mudstone and oil shale; carbonate rocks can be classified as limestone, dolomite, marl, and sandy limestone; the sandstone can be divided into sandstone, gray sandstone and argillaceous sandstone.
And (3) performing statistical analysis on the data of different lithology types, comparing the frequencies of the data, and determining a lithology boundary. Firstly, the lithology is divided into three major categories including carbonate rock, sandstone rock and mudstone, and the DEN value of the carbonate rock in figure 2(a) is more than 2.6g/cm3The DEN value amplitude of sandstone and mudstone is less than 2.6g/cm3And the GR amplitude of sandstone is less than 62API, and the GR amplitude of shale is greater than 62 API. Thus, the DEN amplitude was 2.6g/cm3And GR amplitude 62API as a boundary value for discriminating the three major lithologies of the study area.
And on the basis of finishing the lithology major category identification, further continuously dividing the lithology in each major category into minor categories. For sandstone, including sandstone, limestone and argillaceous sandstone, as shown in fig. 2(b), the DEN and GR amplitudes are used to intersect, and the result shows that the DEN amplitudes of the three are relatively close, and the GR amplitudes are relatively large, specifically, the GR amplitude of the grey sandstone is smaller than 44API, the sandstone is 44-56API, and the argillaceous sandstone is greater than 56 API. Therefore, the GR amplitude values of 44API and 56API can be used as the boundary values of sandstone, gritty rock and argillaceous sandstone.
As shown in fig. 2(c), for the mudstone class including mudstone, gray mudstone, sandy mudstone and oil shale, the GR amplitude and the DEN amplitude are used for intersection, and the results show that the GR amplitudes of the three are close and are not easy to distinguish, and the DEN amplitude difference is obvious, specifically, the DEN amplitude of the gray mudstone is greater than 2.4g/cm3The sandy mudstone is 2.28-2.4g/cm3And mudstone is less than 2.28g/cm3. Thus, the DEN amplitude was set to 2.28g/cm3And 2.4g/cm3As the boundary value of the gray mud rock, the sandy mud rock and the mudstone.
Carbonates include limestone, dolomite, marl and sandstone limestone. The GR amplitude and the DEN amplitude are used for intersection, the result shows that the carbonate rock and the transitional lithology thereof can be distinguished by the DEN amplitude, the dolomite rock and the limestone rock are difficult to distinguish by the DEN amplitude, and the GR curve is used for separating the dolomite rock and the limestone rock. Specifically, DEN amplitudes of both limestone and dolomite are more than 2.6g/cm3When the GR amplitude is larger than 32API, the GR amplitude is dolostone and is smaller than 32API, the limestone amplitude is lower than 2.6g/cm, and the DEN amplitude of impurity limestone is lower than3As shown in fig. 2 (d). Thus, DEN amplitude was 2.6g/cm3As the boundary value of carbonate rock and impurity limestone, the GR amplitude 32API is used as the boundary value of limestone and dolomite.
As shown in FIG. 3, the lithology is recognized on the basis of the method, after the standardized correction value of the lithology is input at GRadd and DENadd, as shown in FIG. 4, the lithology judgment and the top and bottom depths are clicked, the model can correct the well logging curve according to the standardized correction value, the lithology corresponding to the original depth is judged, and the recognized lithology is generated.
As shown in fig. 5, in the section 3066m-3067.5m, the logging result is shown as marl, and the identification result by the lithology identification chart is shown as argillaceous sandstone; in the 3069m-3070m section, the logging result is displayed as grey mud rock, and the identification result of the lithology identification chart is used for displaying argillaceous sandstone; and in the section of 3070.5m-3073.5m, the logging results are displayed as mudstone and gray mudstone, the identification results of the lithology identification chart are used for displaying the argillaceous sandstone, and the identification results are all consistent with the core identification results.
In the embodiment, the lithology of the research area is analyzed, the lithology is rapidly determined by a new method and a new thought, the problems that the cross plot method is difficult to use and cannot meet the standard of fine lithology identification when rock cores, slices and logging information are lacked are solved on the basis that the cross plot method is simple to use, data are converted into graphs, and the lithology can be identified more rapidly and visually.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present 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 solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (6)

1.一种基于交汇图法识别岩性的方法,其特征在于,包括以下步骤:1. a method for identifying lithology based on the intersection diagram method, is characterized in that, comprises the following steps: S01.测井解释系列的选择S01. Selection of logging interpretation series 选取已有钻井资料的井,对其测井、录井和岩性进行研究分析,选择GR、SP、DEN、AC、RILD测井曲线中的任意两条曲线在SPSS中进行聚类分析,选择出岩性识别度高的测井曲线,后续步骤基于选择出的测井曲线进行处理;Select wells with existing drilling data, conduct research and analysis on their logging, logging and lithology, select any two curves in GR, SP, DEN, AC, RILD logging curves for cluster analysis in SPSS, select Logging curves with high lithology identification are produced, and subsequent steps are processed based on the selected logging curves; S02.测井数据归一化处理S02. Logging data normalization processing 确定标准井及标准层,标准层和标准井的选择原则根据中石油测井数据标准化中标准井和标准层的选择原则确定;Determine standard wells and standard layers, and the selection principles of standard layers and standard wells are determined according to the selection principles of standard wells and standard layers in the standardization of CNPC logging data; 以选取的标准井的标准层为基准,采用峰值法对研究区所有井的测井曲线进行校正,计算得到标准井和待测井测井曲线的差值,直至完成工区剩余井的测井数据标准化处理;Based on the standard layer of the selected standard well, the peak value method is used to correct the logging curves of all wells in the study area, and the difference between the standard well and the logging curve to be logged is calculated until the logging data of the remaining wells in the work area are completed. standardized processing; S03.利用两步岩性识别法精细识别岩性S03. Use two-step lithology identification method to finely identify lithology 通过SPSS软件进行Fisher判别分析将岩性划分为三大类,包括碳酸盐岩类、砂岩类和泥岩类;The lithology is divided into three categories by Fisher discriminant analysis by SPSS software, including carbonate rock, sandstone and mudstone; 完成岩性大类判识的基础上,进一步利用Fisher判别模型将每一大类中的岩性继续划分为小类;On the basis of completing the identification of lithology categories, Fisher's discriminant model is further used to further divide the lithology in each major category into sub-categories; S04.确定岩性分界线及岩性类型S04. Determine lithologic boundary and lithologic type 通过对不同岩性类型的数据进行统计分析,对其频率进行对比,确定岩性分界线、岩石类型;Through statistical analysis of data of different lithology types and comparison of their frequencies, the lithology boundary and rock type are determined; S05.判断识别岩性S05. Judging and identifying lithology 根据划分的岩性分界线,输入根据测井数据标准化处理模块得到的校正值,完成岩性判定的识别。According to the divided lithologic boundary, input the correction value obtained by the standardization processing module of logging data to complete the identification of lithologic judgment. 2.根据权利要求1中所述的基于交汇图法识别岩性的方法,其特征在于,S03中,碳酸盐岩类的DEN值大于2.6g/cm3,砂岩类和泥岩类的DEN值幅值小于2.6g/cm3,且砂岩类的GR幅值小于62API,泥岩类的GR幅值大于62API,将DEN幅值2.6g/cm3和GR幅值62API作为判识研究区三大类岩性的分界值。2. The method for identifying lithology based on the intersection diagram method according to claim 1, wherein in S03, the DEN value of carbonate rocks is greater than 2.6g/cm 3 , and the DEN value of sandstone and mudstone The amplitude is less than 2.6g/cm 3 , and the GR amplitude of sandstone is less than 62API, and the GR amplitude of mudstone is greater than 62API. The DEN amplitude of 2.6g/cm 3 and the GR amplitude of 62API are used to identify the three categories of the study area. Lithology cut-off value. 3.根据权利要求1中所述的基于交汇图法识别岩性的方法,其特征在于,S03中,砂岩类包括砂岩、灰质砂岩和泥质砂岩,灰质砂岩的GR幅值小于44API,砂岩为44-56API,而泥质砂岩大于56API,将GR幅值为44API和56API作为砂岩、灰质砂岩和泥质砂岩的分界值。3. The method for identifying lithology based on the intersection diagram method according to claim 1, wherein in S03, the sandstones include sandstone, calcareous sandstone and argillaceous sandstone, the GR amplitude of the calcareous sandstone is less than 44 API, and the sandstone is 44-56API, while the argillaceous sandstone is greater than 56API, the GR amplitudes of 44API and 56API are taken as the demarcation value of sandstone, calcareous sandstone and argillaceous sandstone. 4.根据权利要求1中所述的基于交汇图法识别岩性的方法,其特征在于,泥岩类包括泥岩、灰质泥岩、砂质泥岩和油页岩,灰质泥岩的DEN幅值大于2.4g/cm3,砂质泥岩为2.28-2.4g/cm3,而泥岩小于2.28g/cm3,将DEN幅值为2.28g/cm3和2.4g/cm3作为灰质泥岩、砂质泥岩和泥岩的分界值。4. The method for identifying lithology based on the intersection diagram method according to claim 1, wherein the mudstones include mudstone, calcareous mudstone, sandy mudstone and oil shale, and the DEN amplitude of calcareous mudstone is greater than 2.4g/ cm 3 , the sandy mudstone is 2.28-2.4 g/cm 3 , while the mudstone is less than 2.28 g/cm 3 , the DEN amplitudes are 2.28 g/cm 3 and 2.4 g/cm 3 as the calcareous mudstone, sandy mudstone and mudstone. cutoff value. 5.根据权利要求1中所述的基于交汇图法识别岩性的方法,其特征在于,碳酸盐岩类包括灰岩、白云岩、泥灰岩和砂质灰岩,灰岩和白云岩的DEN幅值均大于2.6g/cm3,当GR幅值大于32API为白云岩,小于32API为灰岩,杂质灰岩DEN幅值小于2.6g/cm3,将DEN幅值2.6g/cm3作为碳酸盐岩和杂质灰岩的分界值,GR幅值32API作为灰岩和白云岩的分界值。5. The method for recognizing lithology based on the intersection diagram method according to claim 1, wherein the carbonate rocks include limestone, dolomite, marl and sandy limestone, and limestone and dolomite The DEN amplitude of the GR is greater than 2.6g/cm 3 , when the GR amplitude is greater than 32API, it is dolomite, less than 32API is limestone, and the DEN amplitude of impurity limestone is less than 2.6g/cm 3 , the DEN amplitude is 2.6g/cm 3 As the boundary value of carbonate rock and impurity limestone, GR amplitude 32API is used as the boundary value of limestone and dolomite. 6.一种基于交汇图法识别岩性的系统,其特征在于,包括以下模块:6. A system for identifying lithology based on intersection diagram method, is characterized in that, comprises following module: 测井数据曲线判别模块,选取已有钻井资料的井,对其测井、录井和岩性进行研究分析,选择GR、SP、DEN、AC、RILD测井曲线中的任意两条曲线在SPSS中进行聚类分析,选择出交汇图识别度高的测井曲线,后续模块基于选择出的测井曲线进行处理;Logging data curve discrimination module, select wells with existing drilling data, conduct research and analysis on their logging, logging and lithology, and select any two curves in GR, SP, DEN, AC, RILD logging curves in SPSS Clustering analysis is performed in the module, and the logging curve with high recognition degree of the intersection graph is selected, and the subsequent modules are processed based on the selected logging curve; 测井数据标准化处理模块,根据中石油测井数据标准化中标准井和标准层的选择原则确定标准井和标准层后,输入剩余各井的数据,采用峰值法进行校正,利用测井数据曲线判别模块筛选出的曲线,计算得到标准井和待测井测井曲线的差值,完成研究区剩余各井的测井数据标准化;Logging data standardization processing module, after the standard wells and standard layers are determined according to the selection principle of standard wells and standard layers in the standardization of CNPC logging data, the data of the remaining wells are input, and the peak value method is used for correction, and the log data curve is used to identify the module The selected curve is calculated to obtain the difference between the standard well and the log curve of the well to be logged, and the standardization of the logging data of the remaining wells in the study area is completed; 岩性分界线生成模块,利用SPSS软件进行Fisher判别分析,对不同岩性的数据进行统计,分析频率分布进行对比得出岩性分界线;The lithological boundary generation module uses SPSS software to carry out Fisher's discriminant analysis, makes statistics on the data of different lithology, analyzes the frequency distribution and compares the lithological boundary; 岩性判定模块,根据划分的岩性分界线,输入原始测井数据和测井数据标准化处理模块得到的校正值,完成岩性判定的识别。The lithology judgment module, according to the divided lithology boundary line, inputs the original logging data and the correction value obtained by the logging data standardization processing module, and completes the identification of lithology judgment.
CN202010281356.2A 2020-04-10 2020-04-10 A method and system for identifying lithology based on intersection map method Expired - Fee Related CN111458767B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010281356.2A CN111458767B (en) 2020-04-10 2020-04-10 A method and system for identifying lithology based on intersection map method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010281356.2A CN111458767B (en) 2020-04-10 2020-04-10 A method and system for identifying lithology based on intersection map method

Publications (2)

Publication Number Publication Date
CN111458767A CN111458767A (en) 2020-07-28
CN111458767B true CN111458767B (en) 2021-06-04

Family

ID=71683729

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010281356.2A Expired - Fee Related CN111458767B (en) 2020-04-10 2020-04-10 A method and system for identifying lithology based on intersection map method

Country Status (1)

Country Link
CN (1) CN111458767B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113325480A (en) * 2021-05-21 2021-08-31 成都理工大学 Seismic lithology identification method based on integrated deep learning
CN114330023B (en) * 2022-01-14 2022-09-23 中国石油大学(北京) Method and device for predicting favorable reservoir distribution in tight sandstone
CN115354992B (en) * 2022-08-31 2024-08-16 成都理工大学 A coal-measure gas reservoir evaluation method based on gas-bearing characteristics of lithology combination
CN115506776A (en) * 2022-09-23 2022-12-23 成都理工大学 A method for identification of easily soluble fabrics in carbonate reservoirs
CN115653580B (en) * 2022-11-07 2024-10-22 成都理工大学 A method for characterizing the internal structure of complex lithology reservoirs
CN115793094B (en) * 2023-02-06 2023-04-25 西北大学 A method and application of curve superposition and reconstruction to identify lithology of complex shale formations
CN117418831B (en) * 2023-12-18 2024-03-22 西安石油大学 Lithology identification method and device for sandstone reservoir using multiple logging parameters

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388531A (en) * 2015-10-19 2016-03-09 成都理工大学 Lithology identification method based on support vector regression machine and Kernel Fisher discriminant analysis
US9995698B2 (en) * 2015-10-20 2018-06-12 W.D. Gonten Laboratories LLC Integrated analysis of pore fluids for characterization of reservoir potential
CN110082840A (en) * 2019-06-05 2019-08-02 中国地质大学(北京) Mud shale lithofacies division methods based on Fisher discriminant analysis
CN110673215A (en) * 2019-09-19 2020-01-10 中国石油大学(华东) Automatic complex lithology distinguishing method based on intersection graph and Fisher distinguishing

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388531A (en) * 2015-10-19 2016-03-09 成都理工大学 Lithology identification method based on support vector regression machine and Kernel Fisher discriminant analysis
US9995698B2 (en) * 2015-10-20 2018-06-12 W.D. Gonten Laboratories LLC Integrated analysis of pore fluids for characterization of reservoir potential
CN110082840A (en) * 2019-06-05 2019-08-02 中国地质大学(北京) Mud shale lithofacies division methods based on Fisher discriminant analysis
CN110673215A (en) * 2019-09-19 2020-01-10 中国石油大学(华东) Automatic complex lithology distinguishing method based on intersection graph and Fisher distinguishing

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
博兴洼陷沙三下段湖相沉积岩性特征及识别;毛振强;《西安石油大学学报(自然科学版)》;20200930;第35卷(第5期);全文 *
测井曲线标准化方法研究;陈熹;《当代化工》;20150228;第44卷(第2期);全文 *

Also Published As

Publication number Publication date
CN111458767A (en) 2020-07-28

Similar Documents

Publication Publication Date Title
CN111458767B (en) A method and system for identifying lithology based on intersection map method
CN107703560B (en) shale lithofacies fine identification method based on triple information
CN104747185B (en) Heterogeneous reservoir reservoir synthetical assortment evaluation method
CN110717301B (en) A method for classifying and identifying flow unit information based on support vector machine algorithm
CN113792936A (en) A method, system, equipment and storage medium for intelligent identification of lithology while drilling
CN107576772A (en) A kind of method using log data quantitative assessment coal body structure type coal
CN110852144A (en) A DTW-based intelligent stratigraphic comparison method and system
CN112363242B (en) Reservoir fluid identification method and device based on logging fusion
CN109031461A (en) Carbonate porosity type oil and gas reservoir quantitative identification method and region quantitative method
CN112946782B (en) Earthquake fine depicting method for dense oil-gas storage seepage body
CN114137610B (en) Low exploration area stratum and lithologic seismic evaluation method combining supervision and unsupervised learning
CN116070147A (en) A method for analyzing shale lithofacies and structure based on grayscale facies
CN116127320A (en) Lost circulation analysis method, system and medium based on interpretable machine learning
CN111155980A (en) Water flow dominant channel identification method and device
CN113640875B (en) Deposition phase identification method, device and system
CN116413831A (en) Multi-well synthetic seismic record automatic calibration and construction interpretation and reservoir prediction method
CN113756796B (en) A Quantitative Evaluation Method for Aqueous Seepage Capacity of Formation Units in Waterflooding Development of Clastic Reservoirs
CN113325474B (en) Method for discriminating biological reef
CN117131650B (en) A method for identifying volcanic rocks and sedimentary rocks based on well logging curves
CN119226957B (en) A method for constructing a relationship model for identifying lithology of clastic rocks
CN118566981B (en) Oil and gas reservoir development well position management system based on data analysis
CN115932967B (en) Crack identification method based on ant body algorithm
CN112505773B (en) Fracture inheritance quantitative calculation and judgment method
CN117687119A (en) An automatic identification method and system for sedimentary microfacies of extremely thick layers of sandy conglomerate body
Burger et al. Geostatistics and the polygonal method: a re-examination

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Teng Juan

Inventor after: Xia Yu

Inventor after: Deng Hucheng

Inventor after: Chen Wenling

Inventor after: Wang Yuanyuan

Inventor after: Xie Xinhui

Inventor before: Xia Yu

Inventor before: Deng Hucheng

Inventor before: Teng Juan

Inventor before: Chen Wenling

Inventor before: Wang Yuanyuan

Inventor before: Xie Xinhui

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210604