Method for lithological analysis of a drilling cuttings sample

Abstract

Method for lithological analysis of a drill cuttings sample. This method (50) comprises the steps of: acquiring (110) an image of the sample composed of grains using a sensor; applying (120) a convolutional neural network to the image (I) to obtain a contour map (C), which associates each pixel of the image with a value equal to one when said pixel belongs to a grain contour and equal to zero otherwise; applying (130) a segmentation algorithm to the contour map (C) to obtain a segmented image (S) comprising a plurality of regions, each region being separated from a neighboring region by a contour with a thickness of one pixel; and applying (140) one or more geometric processing steps to the segmented image (S) to calculate physical quantities for each grain in order to characterize the lithology of the sample (3). Figure for the abstract: Figure 2

Description

Title of the invention: Method for lithological analysis of a sample of drilling debris

[0001] The present invention relates to the field of geological analysis for the optimization of drilling in the energy sector.

[0002] More specifically, the invention relates to methods for lithological characterization of a sample of drilling debris.

[0003] Drilling cuttings analysis is a crucial step in oil and gas drilling processes. It consists of examining samples of drilling mud that rise to the surface during the drilling operation, in order to deduce the geological characteristics of the subsoil.

[0004] These mud samples do indeed contain rock and mineral particles which provide valuable information on the composition and stratigraphy of the subsoil.

[0005] The analysis of sludge makes it possible in particular to determine the presence of oil, gas or hydrogen reservoirs and to estimate their depth, their thickness, as well as the geological and physical properties of the layers crossed.

[0006] It also allows the pressure and temperature in the well to be assessed, which is essential for the safety and design of the installations.

[0007] Drilling mud analysis is carried out using specific instruments and techniques throughout the drilling process.

[0008] Currently, samples are regularly taken at specific intervals during the drilling operation. These samples are taken directly from the drilling mud that rises to the surface through the well.

[0009] The collected samples are sent to a laboratory where a series of analyses are carried out. Commonly used analytical techniques include imaging techniques, mass spectrometry, chromatography, atomic absorption, and other methods for determining the chemical composition of the samples.

[0010] The data obtained from the analysis of the muds are then interpreted by geologists, drilling engineers and other experts, who use this information to assess the geological composition of the subsoil and deduce the properties of the reservoir, the fluid pressure, the presence of contaminants and other important parameters.

[0011] There is a need for an imaging method that can be performed on mud samples, which is fast (almost real-time), can be implemented on or in the immediate vicinity of the drilling site, and uses material means simple and robust, while allowing precise information to be obtained enabling lithological characterization of the geological layer during drilling.

[0012] The invention is therefore intended to meet this need.

[0013] To this end, the invention consists of a lithological analysis method for a sample of drilling debris, which may be wet or dry, and comprising a plurality of grains, said method, implemented by a computer associated with a sensor, being characterized in that it comprises, during an operational phase, the steps of: acquisition of an image of the sample by the sensor;application, on the acquired image, of a first deep learning classification algorithm structured as a convolutional neural network, the first classification algorithm being optimized following the implementation of a training phase of the process, so as to obtain a contour map, the contour map associating with each pixel of the acquired image a probability value, which is equal to one when it is estimated that said pixel of the acquired image belongs to a contour of a grain and equal to zero when it is estimated that said pixel of the acquired image does not belong to the contour of a grain; application, on the contour map, of a segmentation algorithm, so as to obtain a segmented image comprising a plurality of regions, a region being separated from a neighboring region by a contour having a thickness of a single pixel, each region being associated with a probable grain delimited in the acquired image;and, application to the segmented image, of one or more geometric processing(s) to calculate a plurality of physical quantities for each probable grain, the values ​​of the calculated physical quantities allowing a lithological characterization of the sample. ;

[0014] According to particular embodiments, the process comprises one or more of the following characteristics, taken individually or in all technically possible combinations:

[0015] - the first structured deep learning classification algorithm as a convolutional neural network adopts a U-Net structure.

[0016] - the sensor is adapted to acquire a monochrome or polychrome image of the sample.

[0017] - the sensor is a camera with a depth of field extent allowing imaging of multiple layers of the sample.

[0018] - the training phase includes an automatic population step of a training database and a training stage of the first deep learning classification algorithm structured as a convolutional neural network from the data in the training database, the automatic population stage comprising the sub-steps of: acquiring a training image on a training sample using a sensor identical training image to the sensor; application, on the acquired training image, of one or more deterministic image processing algorithm(s) to improve the quality of the acquired training image, to obtain a processed training image; application, on the processed training image, of a second deep learning classification algorithm already trained but in a non-specific way, in order to obtain a segmented training image; and, storage of a new entry in the training database, the new entry associating the acquired training image and the segmented training image estimated for said acquired training image.

[0019] - the application, on the acquired training image, of one or more algorithm(s) deterministic image processing involves the application of one or more algorithms chosen from: an intensity equalization algorithm; a contrast equalization algorithm; and a gamma correction algorithm.

[0020] The invention also relates to an analysis system comprising a computer connected to a sensor, the analysis system being adapted for the implementation of a lithological analysis process of a sample of drilling debris conforming to the previous process.

[0021] The invention also relates to a computer program product comprising software instructions which, when executed by a computer, implement the previous process.

[0022] The invention and its advantages will be better understood upon reading the following detailed description of a particular embodiment, given solely by way of non-limiting example, this description being made with reference to the accompanying drawings in which:

[0023] [Fig-1] The [Fig. 1] is a schematic representation of a system for implementing the analysis method according to the invention;

[0024] [Fig.2] The [Fig.2] is a block representation of a preferred embodiment of the analysis method according to the invention;

[0025] [Fig. 3] Fig. 3 is a representation of the different pixel matrices during the implementation of the process of Figure 2; and,

[0026] [Fig.4] The [Fig.4] is a schematic representation of the architecture of the CNN algorithm trained in the training phase and executed in the operational phase of the process of the [Fig.2].

[0027] The [Fig. 1] is a schematic representation of a system for implementing the analysis method according to the invention.

[0028] The system 1 includes an optical sensor 2, enabling the acquisition of an image of a sludge sample 3 placed in the field of the optical sensor 2.

[0029] The system 1 includes a computer 4. The latter includes computing means, such as a processor 5, storage means, such as a memory 6, and an input / output interface 7, enabling in particular the computer 4 to be connected to the sensor 2.

[0030] The memory 6 stores, among other things, computer program instructions, in particular a program 8 whose execution enables the implementation of the method 50 according to the invention.

[0031] The analysis method 50 according to the invention includes a training phase 200 of a parametric model to obtain a trained model using training data from a training database, and then an operational phase 100 using the trained model to perform the analysis of newly acquired images.

[0032] The computer 4 can be used for both the training phase 200 and the operational phase 100. The program 8 then comprises, on the one hand, a training module 12 whose execution, using the data from a training database 9, allows the completion of the training phase 200, and on the other hand, an operational module 11, whose execution allows the completion of the operational phase 100.

[0033] Alternatively, each phase of the analysis process is carried out on different computer systems, provided that the optical sensors of each of these systems are identical.

[0034] Fig. 2 is a schematic representation, in block form, of a preferred embodiment of the operational phase 100 of the analysis process 50.

[0035] Prior to the implementation of operational phase 100, a sample is prepared.

[0036] A drill head, for example of an oil well, regularly brings mud to the surface.

[0037] These muds are characteristic of the geological layer that the drill head has reached and is drilling.

[0038] Once at the surface, the sludge passes through filter screens to retain only particles (or grains) of a predefined size. For example, the finest and largest grains are removed, so that only grains of an intermediate size remain.

[0039] A sample is taken at the outlet of the filter sieves. It therefore consists of an aggregate of grains, exhibiting low cohesion between them.

[0040] The sample may be wet or dry, after evaporation of all or part of its water content.

[0041] Operational phase 100 begins with a step 110 of acquiring an image I of sample 3.

[0042] The image of the sample is obtained by means of the optical sensor 2.

[0043] Sensor 2 is of the camera type.

[0044] Preferably, it is an improved camera, in the sense that it allows for near-field shooting (macrophotography) of the surface of sample 3, but with an extended depth of field, for example of a few millimeters, preferably between 2 and 3 mm.

[0045] This is advantageous since the sample consists of grains. However, the depth of field in microphotography is shallow, so that only the grains located precisely in the image plane of the camera are sharp and their edges are clearly defined, while the grains located immediately in front of or behind the image plane are blurred and their edges are indistinct. Thus, using a shallow depth of field makes it possible to obtain an image in which all the grains of the sample are sharp and their edges are clearly defined.

[0046] This also allows, when the sample results from the superposition of several layers of grains, to be able to image not only the grains of the upper layer, constituting the surface of the sample, but also the layer(s) immediately below the upper layer.

[0047] In the preferred embodiment, the sensor 2 automatically stacks a series of raw images, obtained from as many different image planes, in order to extend the depth of field. The sensor 2 thus preferably delivers an image in the form of a preprocessed data file. This is preferable, in terms of time savings, to a sensor delivering the series of raw images and a computer programmed to stack the raw images to obtain the acquired image, which is the subject of the subsequent steps in the process.

[0048] The acquired image is a pixel matrix. Typically, the image has on the order of 30,000,000 pixels.

[0049] In one embodiment, the acquired image is a chromatic image, for example red-green-blue - RGB (red-green-blue). Each pixel is then associated with three intensity values, respectively a first value for red, a second value for green, and a third value for blue.

[0050] Figure 3A shows a chromatic image reproduced.

[0051] For the acquisition of a chromatic image, the sample is illuminated by white light.

[0052] Alternatively, the acquired image is a monochromatic image, for example in greyscale. Each pixel is then associated with an intensity value.

[0053] For the acquisition of a monochromatic image, the sample is illuminated by monochromatic light, such as for example ultraviolet - UV light.

[0054] Phase 100 of the process 50 then includes a step 120 of analyzing the acquired image I, in order to obtain a contour map C of the sample 3.

[0055] Step 120 consists of applying to the acquired image I a first deep learning classification algorithm structured as a convolutional neural network - CNN (“Convolutional Neural Network”).

[0056] The CNN algorithm, the structure of which will be presented in more detail below with reference to Figures 5, is a parametric model. The training phase 200 of the process 50 according to the invention calculates the optimal parameter values ​​in order to obtain a trained model. This model is then used in the operational phase 100 of the process 50.

[0057] The CNN algorithm takes as input the image acquired I at the acquisition step 110.

[0058] The CNN algorithm outputs a contour map C. This map has the same pixel resolution as the input image. Each pixel of the map is associated with a unique pixel of the input image. Each pixel of the map is associated with a probability between 0 and 1. A pixel is assigned a value of 1 when the CNN algorithm determines that this pixel belongs to the grain contour in the sample image. A pixel is assigned a value of 0 when the CNN algorithm determines that this pixel does not belong to a grain contour in the sample image.

[0059] It should be noted that on the contour mapping C, two pixels associated with the null value can be separated, along the segment that joins them, by several pixels associated with the unit value or close to the unit.

[0060] In Figure 3B, a contour map is reproduced which is estimated on the acquired image reproduced in Figure 3A.

[0061] Operational phase 100 continues with a step 130 of segmenting the contour mapping C so as to obtain a segmented image S.

[0062] Step 130 consists of applying a segmentation algorithm to the map C, which groups the pixels according to their neighborhood into a plurality of regions, each region corresponding to a grain of the acquired image.

[0063] The segmentation algorithm implemented is for example of the "by watershed" type known to the person skilled in the art.

[0064] The result of this segmentation step is, for example, represented by a segmented image, such as image S in Figure 3B.

[0065] The segmented image S is a partition of the acquired image into a plurality of regions, each of which is closed by a contour line, with two adjacent regions separated by a contour line having a width equal to only one pixel. In other words, the segmented image S is a binary image: the pixels of the different regions are associated with the value zero, while the pixels of the different contours are associated with the value unit.

[0066] Segmentation thus makes it possible to delimit very precisely the outline of a probable grain, as identified by the CNN algorithm in the original acquired image.

[0067] Furthermore, each region of the segmented image receives a unique identifier. Step 130 thus makes it possible to group the pixels of the contour map C according to a neighborhood criterion, by associating them with the same identifier.

[0068] In Figure 3C, a segmented image is reproduced which is obtained from the contour mapping of Figure 3B. False colors are used to identify each delimited region in the image.

[0069] The operational phase 100 advantageously includes a measurement step 140 from the segmented image S and advantageously also from the acquired image I.

[0070] For each region of the segmented image S, different geometric treatments make it possible to measure, for the probable grain corresponding to this region, a set of primary physical quantities, such as its dimensions (surface area, perimeter, average diameter...), its shape (position of the center, principal axis of maximum width, principal axis of minimum width, sphericity...), etc.

[0071] These primary physical quantities are advantageously combined to calculate secondary physical quantities. For example, the roughness of a probable grain can be calculated by the ratio of its perimeter to its surface area.

[0072] The segmented image S is advantageously superimposed, like a mask, on the acquired image I. In the case of an RGB image, for example, for each region of the segmented image S, different colorimetric processing techniques then allow the measurement of a color for the corresponding probable grain. The color is an index of the mineral nature of the associated grain.

[0073] Then, the measurements relating to each region are combined to characterize the sample as a whole.

[0074] For example, a histogram of the sphericity measurements of the grains of the input image is constructed which allows, depending on the distribution it presents, the lithology of the sample to be described.

[0075] For example, a histogram of the surface measurements of the grains of the input image is constructed which allows, depending on the distribution it presents, to say whether the sample is that of a gravel, a sand, a silt or a clay.

[0076] Preferably, phase 100 is iterated to characterize different samples taken at the same time, i.e., corresponding to the same drilling depth. The accumulation of statistics then makes it possible to determine with greater robustness the nature of the geological layer encountered at that drilling depth.

[0077] Advantageously, phase 100 of the process according to the invention is iterated regularly during drilling, so as to obtain geological information for different depths and thus reconstruct the different geological layers traversed by the borehole.

[0078] Since the analysis step of the acquired image implements an artificial intelligence algorithm, the method according to the invention includes a training phase 200 of the CNN algorithm executed in step 120.

[0079] It should be noted that there is no training database available in the field of lithological analysis.

[0080] The training phase 200 of the process 50 includes in particular the steps of automatic population of the training database 9 to then allow the implementation of the actual training steps of the CNN algorithm.

[0081] Thus the training phase 200 of the process 50 includes a step 210 of acquiring a training image.

[0082] This acquisition is carried out under the same conditions as step 110, in particular with the same sensor 2 (or an identical sensor) so as to obtain a training image similar to the acquired image I.

[0083] In the 3D figure, an acquired training image is reproduced.

[0084] Phase 200 continues with a pre-processing step 215 of the training image acquired in step 210. Deterministic image processing is applied to the image to improve its quality and contrast, in order to better highlight the contours of the imaged grains.

[0085] For example, three image processing operations are implemented.

[0086] A first treatment is an image brightness equalization treatment.

[0087] This involves, for example, applying a "histogram equalization" treatment, which allows the redistribution of the intensity of the pixels of the training image to homogenize the intensity throughout this image.

[0088] Alternatively, this involves, for example, applying an "adaptive histogram equalization" process that calculates several histograms, each corresponding to a distinct section of the training image, and uses these different histograms to redistribute the brightness values ​​in the image. This process improves local contrast and edge definition in each section of an image.

[0089] A second treatment is an image contrast equalization treatment.

[0090] This involves, for example, applying a "contrast limited adaptive histogram equalization" (CLAHE) treatment. Equalization). This process divides the training image, preferably the output of the first processing, into contextual tiles. It then creates a histogram for each contextual tile, clips it to a predefined value, and redistributes the clipped amount to the remaining cells of the histogram. This technique makes hidden features of the training image more visible and sharpens edges in each part of the image.

[0091] The third processing is a gamma correction processing, which allows the brightness of the training image to be worked on, preferably at the output of the first or second processing.

[0092] Step 215 produces a processed It training image.

[0093] In Figure 3E, a processed training image is reproduced which is derived from the acquired training image reproduced in Figure 3D.

[0094] The next step, 220, of phase 200 is a segmentation step of the processed training image It.

[0095] A second artificial intelligence segmentation algorithm, already trained, is used to recognize edges in the processed training image. For example, the HED (Holistically-Nested Edge Detection) algorithm is used. This is a convolutional neural network. It is available as open source on the internet. It is trained in a non-specific manner on arbitrary images.

[0096] The result of step 220 is a segmented training image Se.

[0097] It should be noted that the execution of the HED algorithm remains long, on the order of 10 minutes to process an image, even a pre-processed one.

[0098] In Figure 3F, a segmented training image is reproduced which is estimated from the processed training image of Figure 3E.

[0099] In step 250, a new entry is recorded in the training database 9. The segmented training image Se was used as the label of the training image acquired on.

[0100] Once suitably populated, the database 9 is used for training the CNN algorithm in a step 270. The training conforms to the prior art.

[0101] In operation, the CNN algorithm is capable of producing a contour image S from an acquired image I in a few seconds.

[0102] Particularly advantageously, the process according to the invention gives very good results not only on dry samples, but also on wet samples.

[0103] The architecture of the CNN algorithm trained in step 270 and executed in step 120 corresponds to the U-Net structure.

[0104] The U-Net structure is described, for example, in the article Ronneberger, O., Fischer, P., & Brox, T. (2015, October). U-net: Convolutional networks for biomedical image segmentation. In International Conference on Medical image computing and computer-assisted intervention (pp. 234-241). Springer, Cham.

[0105] Thus, the CNN algorithm comprises, in series, an encoder 20, a connection layer 30, a decoder 40 and an output layer 50.

[0106] The encoder 20 takes as input an acquired image I, for example an RGB image.

[0107] The connection layer 30 is located between the output of the encoder 20 and the input of the decoder 40.

[0108] The decoder 40 outputs a contour map C.

[0109] The encoder 20 consists of a plurality of successive convolutional layers, 22; , where i is an index between 1 and N, N denoting the total number of layers in the CNN unit.

[0110] In [Fig.4], N is equal to 5.

[0111] The information flow passing through the encoder 20 goes from the layers with the lowest index (from layer 22i) to the layers with the highest index (towards layer 224).

[0112] Each layer is associated with a specific resolution.

[0113] Each layer 22; typically includes a convolution block in order to extract a feature map from a reduced feature map produced by the previous layer 22^ of the encoder 20.

[0114] Preferably, the convolution block performs a double convolution, so as to successively apply a first 3x3 convolution, a first ReLu function, a second 3x3 convolution and a second ReLu function on the reduced feature map.

[0115] Layer 22; includes a scaling block to calculate the reduced feature map for the next layer, from the feature map.

[0116] Preferably, the scaling block performs a sampling-based discretization process, such as maximum pooling, in particular a 2x2 maximum pooling function.

[0117] The connection layer 30, which links the lower layer 225 of the encoder 20 and the lower layer 425 of the decoder 40, typically includes a convolution block in order to extract a feature map from a reduced-scale feature map emitted by the encoder 20. The reduced-scale feature map is the reduced-scale feature map input to the decoder 40.

[0118] Preferably, the convolution block of the connection layer 30 performs a double convolution, so as to successively apply a first convolution 3x3, a first ReLu function, a second 3x3 convolution and a second ReLu function on the reduced feature map 36.

[0119] The decoder 40 consists of a plurality of successive convolutional layers, 42i, where i is an index between 1 and N, where N denotes the total number of layers. In the present embodiment, N is equal to 5.

[0120] The information flow passing through the decoder 40 goes from the layers with the highest index (from layer 424) to the layers with the lowest index (towards layer 42i).

[0121] Each layer 42; is associated with a specific resolution, identical to the resolution of layer 22; of the encoder 20 of the same index i.

[0122] Each layer 42; typically includes a scaling block which takes the feature map at the output of the previous 42m layer to determine a scaled feature map.

[0123] Layer 42; also includes a concatenation block performing a concatenation operation of the larger-scale feature map and a refined feature map.

[0124] Layer 42; then includes a transposed convolution block, performing a transposed convolution on the output of the concatenation block to obtain a feature map.

[0125] For example, the transposed convolution is a double convolution block which successively applies a first 3x3 convolution, a first ReLu function, a second 3x3 convolution and a second ReLu function.

[0126] According to the U-Net structure, for a particular resolution, the feature map of the encoder side 20 is transmitted to the decoder side 40 via a hop connection 17; on the [Fig.4].

[0127] Alternatively, the input to the first AI algorithm in step 130 takes as input several images of the same sample, for example an RGB image and a UV image. The algorithm must then be trained using training data corresponding to these inputs.

[0128] The present method allows for the automatic population of a training database using a trained but non-specific AI algorithm to label images which must be processed to help this AI algorithm and thus improve the quality of the label.

[0129] The present method makes it possible to train an AI algorithm specific to the application sought, in this case the analysis of sample images.

[0130] The present method allows for near real-time analysis of samples.

Claims

Demands

1. - Method (50) for lithological analysis of a sample (3) of drilling debris, which may be wet or dry, and comprising a plurality of grains, said method being implemented by a computer (4) associated with a sensor (2), characterized in that the method comprises, during an operational phase (100), the steps of: - acquisition (110) of an image of the sample by the sensor;- application (120), on the acquired image (I), of a first classification algorithm by deep learning structured as a convolutional neural network, the first classification algorithm being optimized following the implementation of a training phase (200) of the process, so as to obtain a contour map (C), the contour map associating with each pixel of the acquired image a probability value, which is equal to one when it is estimated that said pixel of the acquired image belongs to a contour of a grain and equal to zero when it is estimated that said pixel of the acquired image does not belong to the contour of a grain;- application (130), on the contour map (C), of a segmentation algorithm, so as to obtain a segmented image (S) comprising a plurality of regions, a region being separated from a neighboring region by a contour having a thickness of a single pixel, each region being associated with a probable grain delimited in the acquired image; and, - application (140) on the segmented image (S), of one or more geometric processing(s) to calculate a plurality of physical quantities for each probable grain, the values ​​of the calculated physical quantities allowing a lithological characterization of the sample (3).;

2. - Method according to claim 1, wherein the first deep learning classification algorithm structured as a convolutional neural network adopts a U-Net structure.

3. - Method according to claim 1 or claim 2, wherein the sensor (2) is adapted to acquire a monochrome or polychrome image of the sample (3).

4. - A method according to any one of claims 1 to 3, wherein the sensor (2) is a camera having a extended depth of field allowing imaging of multiple layers of the sample.

5. - A method according to any one of claims 1 to 4, wherein the training phase (200) comprises a step of automatically populating a training database (9) and a step of learning the first classification algorithm by deep learning structured as a convolutional neural network from the data of the training database, the automatic populating step comprising the substeps of: - acquisition (210) of a training image (le) on a training sample by means of a training sensor identical to the sensor (2); - application (215), on the acquired training image, of one or more deterministic image processing algorithm(s) allowing to improve a quality of the acquired training image, to obtain a processed training image (It);- application (220), on the processed training image, of a second deep learning classification algorithm already trained but in a non-specific way, in order to obtain a segmented training image (Se); and, - storage (250) of a new entry in the training database (9), the new entry associating the acquired training image (le) and the segmented training image (Se) estimated for said acquired training image.;

6. - A method according to claim 5, wherein the application, on the acquired training image, of one or more deterministic image processing algorithm(s) comprises the application of one or more algorithms chosen from: - an intensity equalization algorithm; - a contrast equalization algorithm; and, - a gamma correction algorithm.

7. Analysis system (1) comprising a computer (4) connected to a sensor (2), characterized in that the analysis system is adapted for the implementation of a lithological analysis process of a drill cutting sample according to any one of claims 1 to 6.

8. - Product computer program comprising program code instructions recorded on a computer-readable medium of a system conforming to claim 7, for the execution of the steps of the lithographic analysis process of a drill cutting sample according to any one of claims 1 to 6 when said program is executed on a computer.

Images