Logging lithology identification method based on physical information constraint

Abstract

The invention is suitable for the technical field of logging lithology identification, and provides a logging lithology identification method based on physical information constraint, which comprises the following steps: step 1, data processing; 2, constructing a graph; step 3, constructing a graph attention neural network; 4, constructing a mixed loss function; and 5, lithology prediction. According to the method, through the synergistic effect of data driving and physical information, the performance is obviously superior to that of various pure data driving models. Benefited from the strong regularization effect of physical constraints, the method has stronger generalization ability and prediction stability on new data. The problem that a traditional black box model may generate a result violating physical common knowledge is fundamentally solved, and the geologic rationality of an output result is ensured.

Description

Technical Field

[0001] This invention belongs to the field of well logging lithology identification technology, and particularly relates to a well logging lithology identification method based on physical information constraints. Background Technology

[0002] Well logging lithology identification is a fundamental and core task in oil, gas, geothermal resource exploration, and deep underground engineering. It indirectly infers the type of subsurface rocks by analyzing geophysical logging data continuously collected along the well depth. Accurate lithological profiles are crucial for reservoir evaluation, resource estimation, and development decisions.

[0003] The existing technical routes can be summarized as follows: (1) Traditional geological interpretation methods: such as well logging curve morphology method, cross plot analysis method, etc. These methods rely heavily on the professional experience of the analysts, are highly subjective, inefficient, and have limited discrimination ability when processing multi-dimensional and nonlinear well logging data, making it difficult to meet the needs of large-scale fine exploration. (2) Machine learning methods: such as support vector machine, random forest and XGBoost. These methods classify each depth sampling point as an independent sample, realizing the first step of automation. However, their data point-based processing paradigm ignores the spatial correlation of strata in the vertical direction due to continuous deposition, resulting in the prediction results often showing "noise" in the vertical direction that does not conform to geological laws. (3) Sequence-based deep learning methods: such as one-dimensional convolutional neural network, recurrent neural network and long short-term memory network. These models treat well logging data as a time series, consider the continuity of geology, and have significantly improved performance. However, they typically use fixed-size convolution kernels or windows for information aggregation. This fixed structure inevitably blurs key boundary features when dealing with lithological transition interfaces due to its smoothing effect. (4) Graph neural network-based methods: To model geological sequences more flexibly, one-dimensional well logging sequences are constructed as graph structures, and graph neural networks are used to capture geological context information. Among them, graph attention networks, through self-attention mechanisms, can dynamically assign different attention weights to adjacent strata, showing great potential in preserving geological interface features.

[0004] Despite the continuous advancements in modeling capabilities achieved by the aforementioned methods, they remain essentially data-driven "black box" models. Their core driving force is the statistical regularity within the data, lacking an understanding of the underlying physical laws. This limits the reliability and generalization ability of these models when facing complex geological problems. Summary of the Invention

[0005] The purpose of this invention is to provide a well logging lithology identification method based on physical information constraints, which aims to solve the problems mentioned in the background art.

[0006] The embodiments of the present invention are implemented as follows, such as...Figure 1 As shown, an embodiment of the present invention provides a well logging lithology identification method based on physical information constraints, comprising the following steps: Step 1: Data processing; Raw logging data and known lithological data were collected from multiple wells in the target area. The data were preprocessed to obtain one-dimensional logging sequence data. Step 2: Graph construction; The processed one-dimensional logging sequence data is transformed into a graph structure. Each logging depth sampling point is abstracted as a node, and its characteristics are composed of the multi-dimensional logging response value at that depth point. Edges are established between adjacent nodes at the depth to represent the vertical adjacency relationship of the formation.

[0007] Step 3: Construct a graph attention neural network; The constructed graph attention neural network includes a shared backbone network and two decoupled task branches. The shared backbone network is used to learn general geological context features, which are then fed into the two task branches, which are used for lithology determination and prediction of physical property parameters, respectively. Step 4: Construct a hybrid loss function; The hybrid loss function includes lithological classification loss, multi-scale physical constraint loss, and graph smoothing loss, and the model parameters are optimized by weighting and summing the loss terms through a dynamic weighting strategy. Step 5: Lithological prediction; The data of the new wells to be identified are input into the trained model, and the lithology of each depth point is predicted using the output of the classification branch to generate a lithology profile.

[0008] A further technical solution involves assigning an edge feature to each edge in step 2. This feature consists of the depth difference between two adjacent nodes and the difference in key logging values.

[0009] In a further technical solution, in step 3, the input graph data first passes through several shared graph attention network layers; the general features extracted by the graph attention network layers are then fed into two task branches; the classification branch is used for high-level feature extraction and finally outputs Logits for lithology identification; the physical information branch is used to predict physical property parameters related to physical constraints. Residual connections are introduced between the layers of a graph attention network, allowing the input features of the layer to be directly superimposed onto the output features.

[0010] A further technical solution involves, in step 4, lithological classification loss. The standard cross-entropy loss function is used, and the output of the classification branch is supervised using known lithology labels. The formula is as follows: ; in, Represents the total number of samples in the training set; The total number representing lithological categories; It is an indicator variable, if the sample The true category is the category ,but Otherwise, it is 0; The model predicts the sample. Category The probability is the result of the Logits output by the model being normalized by the Softmax function.

[0011] A further technical solution is that multi-scale physical constraint loss is a set of physical laws that are transformed into differentiable loss functions, which are applied to the output of the physical branch to guide the model to generate solutions that conform to physical laws. Density loss is based on the mass balance equation of rock physics: ; In the formula, Theoretical bulk density, The porosity predicted by the model. The model predicts the density of the rock skeleton. The density of the fluid in the formation; The difference between the theoretical and measured bulk density calculated from the porosity and skeleton density predicted by the loss penalty model is represented by the following loss function: ; in, To account for the overall density loss, For physical fitting density loss, To monitor density loss, This is due to geological prior density loss; in: ; In the formula, Mean square error, This is the measured bulk density; ; In the formula, For the actual measured density of the rock skeleton; ; In the formula, The GR curve value is from the well logging. Pearson correlation coefficient; Porosity-Permeability Relationship Loss: The Kozeny-Carman equation is used to describe the positive correlation between permeability and porosity; by constraining the relationship between permeability and porosity to conform to this physical trend, the loss function is: ; in, The loss is due to the porosity-permeability relationship. This is a vector of predicted penetration rates. Porosity predicted by the model The theoretical porosity term obtained after transformation by the KC equation, where: ; In the formula, For the model to the first Porosity prediction values ​​for each data point The first term of the theoretical porosity term One data point; Darcy's Law Loss: By minimizing the variance of the theoretical vertical velocity profile, the model is constrained to generate a more continuous permeability profile in the vertical direction that better conforms to the laws of fluid transport. ; in, For two adjacent depth points and The vertical fluid velocity calculated between them and Depth points and Predicted penetration rate and They are respectively and The known pressure, and They are respectively and The depth; The loss function is constructed as follows: ; In the formula, It is the variance function. Loss due to Darcy's Law This is the vertical fluid velocity profile vector; In a preferred embodiment of the present invention, the graph smoothing loss is a regularization term based on geological priors. By penalizing the permeability differences between adjacent nodes predicted by the model as belonging to the same lithology, it promotes smoother predictions of physical properties within the lithology. ; in, This represents all pairs of nodes that satisfy the condition. Take the expected value. To smooth the loss, For nodes Predicted lithology For nodes Predicted lithology Let be the edge set of the graph.

[0012] A further technical solution, in step 4, is to adopt a dynamic weighting strategy based on homoscedastic uncertainty to effectively balance the multiple heterogeneous loss terms: Learnable weight parameters: Introduce a learnable log-variance parameter for each loss term; Adaptive optimization: During model training, these parameters, along with the network weights, are optimized through backpropagation. The optimizer automatically adjusts each parameter to minimize the total loss, achieving an adaptive and dynamic balance for multi-task losses. The formula for calculating the total loss is as follows: ; in, It is a learnable logarithmic uncertainty parameter corresponding to each loss. yes Hyperparameters.

[0013] The present invention provides a well logging lithology identification method based on physical information constraints, the beneficial effects of which are as follows: (1) Significantly improve the recognition accuracy: Through the synergistic effect of data-driven and physical information, the model achieved an F1-score of 0.9528 in the prediction of four types of sandstone, which is significantly better than various pure data-driven models.

[0014] (2) Effectively suppresses overfitting and enhances generalization ability: thanks to the strong regularization effect of physical constraints, its validation set loss can continuously converge to a lower and more stable level, proving its stronger generalization ability and prediction stability on new data.

[0015] (3) Ensure the physical authenticity of the prediction results: This fundamentally solves the problem that traditional “black box” models may produce results that violate common sense in physics, and ensures the geological rationality of the output results.

[0016] (4) Achieving more refined geological boundary characterization: The introduction of the graph attention neural network architecture enables the model to process geological sequences more flexibly. Combined with graph smoothing loss, the model can more clearly characterize the contact interface between different lithologies while maintaining the smoothness of the internal physical properties of the lithology. Attached Figure Description

[0017] Figure 1A flowchart illustrating a well logging lithology identification method based on physical information constraints, provided as an embodiment of the present invention; Figure 2 This is the result of lithological prediction; Figure 3 This is the lithology identification confusion matrix of this method on the validation set; Figure 4 This is the predicted fitting result of this method for the rock density of a certain well section. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0019] The specific implementation of the present invention will be described in detail below with reference to specific embodiments.

[0020] like Figure 1 As shown, an embodiment of the present invention provides a well logging lithology identification method based on physical information constraints, comprising the following steps: Step 1: Data processing; First, raw logging data (such as natural gamma, resistivity, sonic transit time, etc.) and known lithological data are collected from multiple wells in the target area. Then, preprocessing operations such as cleaning, alignment, and standardization are performed on the data to obtain one-dimensional logging sequence data.

[0021] Step 2: Graph construction; The processed one-dimensional logging sequence data is transformed into a graph structure to explicitly model the spatial relationships between formations. Each logging depth sampling point is abstracted as a node, and its characteristics are composed of the multi-dimensional logging response values ​​at that depth point. Edges are established between adjacent nodes at depth to represent the vertical adjacency of formations.

[0022] Step 3: Construct a graph attention neural network;

[0023] The constructed graph attention neural network includes a shared backbone network and decoupled task branches. The shared backbone network is responsible for learning general geological context features, which are then fed into two independent branches: a classification branch for final lithology determination and a physical information branch for predicting physical properties such as porosity and permeability. Residual connections are also incorporated into the network to optimize information flow.

[0024] Step 4: Construct a hybrid loss function; The constructed hybrid loss function consists of four parts: Lithology classification loss: Classification loss supervised by lithology labels; Multi-scale physical constraint losses include density loss, porosity-permeability relationship loss, and Darcy's law loss. Graphical smoothing loss: constrains the smoothness of changes in physical properties within the same lithology; Dynamic loss weighting: Automatically balances the contribution of each loss term using learnable parameters.

[0025] The optimizer minimizes the total loss obtained by weighted summation of these parts.

[0026] Step 5: Lithological prediction; The data of the new wells to be identified are input into the trained model. The model uses the output of its classification branch to predict the lithology of each depth point through forward propagation, and finally generates a continuous, high-precision lithological profile that conforms to physical laws.

[0027] In a preferred embodiment of the present invention, in step 2, to enable the model to perceive the gradient of geological changes, each edge is assigned an edge feature, which consists of the depth difference between two adjacent nodes and the difference in key logging values. This allows the attention mechanism to not only focus on connectivity relationships but also perceive the intensity and direction of geological changes.

[0028] In a preferred embodiment of the present invention, in step 3, the input graph data is first passed through several shared graph attention network layers, which are designed to learn a universal rock physical response pattern and geological context features from the input well logging data and formation structure.

[0029] The general features extracted by the graph attention network layer are then fed into two parallel and independent task branches. The classification branch is used for high-level feature extraction, ultimately outputting Logits for lithology identification; the physical information branch is used to predict key physical property parameters related to physical constraints. This decoupled design simulates the qualitative-to-quantitative analysis process in geology and avoids optimization conflicts between different tasks.

[0030] Residual connections are introduced between layers of a graph attention network. By directly superimposing the input features of this layer onto the output features, it effectively solves the gradient vanishing problem in deep graph networks and ensures that the original low-order logging physical information is preserved in the deep network.

[0031] In a preferred embodiment of the present invention, in step 4, lithological classification loss The standard cross-entropy loss function is used, and the output of the classification branch is supervised by known lithology labels, which is the core data-driven term of the model. Its formula is: ; in, Represents the total number of samples in the training set; The total number representing lithological categories; It is an indicator variable, if the sample The true category is the category ,but Otherwise, it is 0; The model predicts the sample. Category The probability is the result of the Logits output by the model being normalized by the Softmax function.

[0032] As a preferred embodiment of the present invention, the multi-scale physical constraint loss is a set of physical laws that are transformed into differentiable loss functions, which are applied to the output of the physical branch to guide the model to generate solutions that conform to physical laws.

[0033] Density loss is based on the mass balance equation of rock physics: ; In the formula, Theoretical bulk density, The porosity predicted by the model. The model predicts the density of the rock skeleton. The density of the fluid in the formation.

[0034] The difference between the theoretical and measured bulk density calculated from the porosity and skeleton density predicted by the loss penalty model is represented by the following loss function: ; in, To account for the overall density loss, For physical fitting density loss, To monitor density loss, This is due to geological prior density loss; in: ; In the formula, Mean square error, This is the measured bulk density.

[0035] ; In the formula, This is the measured density of the rock skeleton.

[0036] ; In the formula, The GR curve value is from the well logging. This is the Pearson correlation coefficient.

[0037] Porosity-Permeability Relationship Loss: The Kozeny-Carman equation is used to describe the positive correlation between permeability and porosity; by constraining the relationship between permeability and porosity to conform to this physical trend, the loss function is: ; in, The loss is due to the porosity-permeability relationship. This is a vector of predicted penetration rates. Porosity predicted by the model The theoretical porosity term obtained after transformation by the KC equation, where: ; In the formula, For the model to the first Porosity prediction values ​​for each data point The first term of the theoretical porosity term Data points.

[0038] Darcy's Law Loss: Based on steady-state Darcy's law, fluid velocity should remain constant within a continuous geological unit. This is achieved by minimizing the variance of the theoretical vertical velocity profile to constrain the model to generate a more continuous permeability profile in the vertical direction that better reflects fluid transport laws. ; in, For two adjacent depth points and The vertical fluid velocity calculated between them and Depth points and Predicted penetration rate and They are respectively and The known pressure, and They are respectively and The depth.

[0039] The loss function is constructed as follows: ; In the formula, It is the variance function. Loss due to Darcy's Law This is the vertical fluid velocity profile vector.

[0040] In a preferred embodiment of the present invention, the graph smoothing loss is a regularization term based on geological priors that penalizes the permeability differences between adjacent nodes predicted by the model as belonging to the same lithology, thereby making the prediction of physical properties within the lithology smoother. ; in, This represents all pairs of nodes that satisfy the condition. Take the expected value. To smooth the loss, For nodes Predicted lithology For nodes Predicted lithology Let be the edge set of the graph.

[0041] In a preferred embodiment of the present invention, in step 4, a dynamic weighting strategy based on homoscedastic uncertainty is adopted to effectively balance the multiple heterogeneous loss terms: Learnable weight parameters: Introduce a learnable log-variance parameter for each loss term; Adaptive optimization: During model training, these parameters, along with the network weights, are optimized through backpropagation. The optimizer automatically adjusts each parameter to minimize the total loss, achieving an adaptive and dynamic balance for multi-task losses. The formula for calculating the total loss is as follows: ; in, It is a learnable logarithmic uncertainty parameter corresponding to each loss. yes Hyperparameters.

[0042] The crucial step occurs during the backpropagation phase. By applying the chain rule, the total loss is... The gradients of all learnable parameters in the model are automatically calculated, including the gradients of each... Parameters. For example, for a specific Its gradient is: ; Subsequently, in the parameter update step, the AdamW optimizer is used to update the values ​​of all parameters in a unified manner based on this gradient signal, including the values ​​of... The optimizer will be adjusted if the gradient is negative. The value of the gradient is used to reduce the weight of that loss term; conversely, if the gradient is positive, it will decrease. The value is used to increase the weight and reduce uncertainty.

[0043] like Figure 2The figure shows a comparison of lithological prediction results using this method in a typical well section. The left side of the figure shows the actual lithological profile of the well section, and the right side shows the lithological profile predicted by this method. The comparison demonstrates a high degree of consistency between the predicted and actual profiles in the distribution, thickness, and stratigraphic contact relationships of different lithologies, intuitively proving the effectiveness of this method in accurately depicting geological boundaries and ensuring prediction accuracy.

[0044] like Figure 3 The image shows the lithology identification confusion matrix of our method on the validation set. This confusion matrix visually demonstrates the model's classification performance on four different lithologies (SiS: silty sandstone, XS: cross-bedding sandstone, SaS: sandy siltstone, SST: sandstone). The diagonal of the matrix represents the proportion of correctly predicted samples, while the off-diagonal represents the proportion of misclassified samples.

[0045] like Figure 4 The figure shows the predicted fitting results of the proposed method for rock density in a certain well section. As can be seen from the figure, the density curve predicted by this method is highly consistent with the fluctuation trend of the actual value, accurately capturing subtle changes in formation properties with depth. This experimental result strongly demonstrates the density loss in the hybrid loss function constructed by this method, successfully integrating the rock physics mass balance equation as a strong constraint into the model training. This ensures that while the model is identifying lithology, its intermediate prediction process for physical parameters such as porosity and skeleton density also strictly adheres to physical laws, fundamentally guaranteeing the physical authenticity and geological reliability of the final lithology prediction results.

[0046] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A well logging lithology identification method based on physical information constraints, characterized in that, Includes the following steps: Step 1: Data processing; Raw logging data and known lithological data were collected from multiple wells in the target area. The data were preprocessed to obtain one-dimensional logging sequence data. Step 2: Graph construction; The processed one-dimensional logging sequence data is transformed into a graph structure. Each logging depth sampling point is abstracted as a node, and its characteristics are composed of the multi-dimensional logging response value of the logging depth sampling point. Edges are established between adjacent nodes in depth to represent the vertical adjacency relationship of the formation. Step 3: Construct a graph attention neural network; The constructed graph attention neural network includes a shared backbone network and two decoupled task branches. The shared backbone network is used to learn general geological context features, which are then fed into the two task branches, which are used for lithology determination and prediction of physical property parameters, respectively. Step 4: Construct a hybrid loss function; The hybrid loss function includes lithological classification loss, multi-scale physical constraint loss, and graph smoothing loss, and the model parameters are optimized by weighting and summing the loss terms through a dynamic weighting strategy. Step 5: Lithological prediction; The data of the new wells to be identified are input into the trained model, and the lithology of each depth point is predicted using the output of the classification branch to generate a lithology profile.

2. The well logging lithology identification method based on physical information constraints according to claim 1, characterized in that, In step 2, each edge is assigned an edge feature, which consists of the depth difference between two adjacent nodes and the difference in key logging values.

3. The well logging lithology identification method based on physical information constraints according to claim 2, characterized in that, In step 3, the input graph data first passes through several shared graph attention network layers; the general features extracted by the graph attention network layers are then fed into two task branches; the classification branch is used for high-level feature extraction and finally outputs Logits for lithology identification; the physical information branch is used to predict physical property parameters related to physical constraints. Residual connections are introduced between the layers of a graph attention network, allowing the input features of the layer to be directly superimposed onto the output features.

4. The well logging lithology identification method based on physical information constraints according to claim 3, characterized in that, In step 4, lithological classification loss The standard cross-entropy loss function is used, and the output of the classification branch is supervised using known lithology labels. The formula is as follows: ； in, Represents the total number of samples in the training set; The total number representing lithological categories; It is an indicator variable, if the sample The true category is the category ,but Otherwise, it is 0; The model predicts the sample. Category The probability is the result of the Logits output by the model being normalized by the Softmax function.

5. The well logging lithology identification method based on physical information constraints according to claim 4, characterized in that, Multiscale physical constraint loss is a set of physical laws that are transformed into differentiable loss functions. It is applied to the output of the physical branch and guides the model to generate solutions that conform to physical laws. Density loss is based on the mass balance equation of rock physics: ； In the formula, Theoretical bulk density, The porosity predicted by the model. The model predicts the density of the rock skeleton. The density of the fluid in the formation; The difference between the theoretical and measured bulk density calculated from the porosity and skeleton density predicted by the loss penalty model is represented by the following loss function: ； in, To account for the overall density loss, For physical fitting density loss, To monitor density loss, This is due to geological prior density loss; in: ； In the formula, Mean square error, This is the measured bulk density; ； In the formula, For the actual measured density of the rock skeleton; ； In the formula, The GR curve value is from the well logging. Pearson correlation coefficient; Porosity-Permeability Relationship Loss: The Kozeny-Carman equation is used to describe the positive correlation between permeability and porosity; by constraining the relationship between permeability and porosity to conform to this physical trend, the loss function is: ； in, The loss is due to the porosity-permeability relationship. This is a vector of predicted penetration rates. Porosity predicted by the model The theoretical porosity term obtained after transformation by the KC equation, where: ； In the formula, For the model to the first Porosity prediction values ​​for each data point The first term of the theoretical porosity term One data point; Darcy's Law Loss: By minimizing the variance of the theoretical vertical velocity profile, the model is constrained to generate a more continuous permeability profile in the vertical direction that better conforms to the laws of fluid transport. ； in, For two adjacent depth points and The vertical fluid velocity calculated between them and Depth points and Predicted penetration rate and They are respectively and The known pressure, and They are respectively and The depth; The loss function is constructed as follows: ； In the formula, It is the variance function. Loss due to Darcy's Law This is the vertical fluid velocity profile vector.

6. The well logging lithology identification method based on physical information constraints according to claim 5, characterized in that, The graph smoothing loss is a regularization term based on geological priors. It penalizes the permeability differences between adjacent nodes predicted by the model as belonging to the same lithology, thereby making the prediction of physical properties within the lithology smoother. ； in, This represents all pairs of nodes that satisfy the condition. Take the expected value. To smooth the loss, For nodes Predicted lithology, For nodes Predicted lithology, Let be the edge set of the graph.

7. The well logging lithology identification method based on physical information constraints according to claim 6, characterized in that, In step 4, to effectively balance the multiple heterogeneous loss terms, a dynamic weighting strategy based on homoscedastic uncertainty is adopted: Learnable weight parameters: Introduce a learnable log-variance parameter for each loss term; Adaptive optimization: During model training, these parameters are optimized together with the network weights through backpropagation; the optimizer automatically adjusts each parameter to minimize the total loss, achieving an adaptive and dynamic balance for multi-task losses; the formula for calculating the total loss is as follows: ； in, It is a learnable logarithmic uncertainty parameter corresponding to each loss. yes Hyperparameters.

Images