CN110966000B - Method and system for evaluating layering property of glutenite reservoir - Google Patents
Method and system for evaluating layering property of glutenite reservoir Download PDFInfo
- Publication number
- CN110966000B CN110966000B CN202010000508.7A CN202010000508A CN110966000B CN 110966000 B CN110966000 B CN 110966000B CN 202010000508 A CN202010000508 A CN 202010000508A CN 110966000 B CN110966000 B CN 110966000B
- Authority
- CN
- China
- Prior art keywords
- index
- reservoir
- conductivity
- row
- column
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000003384 imaging method Methods 0.000 claims abstract description 29
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 37
- 239000011159 matrix material Substances 0.000 claims description 32
- 230000014509 gene expression Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 6
- 238000012512 characterization method Methods 0.000 abstract description 10
- 230000003068 static effect Effects 0.000 description 14
- 238000011160 research Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000013517 stratification Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Abstract
The invention discloses a method and a system for representing a glutenite reservoir layer index, which comprise the following steps: collecting a glutenite reservoir image by using an electrical imaging logging technology; the method comprises the steps that a current image is subjected to gridding division according to preset row intervals and preset column intervals, a current reservoir stratum transverse continuity index is obtained according to the standard deviation of the row conductivity value of each row, and a corresponding longitudinal continuity index is obtained according to the standard deviation of the column conductivity value of each column; and calculating the stratiform index of the current reservoir according to the transverse continuity index and the longitudinal continuity index so as to quantitatively represent the strength of the stratiform of the current reservoir. The invention realizes the quantitative characterization of the reservoir stratum index, determines the stratum index limit of the effective reservoir stratum and provides an effective means for the effective reservoir stratum identification based on the electrical imaging.
Description
Technical Field
The invention relates to the field of petroleum logging engineering, in particular to a method and a system for evaluating the stratifying property of a glutenite reservoir.
Background
In recent years, with the continuous expansion of the field of oil-gas exploration and development in China, the conglomerate oil-gas reservoir becomes one of the objects of deep exploration. The electric imaging well logging can provide visual, clear and high-resolution images of the stratum around the shaft, and the electric imaging images contain geological information such as abundant lithology, rock structures and the like in the stratum. In the glutenite reservoir section, the electrical imaging image can clearly reflect the layered structure characteristics of the reservoir.
At present, qualitative analysis of a large number of examples in a research area finds that the layering property of a glutenite reservoir is closely related to the reservoir capacity, and the strength of the layering property of the reservoir is a key index for judging the effectiveness of the reservoir. The prior people do a lot of research on evaluating the glutenite reservoir structure by utilizing an electrical imaging image, but the research on the glutenite stratifying reservoir structure is less, the research is only limited on the qualitative description of the reservoir stratifying property, an effective means for quantitatively characterizing the reservoir stratifying property is lacked, and the quantitative identification standard of the effective reservoir can not be established on the basis of the reservoir stratifying property.
Disclosure of Invention
In order to solve the above technical problem, an embodiment of the present invention provides a method for evaluating stratifying property of a glutenite reservoir, including: acquiring a conglomerate reservoir image by utilizing an electrical imaging logging technology; step two, gridding and dividing the image according to preset row intervals and column intervals, and obtaining a current reservoir stratum transverse continuity index according to the standard variance of the row conductivity value of each row and a corresponding longitudinal continuity index according to the standard variance of the column conductivity value of each column; and thirdly, calculating the stratiform index of the current reservoir according to the transverse continuity index and the longitudinal continuity index so as to quantitatively represent the strength of the stratiform of the current reservoir.
Preferably, the method further comprises: and determining a stratifying index standard for identifying whether the reservoirs in the block to be researched are effective or not according to the stratifying indexes corresponding to different reservoirs in the block to be researched and combining the productivity information of the different reservoirs, so as to judge the effectiveness of the reservoirs in the block to be researched by utilizing the stratifying index standard.
Preferably, in the second step, a conductivity matrix is constructed for the glutenite reservoir image; calculating the standard deviation of the row conductivity of each row in the image, and calculating the average value of the standard deviations of the row conductivity of each row to obtain the transverse continuity index; and calculating the column conductivity standard deviation of each column in the image, and calculating the average value of the column conductivity standard deviations of each column to obtain the longitudinal continuity index.
Preferably, the longitudinal continuity index and the transverse continuity index are divided to obtain the lamellar index.
Preferably, the method further comprises: scanning the images subjected to gridding division processing according to a preset window to construct a conductivity matrix of each scanning window; calculating the transverse continuity index, the longitudinal continuity index and the corresponding layer index corresponding to each scanning window by executing the second step and the third step; and connecting all the layer indexes according to the depth of the reservoir, and generating a layer index curve aiming at the current reservoir.
In another aspect, there is provided a system for evaluating the stratifying property of a glutenite reservoir, comprising: an image acquisition module configured to acquire a conglomerate reservoir image using an electrical imaging logging technique; the continuity index generation module is configured to grid and divide the image according to preset row intervals and column intervals, obtain a current reservoir transverse continuity index according to the standard deviation of the row conductivity value of each row, and obtain a corresponding longitudinal continuity index according to the standard deviation of the column conductivity value of each column; and the stratiform index generating module is configured to calculate the stratiform index of the current reservoir according to the transverse continuity index and the longitudinal continuity index so as to quantitatively represent the strength of the stratiform of the current reservoir.
Preferably, the system further comprises: and the standard generation module is configured to determine a layered index standard for identifying whether the reservoirs in the block to be researched are effective or not according to the layered indexes corresponding to different reservoirs in the block to be researched and by combining the productivity information of the different reservoirs, so as to judge the effectiveness of the reservoirs in the block to be researched by using the layered index standard.
Preferably, the continuity index generating module includes: a conductivity matrix construction unit configured to construct a conductivity matrix with respect to the glutenite reservoir image; a transverse continuity index calculation unit configured to calculate a row conductivity standard deviation of each row in the image and calculate an average value of the row conductivity standard deviations of the rows to obtain the transverse continuity index; a longitudinal continuity index calculation unit configured to calculate a column conductivity standard deviation of each column in the image and calculate an average of the column conductivity standard deviations of the columns, resulting in the longitudinal continuity index.
Preferably, the layered index generating module is further configured to divide the longitudinal continuity index and the transverse continuity index to obtain the layered index.
Preferably, the system further comprises a curve plotting module, wherein the curve plotting module comprises: the window scanning unit is configured to scan the images subjected to gridding division processing according to preset windows, and a conductivity matrix related to each scanning window is constructed; a window layer index calculation unit configured to execute the continuity index generation module and the layer index generation module, and calculate the transverse continuity index, the longitudinal continuity index and the corresponding layer index corresponding to each scanning window; and the stratiform index curve generating unit is configured to connect all the stratiform indexes according to the reservoir depths and generate the stratiform index curve aiming at the current reservoir.
Compared with the prior art, one or more embodiments in the above scheme can have the following advantages or beneficial effects:
the invention provides a method and a system for representing a glutenite reservoir stratiform index. The method and the system are based on an electrical imaging logging technology, and utilize the conductivity image of the reservoir to calculate the transverse continuity and the longitudinal continuity, so as to further obtain a quantitative index (stratiform index) for representing the stratifying strength of the reservoir. In addition, by utilizing the calculation process of constructing the quantitative index of the reservoir stratum index and combining the capacity conditions of different reservoir stratum samples in the block to be researched, the key stratum index standard for evaluating the reservoir stratum effectiveness in the current block to be researched is determined, and the key stratum index standard is used for automatically judging and identifying the reservoir stratum effectiveness in the block to be researched. The invention realizes the quantitative characterization of the stratifying property of the reservoir stratum, determines the stratiform index limit of the effective reservoir stratum, provides an effective means for the effective reservoir stratum identification based on electrical imaging, and has the effective reservoir stratum identification accuracy rate of 84.2 percent.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
fig. 1 is a step chart of a method for evaluating the stratifying property of a glutenite reservoir according to an embodiment of the present application.
Fig. 2 is a specific flowchart of the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application.
Fig. 3 is a specific flowchart of the method for evaluating the stratifying property of a conglomerate reservoir according to the embodiment of the present application, in which a stratifying index curve is generated.
Fig. 4 shows two examples of glutenite reservoir images in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application.
Fig. 5 is a schematic diagram illustrating calculation of a transverse continuity index and a longitudinal continuity index in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application.
Fig. 6 is a schematic diagram illustrating an example of calculation of a stratiform index in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application.
Fig. 7 is a schematic diagram illustrating an example of determining a stratifying index criterion in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application.
Fig. 8 is a block diagram of a system for evaluating the stratifying property of a glutenite reservoir according to an embodiment of the present invention.
Detailed Description
The following detailed description of the embodiments of the present invention will be provided with reference to the drawings and examples, so that how to apply the technical means to solve the technical problems and achieve the technical effects can be fully understood and implemented. It should be noted that, as long as there is no conflict, the embodiments and the features of the embodiments of the present invention may be combined with each other, and the technical solutions formed are within the scope of the present invention.
In recent years, with the continuous expansion of the field of oil-gas exploration and development in China, the conglomerate oil-gas reservoir becomes one of the objects of deep exploration. The electric imaging well logging can provide visual, clear and high-resolution images of the stratum around the shaft, and the electric imaging images contain geological information such as abundant lithology, rock structures and the like in the stratum. In the glutenite reservoir section, the electrical imaging image can clearly reflect the layered structure characteristics of the reservoir.
At present, qualitative analysis of a large number of examples in a research area finds that the layering property of a glutenite reservoir is closely related to the reservoir capacity, and the strength of the layering property of the reservoir is a key index for judging the effectiveness of the reservoir. The prior people do a lot of research on evaluating the glutenite reservoir structure by utilizing an electrical imaging image, but the research on the glutenite stratifying reservoir structure is less, the research is only limited on the qualitative description of the reservoir stratifying property, an effective means for quantitatively characterizing the reservoir stratifying property is lacked, and the quantitative identification standard of the effective reservoir can not be established on the basis of the reservoir stratifying property.
In order to solve the problems in the prior art, the invention provides a method and a system for evaluating the stratifying property of a glutenite reservoir. The method and the system construct the stratiform index of the reservoir and a corresponding stratiform index continuous curve on the basis of the transverse and longitudinal continuity of a quantitative representation reservoir image based on the electrical imaging logging information of the glutenite reservoir, so that the quantitative representation method of the stratiform strength of the reservoir is realized, and the effectiveness of the glutenite reservoir is predicted or evaluated. In addition, the invention can also establish effective reservoir identification standard (effective reservoir identification threshold parameter) based on the layer index in the area to be researched by utilizing the layer index quantitative characterization method. Therefore, the method not only meets the urgent need of layering quantitative characterization of the glutenite reservoir image, but also forms an effective means for reservoir effectiveness judgment.
Fig. 1 is a step chart of a method for evaluating the stratifying property of a glutenite reservoir according to an embodiment of the present application. Fig. 2 is a specific flowchart of the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application. With reference to fig. 1 and 2, step S110 acquires an image of the glutenite reservoir by using an electrical imaging logging technique. In step S110, an electrical imaging logging construction is performed on the glutenite reservoir in the region to be studied, and a high-resolution static image or dynamic image capable of reflecting geological information such as the lithology of the formation around the wellbore, the rock structure, and the like is obtained. Fig. 4 shows two examples of glutenite reservoir images in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application. FIG. 4 shows two examples of glutenite reservoir images, and FIG. 4a shows an image of a glutenite reservoir with a weak stratification property, the left image being a static image and the right image being a dynamic image; fig. 4b shows an image of a glutenite reservoir with a strong stratifying property, the left image being a static image and the right image being a dynamic image. It should be noted that in the glutenite reservoir image acquired by the electrical imaging logging technology, the formation conductivity value of each pixel point in the image can be calculated, and based on this, the subsequent steps S120 and S130 are used to quantitatively evaluate the stratification of the current reservoir.
Then, step S120 is to perform gridding division on the currently obtained glutenite reservoir image according to preset row intervals and column intervals, and obtain the horizontal continuity index of the current reservoir from the standard deviation of the row conductivity value of each row (grid), and obtain the vertical continuity index of the current reservoir from the standard deviation of the column conductivity value of each column (grid). Specifically, step S1201 is to perform meshing division on the currently obtained gravel reservoir static image or dynamic image according to preset row intervals and preset column intervals, so as to divide the current gravel reservoir image into a plurality of meshes, where the transverse width of each mesh coincides with the length of the column intervals and the longitudinal width of each mesh coincides with the length of the row intervals, and then the process proceeds to step S1202. In the embodiment of the present invention, the length of the row interval and the column interval is not specifically limited, and those skilled in the art can set the length according to the actual calculation accuracy of the layer index, and each grid may include a pixel block formed by a plurality of pixel points, or may be a single pixel point.
Step S1202 constructs a conductivity matrix for a current conglomerate reservoir static image or dynamic image based on the image subjected to the gridding partition processing. Wherein the conductivity matrix is represented by the following matrix expression:
in the above matrix expression, C denotes a conductivity matrix, m denotes the number of rows of the conductivity matrix C, and n denotes the number of columns of the conductivity matrix C. Wherein each matrix element in the conductivity matrix represents a grid, and the value of each matrix element is the conductivity value of the grid. If one grid only contains one pixel point, the conductivity value of the grid is the conductivity value of the pixel point; if the grid is a pixel block, the conductivity value of the grid is preferably the average value of the conductivities of the pixel points in the current block.
After the conductivity matrix construction is completed, the process proceeds to step S1203. Step S1203 is to calculate a row conductivity standard deviation of each row of grids in the static image or the dynamic image of the current glutenite reservoir along the horizontal direction of the formation, and calculate an average value of the row conductivity standard deviations of each row, and use the average value to represent a continuity index of the current reservoir in the horizontal direction, that is, a current reservoir horizontal continuity index. Fig. 5 is a schematic diagram illustrating calculation of a transverse continuity index and a longitudinal continuity index in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application. As shown in FIG. 5, the direction indicated by the right arrow is the direction of the lateral continuity characterization.
Specifically, firstly, carrying out average value operation on conductivity values of all grids in each row, and calculating the row conductivity of each row in the image by using a row conductivity calculation formula; then, according to the row conductivity of each row, calculating the standard deviation of the row conductivity of the corresponding row by using a row conductivity variance calculation formula; and finally, obtaining conductivity continuity characteristics of the current reservoir in the transverse direction reflected by the current image, namely a transverse continuity index, by utilizing a transverse continuity index calculation formula according to the row conductivity standard variance value of each row. Wherein the line conductivity calculation formula is expressed by the following expression:
in the formula (1), i represents a row number, j represents a column number, MEAN _ h (i) represents the row conductivity of the ith row, and C (i, j) represents the conductivity value of the ith row and jth column grid. Further, the line conductivity variance calculation formula is expressed by the following expression:
in the formula (2), STDD _ h (i) represents a line conductivity standard deviation of the i-th line. Further, the lateral continuity index calculation formula is expressed by the following expression:
in equation (3), CONT _ H represents a lateral continuity index of the current image. The lateral continuity index of the image represents the continuity of the reservoir conductivity in the horizontal lateral direction, and the larger the lateral continuity index value, the poorer the continuity of the reservoir in the lateral direction, that is, the weaker the reservoir continuity in the horizontal direction along the formation; conversely, a smaller value for lateral continuity index indicates better reservoir continuity laterally, that is, greater reservoir continuity along the horizontal direction of the formation.
After the calculation of the lateral continuity index is completed, the process proceeds to step S1204. Step S1204 is along the horizontal vertical direction of the stratum, calculate the standard deviation of the column conductivity of each grid in the static image or dynamic image of the current glutenite reservoir, and calculate the mean value of the standard deviation of the column conductivity of each column, use this mean value to represent the continuity of the current reservoir in the longitudinal direction, namely the current reservoir longitudinal continuity index. As shown in FIG. 5, the direction indicated by the upwardly pointing arrow is the direction of the longitudinal continuity characterization.
Specifically, firstly, the conductivity values of all grids in each column are subjected to average value operation, and the column conductivity of each column in the image is calculated by using a column conductivity calculation formula; then, according to the column conductivity of each column, calculating the standard deviation of the column conductivity of the corresponding column by using a column conductivity variance calculation formula; and finally, according to the column conductivity standard variance value of each column, utilizing a longitudinal continuity index calculation formula to obtain a conductivity continuity characteristic, namely a longitudinal continuity index, of the current reservoir in the longitudinal direction reflected by the current image. Wherein the column conductivity calculation formula is represented by the following expression:
in the formula (4), MEAN _ v (i) represents a column conductivity value in the j-th column. Further, the column conductivity variance calculation formula is expressed by the following expression:
in the formula (5), STDD _ v (j) represents the column conductivity standard deviation of the j-th column. Further, the longitudinal continuity index calculation formula is expressed by the following expression:
in equation (6), CONT _ V represents a longitudinal continuity index of the current image. The longitudinal continuity index of the image represents the continuity of the conductivity of the reservoir in the vertical and longitudinal directions, and the larger the longitudinal continuity index value is, the poorer the continuity of the reservoir in the longitudinal direction is, namely, the more the stratifying property of the reservoir in the longitudinal direction of the stratum is possibly; conversely, a smaller value of the longitudinal continuity index indicates a more continuous reservoir continuity in the longitudinal direction, that is, the reservoir stratifying may be weaker in the longitudinal direction of the formation.
After the calculation of the longitudinal continuity index is completed, the above step S120 is ended, and the process proceeds to step S130. Step S130 calculates the stratiform index of the current conglomerate reservoir according to the transverse continuity index and the longitudinal continuity index obtained in step S120, so as to quantitatively characterize the stratifying strength of the current reservoir by using the stratiform index. Further, in step S1301, the longitudinal continuity index and the lateral continuity index are divided by using a layer index calculation formula, and a quantized layer index is expressed by a ratio calculation result of the two. Wherein the calculation formula of the laminar index is represented by the following expression:
in the formula (7), LI represents a lamellarity index. In the practical application process, the smaller the stratiform index value is, the weaker the stratifying property of the reservoir is, and the closer the current reservoir is to the reservoir with the weaker stratifying property; conversely, the larger the value of the stratiform index is, the stronger the stratifying property of the reservoir is, and the closer the current reservoir is to the reservoir with the stronger stratifying property.
Thus, in the embodiment of the present invention, the stratiform index parameter that can be used for quantitatively characterizing the stratifying strength of the current reservoir is calculated by using the steps S110 to S130.
In another embodiment, a stratigraphic index curve for a current reservoir can also be generated based on the methods described above. Fig. 3 is a specific flowchart of the method for evaluating the stratifying property of a conglomerate reservoir according to the embodiment of the present application, in which a stratifying index curve is generated. As shown in fig. 3, first, step S301 acquires a static or dynamic image of a glutenite reservoir using an electrical imaging logging technique. Step S302, according to preset row intervals and column intervals, grid division processing is carried out on the static or dynamic images of the current glutenite reservoir. It should be noted that step S301 is similar to step S110, and step S302 is similar to step S201, and therefore will not be described herein again.
Then, the process proceeds to step S303. Step S303 is to perform window scanning processing on the static or dynamic image of the current glutenite reservoir subjected to gridding division processing according to a preset window and a preset depth step length, and to construct a conductivity matrix of each scanning window. The lateral width of the scanning window corresponds to the width of the current static or dynamic image, and the length of the preset depth step is generally smaller than the longitudinal width of the scanning window. In the current gravel reservoir static or dynamic image, the part framed by the scanning window is the image to be processed. In step S303, the scanning window is moved along the depth direction (up or down) according to a preset scanning window and a preset depth step for scanning, so as to construct a conductivity matrix of each scanning window, and then the process proceeds to step S304. Step S304 calculates a lateral continuity index of each scanning window (each scanning window is taken as an image to be processed), then step S305 calculates a longitudinal continuity index of each scanning window, and then step S306 calculates a layer index of each scanning window. It should be noted that step S304 is similar to step S1203, step S305 is similar to step S1204, and step S306 is similar to step S130, so that the description thereof is omitted here. Finally, step S307 connects all the stratiform indices according to the reservoir depth range corresponding to each scanning window, and generates a stratiform index curve for the current reservoir.
In embodiments of the present invention, it is desirable to acquire images of different reservoirs within the area to be studied using an electrical imaging logging technique. In the glutenite reservoir, the stratifying performance of the reservoir is quantified and represented by utilizing electrical imaging image data, and the stratifying index of the glutenite reservoir is constructed on the basis of the transverse continuity and the longitudinal continuity of the image. The embodiment of the invention for generating the stratiform index parameter is further explained by taking an EL oilfield rising section glutenite reservoir as an example and combining the accompanying drawings, and comprises the following steps:
(1) and selecting a plurality of representative gravel reservoir images with different stratifications in the research block as research samples of the reservoir in the block based on the electrical imaging logging data. Fig. 4 shows two examples of images of glutenite reservoirs with different stratifications.
(2) A dynamic or static image conductivity matrix of the reservoir is obtained.
(3) And calculating the transverse continuity index of the current reservoir by sequentially using the formula (1), the formula (2) and the formula (3).
(4) And calculating the longitudinal continuity index of the current reservoir by sequentially using the formula (4), the formula (5) and the formula (6).
(6) And constructing a reservoir stratum index LI by using the ratio of the longitudinal continuity index CONT _ V to the transverse continuity index CONT _ H. The smaller the stratiform index is, the weaker the stratifying property of the reservoir is, and the closer the stratifying property is to the reservoir with the weaker stratifying property is; conversely, the closer to a reservoir, the more lamellar the reservoir.
Further, the following embodiment of the invention for generating a stratiform index curve is further described by taking an EL oil field for a section of a conglomerate reservoir as an example and combining the accompanying drawings, and comprises the following steps:
(1) collecting an electrical imaging image of a reservoir in a certain depth section;
(2) based on a preset window, carrying out window scanning on the reservoir image from top to bottom according to positions at different depths;
(3) constructing a conductivity matrix of each scanning window;
(4) sequentially calculating a transverse continuity index, a longitudinal continuity index and a layer index according to each conductivity matrix;
(5) and connecting the stratiform indexes according to the depth position (range) corresponding to the scanning window, and drawing a stratiform index curve of the reservoir in the depth section.
Fig. 6 is a schematic diagram illustrating an example of calculation of a stratiform index in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application. As shown in fig. 6, the Tt _ X _ Continuity _ d, Tt _ Y _ Continuity _ d, and Tt _ XY _ LayerIndex _ d curves are the horizontal Continuity, the vertical Continuity, and the stratigraphic index of the reservoir based on the moving images; the Tt _ X _ Continuity, Tt _ Y _ Continuity, and Tt _ XY _ LayerIndex curves are static image-based reservoir lateral Continuity, reservoir longitudinal Continuity, and reservoir stratigraphic indices. In the actual processing, both still images and moving images are used, but moving images are generally used.
In addition, after the calculation of the stratifying index is completed, the stratifying index of the reservoir can reflect the capacity condition of the reservoir, so that whether the reservoirs at different positions in the area to be researched are effective reservoirs can be further conveniently and quickly identified in order to establish the relationship between the stratifying index and the capacity condition.
Further, the quantitative characterization method of the present invention further includes: and determining a layering index standard for identifying whether the reservoir in the block to be researched is effective or not according to the layering indexes corresponding to different reservoirs in the block to be researched and by combining the productivity information corresponding to different reservoirs so as to identify whether the reservoir in the block to be researched is effective or not.
Specifically, firstly, the capacity information of a plurality of representative reservoirs in a block to be researched is collected to be used as a standard to determine a sample set, and each sample in the sample set is subjected to reservoir attribute classification. Wherein the reservoir attributes include: a valid reservoir or an invalid reservoir. The effective reservoir represents that the test oil or production situation reaches the industrial capacity, and the ineffective reservoir represents that the test oil or production situation belongs to the dry layer. Then, the layer index of each sample in the sample set is calculated in the steps S110 to S130. And finally, correlating the reservoir attribute of each sample reservoir with the corresponding layer index, and determining the layer index standard of the reservoir to be researched, so as to judge the effectiveness of the reservoir in the block to be researched by utilizing the layer index standard. The stratifying index standard is an index parameter which can automatically identify whether the reservoir of the block to be researched is effective or not.
Further, when judging whether a certain reservoir (to be evaluated) in the current block to be researched is effective, the layer index standard is used as an effectiveness evaluation threshold, the layer index corresponding to the reservoir to be evaluated is compared with the layer index standard, and the effectiveness of the current reservoir to be evaluated is judged by utilizing the comparison result of the layer index standard and the layer index standard. And judging the current reservoir to be an invalid reservoir when the layer index corresponding to the current reservoir to be evaluated is larger than the current layer index standard. In addition, when the layer index corresponding to the current reservoir to be evaluated is smaller than the current layer index standard, the current reservoir is judged to be an effective reservoir.
In one embodiment, a plurality of representative glutenite reservoir images in the area to be studied are selected as a study sample set of the reservoir in the area based on the electrical imaging logging data. Firstly, determining the reservoir productivity of each reservoir sample according to the test oil and/or test production data of the current region to be researched, and dividing each reservoir sample into an effective reservoir and an ineffective reservoir; secondly, establishing the stratiform indexes of different glutenite reservoir samples by utilizing the steps S110 to S130; and thirdly, rapidly establishing an identification standard of the effective reservoir by comparing the layer indexes of the effective reservoir and the ineffective reservoir and combining the layer indexes corresponding to each reservoir sample. Fig. 7 is a schematic diagram illustrating an example of determining a stratifying index criterion in the method for evaluating the stratifying property of a glutenite reservoir according to the embodiment of the present application. As shown in fig. 7, 19 sample layers are selected from the sample set, wherein the 1 st to 8 th are invalid reservoirs, and the 9 th to 19 th are valid reservoirs, and the stratiform index bound of the valid reservoir and the invalid reservoir in the region is about 1.5 (as shown by the dotted line in fig. 7) as can be seen from the statistical chart. When the stratiform index of a reservoir in the region to be researched is more than 1.5, judging the current reservoir as an invalid reservoir; and when the stratiform index of the reservoir in the region to be researched is less than 1.5, judging that the current reservoir is an effective reservoir, wherein the judgment accuracy of the effectiveness of the whole reservoir is 84.2%.
On the other hand, based on the method for evaluating the stratifying property of the glutenite reservoir, the invention also provides a system for evaluating the stratifying property of the glutenite reservoir. Fig. 8 is a block diagram of a system for evaluating the stratifying property of a glutenite reservoir according to an embodiment of the present invention. As shown in fig. 8, the system includes: an image acquisition module 81, a continuity index generation module 82, and a stratification index generation module 83. The image obtaining module 81 is implemented according to the method described in the step S110, and is configured to collect the gravel reservoir image by using an electrical imaging logging technology. The continuity index generation module 82 is implemented according to the method described in the above step S120, and is configured to perform gridding division on the current image according to preset row intervals and column intervals, obtain a current reservoir transverse continuity index from the standard deviation of the row conductivity value of each row, and obtain a corresponding longitudinal continuity index from the standard deviation of the column conductivity value of each column. The stratifying index generating module 83 is implemented according to the method described in the above step S130, and is configured to calculate the stratifying index of the current reservoir according to the transverse continuity index and the longitudinal continuity index, so as to quantitatively characterize the stratifying strength of the current reservoir.
Further, the continuity index generation module 82 includes: a conductivity matrix construction unit 821, a lateral continuity index calculation unit 822, and a longitudinal continuity index calculation unit 823. Therein, a conductivity matrix construction unit 821 is configured to construct a conductivity matrix with respect to the current conglomerate reservoir image. A transverse continuity index calculation unit 822 configured to calculate a row conductivity standard deviation of each row in the current image and calculate an average value of the row conductivity standard deviations of the rows, resulting in a transverse continuity index. A longitudinal continuity index calculation unit 823 configured to calculate a column conductivity standard deviation of each column in the current image, and calculate an average value of the column conductivity standard deviations of the columns, resulting in a longitudinal continuity index.
Further, the lamellarity index generating module 83 is configured to divide the longitudinal continuity index and the transverse continuity index to obtain a lamellarity index.
The quantitative characterization system of the present invention further comprises: a curve plotting module 84. The curve plotting module 84 includes: a window scanning unit 841, a window slice index calculation unit 842, and a slice index curve generation unit 843. The window scanning unit 841 is configured to scan the current image subjected to the gridding partition processing according to a preset window, and construct a conductivity matrix for each scanning window. A window layer index calculation unit 842 configured to calculate the horizontal continuity index, the vertical continuity index and the corresponding layer index corresponding to each scanning window by executing the continuity index generation module 82 and the layer index generation module 83. And the layer index curve generation unit 843 is configured to connect the layer indexes corresponding to all the scanning windows according to the reservoir depths, and generate a layer index curve for the current reservoir.
In addition, the quantitative characterization system of the present invention further includes: a criteria generation module 85. The standard generation module 85 is configured to determine a stratifying index standard for identifying whether the reservoir in the block to be researched is valid according to the stratifying indexes corresponding to different reservoirs in the block to be researched and by combining the productivity information of different reservoirs, so as to judge the validity of the reservoir in the block to be researched by using the stratifying index standard.
The invention discloses a method and a system for evaluating the stratifying property of a glutenite reservoir. The method and the system are based on an electrical imaging logging technology, and utilize the conductivity image of the reservoir to calculate the transverse continuity and the longitudinal continuity, so as to further obtain a quantitative index (stratiform index) for representing the stratifying strength of the reservoir. In addition, by utilizing the calculation process of constructing the quantitative index of the reservoir stratum index and combining the capacity conditions of different reservoir stratum samples in the block to be researched, the key stratum index standard for evaluating the reservoir stratum effectiveness in the current block to be researched is determined, and the key stratum index standard is used for automatically judging and identifying the reservoir stratum effectiveness in the block to be researched. The invention realizes the quantitative characterization of the stratifying property of the reservoir stratum, determines the stratiform index limit of the effective reservoir stratum, provides an effective means for the effective reservoir stratum identification based on electrical imaging, and has the effective reservoir stratum identification accuracy rate of 84.2 percent.
Although the embodiments of the present invention have been described above, the above descriptions are only for the convenience of understanding the present invention, and are not intended to limit the present invention. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (6)
1. A method for evaluating the stratifying property of a glutenite reservoir is characterized by comprising the following steps:
acquiring a conglomerate reservoir image by utilizing an electrical imaging logging technology;
step two, carrying out gridding division on the image according to preset row intervals and column intervals, obtaining a current reservoir stratum transverse continuity index according to the standard variance of the row conductivity value of each row, and obtaining a corresponding longitudinal continuity index according to the standard variance of the column conductivity value of each column,
constructing a conductivity matrix for the glutenite reservoir image;
calculating the lateral continuity index and the longitudinal continuity index from the conductivity matrix of the image using the following expressions:
wherein CONT _ H denotes a lateral continuity index of the current image, CONT _ V denotes a longitudinal continuity index of the current image, i denotes a row number, j denotes a column number, STDD _ H (i) denotes a row conductivity standard deviation of the ith row, STDD _ V (j) denotes a column conductivity standard deviation of the jth column, C (i, j) denotes a conductivity value of the jth grid of the ith row, MEAN _ H (i) denotes a row conductivity of the ith row, MEAN _ V (i) denotes a column conductivity value of the jth column, wherein the row conductivity value is an average of conductivity values of all grids in each row, and the column conductivity value is an average of conductivity values of all grids in each column;
and thirdly, performing division operation on the longitudinal continuity index and the transverse continuity index to obtain the stratiform index of the current reservoir, and quantitatively representing the strength/weakness of the stratiform property of the current reservoir by using the size/magnitude of the stratiform index.
2. The method of claim 1, further comprising:
and determining a stratifying index standard for identifying whether the reservoirs in the block to be researched are effective or not according to the stratifying indexes corresponding to different reservoirs in the block to be researched and combining the productivity information of the different reservoirs, so as to judge the effectiveness of the reservoirs in the block to be researched by utilizing the stratifying index standard.
3. The method of claim 1, further comprising:
scanning the images subjected to gridding division processing according to a preset window to construct a conductivity matrix of each scanning window;
calculating the transverse continuity index, the longitudinal continuity index and the corresponding layer index corresponding to each scanning window by executing the second step and the third step;
and connecting all the layer indexes according to the depth of the reservoir, and generating a layer index curve aiming at the current reservoir.
4. A system for evaluating the stratifying properties of a conglomerate reservoir, comprising:
an image acquisition module configured to acquire a conglomerate reservoir image using an electrical imaging logging technique;
a continuity index generation module configured to perform gridding division on the image according to preset row intervals and column intervals, obtain a current reservoir transverse continuity index from a standard deviation of a row conductivity value of each row, and obtain a corresponding longitudinal continuity index from a standard deviation of a column conductivity value of each column, where the continuity index generation module includes:
a conductivity matrix construction unit configured to construct a conductivity matrix with respect to the glutenite reservoir image;
a lateral continuity index calculation unit and a longitudinal continuity index calculation unit configured to calculate the lateral continuity index and the longitudinal continuity index, respectively, from the conductivity matrix of the image using the following expressions:
wherein CONT _ H denotes a lateral continuity index of the current image, CONT _ V denotes a longitudinal continuity index of the current image, i denotes a row number, j denotes a column number, STDD _ H (i) denotes a row conductivity standard deviation of the ith row, STDD _ V (j) denotes a column conductivity standard deviation of the jth column, C (i, j) denotes a conductivity value of the jth grid of the ith row, MEAN _ H (i) denotes a row conductivity of the ith row, MEAN _ V (i) denotes a column conductivity value of the jth column, wherein the row conductivity value is an average of conductivity values of all grids in each row, and the column conductivity value is an average of conductivity values of all grids in each column;
and the stratiform index generating module is configured to perform division operation on the longitudinal continuity index and the transverse continuity index to obtain a stratiform index of the current reservoir, and is used for quantitatively representing the strength/weakness of the stratiform of the current reservoir by using the size/magnitude of the stratiform index.
5. The system of claim 4, further comprising:
and the standard generation module is configured to determine a layered index standard for identifying whether the reservoirs in the block to be researched are effective or not according to the layered indexes corresponding to different reservoirs in the block to be researched and by combining the productivity information of the different reservoirs, so as to judge the effectiveness of the reservoirs in the block to be researched by using the layered index standard.
6. The system of claim 4, further comprising a curve plotting module, wherein the curve plotting module comprises:
the window scanning unit is configured to scan the images subjected to gridding division processing according to preset windows, and a conductivity matrix related to each scanning window is constructed;
a window layer index calculation unit configured to execute the continuity index generation module and the layer index generation module, and calculate the transverse continuity index, the longitudinal continuity index and the corresponding layer index corresponding to each scanning window;
and the stratiform index curve generating unit is configured to connect all the stratiform indexes according to the reservoir depths and generate the stratiform index curve aiming at the current reservoir.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010000508.7A CN110966000B (en) | 2020-01-02 | 2020-01-02 | Method and system for evaluating layering property of glutenite reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010000508.7A CN110966000B (en) | 2020-01-02 | 2020-01-02 | Method and system for evaluating layering property of glutenite reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110966000A CN110966000A (en) | 2020-04-07 |
| CN110966000B true CN110966000B (en) | 2021-06-29 |
Family
ID=70037957
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010000508.7A Expired - Fee Related CN110966000B (en) | 2020-01-02 | 2020-01-02 | Method and system for evaluating layering property of glutenite reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110966000B (en) |
Family Cites Families (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7133779B2 (en) * | 2003-05-08 | 2006-11-07 | Schlumberger Technology Corporation | Automated borehole geology and petrophysics interpretation using image logs |
| CN102200008B (en) * | 2010-03-26 | 2013-03-13 | 中国石油天然气股份有限公司 | Reservoir effectiveness identification method based on electrical imaging logging |
| US8725479B2 (en) * | 2011-11-30 | 2014-05-13 | Baker Hughes Incorporated | Use of monte carlo simulations for establishing correlations and their uncertainties |
| EP2948895B1 (en) * | 2013-01-25 | 2019-06-12 | Services Petroliers Schlumberger | Pressure transient testing with sensitivity analysis |
| US9651699B2 (en) * | 2013-02-20 | 2017-05-16 | Apache Corporation | Methods for determining well log attributes for formation characterization |
| CN103115820B (en) * | 2013-03-07 | 2014-08-13 | 北京科技大学 | Method for confirming anisotropism of jointed rock mass |
| GB2512391B (en) * | 2013-03-28 | 2020-08-12 | Reeves Wireline Tech Ltd | Improved borehole log data processing methods |
| US9366135B2 (en) * | 2013-10-08 | 2016-06-14 | Exxonmobil Upstream Research Company | Automatic dip picking from wellbore azimuthal image logs |
| CN104912547A (en) * | 2014-03-11 | 2015-09-16 | 中国石油化工集团公司 | Method for evaluating heterogeneous characteristics of reservoir continuously and quantitatively by applying resistivity imaging logging data |
| US9921329B2 (en) * | 2014-11-20 | 2018-03-20 | Apache Corporation | Automated method for selecting positions within formations from which to extract samples thereof |
| CN104732562A (en) * | 2015-03-24 | 2015-06-24 | 中国石油集团川庆钻探工程有限公司 | Electric imaging logging image similarity degree comparison method |
| CN104808248A (en) * | 2015-04-22 | 2015-07-29 | 长江大学 | Attitude of stratum detecting method based on electric imaging logging full hole image |
| CN104793254B (en) * | 2015-04-28 | 2017-08-25 | 中国科学院地质与地球物理研究所 | A kind of recognition methods of electrical thin layer in underground and device |
| CN104849766A (en) * | 2015-05-19 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司 | Electrical imaging logging image stratum attitude identification method |
| US11307323B2 (en) * | 2015-09-30 | 2022-04-19 | Schlumberger Technology Corporation | Methods and systems to analyze bed boundary detection |
| CN106447670A (en) * | 2016-08-31 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司 | Hole parameter automatic calculation method based on electric imaging logging image |
| CN107219159B (en) * | 2017-04-19 | 2019-07-09 | 中国石油天然气股份有限公司 | A kind of method and apparatus of determining reservoir heterogeneity |
| CN109085663A (en) * | 2018-07-20 | 2018-12-25 | 中国矿业大学 | A kind of tight sandstone reservoir stratification seam recognition methods |
-
2020
- 2020-01-02 CN CN202010000508.7A patent/CN110966000B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN110966000A (en) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107703560B (en) | shale lithofacies fine identification method based on triple information | |
| CN102175832B (en) | Method for determining optimal saturation computing model for typical reservoir | |
| US10228478B2 (en) | Method of searching for oil-gas reservoir based on trap-3D software | |
| CN106526693A (en) | Crack recognition method and device | |
| WO2016070587A1 (en) | Method for looking for oil and gas reservoir on the basis of trap-3d software | |
| US10641923B2 (en) | Method for characterizing and exploiting a subterranean formation comprising a network of fractures | |
| CN103026202A (en) | Method for obtaining consistent and integrated physical properties of porous media | |
| CN111722277B (en) | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration | |
| CN103345566A (en) | Geochemical exploration anomaly recognizing and evaluating method based on contents contained in geology | |
| CN110632659A (en) | Fault space closure analysis method, device and system | |
| Guo et al. | Multiple-point geostatistics-based three-dimensional automatic geological modeling and uncertainty analysis for borehole data | |
| CN115166853A (en) | Shale gas reservoir natural fracture model establishing method and device, electronic equipment and medium | |
| KR101131517B1 (en) | Processing system and method for subsurface information in oil sands reservoir using statistical analysis of well logging data | |
| CN110966000B (en) | Method and system for evaluating layering property of glutenite reservoir | |
| Lu et al. | Identifying flow units by FA-assisted SSOM—An example from the Eocene basin-floor-fan turbidite reservoirs in the Daluhu Oilfield, Dongying Depression, Bohai Bay Basin, China | |
| Zhou et al. | Fracture prediction of tight sandstone reservoirs using outcrops and log curve-based extremum method: a case study of the Chang 7 member of the Yanchang Formation in block X, Ordos Basin | |
| CN109738955B (en) | Metamorphic rock lithology comprehensive judgment method based on component-structure classification | |
| Sarkheil et al. | The fracture network modeling in naturally fractured reservoirs using artificial neural network based on image loges and core measurements | |
| CN115437013A (en) | Oil reservoir small slip distance slip fracture analysis method, medium and equipment | |
| CN108804848B (en) | Method for calculating box dimension of logging curve | |
| CN114153002A (en) | Three-dimensional geological modeling method and device for natural fracture of reservoir, electronic equipment and medium | |
| CN114460639A (en) | Method and device for predicting permeability of shale oil reservoir | |
| CN113627069A (en) | Well testing dynamic yield evaluation method and system for fracture-cavity type oil reservoir oil and gas well | |
| CN117687096B (en) | Proxy model construction method for predicting small-scale fracture-cavity distribution | |
| CN112906242B (en) | Geophysical modeling method based on combination of naive Bayes method and proximity classification method |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200917 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Applicant after: CHINA NATIONAL PETROLEUM Corp. Applicant after: CHINA PETROLEUM LOGGING Co.,Ltd. Applicant after: North China branch of CNPC logging Co.,Ltd. Address before: 062550 No. 002, Bohai East Road, Renqiu City, Cangzhou City, Hebei Province Applicant before: North China branch of CNPC logging Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210629 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |