CN108318534B - Core-constrained electrical imaging logging image processing method and device - Google Patents

Core-constrained electrical imaging logging image processing method and device Download PDF

Info

Publication number
CN108318534B
CN108318534B CN201711362904.9A CN201711362904A CN108318534B CN 108318534 B CN108318534 B CN 108318534B CN 201711362904 A CN201711362904 A CN 201711362904A CN 108318534 B CN108318534 B CN 108318534B
Authority
CN
China
Prior art keywords
resistivity
color
component value
current depth
depth window
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
Application number
CN201711362904.9A
Other languages
Chinese (zh)
Other versions
CN108318534A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201711362904.9A priority Critical patent/CN108318534B/en
Publication of CN108318534A publication Critical patent/CN108318534A/en
Application granted granted Critical
Publication of CN108318534B publication Critical patent/CN108318534B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/001—Texturing; Colouring; Generation of texture or colour
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs

Abstract

The embodiment of the application provides a method and a device for processing a core-constrained electrical imaging logging image, wherein the method comprises the following steps: acquiring electrical imaging logging data and logging data; determining an electrical imaging logging image resistivity distribution diagram of the target area according to the electrical imaging logging information; establishing a conversion relation between the resistivity and the rock color according to the resistivity distribution map and logging information of the electrical imaging logging image; and converting the resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to the conversion relation between the resistivity and the rock color. According to the scheme, the conversion relation between the resistivity and the rock color is established firstly through the characteristics of the electric imaging logging data and the logging data; and then, converting the resistivity distribution map of the electrical imaging logging image into an optical logging image by utilizing the conversion relation so as to replace the optical imaging image. Therefore, the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor in the existing method are solved.

Description

Core-constrained electrical imaging logging image processing method and device
Technical Field
The application relates to the technical field of oil and gas exploration, in particular to a method and a device for processing a core-constrained electrical imaging logging image.
Background
In the course of oil and gas exploration, it is often necessary to acquire optical imaging images of the core in the target area. The optical imaging image is usually obtained through continuous optical imaging, and is different from a common electrical logging imaging image, the optical imaging image also carries characteristic information such as reservoir color, and the like, and reservoir evaluation can be performed on a target area more intuitively and accurately by using the characteristic information.
Currently, in order to obtain an optical imaging image, drilling and coring are generally performed first, and then a core scanning system is used to perform continuous optical scanning imaging on a cored rock sample so as to obtain a corresponding optical imaging image (i.e., a core scanning image). However, in practice, the cost of coring while drilling is generally high; in addition, in specific implementation, limited by cost and other objective factors, only the key interval in the target region can be cored, and therefore, the continuity of the optical imaging image obtained by the existing method is relatively poor. In addition, for a fracture development well section, it is often difficult to obtain a relatively complete core, resulting in a poor characterization effect of obtaining an optical imaging image. In summary, when the existing method is implemented, the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor often exist.
In view of the above problems, no effective solution has been proposed.
Disclosure of Invention
The embodiment of the application provides a core-constrained electrical imaging logging image processing method and device, aims to solve the technical problems that the cost of determining an optical imaging image is high and the continuity of the determined optical imaging image is poor in the existing method, and achieves the technical effect that the optical logging image capable of well representing the integral geological characteristics of a reservoir can be obtained at low cost and high efficiency.
The embodiment of the application provides a rock core constrained electrical imaging logging image processing method, which comprises the following steps:
acquiring electrical imaging logging information and logging information of a target area;
determining an electrical imaging logging image resistivity distribution diagram of a target area according to the electrical imaging logging information;
establishing a conversion relation between the resistivity and the rock color according to the electrical imaging logging image resistivity distribution map and the logging information;
and converting the resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to the conversion relation between the resistivity and the rock color.
In one embodiment, determining an electrical imaging log resistivity profile of the target region from the electrical imaging log data comprises:
and carrying out shallow resistivity scale and whole borehole image generation processing on the electrical imaging logging information to obtain the electrical imaging logging image resistivity distribution diagram.
In one embodiment, establishing a conversion relationship between resistivity and rock color according to the electrical imaging logging image resistivity profile and the logging information comprises:
determining the color and lithology information of continuous stratums in the target area according to the logging information;
dividing the electrical imaging logging image resistivity distribution map into a plurality of depth windows according to the depth;
and respectively establishing a conversion relation between the resistivity and the rock color corresponding to each depth window in the plurality of depth windows according to the color and lithology information of the continuous stratum.
In one embodiment, respectively establishing a conversion relationship between the resistivity and the rock color corresponding to each depth window of the plurality of depth windows according to the color and lithology information of the continuous stratum includes:
establishing a conversion relation between the resistivity corresponding to the current depth window and the rock color according to the following mode:
determining a first component value corresponding to a current depth window, a second component value corresponding to the current depth window and a third component value corresponding to the current depth window according to the color and lithology information of the continuous stratum, wherein the first component value, the second component value and the third component value correspond to different color component values respectively;
dividing the current depth window into a plurality of resistivity intervals according to the resistivity;
and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
In one embodiment, the first component value corresponding to the current depth window is a red component value corresponding to the current depth window, the second component value corresponding to the current depth window is a green component value corresponding to the current depth window, and the third component value corresponding to the current depth window is a blue component value corresponding to the current depth window.
In one embodiment, the determining the rock color corresponding to each of the resistivity intervals comprises:
determining the rock color corresponding to the current resistivity interval in the resistivity intervals according to the following formula:
Figure BDA0001512214190000031
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrFor the red component value, R, corresponding to the current depth windowgFor the green component value, R, corresponding to the current depth windowbThe blue component value corresponding to the current depth window.
In one embodiment, converting the electrical imaging log resistivity profile into a corresponding optical log image according to the relationship between the resistivity and the rock color comprises:
and respectively converting pixel points in the resistivity interval in the electrical imaging logging image resistivity distribution diagram into corresponding rock colors according to the rock colors corresponding to the resistivity interval.
In one embodiment, after converting the electrical imaging log resistivity profile to a corresponding optical log image, the method further comprises:
and according to the optical logging image, performing reservoir evaluation on a target area.
In one embodiment, the determining, according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window, the rock color corresponding to each of the resistivity intervals further includes:
acquiring coring data of a target area;
determining the color of the slot hole filler and the color information of the massive breccia according to the coring data;
and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the color of the slot hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
In one embodiment, determining, according to the color of the hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window, the rock color corresponding to each of the resistivity intervals in the plurality of resistivity intervals respectively includes:
determining the rock color corresponding to the current resistivity interval in the resistivity intervals according to the following formula:
Figure BDA0001512214190000041
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrIs a first component value, R, corresponding to the current depth windowgA second component value, R, corresponding to the current depth windowbA third component value corresponding to the current depth window, cr1A fourth component value corresponding to the current depth window, cg1A fifth value of the component corresponding to the current depth window, cb1A sixth component value, c, corresponding to the current depth windowr2A seventh component value, c, corresponding to the current depth windowg2For an eighth component value, c, corresponding to the current depth windowb2Is the ninth component value corresponding to the current depth window, wherein,cr1、cg1、cb1、cr2、cg2、cb2is determined according to the color of the slot hole filling material and the color information of the massive breccia.
The application also provides an electrical imaging logging image processing apparatus of rock core restraint, includes:
the acquisition module is used for acquiring electrical imaging logging data and logging data of a target area;
the determining module is used for determining an electrical imaging logging image resistivity distribution map of the target area according to the electrical imaging logging information;
the establishing module is used for establishing a conversion relation between the resistivity and the rock color according to the electrical imaging logging image resistivity distribution map and the logging information;
and the conversion module is used for converting the resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to the conversion relation between the resistivity and the rock color.
In one embodiment, the establishing module comprises:
the determining unit is used for determining the color and lithology information of the continuous stratum in the target area according to the logging information;
the dividing unit is used for dividing the electrical imaging logging image resistivity distribution map into a plurality of depth windows according to the depth;
and the establishing unit is used for respectively establishing the conversion relation between the resistivity and the rock color corresponding to each depth window in the plurality of depth windows according to the color and lithology information of the continuous stratum.
In the embodiment of the application, the conversion relation between the resistivity and the rock color is established firstly by integrating the characteristics of the electrical imaging logging data and the logging data; and then, converting the resistivity distribution map of the electrical imaging logging image into an optical logging image by utilizing the conversion relation so as to replace the optical imaging image. Therefore, the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor in the existing method are solved, and the technical effect that the optical logging image capable of well representing the integral geological characteristics of the reservoir can be obtained at low cost and high efficiency is achieved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort.
FIG. 1 is a process flow diagram of a method for processing a core-constrained electrical imaging log image provided in accordance with an embodiment of the present application;
FIG. 2 is a block diagram of a core-constrained electrographic logging image processing apparatus provided in accordance with an embodiment of the present application;
FIG. 3 is a schematic diagram of results obtained after resistivity scale and whole borehole image processing by applying the method and apparatus for electrical imaging log image processing for core constraint provided by embodiments of the present application in one example scenario;
FIG. 4 is a schematic diagram of a resistivity profile and a resistivity interval obtained by applying the method and apparatus for processing electrical imaging log images for core constraint provided by the embodiments of the present application in one example scenario;
FIG. 5 is a graphical representation of a comparison of core scan images simulated using optical electrical logging images obtained using the method and apparatus for electrical imaging logging image processing for core constraint provided by embodiments of the present application in one example scenario.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. 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 application.
In consideration of the fact that the drilling and coring cost is relatively high when the existing optical imaging image acquisition method is implemented, and the existing method is limited by coring rock samples, the continuity of the optical imaging image obtained by the existing method is relatively poor, and for a fracture development well section, the characterization effect of the obtained optical imaging image is poor due to the fact that relatively complete rock cores are difficult to obtain. In summary, when the existing method is implemented, the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor often exist. Aiming at the root cause of the technical problems, the method considers that the characteristics of the imaging logging information and the logging information can be effectively utilized to establish the conversion relation between the resistivity and the rock color, and further can determine the optical logging image capable of representing the color, and the optical logging image is utilized to replace the optical imaging image, so that the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor in the existing method are solved, and the technical effect of efficiently obtaining the optical logging image capable of well representing the integral geological features of the reservoir with low cost is achieved.
Based on the thought, the embodiment of the application provides a rock core constrained electrical imaging logging image processing method. Specifically, refer to fig. 1, which is a flow chart illustrating a processing method of a core-constrained electrical imaging log image according to an embodiment of the present disclosure. The method for processing the core-constrained electrical imaging logging image provided by the embodiment of the application can be implemented specifically by the following steps.
S11: and acquiring the electrical imaging logging data and logging data of the target area.
In this embodiment, the electrical imaging logging information may be logging information obtained by entering a corresponding device to perform electrical imaging logging acquisition in a drilling well of the target area. Compared with other conventional logging information, the electrical imaging logging information usually has higher longitudinal resolution and experience coverage rate, and is suitable for specific geological research and reservoir evaluation. Specifically, a resistivity change image (namely, a resistivity distribution map of the electrical imaging logging image) of the well wall can be measured and obtained through electrical imaging logging, wherein the resistivity change image can generally have good response to various fracture and cavity bodies, deposition, structural characteristics and the like, has good continuity, and can well represent the integral geological characteristics of the reservoir. However, the resistivity change image is different from an optical imaging image based on core scanning, and cannot reflect the formation rock characteristics in the well based on optical principles, for example, the specific rock of the formation rock cannot be visually characterized. Therefore, it is generally not possible to directly replace the optical imaging image with the electrical imaging image.
In this embodiment, the logging information may be logging cuttings obtained during a drilling process. The logging rock debris can be used as a rock sample in logging, and can accurately reflect some basic information of rocks in the well. Such as the color, lithology, etc. of the reservoir rock in the log.
S12: and determining the resistivity distribution map of the electrical imaging logging image of the target area according to the electrical imaging logging information.
In an embodiment, the determining the electrical imaging log resistivity profile of the target region according to the electrical imaging log data may specifically include: and carrying out shallow resistivity scale and whole borehole image generation processing on the electrical imaging logging information to obtain the electrical imaging logging image resistivity distribution diagram. It should be noted that the electrical imaging log resistivity profile of the target region determined above can be considered as an electrical imaging log. However, the above-mentioned electrical imaging log images can only characterize the correlation of the resistivity in the target region, and cannot characterize the condition reflected by the corresponding optical imaging image.
In this embodiment, the performing the shallow resistivity calibration process on the electrical imaging logging data may specifically include: and correspondingly processing the resistivity and the depth of the resistivity image in the electrical imaging logging data.
In the embodiment, it is considered that the conventional electrical imaging logging is mostly measured in a manner that the electrode plate is attached to the borehole wall, and therefore, the obtained data cannot completely cover the whole borehole (usually, the coverage rate can only reach 60% to 80%), and therefore, the whole borehole image generation processing needs to be performed on the electrical imaging logging data. The above-mentioned processing of producing the whole borehole image of the electrical imaging logging data may specifically include: and performing probability difference on uncovered parts (namely parts without numerical values) in the resistivity image after calibration so as to obtain a resistivity distribution map covering the whole borehole, namely the resistivity distribution map of the electrical imaging logging image.
In one embodiment, before determining the electrical imaging log resistivity profile of the target region according to the electrical imaging log data, the electrical imaging log data may be preprocessed to obtain a higher-precision electrical imaging log resistivity profile. Wherein, the pretreatment specifically comprises: acceleration correction processing, and/or equalization processing, and the like.
S13: and establishing a conversion relation between the resistivity and the rock color according to the electrical imaging logging image resistivity distribution map and the logging information.
In one embodiment, in order to integrate the advantages of the electrical imaging log image (corresponding to the resistivity distribution map of the electrical imaging log image in this embodiment) and the optical imaging image, both the electrical imaging log image has better continuity and the optical imaging image can intuitively and accurately reflect the rock characteristics of the reservoir, such as the color of the rock, so that the conversion relationship between the resistivity and the color of the rock can be established first; and then the conversion relation between the resistivity and the rock color can be utilized, the electric imaging logging image is taken as a basis, and the electric imaging logging image is converted into an optical logging image which can accurately and intuitively reflect the characteristics of the reservoir rock color and the like. In specific implementation, in order to establish the above-mentioned conversion relationship between resistivity and rock color, the following steps can be performed.
S13-1: and determining the color and lithology information of the continuous stratum in the target area according to the logging information.
In this embodiment, in a specific implementation, the logging rock debris in the logging information may be described to determine characteristic information such as rock color, lithology, oil-gas content, and the like of various rocks in the formation in the target region; and performing classified statistics on the characteristic information to obtain the continuous stratum rock color, lithology and hydrocarbon-containing information of the whole well section, namely the color and lithology information of the continuous stratum in the target area.
S13-2: and dividing the electrical imaging logging image resistivity distribution graph into a plurality of depth windows according to the depth.
In this embodiment, in specific implementation, the imaging resistivity distribution map may be divided into a plurality of depth windows corresponding to different depths according to the distribution of the depths of the formation, so that a single depth window may be used as a specific processing unit for different depth windows in the following, so as to determine the conversion relationship between the resistivity and the rock color. Specifically, for example, the following may be mentioned: 50 measuring points are used as the vertical length, and 200 measuring points are used as the division standard of the horizontal length which is a unit depth window. Of course, the above-mentioned division criteria of the depth window are only for better explaining the embodiment of the present application, and in the specific implementation, the longitudinal length and the transverse length of the depth window may be determined according to specific situations and construction requirements. The present application is not limited thereto.
S13-3: and respectively establishing a conversion relation between the resistivity and the rock color corresponding to each depth window in the plurality of depth windows according to the color and lithology information of the continuous stratum.
In one embodiment, in implementation, a single depth window may be used as a specific processing unit, and the resistivity and the conversion relation of the rock corresponding to the depth window may be determined in the single depth window. Specifically, the above-mentioned respectively establishing a conversion relationship between the resistivity and the rock color corresponding to each depth window of the plurality of depth windows according to the color and lithology information of the continuous stratum may include establishing a conversion relationship between the resistivity and the rock color corresponding to a current depth window in the following manner:
s1: and determining a first component value corresponding to the current depth window, a second component value corresponding to the current depth window and a third component value corresponding to the current depth window according to the color and lithology information of the continuous stratum, wherein the first component value, the second component value and the third component value correspond to different color component values respectively.
In an embodiment, the first component value corresponding to the current depth window may specifically be a red component value corresponding to the current depth window, the second component value corresponding to the current depth window may specifically be a green component value corresponding to the current depth window, and the third component value corresponding to the current depth window may specifically be a blue component value corresponding to the current depth window. Therefore, the three data, namely the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window, can be subsequently combined to accurately represent the color corresponding to the specific resistivity.
S2: and dividing the current depth window into a plurality of resistivity intervals according to the resistivity.
In this embodiment, in order to more precisely characterize the colors corresponding to the different resistivity ranges in the depth window, the current depth window may be further subdivided into a plurality of resistivity intervals corresponding to the different resistivity ranges according to the resistivity ranges. Each resistivity interval in the plurality of resistivity intervals may include a plurality of pixels, and the plurality of pixels in one resistivity interval correspond to the same resistivity range.
S3: and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
In one embodiment, the rock color corresponding to each resistivity interval in the current depth window may be determined; and then the rock color corresponding to each resistivity interval in the current depth window can be used as the conversion relation between the resistivity corresponding to the current depth window and the rock color. In specific implementation, the rock color corresponding to the current resistivity interval in the resistivity intervals can be determined according to the following formula:
Figure BDA0001512214190000091
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrFor the red component value, R, corresponding to the current depth windowgFor the green component value, R, corresponding to the current depth windowbThe blue component value corresponding to the current depth window.
In this way, the resistivity intervals corresponding to different resistivity ranges in the current depth window can be taken as basic processing units, and the colors corresponding to the resistivity intervals in the current depth window can be respectively determined; and then, the color corresponding to each resistivity interval in the current depth window is used as the conversion relation between the resistivity corresponding to the current depth window and the rock color.
S14: and converting the resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to the conversion relation between the resistivity and the rock color.
In one embodiment, in order to integrate the advantages of the electrical imaging logging image resistivity distribution map, the electrical imaging logging image resistivity distribution map is converted into an optical logging image capable of visually representing the color of reservoir rock, and in specific implementation, according to the rock color corresponding to the resistivity interval, pixel points in the resistivity interval in the electrical imaging logging image resistivity distribution map can be converted into the corresponding rock color respectively. Therefore, the optical logging image based on the resistivity distribution map of the electrical imaging logging image can be obtained, the image has the advantages of electrical imaging logging information, has better continuity and can more accurately reflect the integral geological characteristics of the stratum; meanwhile, the advantages of optical imaging images are considered, certain optical characteristics are achieved, the conditions such as the real color of stratum rocks can be reflected visually and accurately, and further specific research on a target area is facilitated.
Compared with the prior art, the method has the advantages that the conversion relation between the resistivity and the rock color is established firstly by integrating the characteristics of the electrical imaging logging data and the logging data; and then, converting the resistivity distribution map of the electrical imaging logging image into an optical logging image by utilizing the conversion relation so as to replace the optical imaging image. Therefore, the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor in the existing method are solved, and the technical effect that the optical logging image capable of well representing the integral geological characteristics of the reservoir can be obtained at low cost and high efficiency is achieved.
In one embodiment, after converting the electrical imaging log resistivity profile into a corresponding optical log, the method may further include: and according to the optical logging image, performing reservoir evaluation on a target area. Therefore, the optical logging image can replace the optical imaging image with higher cost, the respective advantages of the electrical imaging logging information and the optical imaging image can be simultaneously integrated, and the corresponding reservoir evaluation can be more accurately carried out on the target area, so that the oil and gas distribution condition of the target area can be more accurately analyzed and researched in the following process.
In an embodiment, in order to determine the color corresponding to the resistivity interval more accurately, when the rock color corresponding to each of the resistivity intervals is determined according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window, the following may be further included:
s1: acquiring coring data of a target area;
s2: determining the color of the slot hole filler and the color information of the massive breccia according to the coring data;
s3: and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the color of the slot hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
Therefore, the color of the slot hole filler and the color information of the massive cobbles determined according to the coring data can be used as reference data, so that the rock color corresponding to each resistivity interval can be determined more accurately in the following process.
In this embodiment, the first component value corresponding to the current depth window may specifically be a red component value corresponding to the current depth window, the second component value corresponding to the current depth window may specifically be a green component value corresponding to the current depth window, and the third component value corresponding to the current depth window may specifically be a blue component value corresponding to the current depth window. Therefore, the three data, namely the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window, can be subsequently combined to accurately represent the color corresponding to the specific resistivity.
In an embodiment, in the above determining, according to the color of the hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window, the rock color corresponding to each resistivity interval in the multiple resistivity intervals is determined, in a specific implementation, the rock color corresponding to the current resistivity interval in the multiple resistivity intervals may be determined according to the following formula:
Figure BDA0001512214190000111
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrIs a first component value, R, corresponding to the current depth windowgA second component value, R, corresponding to the current depth windowbA third component value corresponding to the current depth window, cr1A fourth component value corresponding to the current depth window, cg1A fifth component corresponding to the current depth windowValue, cb1A sixth component value, c, corresponding to the current depth windowr2A seventh component value, c, corresponding to the current depth windowg2For an eighth component value, c, corresponding to the current depth windowb2A ninth component value corresponding to the current depth window, wherein cr1、cg1、cb1、cr2、cg2、cb2Is determined according to the color of the slot hole filling material and the color information of the massive breccia.
In the present embodiment, it should be added that c isr1(fourth component value corresponding to current depth window), cg1(fifth component value corresponding to current depth window), cb1The sixth component value (corresponding to the current depth window) may be specifically determined according to the color of the hole filler with relatively low resistivity and the color of the massive boulder in the current depth window; c abover2(seventh component value corresponding to current depth window), cg2(eighth component value corresponding to current depth window), cb2The (ninth component value corresponding to the current depth window) may specifically be determined according to the color of the hole filler with relatively large resistivity and the color of the massive boulder in the current depth window. Therefore, the rock color corresponding to each resistivity interval can be determined more accurately through more component values.
In one embodiment, in order to further accurately establish the conversion relation between the resistivity and the rock color, while acquiring the electrical imaging logging data and the logging data of the target area, reference data such as conventional logging data and core scanning data of the target area can be acquired. Therefore, the information such as the core color, the structure, the fracture-cavity filling materials and the like of the stratum of the target area can be preliminarily analyzed according to the reference data such as the conventional logging data, the core scanning data and the like so as to determine the basic storage space type of the stratum of the target area, and the basic storage space type is used as a reference to guide the establishment of the subsequent resistivity and rock color conversion relation.
In this embodiment, the conventional logging data may specifically refer to other types of logging data different from the electrical imaging logging data and the logging data. The specific data type of the conventional logging data is not limited in the application.
From the above description, it can be seen that the rock core constrained electrical imaging logging image processing method provided by the embodiment of the present application establishes a conversion relationship between resistivity and rock color first by integrating the characteristics of electrical imaging logging data and logging data; then, the conversion relation is utilized to convert the resistivity distribution map of the electrical imaging logging image into an optical logging image to replace the optical imaging image, so that the technical problems that the cost for determining the optical imaging image is high and the continuity of the determined optical imaging image is poor in the existing method are solved, and the technical effect that the optical logging image capable of well representing the integral geological characteristics of the reservoir can be obtained at low cost and high efficiency is achieved; and determining the color of the slot hole filler and the color information of the blocky breccia by combining the coring data, and determining the rock color corresponding to the more accurate resistivity interval according to the color of the slot hole filler and the color information of the blocky breccia, so that the accuracy of the determined optical logging image is improved.
Based on the same inventive concept, the embodiment of the invention also provides a core-constrained electrical imaging logging image processing device, which is described in the following embodiment. Because the principle of solving the problems of the device is similar to the core-constrained electrical imaging logging image processing method, the implementation of the device can refer to the implementation of the core-constrained electrical imaging logging image processing method, and repeated parts are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated. Referring to fig. 2, a block diagram of an apparatus for processing a core-constrained electrographic log image according to an embodiment of the present disclosure may include: the system comprises an acquisition module 21, a determination module 22, a creation module 23, and a conversion module 24, and the structure is described in detail below.
The obtaining module 21 may be specifically configured to obtain electrical imaging logging data and logging data of the target area.
The determining module 22 may be specifically configured to determine an electrical imaging log resistivity profile of the target region according to the electrical imaging log data.
The establishing module 23 may be specifically configured to establish a conversion relationship between the resistivity and the rock color according to the electrical imaging logging image resistivity distribution map and the logging information.
The conversion module 24 may be specifically configured to convert the electrical resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to a conversion relationship between the electrical resistivity and the rock color.
In one embodiment, in order to determine an electrical imaging log resistivity profile of the target region based on the electrical imaging log data, the determining module 22 may be implemented according to the following procedures: and carrying out shallow resistivity scale and whole borehole image generation processing on the electrical imaging logging information to obtain the electrical imaging logging image resistivity distribution diagram.
In one embodiment, in order to establish a conversion relationship between resistivity and rock color according to the electrical imaging log resistivity profile and the logging data, the establishing module 23 may specifically include the following structural units:
the determining unit is specifically used for determining the color and lithology information of the continuous stratum in the target area according to the logging information;
the dividing unit is specifically used for dividing the electrical imaging logging image resistivity distribution map into a plurality of depth windows according to the depth;
and the establishing unit is specifically configured to respectively establish a conversion relationship between the resistivity and the rock color corresponding to each depth window in the multiple depth windows according to the color and lithology information of the continuous stratum.
In an embodiment, in order to respectively establish a conversion relationship between the resistivity and the rock color corresponding to each of the multiple depth windows according to the color and lithology information of the continuous formation, a specific implementation of the establishing unit may include the following structural sub-units:
the first determining subunit is specifically configured to determine, according to the color and lithology information of the continuous formation, a first component value corresponding to the current depth window, a second component value corresponding to the current depth window, and a third component value corresponding to the current depth window;
the dividing subunit is specifically configured to divide the current depth window into a plurality of resistivity intervals according to the resistivity;
the second determining subunit may be specifically configured to determine, according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window, a rock color corresponding to each of the resistivity intervals.
In one embodiment, the first component value corresponding to the current depth window is a red component value corresponding to the current depth window, the second component value corresponding to the current depth window is a green component value corresponding to the current depth window, and the third component value corresponding to the current depth window is a blue component value corresponding to the current depth window.
In one embodiment, in order to be able to determine the rock color corresponding to each of the resistivity intervals, the second determining subunit may determine the rock color corresponding to the current resistivity interval in the resistivity intervals according to the following formula:
Figure BDA0001512214190000141
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrFor the red component value, R, corresponding to the current depth windowgFor the green component value, R, corresponding to the current depth windowbThe blue component value corresponding to the current depth window.
In one embodiment, in order to convert the electrical imaging log resistivity profile into a corresponding optical log image according to the relationship between the resistivity and the rock color, the conversion module 24 may be implemented according to the following procedures: and respectively converting pixel points in the resistivity interval in the electrical imaging logging image resistivity distribution diagram into corresponding rock colors according to the rock colors corresponding to the resistivity interval.
In one embodiment, the apparatus may further include an evaluation module, wherein the evaluation module, when implemented, may perform reservoir evaluation on a target region according to the optical log image.
In one embodiment, in order to determine the rock color corresponding to each resistivity interval more accurately, the second determining subunit may further specifically perform the following steps: acquiring coring data of a target area; determining the color of the slot hole filler and the color information of the massive breccia according to the coring data; and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the color of the slot hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
In one embodiment, in order to determine the rock color corresponding to the resistivity interval more accurately, in a specific implementation, the second determining subunit may determine the rock color corresponding to the current resistivity interval in the multiple resistivity intervals according to the following formula:
Figure BDA0001512214190000142
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrIs a first component value, R, corresponding to the current depth windowgA second component value, R, corresponding to the current depth windowbA third component value corresponding to the current depth window, cr1A fourth component value corresponding to the current depth window, cg1A fifth value of the component corresponding to the current depth window, cb1A sixth component value, c, corresponding to the current depth windowr2A seventh component value, c, corresponding to the current depth windowg2For an eighth component value, c, corresponding to the current depth windowb2A ninth component value corresponding to the current depth window, wherein cr1、cg1、cb1、cr2、cg2、cb2Is determined according to the color of the slot hole filling material and the color information of the massive breccia.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
It should be noted that, the systems, devices, modules or units described in the above embodiments may be implemented by a computer chip or an entity, or implemented by a product with certain functions. For convenience of description, in the present specification, the above devices are described as being divided into various units by functions, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
Moreover, in the subject specification, adjectives such as first and second may only be used to distinguish one element or action from another element or action without necessarily requiring or implying any actual such relationship or order. References to an element or component or step (etc.) should not be construed as limited to only one of the element, component, or step, but rather to one or more of the element, component, or step, etc., where the context permits.
From the above description, it can be seen that the core-constrained electrical imaging logging image processing apparatus provided by the embodiment of the present application integrates the characteristics of electrical imaging logging information and logging information, and establishes a conversion relationship between resistivity and rock color through the establishing module; the conversion relation is utilized by the conversion module to convert the resistivity distribution map of the electrical imaging logging image into an optical logging image to replace the optical imaging image, so that the technical problems of high cost for determining the optical imaging image and poor continuity of the determined optical imaging image in the existing method are solved, and the technical effect of obtaining the optical logging image which can well represent the integral geological characteristics of the reservoir at low cost and high efficiency is achieved; and determining the color of the slot hole filler and the color information of the blocky breccia by combining the establishing module with the coring data, and determining the rock color corresponding to the more accurate resistivity interval according to the color of the slot hole filler and the color information of the blocky breccia, thereby improving the accuracy of the determined optical logging image.
In a specific implementation scenario, the method and the device for processing an electrical imaging log image of core constraint provided by the present application are specifically applied to produce an optical log image capable of representing a reservoir color according to electrical imaging log data and log data of a certain work area, so as to replace a conventional optical imaging image. For a specific implementation process, the following contents can be referred to.
S1: collecting and collating well logging data of a work area, which may specifically include: logging data, electrical imaging logging data, core data, and the like.
In the embodiment, during specific implementation, conventional logging information, electrical imaging logging information, logging rock debris color and lithology description information of the interval to be researched in the work area can be collected, sorted and sorted; meanwhile, the related data such as the core color, the structure, the fracture-cavity filler description data, the core scanning image and the like are analyzed, so that the integral knowledge of the formation lithology and the main reservoir space type is established, and the subsequent establishment of the specific resistivity and rock color conversion relationship is guided.
S2: and obtaining information such as the color of the stratum rock, the lithology description and the like (namely the color and lithology information of the continuous stratum in the target area) according to the description of the logging rock debris (namely logging information).
In the embodiment, the logging cuttings may be rock cuttings drilled at the bottom of the well by the drill bit during the drilling process, and then the logging cuttings are returned to the ground after circulating the drilling fluid, so that the logging cuttings are visual data capable of better reflecting formation lithology and oil-gas-containing information. The cuttings description data of the logging record may specifically include: formation rock color, lithology, oil and gas containing property and the like. Taking a certain well in a certain work area as an example, the logging depth is 5731.0-5735.0m, and the rock debris is displayed as gray fluorescent sand debris limestone, so that the lithology of the stratum at the section can be judged as sand debris limestone, and the color is gray. And then the rock debris information can be classified and counted to obtain information such as the color, lithology and the like of the continuous stratum rock in the whole well section.
S3: and processing the imaging logging data to obtain well wall resistivity image information covered by the whole borehole (namely determining an electrical imaging logging image resistivity distribution diagram of the target area according to the electrical imaging logging data).
In the embodiment, during specific processing, preprocessing such as acceleration correction processing, equalization processing and the like can be performed on the electrical imaging logging data to obtain original electrical imaging logging image information; and then, the information of the original electrical imaging logging image is scaled window by adopting shallow resistivity so as to obtain a well wall resistivity change image of the whole well section.
In the embodiment, it is considered that the conventional method for performing electrical imaging logging mostly adopts a polar plate-borehole wall-to-wall mode to perform corresponding measurement, and the whole borehole cannot be covered completely (the coverage rate can only reach 60% -80%). Therefore, a whole borehole image processing method is required to perform probability interpolation processing on uncovered parts in the resistivity image after calibration so as to obtain resistivity image information covered by the whole borehole. Blank parts in the original image can be effectively made up through the processing of the full-borehole image, and the display effect of the later-stage simulated rock core scanning image can be improved.
In this embodiment, taking a certain well in a certain working area as an example, specific situations of a resistivity image obtained after scaling an original imaging image of a well section with a depth of 5730.0-5734.0m and a borehole wall resistivity image obtained after full borehole processing can be referred to as a schematic diagram of a result obtained after processing the resistivity scale and the full borehole image by applying the electrical imaging logging image processing method and device for rock core constraint provided by the embodiment of the present application in one scene example shown in fig. 3.
S4: and establishing a resistivity-rock color conversion relation (namely the conversion relation between the resistivity and the rock color).
In this embodiment, in a specific implementation, statistics may be performed on the borehole wall resistivity images covered by the whole borehole one by one according to the depth windows obtained in S3, so as to determine the resistivity distribution map of each depth window image. The resistivity distribution map of each depth window image is divided into N intervals (namely resistivity distribution intervals) according to a given classification standard (namely according to the resistivity range). Specifically, all the pixel points in the depth window can be divided into N equal-pixel-number intervals according to the order of the resistivity values of the pixel points in the depth window, so as to obtain N intervals. For a typical resistivity distribution and interval division, a resistivity distribution graph and a resistivity interval division schematic diagram obtained by applying the method and the device for processing the electrical imaging log image of the core constraint provided by the embodiment of the present application in one scenario example shown in fig. 4 can be referred to.
In this embodiment, in specific implementation, the rock color corresponding to the current depth window may be determined according to characteristic information such as formation rock color information obtained from logging rock debris data. Specifically, the RGB color values of each resistivity interval in the depth window (i.e., the rock color corresponding to the resistivity interval) may be determined according to the following formula:
Figure BDA0001512214190000171
in the above formula, RankRGB value representing the color corresponding to the k resistivity interval; rb、Rg、RrRespectively representing the RGB component values of the rock color corresponding to the current depth window.
In this embodiment, in order to determine the RGB color values of each resistivity interval more accurately, the coring data of the well section to be analyzed may be acquired, and the core observation may be performed on the coring data, so as to obtain the color information of the fracture-cavity filler and the massive gravel (i.e., the color of the fracture-cavity filler and the color information of the massive gravel), synthesize the core data observation to obtain information, and determine the RGB color values of each resistivity interval according to the following formula:
Figure BDA0001512214190000172
in the above formula, Cb1、Cg1、Cr1Respectively representing RGB components of low-resistance type slot hole fillers or block-shaped breccia colors observed on the corresponding rock core of the current depth window; cb2、Cg2、Cr2Respectively representing RGB components of high-resistance type slot hole filling materials or massive breccia colors observed on the corresponding rock core of the current depth window.
In this embodiment, for example, when a well in a work area is used and cuttings with a logging depth of 5731.0-5735.0m are displayed as gray fluorescent sand-debris limestone, the color of the rock in the depth window is set to gray, and the RGB components are Rb＝142、Rg＝149、Rr142. Furthermore, a clear white calcite filling gap is visible from the core, of the type with high resistance filling, Cb2＝255、Cg2＝255、Cr22255; the low-resistance crack is not filled, and the color component is set as Cb1＝0、Cg1＝0、C r10. The RGB color values of each resistivity interval can be calculated and determined by substituting the formula.
S5: the borehole wall resistivity image of the full borehole coverage is converted to a rock color image (i.e., an optical log image).
In this embodiment, in specific implementation, for each pixel point in the depth window, the resistivity interval number where the pixel point is located is determined according to the resistivity-rock color relationship determined in S4, and then the resistivity value may be converted into a corresponding color value by using formula (1) or (2). Therefore, the conversion from the borehole wall resistivity image to the rock color image can be carried out window by window on the whole borehole section. Specifically, the simulated core scanned image (i.e., the optical logging image) obtained by performing the above processing on a well section with a well depth of 5732.0-5733.5m in a certain work area of the tower is compared with the actual core scanned image (i.e., the optical imaging image), so that the obtained simulated core scanned image really has a good characterization effect, and can actually replace the core scanned image with high cost to a certain extent. Specifically, a schematic diagram of a comparison result of a core scan image simulated by an optical electrical logging image obtained by applying the method and the device for processing the electrical imaging logging image of the core constraint in one scenario example shown in fig. 5 can be provided.
Through the scene example, the rock core constraint electrical imaging logging image processing device provided by the embodiment of the application is verified, and the conversion relation between the resistivity and the rock color is established firstly through the characteristics of the electrical imaging logging data and the logging data; and then, the conversion relation is utilized to convert the resistivity distribution map of the electrical imaging logging image into an optical logging image to replace the optical imaging image, so that the technical problems of high cost for determining the optical imaging image and poor continuity of the determined optical imaging image in the existing method are really solved, and the technical effect of obtaining the optical logging image which can well represent the integral geological characteristics of the reservoir with low cost and high efficiency is really achieved.
Although various specific embodiments are mentioned in the present application, the present application is not limited to the cases described in the industry standards or examples, and the like, and some industry standards or the embodiments slightly modified based on the implementation described in the custom manner or examples can also achieve the same, equivalent or similar implementation effects as those of the above embodiments or the implementation effects expected after the modifications. Embodiments employing such modified or transformed data acquisition, processing, output, determination, etc., may still fall within the scope of alternative embodiments of the present application.
Although the present application provides method steps as described in an embodiment or flowchart, more or fewer steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an apparatus or client product in practice executes, it may execute sequentially or in parallel (e.g., in a parallel processor or multithreaded processing environment, or even in a distributed data processing environment) according to the embodiments or methods shown in the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded.
The devices or modules and the like explained in the above embodiments may be specifically implemented by a computer chip or an entity, or implemented by a product with certain functions. For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, in implementing the present application, the functions of each module may be implemented in one or more pieces of software and/or hardware, or a module that implements the same function may be implemented by a combination of a plurality of sub-modules, and the like. The above-described apparatus embodiments are merely illustrative, and for example, the division of the modules is merely a logical division, and other divisions may be realized in practice, for example, a plurality of modules or components may be combined or integrated into another system, or some features may be omitted, or not executed.
Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may therefore be considered as a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which may be stored in a storage medium, such as a ROM/RAM, a magnetic disk, an optical disk, or the like, and includes several instructions for enabling a computer device (which may be a personal computer, a mobile terminal, a server, or a network device) to execute the method according to the embodiments or some parts of the embodiments of the present application.
The embodiments in the present specification are described in a progressive manner, and the same or similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet-type devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
While the present application has been described by way of examples, those of ordinary skill in the art will appreciate that there are numerous variations and permutations of the present application that do not depart from the spirit of the present application and that the appended embodiments are intended to include such variations and permutations without departing from the present application.

Claims (10)

1. A method for processing a core-constrained electrical imaging log image is characterized by comprising the following steps:
acquiring electrical imaging logging information and logging information of a target area;
determining an electrical imaging logging image resistivity distribution diagram of a target area according to the electrical imaging logging information;
establishing a conversion relation between the resistivity and the rock color according to the electrical imaging logging image resistivity distribution map and the logging information;
converting the resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to the conversion relation between the resistivity and the rock color;
wherein, according to the resistivity distribution map of the electric imaging logging image and the logging information, establishing a conversion relation between the resistivity and the rock color, and the method comprises the following steps: determining the color and lithology information of continuous stratums in the target area according to the logging information; dividing the electrical imaging logging image resistivity distribution map into a plurality of depth windows according to the depth; and respectively establishing a conversion relation between the resistivity and the rock color corresponding to each depth window in the plurality of depth windows according to the color and lithology information of the continuous stratum.
2. The method of claim 1, wherein determining an electrical imaging log resistivity profile of the target region from the electrical imaging log data comprises:
and carrying out shallow resistivity scale and whole borehole image generation processing on the electrical imaging logging information to obtain the electrical imaging logging image resistivity distribution diagram.
3. The method of claim 1, wherein the step of respectively establishing a conversion relationship between the resistivity and the rock color corresponding to each of the plurality of depth windows according to the color and lithology information of the continuous stratum comprises:
establishing a conversion relation between the resistivity corresponding to the current depth window and the rock color according to the following mode:
determining a first component value corresponding to a current depth window, a second component value corresponding to the current depth window and a third component value corresponding to the current depth window according to the color and lithology information of the continuous stratum, wherein the first component value, the second component value and the third component value correspond to different color component values respectively;
dividing the current depth window into a plurality of resistivity intervals according to the resistivity;
and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
4. The method according to claim 3, wherein the first component value corresponding to the current depth window is a red component value corresponding to the current depth window, the second component value corresponding to the current depth window is a green component value corresponding to the current depth window, and the third component value corresponding to the current depth window is a blue component value corresponding to the current depth window.
5. The method of claim 4, wherein separately determining a rock color for each of the plurality of resistivity intervals comprises:
determining the rock color corresponding to the current resistivity interval in the resistivity intervals according to the following formula:
Figure FDA0002438361450000021
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrFor the red component value, R, corresponding to the current depth windowgIs the current depthGreen component value, R, corresponding to the window of degreesbThe blue component value corresponding to the current depth window.
6. The method of claim 3, wherein converting the electrical imaging log resistivity profile into a corresponding optical log based on the resistivity and rock color conversion comprises:
and respectively converting pixel points in the resistivity interval in the electrical imaging logging image resistivity distribution diagram into corresponding rock colors according to the rock colors corresponding to the resistivity interval.
7. The method of claim 1, wherein after converting the electrical imaging log resistivity profile to a corresponding optical log image, the method further comprises:
and according to the optical logging image, performing reservoir evaluation on a target area.
8. The method of claim 3, wherein determining the rock color corresponding to each of the plurality of resistivity intervals according to the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window comprises:
acquiring coring data of a target area;
determining the color of the slot hole filler and the color information of the massive breccia according to the coring data;
and respectively determining the rock color corresponding to each resistivity interval in the multiple resistivity intervals according to the color of the slot hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window and the third component value corresponding to the current depth window.
9. The method according to claim 8, wherein determining a rock color corresponding to each of the plurality of resistivity intervals according to the color of the hole filler and the color information of the massive cobbles, the first component value corresponding to the current depth window, the second component value corresponding to the current depth window, and the third component value corresponding to the current depth window comprises:
determining the rock color corresponding to the current resistivity interval in the resistivity intervals according to the following formula:
Figure FDA0002438361450000031
in the above formula, RankThe rock color value corresponding to the current resistivity interval, k is the number of the current resistivity interval, N is the total number of the resistivity intervals in the current depth window where the current resistivity interval is located, and R is the total number of the resistivity intervals in the current depth window where the current resistivity interval is locatedrIs a first component value, R, corresponding to the current depth windowgA second component value, R, corresponding to the current depth windowbA third component value corresponding to the current depth window, cr1A fourth component value corresponding to the current depth window, cg1A fifth value of the component corresponding to the current depth window, cb1A sixth component value, c, corresponding to the current depth windowr2A seventh component value, c, corresponding to the current depth windowg2For an eighth component value, c, corresponding to the current depth windowb2A ninth component value corresponding to the current depth window, wherein cr1、cg1、cb1、cr2、cg2、cb2Is determined according to the color of the slot hole filling material and the color information of the massive breccia.
10. A core-constrained electrical imaging logging image processing apparatus, comprising:
the acquisition module is used for acquiring electrical imaging logging data and logging data of a target area;
the determining module is used for determining an electrical imaging logging image resistivity distribution map of the target area according to the electrical imaging logging information;
the establishing module is used for establishing a conversion relation between the resistivity and the rock color according to the electrical imaging logging image resistivity distribution map and the logging information;
the conversion module is used for converting the resistivity distribution map of the electrical imaging logging image into a corresponding optical logging image according to the conversion relation between the resistivity and the rock color;
wherein the establishing module comprises: the determining unit is used for determining the color and lithology information of the continuous stratum in the target area according to the logging information; the dividing unit is used for dividing the electrical imaging logging image resistivity distribution map into a plurality of depth windows according to the depth; and the establishing unit is used for respectively establishing the conversion relation between the resistivity and the rock color corresponding to each depth window in the plurality of depth windows according to the color and lithology information of the continuous stratum.
CN201711362904.9A 2017-12-18 2017-12-18 Core-constrained electrical imaging logging image processing method and device Active CN108318534B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711362904.9A CN108318534B (en) 2017-12-18 2017-12-18 Core-constrained electrical imaging logging image processing method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711362904.9A CN108318534B (en) 2017-12-18 2017-12-18 Core-constrained electrical imaging logging image processing method and device

Publications (2)

Publication Number Publication Date
CN108318534A CN108318534A (en) 2018-07-24
CN108318534B true CN108318534B (en) 2020-07-10

Family

ID=62892340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711362904.9A Active CN108318534B (en) 2017-12-18 2017-12-18 Core-constrained electrical imaging logging image processing method and device

Country Status (1)

Country Link
CN (1) CN108318534B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109577947A (en) * 2018-11-23 2019-04-05 长江大学 A kind of well logging and rock core information linkage comprehensive analysis method and system
CN112067663B (en) * 2020-08-05 2024-01-26 山东天岳先进科技股份有限公司 Method and device for detecting resistivity of high-purity silicon carbide crystal
CN111983693A (en) * 2020-08-18 2020-11-24 东华理工大学 Geological disaster multi-mode monitoring data fusion imaging method based on resistivity method
CN116184527B (en) * 2023-04-25 2023-07-11 西南石油大学 Magma rock lithology discrimination method and device based on element logging and imaging logging

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3213105B2 (en) * 1993-02-23 2001-10-02 日本電信電話株式会社 Video composition automatic decision processing method
US9047689B2 (en) * 2009-06-04 2015-06-02 Schlumberger Technology Corporation Method and apparatus for visualizing multi-dimensional well logging data with shapelets
CN101799559B (en) * 2010-02-09 2011-12-07 康志勇 Formation data processing method based on resistivity of sandstone reservoir oil layer
CN102200008B (en) * 2010-03-26 2013-03-13 中国石油天然气股份有限公司 Reservoir effectiveness identification method based on electrical imaging logging
CN102373923B (en) * 2010-08-20 2013-03-13 中国石油天然气股份有限公司 Reservoir stratum identification method
CN103529475B (en) * 2013-04-19 2016-08-03 中国石油大学(华东) A kind of method identified and explain carbonate paleokarst reservoir three dimensional structure
CN103592690B (en) * 2013-10-24 2016-09-14 长江大学 Method based on electric imaging logging porosity spectrum information automatic Reservoir Recognition crack
CN104808248A (en) * 2015-04-22 2015-07-29 长江大学 Attitude of stratum detecting method based on electric imaging logging full hole image
CN105891905A (en) * 2016-04-13 2016-08-24 成都创源油气技术开发有限公司 Shale lithofacies well logging fast identification method
CN107315203A (en) * 2017-08-04 2017-11-03 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 A kind of tight sandstone reservoir dessert identification and integrated evaluating method

Also Published As

Publication number Publication date
CN108318534A (en) 2018-07-24

Similar Documents

Publication Publication Date Title
CN108318534B (en) Core-constrained electrical imaging logging image processing method and device
CN106526693B (en) Crack identification method and apparatus
Casini et al. Fracture characterization and modeling from virtual outcrops
CN107402176B (en) method and device for determining porosity of crack
Wilson et al. From outcrop to flow simulation: Constructing discrete fracture models from a LIDAR survey
EP3217192A1 (en) Method for looking for oil and gas reservoir on the basis of trap-3d software
CN106370817A (en) Quantitative identification of fractures and caverns based on core analysis and electrical imaging logging
CN108363110A (en) Imaging logging calculates shale reservoir mineral content and the spectral analysis method of brittleness index
Harraden et al. Automated core logging technology for geotechnical assessment: A study on core from the Cadia East porphyry deposit
EA016499B1 (en) Automated borehole image interpretation
Minisini et al. Geological data extraction from lidar 3-D photorealistic models: A case study in an organic-rich mudstone, Eagle Ford Formation, Texas
Rivard et al. Inferring total organic carbon and major element geochemical and mineralogical characteristics of shale core from hyperspectral imagery
BR112015021666B1 (en) computer-implemented method for updating an estimate for a material property of a volume and computer system for updating an estimate for a material property of a volume
Harvey et al. Measuring the grain‐size distributions of mass movement deposits
Hansberry et al. How wide is a fault damage zone? Using network topology to examine how fault-damage zones overprint regional fracture networks
Hood et al. Assessing continuous resources–building the bridge between static and dynamic analyses
Holland et al. Evolution of fractures in a highly dynamic thermal, hydraulic, and mechanical system-(II) Remote sensing fracture analysis, Jabal Shams, Oman Mountains
Harraden et al. Proposed methodology for utilising automated core logging technology to extract geotechnical index parameters
CN109738955B (en) Metamorphic rock lithology comprehensive judgment method based on component-structure classification
CN115690292B (en) Working face stratum modeling rendering optimization method and system
CN113592823B (en) Rock stratum porosity distribution calculation method
Koch et al. Does fluvial channel-belt clustering predict net sand to gross rock volume? Architectural metrics and point-pattern analysis of a digital outcrop model
CN109598049B (en) Method for drilling rock fracture development degree and regional rock fracture development rule
Aliverti et al. Data analysis, processing and 3D fracture network simulation at wellbore scale for fractured reservoir description
Farmakis et al. Assessment of the GSI along rock slopes based on LiDAR and photogrammetry point clouds

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