Method and device for dynamically predicting yield of multi-section fractured horizontal well

Abstract

The invention discloses a multi-section fractured horizontal well yield dynamic prediction method and device, and relates to the technical field of yield dynamic prediction.The method comprises the steps that feature extraction and normalization processing are conducted on multi-section crack information of a target well, and crack statistical features and crack spatial features are obtained through extraction; wherein spatial features of multi-segment crack information are extracted through a convolutional neural network in combination with an attention mechanism; performing feature splicing on the reservoir geological parameter information, the fracture statistical features and the fracture spatial features of the target well to obtain a static feature vector, and inputting the static feature vector into the trained yield dynamic prediction model to obtain a yield dynamic prediction result of the target well; according to the scheme, the yield data of the target well at a plurality of time points in the whole life cycle from production putting to final production stopping can be quickly and accurately obtained, and the problems that a traditional numerical simulation method is large in calculation resource consumption and high in dependence on historical production data can be well solved.

Description

Technical Field

[0001] The present application relates to the technical field of production dynamic prediction, and particularly to a method and device for predicting the production dynamics of a multi-stage fractured horizontal well. Background Art

[0002] Unconventional oil and gas reservoirs constitute an important energy resource base. In particular, the multi-stage fractured horizontal well technology has become a key technology for developing these oil and gas reservoirs due to its excellent stimulation effect and high economy. During the development process, accurately predicting the production dynamics of multi-stage fractured horizontal wells is crucial. This not only provides a scientific basis for well placement, fracturing parameter optimization, and production drainage system formulation, supporting the optimization of the development plan; at the same time, it also provides decision-making support for reserve dynamic adjustment, economic benefit evaluation, and production management. This helps to improve resource recovery rate, reduce development costs, and effectively control economic risks.

[0003] Currently, the methods for predicting the production dynamics of multi-stage fractured horizontal wells in unconventional oil and gas reservoirs mainly include empirical formula method, numerical simulation method, and machine learning method. The empirical formula method uses statistical regression and historical data to construct a simple model. Its advantages are fast calculation speed and easy application, which are suitable for preliminary estimation. However, its disadvantages are that it is difficult to accurately capture complex non-linear relationships, highly dependent on data, and its applicability is limited by the specific conditions at the time of formula development. Although the numerical simulation method has high prediction accuracy and wide applicability, it requires a large amount of computing resources and professional software support. In contrast, the application of the machine learning method in the field of unconventional oil and gas reservoirs is increasing. However, most current studies mainly predict future production changes based on historical production data and related influencing factors.

[0004] Therefore, solving the problems of large computing resource consumption and strong dependence on historical production data of the numerical simulation method is of great significance for realizing the rapid and accurate prediction of the production dynamics of multi-stage fractured horizontal wells in unconventional oil and gas reservoirs throughout the entire life cycle from production input to final production cessation. Summary of the Invention

[0005] The purpose of the present application is to provide a method and device for predicting the production dynamics of a multi-stage fractured horizontal well, which can quickly and accurately predict the production dynamic data of the target well.

[0006] To achieve the above purpose, the present application provides the following solutions: In the first aspect, the present application provides a method for predicting the production dynamics of a multi-stage fractured horizontal well, including: Obtaining the reservoir geological parameter information and multi-stage fracture information of the target well.

[0007] Extract the characteristics and normalize the multi-segment fracture information of the target well to obtain fracture statistical characteristics and fracture spatial characteristics; extract the spatial characteristics of the multi-segment fracture information through a convolutional neural network combined with an attention mechanism.

[0008] Concatenate the reservoir geological parameter information, the fracture statistical characteristics, and the fracture spatial characteristics of the target well to obtain a static feature vector.

[0009] Input the static feature vector into the trained production dynamic prediction model to obtain the production dynamic prediction result of the target well; the production dynamic prediction model is a model constructed based on the BiLSTM time series prediction model; the production dynamic prediction result is the production data at several time points during the entire life cycle of the target well from the start of production to the final shutdown.

[0010] Optionally, extracting the characteristics and normalizing the multi-segment fracture information of the target well to obtain fracture statistical characteristics and fracture spatial characteristics specifically includes: Parse each row of data of the multi-segment fracture information of the target well to obtain initial fracture statistical characteristics.

[0011] Convert the multi-segment fracture information of the target well into a picture form, and extract initial fracture spatial characteristics through a convolutional neural network combined with an attention mechanism.

[0012] Normalize the initial fracture statistical characteristics and the initial fracture spatial characteristics respectively to obtain fracture statistical characteristics and fracture spatial characteristics.

[0013] Optionally, it further includes: constructing an initial production dynamic prediction model, and training it according to the reservoir geological parameters and multi-segment fracture information of several well groups, and optimizing the initial production dynamic prediction model to obtain a trained production dynamic prediction model; specifically including: Construct an initial production dynamic prediction model based on the BiLSTM time series prediction model.

[0014] Obtain the reservoir geological parameters and multi-segment fracture information of several historical well groups.

[0015] For any historical well group, perform reservoir numerical simulation according to the reservoir geological parameters and multi-segment fracture information of the historical well group to obtain the simulated production dynamic data of the historical well group.

[0016] Extract the characteristics and normalize the multi-segment fracture information of the historical well group to obtain the fracture statistical characteristics and fracture spatial characteristics of the historical well group.

[0017] Taking the reservoir geological parameters, fracture statistical characteristics, and fracture spatial characteristics of historical well groups as the input of the initial production dynamic prediction model, and taking the simulated production dynamic data of historical well groups as the target output of the initial production dynamic prediction model, training and optimizing the initial production dynamic prediction model to obtain a trained production dynamic prediction model.

[0018] Optionally, performing reservoir numerical simulation based on the reservoir geological parameters and multi-stage fracture information of the historical well groups to obtain the simulated production dynamic data of the historical well groups, specifically including: For any historical well group, designing a reservoir geological parameter plan based on the reservoir geological parameters of the historical well group, and designing a fracture design plan based on the multi-stage fracture information of the historical well group.

[0019] Establishing a regular grid-based numerical simulation model based on the multi-stage fracture information of the historical well group.

[0020] According to the reservoir geological parameter plan and fracture design plan of the historical well group, calling the numerical simulation model to perform reservoir simulation to obtain the simulated production dynamic data of the historical well group.

[0021] Optionally, it further includes: keeping the reservoir geological parameters, fracture statistical characteristics, fracture spatial characteristics, and simulated production dynamic data of several historical well groups corresponding, and dividing them into a training set and a test set; using the data in the training set to train and optimize the initial production dynamic prediction model, and using the data in the test set to test the performance of the trained initial production dynamic prediction model.

[0022] Optionally, the production dynamic prediction model includes two layers of bidirectional LSTM. The number of hidden units in the first layer of bidirectional LSTM is 208, the number of hidden units in the second layer of bidirectional LSTM is 176, the learning rate is 0.0004, and the batch_size is 12.

[0023] In a second aspect, the present application provides a multi-stage fractured horizontal well production dynamic prediction device, including: A target well parameter acquisition module, configured to acquire the reservoir geological parameter information and multi-stage fracture information of a target well.

[0024] A fracture parameter feature extraction module, configured to perform feature extraction and normalization processing on the multi-stage fracture information of the target well, and extract fracture statistical characteristics and fracture spatial characteristics; fully extract the spatial characteristics of the multi-stage fracture information through a convolutional neural network combined with an attention mechanism.

[0025] A static feature vector splicing module, configured to splice the reservoir geological parameter information, the fracture statistical characteristics, and the fracture spatial characteristics of the target well to obtain a static feature vector.

[0026] The production dynamic prediction module is used to input the static feature vector into the trained production dynamic prediction model to obtain the production dynamic prediction result of the target well; the production dynamic prediction model is a model constructed based on the BiLSTM time series prediction model; the production dynamic prediction result is the production data at several time points during the entire life cycle of the target well from the start of production to the final shutdown.

[0027] According to the specific embodiments provided by this application, the following technical effects are disclosed in this application: This application provides a method and device for predicting the production dynamics of a multi-stage fractured horizontal well. In this method, feature extraction and normalization processing are performed on the multi-stage fracture information of the target well, and fracture statistical features and fracture spatial features are extracted. Among them, the spatial features of the multi-stage fracture information are extracted through a convolutional neural network combined with an attention mechanism, and the fracture distribution and geometric characteristics are represented by the spatial features and statistical features of the multi-stage fracture information, so as to solve the problem that the data lengths of the multi-stage fracture information of different horizontal wells are inconsistent due to different numbers of fractures; and the reservoir geological parameter information, fracture statistical features, and fracture spatial features of the target well are feature-stitched to obtain a static feature vector and input it into the trained production dynamic prediction model to obtain the production dynamic prediction result of the target well; the production dynamic prediction model used is a model constructed based on the BiLSTM time series prediction model, which can quickly and accurately obtain the production data at several time points during the entire life cycle of the target well from the start of production to the final shutdown, and can well overcome the problems of large consumption of computing resources and strong dependence on historical production data in the traditional numerical simulation method, which is of great significance for realizing the rapid and accurate prediction of the production dynamics of a multi-stage fractured horizontal well during the entire life cycle from the start of production to the final shutdown in unconventional oil and gas reservoirs. Description of the Drawings

[0028] In order to more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the following will briefly introduce the drawings required to be used in the embodiments. Obviously, the drawings in the following description are only some embodiments of this application. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0029] Figure 1 It is a flowchart of a method for predicting the production dynamics of a multi-stage fractured horizontal well provided by an embodiment of this application.

[0030] Figure 2 It is a schematic diagram of the multi-stage fracture design of a typical horizontal well in a method for predicting the production dynamics of a multi-stage fractured horizontal well provided by an embodiment of this application.

[0031] Figure 3It is a flowchart of step A2 in a method for dynamically predicting the production of a multi-stage fractured horizontal well provided by an embodiment of the present application.

[0032] Figure 4 It is a schematic structural diagram of a traditional LSTM.

[0033] Figure 5 It is a schematic structural diagram of the BiLSTM adopted by the present application.

[0034] Figure 6 It is a schematic diagram of the visualization of the production dynamics of multiple horizontal wells in a method for dynamically predicting the production of a multi-stage fractured horizontal well provided by an embodiment of the present application.

[0035] Figure 7 It is a technical roadmap of a method for dynamically predicting the production of a multi-stage fractured horizontal well provided by an embodiment of the present application.

[0036] Figure 8 It is a schematic diagram of the functional modules of a device for dynamically predicting the production of a multi-stage fractured horizontal well provided by an embodiment of the present application.

[0037] Figure 9 It is a schematic structural diagram of a computer device provided by an embodiment of the present application. Detailed implementation manners

[0038] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present application.

[0039] To make the above objects, features, and advantages of the present application more obvious and understandable, the present application will be further described in detail below in conjunction with the drawings and specific implementation manners.

[0040] A method for dynamically predicting the production of a multi-stage fractured horizontal well provided by an embodiment of the present application, in an exemplary embodiment, as Figure 1 shown, includes the following steps: A1. Obtain the reservoir geological parameter information and multi-stage fracture information of the target well. Specifically, in this embodiment, the reservoir geological parameter information includes key parameters such as porosity, permeability, thickness, Young's modulus, relative permeability, reservoir pressure, ratio of desorption pressure to reservoir pressure, Langmuir volume, Langmuir pressure, fracture permeability, well control area, etc. Subsequently, according to the reservoir geological parameter information of different types of unconventional oil and gas reservoirs (assuming predicting the production dynamics of a coalbed methane reservoir, only the reservoir geological parameters related to the coalbed methane reservoir can be used, and they cannot be confused with the reservoir geological parameters of other unconventional oil and gas reservoirs), and combined with the number of wells participating in the training, a reasonable reservoir geological parameter scheme is formulated to ensure the adaptability and accuracy of the prediction model.

[0041] The multi-stage fracture information includes the ratio of the well control area to the well-controlled reservoir (the well control area is also called the effective area), the number of fractures, the fracture spacing, and the fracture length. The specific design range is adjusted according to the reservoir type and geological conditions to optimize the fracture distribution and conductivity. For example, the number of fractures is set according to the reservoir permeability and stress field characteristics, the fracture spacing is flexibly adjusted in combination with the reservoir stimulation requirements, and the fracture length is restricted by the width of the well control area (the effective area represents its length along the horizontal direction and its width along the vertical direction; the fracture stimulation area calculates the fracture spacing along the horizontal direction and the fracture length along the vertical direction) to ensure the effectiveness of the fracturing stimulation. Specifically, the Latin hypercube sampling method can be used to design the positions of multi-stage fractures in horizontal wells to form fracture design schemes applicable to different types of unconventional oil and gas reservoirs. A schematic diagram of the multi-stage fracture design of a typical horizontal well is shown as Figure 2 shown.

[0042] A2. Extract features and perform normalization processing on the multi-stage fracture information of the target well, and extract fracture statistical features and fracture spatial features; extract the spatial features of the multi-stage fracture information through a convolutional neural network combined with an attention mechanism. In this embodiment, as shown in Figure 3 shown, step A2 specifically includes the following steps: A21. Parse each row of data of the multi-stage fracture information of the target well to extract the initial fracture statistical features.

[0043] A22. Convert the multi-stage fracture information of the target well into a picture form, and extract the initial fracture spatial features through a convolutional neural network combined with an attention mechanism.

[0044] A23. Perform normalization processing on the initial fracture statistical features and the initial fracture spatial features respectively to obtain the fracture statistical features and the fracture spatial features.

[0045] A3. Perform feature splicing on the reservoir geological parameter information, fracture statistical features, and fracture spatial features of the target well to obtain a static feature vector.

[0046] A4. Input the static feature vector into the trained production dynamic prediction model to obtain the production dynamic prediction result of the target well; the production dynamic prediction model is a model constructed based on the BiLSTM time series prediction model; the production dynamic prediction result is the production data at several time points during the entire life cycle of the target well from the start of production to the final shutdown.

[0047] In another exemplary embodiment of the present application, the production dynamic prediction method for multi-stage fractured horizontal wells further includes: B1. Construct an initial production dynamic prediction model, and train it according to the reservoir geological parameters and multi-stage fracture information of several well groups, and optimize the initial production dynamic prediction model to obtain the trained production dynamic prediction model; specifically including: B11. Based on the BiLSTM time series prediction model, construct an initial production dynamic prediction model. Since the predicted production dynamics have time series characteristics, BiLSTM is selected to predict the production dynamics of unconventional oil and gas. At the same time, compared with general time series models such as LSTM (the structural schematic diagram of LSTM is as Figure 4 shown), BiLSTM introduces a bidirectional structure in the network, has stronger time series feature extraction ability, and is more suitable for long-term prediction. The BiLSTM structure schematic diagram adopted in the present application is as Figure 5 shown.

[0048] B12. Obtain the reservoir geological parameters and multi-stage fracture information of several historical well groups. In this embodiment, 1000 sets of reservoir geological parameters and multi-stage fracture information are collected to form a reservoir geological parameter plan and a fracture design plan respectively.

[0049] B13. For any historical well group, perform reservoir numerical simulation according to the reservoir geological parameters and multi-stage fracture information of the historical well group to obtain the simulated production dynamic data of the historical well group. After obtaining the reservoir geological parameter plan and the fracture design plan, use numerical simulation software MATLAB and CMG to simulate the production process of the multi-stage fractured horizontal well for 20 years, and finally extract the production dynamic data including daily gas production and cumulative gas production. For each horizontal well, with a step size of 1 month, a total of 240 time steps of corresponding production data are obtained. In this embodiment, step B13 specifically includes: B131. For any historical well group, design a reservoir geological parameter plan according to the reservoir geological parameters of the historical well group, and design a fracture design plan according to the multi-stage fracture information of the historical well group.

[0050] B132. Establish a numerical simulation model based on a regular grid using the multi-segment fracture information of the historical well group. In this embodiment, a numerical simulation model based on a regular grid is established according to the multi-segment fracture information of the historical well group, aiming to improve the calculation accuracy of the numerical simulation process and reasonably characterize the fracture distribution. The model is divided into several grid units in both the horizontal and vertical directions, and the size of each grid unit is set according to the actual geological conditions to balance the calculation efficiency and simulation accuracy. This grid division scheme ensures the accurate characterization of multi-segment fractures and provides a reliable data basis for subsequent production performance prediction.

[0051] B133. According to the reservoir geological parameter scheme and fracture design scheme of the historical well group, call the numerical simulation model to perform reservoir simulation and obtain the simulated production performance data of the historical well group. Package the reservoir geological parameter scheme and fracture design scheme into files (such as CSV or MAT files) respectively, and use the numerical calculation tool MATLAB and the reservoir numerical simulation software CMG to perform reservoir simulation; for example: First, write a script file for performing reservoir simulation using MATLAB. Read the reservoir geological parameter scheme and fracture design scheme in this script file, and then implement the numerical simulation model based on a regular grid in this script file. Subsequently, call the GEM tool in CMG in the script file to batch implement the reservoir simulation process. Finally, execute the script file, extract the required production performance data (such as gas production, oil production, etc.) from the generated simulation results, and save them in an Excel file.

[0052] B14. Extract features and perform normalization processing on the multi-segment fracture information of the historical well group to obtain the fracture statistical features and fracture spatial features of the historical well group.

[0053] After step B13 processing, the shape of the reservoir geological parameters is (1000, 15), indicating 1000 samples, and each sample includes 15 geological parameters. The shape of the statistical features extracted from the multi-segment fractures is (1000, 6), and each sample includes 6 kinds of information: the proportion of the effective fracture area, the number of fractures, the mean and variance of the fracture length, and the mean and variance of the fracture spacing. Then convert the multi-segment fracture information into a picture with a shape of 128 * 256, input it into the CNN-AM model for spatial feature extraction, and the shape of the spatial features after screening is (1000, 20). Finally, splice the reservoir geological parameters, the statistical features and spatial features of the multi-segment fracture information together to form a static feature vector with a shape of (1000, 41).

[0054] For reservoir geological parameters, whether to perform data normalization should be determined in combination with the distribution of the actually collected reservoir geological parameters. According to observations, the values of all the collected data are between 0 and 1, so data normalization is not required. At the same time, through detection, it is found that there are no missing values and outliers, and the lengths of the reservoir geological parameters of each horizontal well are the same. For multi-segment fracture information, due to the fact that the number of fractures in each horizontal well is not exactly the same, the data lengths of their fracture information are inconsistent.

[0055] Specifically, for multi-segment fracture information with inconsistent input data lengths, by parsing each row of the input data, statistical features such as the proportion of the well-controlled area, the number of fractures, the length, and the spacing are extracted. First, calculate the ratio of the well-controlled area to the well-controlled reservoir according to the boundary values of the given area; then, parse the fracture center position and half-length information, calculate the length of each fracture, the spacing between adjacent fractures, and further calculate their mean and variance; finally, summarize the calculation results and perform min-max normalization on the eigenvalue of each column to standardize the data for subsequent analysis.

[0056] Fracture statistical features include the proportion of the well-controlled area, the number of fractures, the mean and variance of the fracture length, and the mean and variance of the fracture spacing. The proportion of the well-controlled area is the ratio of the well-controlled area to the well-controlled reservoir; the number of fractures will be processed by min-max normalization, as shown in formula (1): (1).

[0057] Among them, x is the number of fractures in a certain horizontal well, min is the minimum value of the number of fractures among all the horizontal wells, max is the maximum value of the number of fractures, x' is the number of fractures after normalization, and the result is between [0, 1].

[0058] The calculation of the mean and variance of the fracture length and spacing is shown in formulas (2) to (5): (2).

[0059] (3).

[0060] (4).

[0061] (5).

[0062] Among them, μ L is the mean of the fracture length in a certain horizontal well, is the variance of the fracture length in this well, N Lis the number of fractures in the well, x L,i is the length of one of the fractures in the well, μ D is the mean of the spacings in a horizontal well, is the variance of the spacings in the well, N D is the number of spacings in the well, x D,i is the spacing between two adjacent fractures in the well.

[0063] After that, convert the multi - fracture information into a picture format, and use CNN and AM to extract its spatial features. CNN and AM are mainly used to extract the spatial features of multi - fracture information. Using spatial and channel attention mechanisms can enhance the attention to fractures. Through the convolution operation of CNN, gradually extract the local and global spatial features of fractures. At the same time, combine the channel attention and spatial attention mechanisms, which respectively focus on important channels and significant spatial regions in the feature map, highlight key features and suppress redundant information, so as to achieve more efficient and accurate feature extraction of fracture distribution, morphology and characteristics. After extracting the spatial features of multi - fractures, it is necessary to screen the blank column data caused by the sparse distribution of fractures, so as to reduce the data dimension to improve the accuracy of subsequent prediction. Finally, also perform maximum - minimum normalization processing on the eigenvalue of each column of spatial features.

[0064] Specifically in this embodiment, the constructed CNN - AM model architecture includes an input layer, a convolution and pooling layer, a global average pooling layer, a fully - connected layer and an output layer. Among them, the input layer receives two - dimensional input data with a shape of (128, 256, 1). The convolution and pooling layer contains 3 layers: The first - layer convolution extracts low - level features, uses 32 3x3 convolution kernels, and the activation function is ReLU. Subsequently, reduce the feature map resolution through 2x2 max - pooling; The second - layer convolution uses 64 3x3 convolution kernels and introduces a spatial attention module to enhance the model's focus on key regions. Subsequently, reduce the dimension again through 2x2 max - pooling. The third - layer convolution uses 128 3x3 convolution kernels and introduces a channel attention module to enhance the feature selection ability. Global average pooling layer: Convert the feature map into a fixed - size feature vector to reduce the number of parameters. Fully - connected layer: Generate the final feature vector, with the dimension of num_features = 32, and the activation function is ReLU. The output layer outputs a 32 - dimensional feature vector.

[0065] For the simulated production dynamic data, through such as Figure 6It is found from the visualization diagram of the production performance of the multi-lateral horizontal wells shown that there are significant differences in the cumulative production (cumulative gas production or cumulative oil production) values among different horizontal wells. Therefore, during the training of the model, the min-max normalization process is also performed on it to enhance the stability of the model and improve the model training efficiency.

[0066] As an optional implementation manner, to ensure a better production performance prediction model for training, the method further includes: keeping the reservoir geological parameters, fracture statistical characteristics, fracture spatial characteristics, and simulated production performance data of several historical well groups corresponding, and dividing them into a training set and a test set; the training set is used for the training of the model, and the test set is used to verify the prediction performance of the model, and the repeatability of data division is ensured by fixing the random seed.

[0067] B15. Use the reservoir geological parameters, fracture statistical characteristics, and fracture spatial characteristics of the historical well group as the input of the initial production performance prediction model, and use the simulated production performance data of the historical well group as the target output of the initial production performance prediction model, and perform training and optimization on the initial production performance prediction model to obtain a trained production performance prediction model.

[0068] The above model training method of this application is implemented based on TensorFlow in Python. The static feature vectors and production performance data input into the BiLSTM model need to be reshaped into three-dimensional tensors with a shape of (number of samples, time steps, feature dimension) to adapt to the input requirements of the LSTM before the model training can be carried out. The BiLSTM model architecture constructed in this embodiment mainly consists of three parts: First, the first layer of the model is a bidirectional LSTM (Bidirectional LSTM), which consists of a forward and a backward LSTM network. Each LSTM unit contains 208 hidden units and uses the ReLU activation function to capture the non-linear features in the input data. By setting return_sequences=True, this layer outputs the hidden state at each time step to provide complete time series information for the subsequent LSTM layer. The second layer is another bidirectional LSTM, which contains 176 hidden units but does not return the sequence (return_sequences=False). It extracts global time features and further compresses the time dimension information. Finally, the model outputs the time step prediction values of the target sequence through a fully connected layer (Dense), and the output dimension is the same as that of the target sequence.

[0069] The BiLSTM model can comprehensively extract context information by capturing both forward and backward temporal dependencies in time series data. It combines the gate mechanism of LSTM to effectively solve the problems of gradient vanishing and gradient explosion in ordinary recurrent neural networks. At the same time, it extracts historical and future information through two LSTM layers in the forward and backward directions respectively, and then fuses them to generate a more accurate feature representation.

[0070] The calculation of the forward LSTM is the same as that of the standard LSTM, starting from the first time step of the sequence. The following calculation formulas (6) - (11) are, from top to bottom, the forget gate, input gate, candidate memory cell, memory cell update, output gate, and hidden state respectively.

[0071] (6).

[0072] (7).

[0073] (8).

[0074] (9).

[0075] (10).

[0076] (11).

[0077] Among them, the superscript ( f ) represents the calculation of the forward LSTM, x t is the cell input vector, , , represent the output of the forget gate, input gate, and output gate of the forward LSTM at time t respectively; , , represent the memory cell state, candidate memory cell, and hidden state of the forward LSTM at time t respectively; , , , represent the weight matrices of the forget gate, input gate, memory cell, and output gate of the forward LSTM respectively; , , , represent the bias terms of the forget gate, input gate, memory cell, and output gate of the forward LSTM respectively; is the Sigmoid activation function, ranging from 0 to 1; tanh is the hyperbolic tangent activation function, ranging from -1 to 1.

[0078] The calculation of the reverse LSTM starts from the last time step of the sequence and performs backpropagation. The calculation formulas are as shown in Eqs. (12) to (18): (12).

[0079] (13).

[0080] (14).

[0081] (15).

[0082] (16).

[0083] (17).

[0084] Among them, the superscript ( b ) represents the calculation of the reverse LSTM, and the meanings of the other symbols are the same as those of the forward LSTM.

[0085] The final output of the BiLSTM is the concatenation (or weighted average) of the hidden states of the forward and reverse LSTMs. Assuming the outputs of the forward and reverse LSTMs are and , then the final output is:[[]] (18).

[0086] The model is trained using the Adam optimizer (setting the learning rate) and the mean squared error (MSE) as the loss function to compile the model. In the training stage, the model is trained with a batch size of 8, iterated at most 300 times, and 10% of the training data is used as the validation set. To prevent overfitting, an early stopping mechanism is set, and the training is terminated early when the validation loss does not improve for 20 consecutive epochs.

[0087] The model optimization process used Keras Tuner for hyperparameter tuning. First, a hyperparameter model construction function was defined, which included the number of BiLSTM layers, the number of units in each BiLSTM layer, and the hyperparameter range of the learning rate. Hyperparameter search was performed using the random search method, with the maximum number of trials set to 20, and each trial was executed twice to find the best hyperparameter combination. During the search process, an early stopping strategy was used to prevent overfitting, and the model was evaluated on the validation set. Finally, the final model was constructed and trained using the best hyperparameters, and its performance was evaluated on the test set. In this embodiment, the production dynamic prediction model contains two layers of bidirectional LSTM. The number of hidden units in the first layer of bidirectional LSTM is 208, the number of hidden units in the second layer of bidirectional LSTM is 176, the learning rate is 0.0004, and the batch_size is 12.

[0088] The prediction and evaluation process of the model first makes predictions on the training set and the test set through the final model. Then, a custom evaluation function is used to calculate and return the model evaluation metrics for the training set and the test set, including the correlation coefficient ( R 2 ), mean squared error ( MSE ), root mean squared error ( RMSE ), and mean absolute error ( MAE ), as shown in formulas (19) - (22) respectively. Finally, the evaluation results of the training set and the test set are output, and some prediction images are displayed to help further analyze the model performance.

[0089] (19).

[0090] (20).

[0091] (21).

[0092] (22).

[0093] Among them, y i is the true value in the production data, is the predicted value in the production data, n is the number of time steps of the production dynamics in the multi-stage fractured horizontal well. Among these evaluation metrics, R 2 the closer the value of MSE , RMSE and MAE is to 1, and the smaller the values of Figure 7As shown, it includes five steps: data acquisition, data preprocessing, data preparation and partitioning, model construction and training, and model optimization and prediction. After the above training and optimization process, the prediction effects of the trained model on the training set and the test set are shown in Table 1.

[0094] Table 1 Prediction Result Evaluation

[0095] The above experimental cases prove that the above method proposed in the above embodiments of the present application performs well in predicting the production dynamics of multi-stage fractured horizontal wells in coalbed methane reservoirs using static reservoir geological parameters and multi-stage fracture information. It not only has high accuracy but also has a faster prediction speed compared with empirical formulas and numerical simulation methods.

[0096] Based on the same inventive concept, the embodiments of the present application also provide a device for implementing the above-mentioned method for predicting the production dynamics of multi-stage fractured horizontal wells. The solution provided by the device to solve the problem is similar to the solution described in the above method. In an exemplary embodiment, as Figure 8 shown, a device for predicting the production dynamics of multi-stage fractured horizontal wells is provided, including: A target well parameter acquisition module for acquiring reservoir geological parameter information and multi-stage fracture information of the target well.

[0097] A fracture parameter feature extraction module for performing feature extraction and normalization processing on the multi-stage fracture information of the target well, and extracting fracture statistical features and fracture spatial features; fully extracting the spatial features of the multi-stage fracture information through a convolutional neural network combined with an attention mechanism.

[0098] A static feature vector splicing module for splicing the reservoir geological parameter information, the fracture statistical features, and the fracture spatial features of the target well to obtain a static feature vector.

[0099] A production dynamics prediction module for inputting the static feature vector into a trained production dynamics prediction model to obtain a production dynamics prediction result of the target well; the production dynamics prediction model is a model constructed based on a BiLSTM time series prediction model; the production dynamics prediction result is the production data of the target well at several time points during the entire life cycle from the start of production to the final shutdown.

[0100] Of course, Figure 8 the architecture shown is only exemplary. When implementing different functions, one or at least two components in the device shown may be omitted according to actual needs. Figure 8 above In an exemplary embodiment, a computer device is provided. The computer device can be a server or a terminal, and its internal structure diagram can be as shown in Figure 9 . The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O), and a communication interface. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals through a network connection. When the computer program is executed by the processor, the multi-stage fractured horizontal well production dynamic prediction method provided in the above embodiment can be implemented.

[0101] Those skilled in the art can understand that Figure 9 the structure shown in is only a block diagram of some structures related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have a different component layout.

[0102] In an exemplary embodiment, a computer device is further provided, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

[0103] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by the processor, the steps in the above method embodiments are implemented.

[0104] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by the processor, the steps in the above method embodiments are implemented.

[0105] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or fully authorized by all parties, and the collection, use, and processing of relevant data need to comply with relevant regulations.

[0106] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, database, or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memories. Non-volatile memories can include read-only memory (ROM), magnetic tapes, floppy disks, flash memories, optical memories, high-density embedded non-volatile memories, resistive random access memories (ReRAM), magnetoresistive random access memories (MRAM), ferroelectric random access memories (FRAM), phase change memories (PCM), graphene memories, etc. Volatile memories can include random access memory (RAM) or external cache memories, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0107] The databases involved in the embodiments provided in the present application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., without limitation. The processors involved in the embodiments provided in the present application can be general-purpose processors, central processors, graphics processors, digital signal processors, programmable logics, data processing logics based on quantum computing, etc., without limitation.

[0108] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope described in this specification.

[0109] Specific examples are used in this article to elaborate on the principles and implementation manners of the present application. The descriptions of the above embodiments are only used to help understand the method and its core idea of the present application; at the same time, for those of ordinary skill in the art, according to the idea of the present application, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present application.

Claims

1. A method for dynamically predicting the production of a multi-stage fractured horizontal well, characterized in that: include: Obtain reservoir geological parameter information and multi-segment fracture information of the target well; Performing feature extraction and normalization processing on the multiple fracture information of the target well to extract fracture statistical features and fracture spatial features; Extracting spatial features of the multiple crack information by using a convolutional neural network combined with an attention mechanism; Performing feature splicing on the reservoir geological parameter information, the fracture statistical characteristics and the fracture spatial characteristics of the target well to obtain a static feature vector; The static feature vector is input into the trained dynamic production prediction model to obtain the dynamic production prediction result of the target well; the dynamic production prediction model is a model built based on the BiLSTM time series prediction model; the dynamic production prediction result is the production data of the target well at several time points in the entire life cycle from the start of production to the final shutdown.

2. The method for dynamic prediction of production of multi-stage fractured horizontal wells according to claim 1, characterized in that: Feature extraction and normalization processing are performed on the multiple fracture information of the target well to extract fracture statistical features and fracture spatial features, specifically including: Analyze each row of data of the multiple sections of fracture information of the target well to extract initial fracture statistical features; The multiple fracture information of the target well is converted into a picture form, and the initial fracture space features are extracted by combining a convolutional neural network with an attention mechanism; The initial fracture statistical characteristics and the initial fracture spatial characteristics are respectively normalized to obtain fracture statistical characteristics and fracture spatial characteristics.

3. The method for dynamic prediction of production of multi-stage fractured horizontal wells according to claim 1, characterized in that: Also includes: Constructing an initial production dynamic prediction model, and training it according to reservoir geological parameters and multi-segment fracture information of several well groups, optimizing the initial production dynamic prediction model, and obtaining a trained production dynamic prediction model; specifically including: Based on the BiLSTM time series prediction model, the initial production dynamic prediction model is constructed; Obtain reservoir geological parameters and multi-segment fracture information of several historical well groups; For any historical well group, numerical simulation of the reservoir is performed according to the reservoir geological parameters and multi-segment fracture information of the historical well group to obtain simulated production dynamic data of the historical well group; Performing feature extraction and normalization processing on the multi-segment fracture information of the historical well group, and extracting fracture statistical features and fracture spatial features of the historical well group; The reservoir geological parameters, fracture statistical characteristics and fracture spatial characteristics of the historical well group are used as the input of the initial production dynamic prediction model, and the simulated production dynamic data of the historical well group is used as the target output of the initial production dynamic prediction model. The initial production dynamic prediction model is trained and optimized to obtain a trained production dynamic prediction model.

4. The method for dynamic prediction of production of multi-stage fractured horizontal wells according to claim 3, characterized in that: According to the reservoir geological parameters and multi-segment fracture information of the historical well group, a numerical simulation of the oil reservoir is performed to obtain the simulated production dynamic data of the historical well group, specifically including: For any historical well group, a reservoir geological parameter scheme is obtained according to the reservoir geological parameter design of the historical well group, and a fracture design scheme is obtained according to the multi-segment fracture information design of the historical well group; Establishing a numerical simulation model based on a regular grid based on the multi-segment fracture information of the historical well group; According to the reservoir geological parameter scheme and fracture design scheme of the historical well group, the numerical simulation model is called to perform reservoir simulation to obtain simulated production dynamic data of the historical well group.

5. The method for dynamic prediction of production of multi-stage fractured horizontal wells according to claim 3, characterized in that: Also includes: The reservoir geological parameters, fracture statistical characteristics, fracture spatial characteristics and simulated production dynamic data of several historical well groups are kept corresponding and divided into a training set and a test set; the initial production dynamic prediction model is trained and optimized using the data in the training set, and the performance of the trained initial production dynamic prediction model is tested using the data in the test set.

6. The method for dynamic prediction of production of multi-stage fractured horizontal wells according to claim 1, characterized in that: The dynamic output prediction model includes two layers of bidirectional LSTM, the number of hidden units of the first layer of bidirectional LSTM is 208, the number of hidden units of the second layer of bidirectional LSTM is 176, the learning rate is 0.0004, and the batch_size is 12.

7. A dynamic prediction device for production of multi-stage fracturing horizontal wells, characterized in that: include: A target well parameter acquisition module is used to obtain reservoir geological parameter information and multi-segment fracture information of the target well; A fracture parameter feature extraction module is used to extract and normalize the multiple fracture information of the target well to obtain fracture statistical features and fracture spatial features; a convolutional neural network is combined with an attention mechanism to fully extract the spatial features of the multiple fracture information; A static feature vector splicing module is used to perform feature splicing on the reservoir geological parameter information, the fracture statistical characteristics and the fracture spatial characteristics of the target well to obtain a static feature vector; The production dynamic prediction module is used to input the static feature vector into the trained production dynamic prediction model to obtain the production dynamic prediction result of the target well; the production dynamic prediction model is a model built based on the BiLSTM time series prediction model; the production dynamic prediction result is the production data of the target well at several time points in the entire life cycle from the start of production to the final shutdown.

Images