CN113553546B - Method, system and computer readable storage medium for extracting rock continuous cementation index from electric imaging data - Google Patents
Method, system and computer readable storage medium for extracting rock continuous cementation index from electric imaging data Download PDFInfo
- Publication number
- CN113553546B CN113553546B CN202110795999.3A CN202110795999A CN113553546B CN 113553546 B CN113553546 B CN 113553546B CN 202110795999 A CN202110795999 A CN 202110795999A CN 113553546 B CN113553546 B CN 113553546B
- Authority
- CN
- China
- Prior art keywords
- index
- prh
- electric imaging
- imaging data
- proportion
- 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
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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/15—Correlation function computation including computation of convolution operations
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the technical field of oilfield development, in particular to a method, a system and a computer readable storage medium for extracting rock continuous cementation index from electric imaging data. A method of extracting rock continuity bond index from electrical imaging data, comprising: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The value and the electric imaging original data are preprocessed to obtain processed electric imaging data, and imaging is synthesized; counting and analyzing the preprocessed electric imaging data, and describing the continuous resistance patches and conductive spots of the whole well section; carrying out numerical quantification on the characteristics of the resistor patches and the conductive spots to obtain the resistor inclusion proportion PRH of each depth point of the whole well section; cementing index m of resistor inclusion proportion PRH and core sampling point experimental result 0 Obtaining a functional relation between the cementing index m and the resistor inclusion proportion PRH through value fitting; and substituting the PRH values of the depth points to calculate the cementing index m of the whole well section.
Description
Technical Field
The invention relates to the technical field of oilfield development, in particular to a method, a system and a computer readable storage medium for extracting rock continuous cementation index from electric imaging data.
Background
Along with the large-scale exploitation of unconventional tight oil and gas reservoirs, the progress of unconventional energy exploration and development is gradually accelerated, and meanwhile, great challenges are brought to well logging reservoir evaluation. Unconventional tight oil reservoirs often have the characteristics of strong diagenetic effect, low pore hypotonic property, complex pore structure and the like, so that fluid identification is particularly difficult. In the existing industry, a conventional saturation formula and a fixed cementation index are mostly adopted to calculate a fluid identification and saturation evaluation model, and as the reservoir heterogeneity of a tight reservoir is strong, the cementation index reflecting the reservoir heterogeneity is dynamically changed, the saturation calculated by the conventional Alqi formula can not accurately reflect the property of the reservoir fluid due to the fixed cementation index. Chinese patent application, publication No.: CN101649738A discloses a method for determining water saturation of a stratum, which mainly comprises the following steps: selecting a series of cores in a measurement area, which can represent geological features of the local area, and performing experiments and calculation to obtain the porosity phi of the cores; obtaining a stratum water resistivity R value saturated by the rock core through experiments and calculation; and obtaining a cementing index and a saturation index by adopting a public expression, and finally obtaining the water saturation of the stratum according to the cementing index and the saturation index. When the disclosed technical scheme solves the cementation index, a large amount of core sampling analysis is needed, so that more time and funds are consumed, and the rock cementation condition of all well sections of each well in a research area cannot be finely described.
Disclosure of Invention
The invention provides a method for extracting rock continuous cementation indexes by using electric imaging data, and also provides a system for extracting rock continuous cementation indexes by using electric imaging data and a computer readable storage medium, which aims to solve the problems that a large amount of core sampling analysis is needed and time and funds are wasted in the prior art.
In order to solve the technical problems, the invention adopts the following technical scheme: a method for extracting rock continuous bond index from electrical imaging data, comprising the steps of: s1: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The method comprises the steps of obtaining processed electric imaging data by preprocessing value and electric imaging original data, and synthesizing imaging according to the processed electric imaging data;
s2: counting and analyzing the electric imaging data preprocessed in the step S1, and describing the characteristics of continuous resistance patches and conductive spots of the whole well section;
s3: performing numerical quantification on the resistor patches and the conductive spot characteristics in the step S2 to obtain the resistor inclusion proportion PRH of each depth point of the whole well section;
s4: the resistor inclusion proportion PRH of the step S3 and the cementing index m of the core sampling point experimental result of the step S1 0 Obtaining a functional relation between the cementing index m and the resistance inclusion proportion PRH applicable to the whole well section through value fitting;
s5: substituting the PRH values of the depth points into the functional relation in the step S4 to calculate the cementing index m of the whole well section;
in step S3, the specific process of obtaining the resistance inclusion ratio PRH of each depth point of the full well section is as follows:
selecting a certain number of depth points as a window length, and properly cutting off the resistor patches and the conductive spots of the potential values of the transverse polar plates; maximum normalization is carried out on the values of the resistance patches and the conductive spots obtained by calculation of all the depth points in each window length, and the proportion of the resistance patches and the conductive spots is obtained, wherein the proportion of the resistance patches is the resistance inclusion proportion PRH of the points;
the functional relationship between the cementing index m and the resistance inclusion ratio PRH is:
m＝a*PRH+b
where a, b represents a coefficient, which can be obtained by the least squares method, a=0.2669 and b= 1.4766 in this embodiment.
In the technical scheme, the electric imaging data are utilized to describe continuous resistance patches and conductive spot characteristics of the whole well section, then the resistance patches and the conductive spot characteristics are quantified in numerical value, the parameter resistance containing proportion capable of more intuitively reflecting the heterogeneity of the reservoir is obtained, the cementing index suitable for the whole well section is obtained through the cementing index of the resistance containing proportion and the sampling point experimental result, and the cementing index of the whole well section is obtained through a functional relation formula of the cementing index and the resistance containing proportion. According to the technical scheme, the time and fund waste is reduced, a bridge of electric imaging data and reservoir heterogeneity and cementation indexes is established, and a more accurate saturation model is provided for subsequent reservoir fluid identification and quantitative calculation.
Preferably, in step S1, the preprocessing process includes bad data rejection, missing value processing, feature normalization and denoising processing on the button electrode potential, and shallow resistivity curves and scale button electrode measurement information are obtained according to the preprocessed electric imaging data synthesis.
Preferably, in step S2, the specific process of counting and analyzing the electrical imaging data preprocessed in step S1 is: calculating potential values of all the transverse polar plates in the electric imaging data, and longitudinally extending all the depth points to generate a characteristic curve of the potential values; and (3) describing the electrode plate potential by combining the composite imaging graph of the electric imaging and the potential value characteristic curve, and selecting a curve which can best show the potential change characteristic for normalization processing.
Preferably, the potential value includes a potential value average value, a potential value median and a potential value variance, the characteristic curve includes a potential value average value characteristic curve, a potential value median characteristic curve and a potential value variance characteristic curve, and the characteristic curve which best represents the potential change characteristic is selected for normalization processing by combining the electric imaging synthetic imaging graph and the potential value average value characteristic curve, the potential value median characteristic curve and the potential value variance characteristic curve.
Preferably, in step S2, the specific process of characterizing the resistive patches and conductive spots that are continuous throughout the well is: dividing intervals of effective current values by counting and analyzing the electric imaging data, and dividing the intervals of the potential values of each polar plate based on the potential value distribution condition of each transverse polar plate of all depth points of the whole well section and the curve of normalization processing; and according to the distribution result of the intervals of the effective current values on the polar plates, the interval of the current value with the largest proportion is defined, the proportion of the minimum current value in the remaining interval is calculated, the proportion of the minimum current value is the characteristic parameter of the resistor patch, and the proportion of the maximum current value interval in the interval of the current value is the characteristic parameter of the conductive spot.
Preferably, in step S3, the specific process of obtaining the resistance inclusion ratio PRH of each depth point of the full interval is: selecting a certain number of depth points as a window length, and properly cutting off the resistor patches and the conductive spots of the potential values of the transverse polar plates; and carrying out maximum normalization on the values of the resistance patches and the conductive spots obtained by calculating all the depth points in each window length to obtain the proportion of the resistance patches and the conductive spots, wherein the proportion of the resistance patches is the resistance inclusion proportion PRH of the points.
Preferably, a continuously varying curve is constructed based on all of the cementation indexes m of step S5.
Preferably, the well logging data further includes a formation factor versus porosity plot and a log interpretation result plot.
The invention also provides a system for extracting rock continuous cementation index from electric imaging data, which comprises a memory and a processor, wherein the memory comprises a method program for extracting rock continuous cementation index from electric imaging data, and the memory comprises the following steps when the method program for extracting rock continuous cementation index from electric imaging data is executed by the processor:
s1: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The method comprises the steps of obtaining processed electric imaging data by preprocessing value and electric imaging original data, and synthesizing imaging according to the processed electric imaging data;
s2: counting and analyzing the electric imaging data preprocessed in the step S1, and describing the characteristics of continuous resistance patches and conductive spots of the whole well section;
s3: performing numerical quantification on the resistor patches and the conductive spot characteristics in the step S2 to obtain the resistor inclusion proportion PRH of each depth point of the whole well section;
s4: the cementing index m of the core sampling point experimental result is obtained by the resistor inclusion ratio PRH described in the step S3 and the core sampling point experimental result described in the step S1 0 Obtaining a functional relation between the cementing index m and the resistor inclusion proportion PRH applicable to the whole block through value fitting;
s5: substituting the PRH values of the depth points into the functional relation in the step S4 to calculate the cementing index m of the whole well section;
in step S3, the specific process of obtaining the resistance inclusion ratio PRH of each depth point of the full well section is as follows:
selecting a certain number of depth points as a window length, and properly cutting off the resistor patches and the conductive spots of the potential values of the transverse polar plates; maximum normalization is carried out on the values of the resistance patches and the conductive spots obtained by calculation of all the depth points in each window length, and the proportion of the resistance patches and the conductive spots is obtained, wherein the proportion of the resistance patches is the resistance inclusion proportion PRH of the points;
the functional relationship between the cementing index m and the resistance inclusion ratio PRH is:
m＝a*PRH+b
where a, b represents a coefficient, which can be obtained by the least squares method, a=0.2669 and b= 1.4766 in this embodiment.
The invention also provides a computer readable storage medium, which comprises the method program for extracting rock continuous cementation index from the electric imaging data.
Compared with the prior art, the invention has the beneficial effects that:
according to the invention, the electric imaging data are utilized to describe the continuous resistance patch and conductive spot characteristics of the whole well section, then the numerical quantization is carried out on the resistance patch and the conductive spot characteristics, the parameter resistance inclusion proportion capable of more intuitively reflecting the heterogeneity of the reservoir is obtained, and the cementing index applicable to the whole well section is further obtained through the functional relation formula of the cementing index and the resistance inclusion proportion of the experimental result of the resistance inclusion proportion and the sampling point, so that the continuous cementing index curve is further obtained, and the cementing index of the whole well section is obtained. The invention reduces the time and fund waste, establishes a bridge of electric imaging data and reservoir heterogeneity and cementation index, and provides a more accurate saturation model for subsequent reservoir fluid identification and quantitative calculation.
Drawings
FIG. 1 is a schematic flow chart of a method for extracting rock continuous bond index from electrical imaging data of example 1 of the present invention;
FIG. 2 is a diagram of the pre-processed electrical imaging data and composite imaging of example 1 of the present invention;
FIG. 3 is a graph showing the potential value distribution of the lateral plate according to example 1 of the present invention;
FIG. 4 is a graph of the PRH value of the resistor inclusion ratio in example 1 of the present invention;
FIG. 5 is a graph of m values of the continuous bond index in example 1 of the present invention.
Detailed Description
The drawings are for illustrative purposes only and are not to be construed as limiting the present patent; for the purpose of better illustrating the embodiments, certain elements of the drawings may be omitted, enlarged or reduced and do not represent the actual product dimensions; it will be appreciated by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted. The positional relationship depicted in the drawings is for illustrative purposes only and is not to be construed as limiting the present patent.
The same or similar reference numbers in the drawings of embodiments of the invention correspond to the same or similar components; in the description of the present invention, it should be understood that, if there are orientations or positional relationships indicated by terms "upper", "lower", "left", "right", "long", "short", etc., based on the orientations or positional relationships shown in the drawings, this is merely for convenience in describing the present invention and simplifying the description, and is not an indication or suggestion that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, so that the terms describing the positional relationships in the drawings are merely for exemplary illustration and are not to be construed as limitations of the present patent, and that it is possible for those of ordinary skill in the art to understand the specific meaning of the terms described above according to specific circumstances.
The technical scheme of the invention is further specifically described by the following specific embodiments with reference to the accompanying drawings:
example 1
As shown in fig. 1, a method for extracting rock continuous cementation index from electric imaging data comprises the following steps:
s1: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The method comprises the steps of obtaining processed electric imaging data by preprocessing value and electric imaging original data, and synthesizing imaging according to the processed electric imaging data;
s2: counting and analyzing the electric imaging data preprocessed in the step S1, and describing the characteristics of continuous resistance patches and conductive spots of the whole well section;
s3: performing numerical quantification on the resistor patches and the conductive spot characteristics in the step S2 to obtain the resistor inclusion proportion PRH of each depth point of the whole well section;
s4: the resistance inclusion proportion PRH of the step S3 and the cementing index m of the core sampling point experimental result of the step S1 0 Obtaining a functional relation between the cementing index m and the resistance inclusion proportion PRH applicable to the whole well section through value fitting;
s5: and (3) substituting the PRH values of the depth points according to the functional relation in the step S4 to calculate the cementing index m of the whole well section. In this embodiment, the electrical imaging data is used to characterize the continuous resistive patch and conductive spot feature of the whole well section, then the resistive patch and conductive spot feature are quantified numerically to obtain the parameter resistance inclusion proportion capable of more intuitively reflecting the reservoir heterogeneity, and further the cementing index is obtained by using the resistance inclusion proportion and the sampling point experimental result cementing index to obtain a functional relation between the cementing index and the resistance inclusion proportion suitable for the whole well section, and the cementing index of the whole well section is obtained. The invention reduces the waste of time and funds, establishes a bridge of electric imaging data and reservoir heterogeneity and cementation index, provides a more accurate saturation model for the subsequent reservoir fluid identification and quantitative calculation,
meanwhile, in this embodiment, the functional relation between the cementing index m and the resistance inclusion ratio PRH is:
m＝a*PRH+b
where a, b represents a coefficient, which can be obtained by the least squares method, a=0.2669 and b= 1.4766 in this embodiment.
In step S1, the preprocessing process includes removing bad data of button electrode potential, processing missing value, normalizing features, denoising, and synthesizing to obtain shallow resistivity curve and scale button electrode measurement information according to preprocessed electric imaging data.
In addition, as shown in fig. 2, where X represents a circumferential plate value and Y represents a longitudinal depth, in step S2, a specific procedure for counting and analyzing the electrical imaging data preprocessed in step S1 is as follows: calculating potential values of all the transverse polar plates in the electric imaging data, and longitudinally extending all the depth points to generate a characteristic curve of the potential values; and (3) describing the electrode plate potential by combining the composite imaging graph of the electric imaging and the potential value characteristic curve, and selecting a curve which can best show the potential change characteristic for normalization processing.
The characteristic curves comprise a potential value average value characteristic curve, a potential value median characteristic curve and a potential value variance characteristic curve, and the characteristic curves which can represent the potential change characteristic most are selected for normalization processing by combining the synthesized imaging graph of the electric imaging with the potential value average value characteristic curve, the potential value median characteristic curve and the potential value variance characteristic curve.
In addition, as shown in fig. 3, four conditions of the polar plate potential distribution diagram are shown, the abscissa is the potential value, and the ordinate is the frequency, and in step S2, the specific process of characterizing the continuous resistance patches and conductive spots of the whole well section is as follows: dividing the interval of the effective current value by using statistical and analyzed electric imaging data, and dividing the interval of the potential value of each polar plate based on the potential value distribution condition of each polar plate in the transverse direction of all depth points of the whole well section and a curve processed by normalization; and according to the distribution result of the intervals of the effective current values on the polar plates, the current value interval with the largest occupation ratio is defined, the proportion of the minimum current value in the residual interval is calculated, the proportion of the minimum current value is the characteristic parameter of the resistor patch, and the proportion of the maximum current value interval in the current value interval is the characteristic parameter of the conductive spot.
As shown in fig. 4, in step S3, the specific process of obtaining the resistance inclusion ratio PRH of each depth point of the whole well section is as follows: selecting a certain number of depth points as a window length, and properly cutting off the resistor patches and the conductive spots of the potential values of the transverse polar plates; and carrying out maximum normalization on the values of the resistance patches and the conductive spots obtained by calculating all depth points in each window length to obtain the proportion of the resistance patches and the conductive spots, wherein the proportion of the resistance patches is the resistance inclusion proportion PRH of the points.
In addition, as shown in fig. 5, a continuously varying curve is constructed based on all the cement indexes m of step S5. In this embodiment, the continuous consolidated index profile contains more information about formation heterogeneity.
In addition, the logging data also comprises a stratum factor and porosity relation chart and a logging interpretation result chart.
The working principle is as follows;
according to the invention, the electric imaging data are utilized to describe the continuous resistance patch and conductive spot characteristics of the whole well section, then the numerical quantization is carried out on the resistance patch and the conductive spot characteristics, the parameter resistance inclusion proportion capable of more intuitively reflecting the heterogeneity of the reservoir is obtained, and the cementing index applicable to the whole well section is further obtained through the functional relation formula of the cementing index and the resistance inclusion proportion of the experimental result of the resistance inclusion proportion and the sampling point, so that the continuous cementing index curve is further obtained, and the cementing index of the whole well section is obtained. The invention reduces the time and fund waste, establishes a bridge of electric imaging data and reservoir heterogeneity and cementation index, and provides a more accurate saturation model for subsequent reservoir fluid identification and quantitative calculation.
Example 2
A system for extracting rock continuity bond index from electroimaging data, comprising a memory and a processor, the memory comprising a method program for extracting rock continuity bond index from electroimaging data, when executed by the processor, effecting the steps of:
s1: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The method comprises the steps of obtaining processed electric imaging data by preprocessing value and electric imaging original data, and synthesizing imaging according to the processed electric imaging data;
s2: counting and analyzing the electric imaging data preprocessed in the step S1, and describing the characteristics of continuous resistance patches and conductive spots of the whole well section;
s3: performing numerical quantification on the resistor patches and the conductive spot characteristics in the step S2 to obtain the resistor inclusion proportion PRH of each depth point of the whole well section;
s4: the cementing index m is obtained by the core sampling point experimental result of the resistor inclusion ratio PRH and the core sampling point experimental result of the step S3 0 Obtaining a functional relation between the cementing index m and the resistor inclusion proportion PRH applicable to the whole block through value fitting;
s5: and (3) substituting the PRH values of the depth points according to the functional relation in the step S4 to calculate the cementing index m of the whole well section.
Example 3
A computer readable storage medium comprising a method program for extracting rock continuity indicator from electrical imaging data.
It is to be understood that the above examples of the present invention are provided by way of illustration only and not by way of limitation of the embodiments of the present invention. Other variations or modifications of the above teachings will be apparent to those of ordinary skill in the art. It is not necessary here nor is it exhaustive of all embodiments. Any modification, equivalent replacement, improvement, etc. which come within the spirit and principles of the invention are desired to be protected by the following claims.
Claims (9)
1. A method for extracting rock continuous cementation index from electric imaging data, which is mainly characterized by comprising the following steps:
s1: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The method comprises the steps of obtaining processed electric imaging data by preprocessing value and electric imaging original data, and synthesizing imaging according to the processed electric imaging data;
s2: counting and analyzing the electric imaging data preprocessed in the step S1, and describing the characteristics of continuous resistance patches and conductive spots of the whole well section;
s3: performing numerical quantification on the resistor patches and the conductive spot characteristics in the step S2 to obtain the resistor inclusion proportion PRH of each depth point of the whole well section;
s4: the resistor inclusion proportion PRH of the step S3 and the cementing index m of the core sampling point experimental result of the step S1 0 Obtaining a functional relation between the cementing index m and the resistance inclusion proportion PRH applicable to the whole well section through value fitting;
s5: substituting the PRH values of the depth points into the functional relation in the step S4 to calculate the cementing index m of the whole well section;
in step S3, the specific process of obtaining the resistance inclusion ratio PRH of each depth point of the full well section is as follows: selecting a certain number of depth points as a window length, and properly cutting off the resistor patches and the conductive spots of the potential values of the transverse polar plates; maximum normalization is carried out on the values of the resistance patches and the conductive spots obtained by calculation of all the depth points in each window length, and the proportion of the resistance patches and the conductive spots is obtained, wherein the proportion of the resistance patches is the resistance inclusion proportion PRH of the points;
the functional relationship between the cementing index m and the resistance inclusion ratio PRH is:
m＝a*PRH+b
where a, b represents a coefficient, which can be obtained by the least squares method, a=0.2669 and b= 1.4766 in this embodiment.
2. The method for extracting rock continuous bond index from electric imaging data according to claim 1, mainly characterized in that: in step S1, the preprocessing process includes removing bad data of the button electrode potential, processing missing values, normalizing features, denoising, and synthesizing according to the preprocessed electric imaging data to obtain a shallow resistivity curve and scale button electrode measurement information.
3. The method for extracting rock continuous bond index from electric imaging data according to claim 1, wherein in step S2, the specific process of counting and analyzing the electric imaging data pretreated in step S1 is as follows: calculating potential values of all the transverse polar plates in the electric imaging data, and longitudinally extending all the depth points to generate a characteristic curve of the potential values; and (3) describing the electrode plate potential by combining the composite imaging graph of the electric imaging and the potential value characteristic curve, and selecting a curve which can best show the potential change characteristic for normalization processing.
4. A method of extracting rock continuous bond index from electrical imaging data according to claim 3, characterized in that: the characteristic curves comprise a potential value average value characteristic curve, a potential value median characteristic curve and a potential value variance characteristic curve, and the characteristic curves which can best represent potential change characteristics are selected for normalization processing by combining a synthetic imaging diagram of electric imaging with the potential value average value characteristic curve, the potential value median characteristic curve and the potential value variance characteristic curve.
5. A method for extracting rock continuous cementation index from electric imaging data according to claim 3, wherein in step S2, the specific process of characterizing the continuous resistive patches and conductive spots of the whole well section is as follows: dividing intervals of effective current values by counting and analyzing the electric imaging data, and dividing the intervals of the potential values of each polar plate based on the potential value distribution condition of each transverse polar plate of all depth points of the whole well section and the curve of normalization processing; and according to the distribution result of the intervals of the effective current values on the polar plates, the interval of the current value with the largest proportion is defined, the proportion of the minimum current value in the remaining interval is calculated, the proportion of the minimum current value is the characteristic parameter of the resistor patch, and the proportion of the maximum current value interval in the interval of the current value is the characteristic parameter of the conductive spot.
6. The method for extracting rock continuous bond index from electric imaging data according to claim 1, mainly characterized in that: a continuously varying curve is constructed based on all of the cementation indexes m of step S5.
7. The method for extracting rock continuous bond index from electric imaging data according to claim 1, mainly characterized in that: the logging data also includes a map of formation factors versus porosity and a log interpretation result chart.
8. A system for extracting rock continuity bond index from electroimaging data, comprising a memory and a processor, the memory comprising a method program for extracting rock continuity bond index from electroimaging data, when executed by the processor, effecting the steps of:
s1: collecting logging data, wherein the logging data comprises a core sampling point experimental result cementing index m 0 The method comprises the steps of obtaining processed electric imaging data by preprocessing value and electric imaging original data, and synthesizing imaging according to the processed electric imaging data;
s2: counting and analyzing the electric imaging data preprocessed in the step S1, and describing the characteristics of continuous resistance patches and conductive spots of the whole well section;
s3: performing numerical quantification on the resistor patches and the conductive spot characteristics in the step S2 to obtain the resistor inclusion proportion PRH of each depth point of the whole well section;
s4: the cementing index m of the core sampling point experimental result is obtained by the resistor inclusion ratio PRH described in the step S3 and the core sampling point experimental result described in the step S1 0 Obtaining a functional relation between the cementing index m and the resistor inclusion proportion PRH applicable to the whole block through value fitting;
s5: substituting the PRH values of the depth points into the functional relation in the step S4 to calculate the cementing index m of the whole well section;
in step S3, the specific process of obtaining the resistance inclusion ratio PRH of each depth point of the full well section is as follows: selecting a certain number of depth points as a window length, and properly cutting off the resistor patches and the conductive spots of the potential values of the transverse polar plates; maximum normalization is carried out on the values of the resistance patches and the conductive spots obtained by calculation of all the depth points in each window length, and the proportion of the resistance patches and the conductive spots is obtained, wherein the proportion of the resistance patches is the resistance inclusion proportion PRH of the points;
the functional relationship between the cementing index m and the resistance inclusion ratio PRH is:
m＝a*PRH+b
where a, b represents a coefficient, which can be obtained by the least squares method, a=0.2669 and b= 1.4766 in this embodiment.
9. A computer readable storage medium, comprising a method program for implementing an electroimaging data extraction rock continuous bond index according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110795999.3A CN113553546B (en) | 2021-07-14 | 2021-07-14 | Method, system and computer readable storage medium for extracting rock continuous cementation index from electric imaging data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110795999.3A CN113553546B (en) | 2021-07-14 | 2021-07-14 | Method, system and computer readable storage medium for extracting rock continuous cementation index from electric imaging data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113553546A CN113553546A (en) | 2021-10-26 |
| CN113553546B true CN113553546B (en) | 2023-08-18 |
Family
ID=78131807
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110795999.3A Active CN113553546B (en) | 2021-07-14 | 2021-07-14 | Method, system and computer readable storage medium for extracting rock continuous cementation index from electric imaging data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113553546B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009126881A2 (en) * | 2008-04-10 | 2009-10-15 | Services Petroliers Schlumberger | Method to generate numerical pseudocores using borehole images, digital rock samples, and multi-point statistics |
| CN101929973A (en) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | Quantitative calculation method for hydrocarbon saturation of fractured reservoir |
| CN104453874A (en) * | 2014-10-23 | 2015-03-25 | 中国石油天然气集团公司 | Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance |
| CN104730095A (en) * | 2015-02-12 | 2015-06-24 | 中国石油大学(华东) | Complex formation clay content measuring method based on electric imaging logging information |
-
2021
- 2021-07-14 CN CN202110795999.3A patent/CN113553546B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009126881A2 (en) * | 2008-04-10 | 2009-10-15 | Services Petroliers Schlumberger | Method to generate numerical pseudocores using borehole images, digital rock samples, and multi-point statistics |
| CN101929973A (en) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | Quantitative calculation method for hydrocarbon saturation of fractured reservoir |
| CN104453874A (en) * | 2014-10-23 | 2015-03-25 | 中国石油天然气集团公司 | Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance |
| CN104730095A (en) * | 2015-02-12 | 2015-06-24 | 中国石油大学(华东) | Complex formation clay content measuring method based on electric imaging logging information |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113553546A (en) | 2021-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8583410B2 (en) | Method for obtaining consistent and integrated physical properties of porous media | |
| CN106761729B (en) | Multi-parameter interpretation method for low-permeability conglomerate logging | |
| CN108694264B (en) | Method for determining permeability of shale gas reservoir | |
| CN103670384A (en) | Reservoir classification method and system | |
| CN104863574A (en) | Fluid identification method applicable to compact sandstone reservoir | |
| CN112145165B (en) | Microcrack-pore type reservoir dynamic and static permeability conversion method | |
| CN106777649B (en) | Quantitative evaluation method for pore structure of fractured reservoir | |
| CN112963145B (en) | Method for predicting capacity of carbonate reservoir gas well | |
| CN113553546B (en) | Method, system and computer readable storage medium for extracting rock continuous cementation index from electric imaging data | |
| CN110017136B (en) | Water flooded layer identification and water production rate prediction method based on apparent water layer resistivity | |
| CN114086938B (en) | Method for predicting gas saturation of heterogeneous sandstone reservoir | |
| CN112505154B (en) | Shale reservoir mineral component content analysis and lithofacies identification characterization method | |
| CN114060015B (en) | Method and device for evaluating gas content of compact sandstone | |
| CN110873904B (en) | Fluid identification method and device | |
| CN114428365A (en) | Gas-containing characteristic-based tight sandstone gas reservoir logging identification method and device | |
| CN113720745A (en) | Method for calculating porosity of reservoir stratum containing carbon debris by geophysical logging | |
| CN114991745B (en) | Shale oil reservoir dessert identification method and device | |
| CN112922587B (en) | Method and device for identifying salt frost position | |
| CN111827989B (en) | Method for identifying asphalt-containing conglomerate oil layer by asphalt cutting index | |
| CN114428049B (en) | Method for calculating asphalt content of ancient carbonate reservoir | |
| CN110735634B (en) | Method and device for determining dynamic permeability of limestone oil reservoir | |
| CN116148942A (en) | Deep high-quality sandstone reservoir identification method based on logging curve fusion | |
| CN112668814A (en) | Cementing grade determining method and device | |
| Filiptsova et al. | Determining petrophysical and hydrogeological parameters from historical bore logs for the Leederville-Parmelia aquifer, northern Perth Basin, using regression methods | |
| Rushton | Impact of aquitard storage on leaky aquifer pumping test analysis |
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 |