An automatic logging operation system for plateau areas

Abstract

The present invention relates to the field of automated logging technology, specifically to an automated logging operation system in plateau areas, the system comprising: a rock formation characteristic analysis module performs pixel point detection and regional labeling based on rock formation texture and fault line image data in plateau areas, performs feature extraction and node association analysis on each node, and establishes a structured geological map. In the present invention, by accurately analyzing rock formation texture and fault line data in plateau areas, the detailed and accurate geological data are ensured, and these high-quality data are used for real-time feedback and parameter adjustment, effectively reducing the drilling risk caused by rapid changes in geological conditions, dynamically adjusting drilling parameters to reduce the possibility of equipment failure and operation interruption, optimizing the continuity and safety of drilling operations, and enhancing the creation of geological maps. By improving the relationship between data points, the data has more predictive value and operational guidance significance, thereby improving the accuracy of decision support and the overall efficiency of oilfield development.

Description

Technical Field

[0001] The invention relates to the technical field of automated well logging, and in particular to an automated well logging operation system in plateau areas. Background Art

[0002] The field of automated logging technology involves the use of highly automated equipment and systems to perform logging tasks in drilling operations in underground oil and gas fields. Logging refers to the measurement of geological, geophysical, and geochemical characteristics using specialized tools and equipment during the drilling process. Automation technology makes the logging process more efficient and accurate, and can analyze underground structures in real time, providing key data to support oilfield development decisions. This includes automatic control of the operation of drilling equipment, real-time data collection, processing and interpretation, and instant adjustments based on data feedback. Automated logging systems typically integrate advanced software and hardware, such as artificial intelligence, machine learning, and remote operation technologies to improve operational safety and efficiency.

[0003] Among them, the automation system for logging operations in plateau areas involves the development and application of an automated logging system designed specifically for plateau areas. This system takes into account the unique environmental conditions in plateau areas, such as low oxygen, large temperature changes, and inconvenient transportation, and uses special equipment and highly adaptable control strategies. Its main purpose is to effectively perform logging operations under these complex and extreme conditions to ensure data accuracy and operational safety. Through automated control, the system can reduce personnel exposure to harsh environments, while improving the efficiency and accuracy of data collection, providing a scientific basis for oil and gas development.

[0004] Although existing automated logging technologies show high efficiency and accuracy in most cases, their application performance in special environments such as plateaus is not ideal. These technologies often lack the necessary adaptability, especially under conditions of low oxygen and drastic temperature fluctuations. The system may not be able to adjust to the optimal drilling parameters in a timely manner, resulting in inaccurate data collection and reduced operational efficiency. In addition, traditional systems often lack flexibility when processing highly variable geological data, making it difficult to effectively distinguish and filter noise or erroneous data, which in turn affects the accuracy of data interpretation. This lack of highly adaptable and optimized data processing capabilities is particularly evident in extreme conditions such as plateaus, often leading to major decision-making errors and increasing the risks and costs of oil and gas development. Summary of the invention

[0005] The purpose of the present invention is to solve the shortcomings of the prior art and to propose an automated well logging operation system for plateau areas.

[0006] In order to achieve the above-mentioned purpose, the present invention adopts the following technical scheme: A well logging operation automation system in plateau areas comprises:

[0007] The rock formation characteristic analysis module performs pixel detection and region annotation based on the rock formation texture and fault line image data of the plateau area, and performs feature extraction and node association analysis on each node to establish a structured geological map;

[0008] The data purification module identifies and filters noise on the structured geological map, performs quality assessment on each data point and removes erroneous data, obtains a clean geological map, performs data enhancement processing on the basis of the clean geological map, adjusts the relationship strength between data points, and creates an enhanced geological map;

[0009] The parameter optimization module extracts rock hardness and wear rate from the enhanced geological map, performs parameter impact analysis and sensitivity test in combination with current drilling parameters, formulates a preliminary adjustment plan, and dynamically adjusts the learning rate and micro-batch data in the preliminary adjustment plan according to real-time drilling feedback to obtain optimized drilling parameters;

[0010] The efficiency evaluation module uses the optimized drilling parameters to perform drilling operations, collects drilling speed and drill bit wear data in real time, performs correlation analysis between speed and wear, generates drilling efficiency results, and evaluates the effectiveness of drilling parameters by comparing historical data and real-time data to obtain parameter efficiency evaluation.

[0011] As a further solution of the present invention, the steps of obtaining the structured geological map are specifically as follows:

[0012] Load the rock texture and fault line image data of the plateau area, and use the edge enhancement algorithm to process the image using the formula:

[0013]

[0014] Calculate edge strength to extract rock texture and fault lines and obtain characteristic boundary map;

[0015] in, Indicates location The edge strength, , is a parameter to adjust the sensitivity of edge detection, represents the image intensity, , Respectively represent the images in and Partial derivatives of direction;

[0016] The feature boundary map is used to perform pixel detection and region annotation, and an image segmentation algorithm based on brightness and texture is used to refine the region division, using the formula:

[0017]

[0018] Perform weighted summation on the pixels in each area to obtain a regional map of geological characteristics;

[0019] in, represents the score of the region, It is The importance weight of each pixel, It is The attribute value of each pixel, is the brightness reference value, is the brightness adjustment factor;

[0020] Based on the geological feature area map, node association analysis is performed to construct a connection network between nodes using the formula:

[0021]

[0022] Evaluate the strength of associations between nodes to obtain a structured geological map;

[0023] in, Representation Node For nodes The strength of association, Is a node and The relationship weight between is the stability parameter.

[0024] As a further solution of the present invention, the steps of obtaining the clean geological map are specifically as follows:

[0025] The structured geological map is subjected to noise identification, and the image data is analyzed using the frequency domain filtering method, using the formula:

[0026]

[0027] The noise pattern was identified by Fourier transform and a Gaussian decay function was used to generate a preliminary filtered geological map;

[0028] in, Indicates location The image intensity after noise attenuation, represents the Fourier transform, is the image at position The original strength, is the attenuation parameter, , The coordinates representing the pixel location;

[0029] The data point quality assessment was performed using the geological map after the preliminary filtering to exclude erroneous data using the formula:

[0030]

[0031] Calculate the consistency of each data point with its neighboring points, exclude abnormal data points by setting thresholds, and generate a clean geological map;

[0032] in, Indicates The quality score of the data point, , Represents data points and neighboring points The attribute value of Yes The set of neighboring points of is a small amount, avoiding the denominator being zero.

[0033] As a further solution of the present invention, the steps of acquiring the enhanced geological map are specifically as follows:

[0034] Starting from the clean geological map, the relationship between data points is strengthened using the formula:

[0035]

[0036] By adjusting the weights between data points, a data relationship reinforcement graph is generated;

[0037] in, Represents data points and The strength of the relationship between , is the attribute value of the corresponding data point, is the stability coefficient, avoiding the denominator being zero, is the spatial attenuation coefficient, is the distance between points;

[0038] Using the data relationship enhancement diagram, spatial continuity analysis is performed using the formula:

[0039]

[0040] By referring to the effect of geometric distance on the strength of relationship, a spatial continuity enhancement graph is generated;

[0041] in, Represents data points The spatial continuity score of is the adjustment parameter of spatial attenuation;

[0042] Integrate the spatial continuity enhancement map and perform image rendering using the formula:

[0043]

[0044] Enhance the readability and detail of geological maps to obtain enhanced geological maps;

[0045] in, represents an enhanced geological map, It is a data point The original attribute value of is the logarithmic gain parameter.

[0046] As a further solution of the present invention, the step of obtaining the optimized drilling parameters is specifically as follows:

[0047] From the enhanced geological map, the data of rock hardness and wear rate are extracted, and the rock formation characteristics are analyzed using image processing technology, using the formula:

[0048]

[0049]

[0050] in, Represents the light intensity value of the rock layer, , Represent the texture and edge strength information of the rock formation, Indicates location The hardness rating of is the weight coefficient of multiple attributes;

[0051] Based on the rock hardness and wear data, combined with the current drilling parameters, parameter impact analysis and sensitivity testing were performed using the formula:

[0052]

[0053] Determine key adjustment points and generate preliminary drilling adjustment plans;

[0054] in, represents the drilling depth, represents the drilling performance prediction, represent the depth and hardness parameters respectively, is a parameter , The influence function of , are the adjustment factors for impact and depth attenuation, respectively;

[0055] The learning rate and micro-batch data are dynamically adjusted according to the preliminary drilling adjustment plan and real-time drilling feedback, using the formula:

[0056]

[0057] After adjustment, the optimized drilling parameters are obtained;

[0058] in, is the learning rate, is a fine-tuning parameter, represents the adjustment amount from drilling feedback, is with Associated drilling performance indicators, Represents the optimization results of drilling parameters.

[0059] As a further solution of the present invention, the steps for obtaining the drilling efficiency result are specifically as follows:

[0060] Start the drilling operation, use the optimized drilling parameters, monitor the drilling speed and drill bit wear data in real time, and use the formula:

[0061]

[0062] generating adjusted drilling speed data;

[0063] in, is the adjusted drilling speed data, is the initial drilling speed, , Adjust the drilling depth and the influence of drill bit wear respectively, Indicates the current drilling depth. is the drill wear index;

[0064] The adjusted drilling speed data was used to perform a speed-wear correlation analysis using the formula:

[0065]

[0066] Calculate the correlation between speed and wear, and generate correlation analysis results between drilling speed and wear data;

[0067] in, is the result of correlation analysis between drilling speed and wear data, is the drilling speed data, is the corresponding drill wear rate;

[0068] According to the correlation analysis results of the drilling speed and wear data, the drilling efficiency is calculated using the formula:

[0069]

[0070] Get drilling efficiency results;

[0071] in, Drilling efficiency results.

[0072] As a further solution of the present invention, the step of obtaining the parameter efficiency evaluation is specifically as follows:

[0073] Based on the drilling efficiency results, the data of drilling speed and drill bit wear rate were collected from the historical database and real-time monitoring system, and the formula was used:

[0074]

[0075] Calculate and generate real-time drilling parameter data;

[0076] in, For real-time drilling parameter data, represents the weight coefficient of speed, represents the weight coefficient of wear, represents the square root of the real-time drilling speed, Represents the natural logarithmic increment of the real-time drill wear rate;

[0077] The real-time drilling parameter data is compared with the stored historical drilling parameter data using the formula:

[0078]

[0079] Calculate the difference between historical data and real-time data to obtain parameter difference analysis results;

[0080] in, is the result of parameter difference analysis. represents the historical drilling parameters, Represents the number of data points;

[0081] According to the results of the parameter difference analysis, the formula is used:

[0082]

[0083] The average difference value was calculated to obtain parameter efficiency evaluation;

[0084] in, is the parameter efficiency evaluation, Represents the total number of data points involved in the calculation.

[0085] Compared with the prior art, the advantages and positive effects of the present invention are:

[0086] In the present invention, the geological data is ensured to be detailed and accurate through precise analysis of rock texture and fault line data in plateau areas, and these high-quality data are used for real-time feedback and parameter adjustment, effectively reducing the drilling risk caused by rapid changes in geological conditions, dynamically adjusting drilling parameters to reduce the possibility of equipment failure and operation interruption, optimizing the continuity and safety of drilling operations, and enhancing the creation of geological maps by improving the relationship between data points, making the data more predictive and operationally guiding, thereby improving the accuracy of decision support and the overall efficiency of oilfield development. BRIEF DESCRIPTION OF THE DRAWINGS

[0087] Figure 1 is a system flow chart of the present invention;

[0088] Figure 2 A flow chart of the steps for obtaining a structured geological map of the present invention;

[0089] Figure 3 A flow chart of the steps for obtaining a clean geological map of the present invention;

[0090] Figure 4 A flow chart of the steps for obtaining an enhanced geological map of the present invention;

[0091] Figure 5 A flow chart of the steps for obtaining the optimized drilling parameters of the present invention;

[0092] Figure 6 A flowchart of the steps for obtaining the drilling efficiency results of the present invention;

[0093] Figure 7 Flow chart of the steps for obtaining parameter efficiency evaluation of the present invention. DETAILED DESCRIPTION

[0094] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0095] In the description of the present invention, it should be understood that the terms "length", "width", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate positions or positional relationships based on the positions or positional relationships shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, in the description of the present invention, "multiple" means two or more, unless otherwise clearly and specifically defined. Example

[0096] See also Figure 1 , a well logging operation automation system in plateau areas includes:

[0097] The rock formation characteristic analysis module performs pixel detection and region annotation based on the rock formation texture and fault line image data of the plateau area, and performs feature extraction and node association analysis on each node to establish a structured geological map;

[0098] The data purification module identifies and filters noise in structured geological maps, performs quality assessment and removes erroneous data for each data point, obtains clean geological maps, performs data enhancement processing on the basis of clean geological maps, adjusts the relationship strength between data points, and creates enhanced geological maps;

[0099] The parameter optimization module extracts rock hardness and wear rate from the enhanced geological map, combines the current drilling parameters, performs parameter impact analysis and sensitivity testing, formulates a preliminary adjustment plan, and dynamically adjusts the learning rate and micro-batch data in the preliminary adjustment plan based on real-time drilling feedback to obtain optimized drilling parameters;

[0100] The efficiency evaluation module uses optimized drilling parameters to perform drilling operations, collects drilling speed and drill wear data in real time, performs correlation analysis between speed and wear, generates drilling efficiency results, and evaluates the effectiveness of drilling parameters by comparing historical data with real-time data to obtain parameter efficiency evaluation.

[0101] The structured geological map includes geological node characteristics, connection relationship indicators, and terrain mapping data. The enhanced geological map specifically adjusts the relationship strength of data points. The optimized drilling parameters include hardness grade, wear rate range, and adjustment frequency. The parameter efficiency evaluation specifically includes efficiency improvement rate, parameter stability, and adjustment response time.

[0102] See also Figure 2 , the steps to obtain the structured geological map are as follows:

[0103] Load the rock texture and fault line image data of the plateau area, and use the edge enhancement algorithm to process the image using the formula:

[0104]

[0105] Calculate edge strength to extract rock texture and fault lines and obtain characteristic boundary map;

[0106] in, Indicates location The edge strength, , is a parameter to adjust the sensitivity of edge detection, represents the image intensity, , Respectively represent the images in and Partial derivatives of direction;

[0107] The feature boundary map is used for pixel detection and region annotation, and the image segmentation algorithm based on brightness and texture is used to refine the region division. The formula is:

[0108]

[0109] Perform weighted summation on the pixels in each area to obtain a regional map of geological characteristics;

[0110] in, represents the score of the region, It is The importance weight of each pixel, It is The attribute value of each pixel, is the brightness reference value, is the brightness adjustment factor;

[0111] Based on the geological feature area map, node association analysis is performed to build a connection network between nodes using the formula:

[0112]

[0113] Evaluate the strength of associations between nodes to obtain a structured geological map;

[0114] in, Representation Node For nodes The strength of association, Is a node and The relationship weight between is the stability parameter.

[0115] Formula 1: Edge detection formula

[0116]

[0117] Parameter explanation and derivation:

[0118] : Image intensity

[0119] , : The image is and The first-order derivative in the direction represents the gradient of the image in the horizontal and vertical directions.

[0120] , : The image is and The second derivative of the direction is used to detect the curvature of the edge.

[0121] , : Adjustment parameters used to balance the effects of first-order and second-order derivatives. Usually, they are adjusted according to image characteristics in experiments. Here we set and .

[0122] Example derivation:

[0123] Assume that at some point At, image intensity ,exist First derivative of direction , First derivative of direction , Second derivative of direction , Second derivative of direction .

[0124] Compute the contribution of the first-order derivative part:

[0125] Compute the contribution of the second-order derivative part:

[0126] Final edge strength

[0127] Formula 2: Weighted summation formula

[0128]

[0129] Parameter explanation and derivation:

[0130] : No. The importance weight of the pixel, assuming the weight is 1

[0131] : No. The attribute value of each pixel is assumed to be 255 (the maximum brightness value of the image)

[0132] : Brightness reference value, assumed to be 100

[0133] : Brightness adjustment coefficient, assuming it is 10

[0134] Example derivation:

[0135] Assume there is only one pixel :

[0136] Calculate the ratio:

[0137] Weighted sum:

[0138] Formula 3: Node association strength formula

[0139]

[0140] Parameter explanation and derivation:

[0141] :node and nodes The relationship weight between is assumed to be 0.8

[0142] : Stability parameter, assumed to be 0.01

[0143] Example derivation:

[0144] Assume Node Total weight with other nodes :

[0145] Calculate the strength of association:

[0146] See also Figure 3 , the specific steps for obtaining a clean geological map are:

[0147] To identify the noise of structured geological maps, the frequency domain filtering method is used to analyze the image data, using the formula:

[0148]

[0149] The noise pattern was identified by Fourier transform and a Gaussian decay function was used to generate a preliminary filtered geological map;

[0150] in, Indicates location The image intensity after noise attenuation, represents the Fourier transform, is the image at position The original strength, is the attenuation parameter, , The coordinates representing the pixel location;

[0151] The data point quality was evaluated using the geological map after preliminary filtering to exclude erroneous data using the formula:

[0152]

[0153] Calculate the consistency of each data point with its neighboring points, exclude abnormal data points by setting thresholds, and generate a clean geological map;

[0154] in, Indicates The quality score of the data point, , Represents data points and neighboring points The attribute value of Yes The set of neighboring points of is a small amount, avoiding the denominator being zero.

[0155] Formula 1: Noise Identification and Attenuation

[0156] formula:

[0157]

[0158] Parameter explanation and calculation process:

[0159] :This is for the image Perform Fourier transform. Fourier transform is used to convert an image from the time domain (spatial domain) to the frequency domain, and is mainly used to identify the frequency components in the image, that is, the noise and signal parts.

[0160] : This represents the absolute value of the Fourier transform result and is used to obtain the amplitude of each point in the frequency domain, that is, the noise intensity.

[0161] : Logarithmic transformation of the amplitude of the Fourier transform can help expand the dynamic range of low-intensity values ​​while compressing the range of high-intensity values, making the noise characteristics of the image more obvious.

[0162] : Attenuation coefficient, which controls the attenuation rate of the Gaussian function. In actual use, this coefficient can be adjusted according to the distribution and density of the noise.

[0163] : Gaussian attenuation function, used to reduce the impact of noise, where and Represents the coordinates of the pixel position. This function decays faster as the distance from the center increases, effectively suppressing edge noise.

[0164] Calculation example:

[0165] Suppose there is a pixel , the intensity of the image at that point is 100. The modulus of the Fourier transform result is assumed to be 50, and the attenuation coefficient Set to 0.001.

[0166]

[0167] here, The value of indicates the intensity of the noise at the position (10, 10) after the noise attenuation process. This value is low, indicating that the noise has been effectively attenuated.

[0168] Formula 2: Data point quality assessment and error data removal

[0169] formula:

[0170]

[0171] Parameter explanation and calculation process:

[0172] : No. The attribute value of a data point.

[0173] : No. Neighboring points of data points The attribute value of .

[0174] :point The set of neighboring points.

[0175] : Used to avoid small quantities with a denominator of zero, set to a small positive value such as 0.01.

[0176] Calculation example:

[0177] Assumptions Point Value , its neighboring points Value and 90.

[0178]

[0179] here, A high value indicates that The consistency with its neighboring points is good and the data quality is high.

[0180] See also Figure 4 , the specific steps for obtaining the enhanced geological map are:

[0181] Starting from the clean geological map, the relationship between data points is strengthened using the formula:

[0182]

[0183] By adjusting the weights between data points, a data relationship reinforcement graph is generated;

[0184] in, Represents data points and The strength of the relationship between , is the attribute value of the corresponding data point, is the stability coefficient, avoiding the denominator being zero, is the spatial attenuation coefficient, is the distance between points;

[0185] Use the data relationship enhancement diagram to perform spatial continuity analysis using the formula:

[0186]

[0187] Enhance the readability and detail of geological maps to obtain enhanced geological maps;

[0188] in, represents an enhanced geological map, It is a data point The original attribute value of is the logarithmic gain parameter.

[0189] Formula analysis for step 1

[0190] formula:

[0191]

[0192] Parameter explanation and examples:

[0193] : Data points and For example, .

[0194] : Stability coefficient, used to avoid the denominator being zero, set to 1.

[0195] : Spatial attenuation coefficient, assumed to be 0.1.

[0196] :point and The distance between them is assumed to be 2.

[0197] Calculation process:

[0198] Calculate the numerator:

[0199] Calculate the denominator:

[0200] Apply a decay factor:

[0201] Final result:

[0202] Result interpretation:

[0203] Represents data points and The relationship strength is the attenuated value, and the closer it is to 1, the stronger the relationship.

[0204] Formula analysis for step 2

[0205] formula:

[0206]

[0207] Parameter explanation and examples:

[0208] : Calculated from step 1, assuming that for all , .

[0209] : The adjustment parameter of spatial attenuation, assumed to be 0.05.

[0210] : Assume that for all , .

[0211] Calculation process:

[0212] For each exist In the neighborhood of .

[0213] Assumptions There are 3 neighbors ,but .

[0214] Result interpretation:

[0215] It means that after spatial attenuation adjustment, point The higher the value, the better the spatial continuity.

[0216] Formula analysis for step 3

[0217] formula:

[0218]

[0219] Parameter explanation and examples:

[0220] : Obtained from step 2, assuming that for all , .

[0221] : Attribute value, assumed to be 20.

[0222] : Logarithmic gain parameter, set to 0.1.

[0223] Calculation process:

[0224] calculate

[0225] Assume there are 10 points in total ,but .

[0226] Result interpretation:

[0227] Indicates the overall strength of the enhanced geological map. The higher the value, the better the readability and detail of the image.

[0228] See also Figure 5 , the specific steps for obtaining optimized drilling parameters are:

[0229] From the enhanced geological map, the rock hardness and wear rate data were extracted, and the rock formation characteristics were analyzed using image processing technology, using the formula:

[0230]

[0231]

[0232] in, Represents the light intensity value of the rock layer, , Represent the texture and edge strength information of the rock formation, Indicates location The hardness rating of is the weight coefficient of multiple attributes;

[0233] Based on the rock hardness and wear data, combined with the current drilling parameters, parameter impact analysis and sensitivity test were performed using the formula:

[0234]

[0235] Determine key adjustment points and generate preliminary drilling adjustment plans;

[0236] in, represents the drilling depth, represents the drilling performance prediction, represent the depth and hardness parameters respectively, is a parameter , The influence function of , are the adjustment factors for impact and depth attenuation, respectively;

[0237] Dynamically adjust the learning rate and micro-batch data according to the preliminary drilling adjustment plan and real-time drilling feedback, using the formula:

[0238]

[0239] After adjustment, the optimized drilling parameters are obtained;

[0240] in, is the learning rate, is a fine-tuning parameter, represents the adjustment amount from drilling feedback, is with Associated drilling performance indicators, Represents the optimization results of drilling parameters.

[0241] Step 1: Rock hardness and wear data extraction

[0242] formula:

[0243]

[0244] Detailed explanation and derivation of formula parameters:

[0245] : This is the rock image at pixel points Assume that the light intensity value of a pixel is 120.

[0246] : This means pixel Texture analysis results in the image. Assume the texture value is 0.5.

[0247] : This means pixel The edge strength in the image. Assume the edge strength is 0.3.

[0248] : These are weight coefficients used to balance the influence of various physical properties. Assume that the weight coefficients .

[0249] Calculation process:

[0250] Calculate the squared contribution of the light intensity value:

[0251]

[0252] Calculate the product contribution of texture and edge:

[0253]

[0254] Calculate the square root contribution of edge strength:

[0255]

[0256] Total hardness score:

[0257]

[0258] Explanation of the example results:

[0259] The final hardness score is approximately 11520.29, which is dominated by the contribution of the light intensity value, indicating that the rock has a higher hardness at this point and is suitable for predicting the resistance encountered in drilling.

[0260] Step 2: Parameter impact analysis and sensitivity testing

[0261] formula:

[0262]

[0263] Detailed explanation and derivation of formula parameters:

[0264] : Drilling depth, assumed to be 300 meters.

[0265] : The hardness score obtained from step 1, i.e. 11520.29.

[0266] : Influence function, assumed to be .

[0267] : Sensitivity coefficient and attenuation rate, assuming .

[0268] Calculation process:

[0269] Calculate the influence function:

[0270]

[0271] Calculate the effect of exponential decay:

[0272]

[0273] The resulting performance forecast is:

[0274]

[0275] Explanation of the example results:

[0276] The performance prediction value is 3.825, which reflects the drilling performance at a given depth and rock hardness, where a lower performance prediction value indicates more difficult drilling conditions.

[0277] Step 3: Dynamic Adjustment of Real-Time Drilling Feedback

[0278] formula:

[0279]

[0280] Detailed explanation and derivation of formula parameters:

[0281] : Learning rate and fine-tuning parameters, assuming .

[0282] : Adjustment amount, assumed to be 0.02.

[0283] : Performance indicators, using the previous steps .

[0284] Calculation process:

[0285] Compute the time derivative contribution: (Here we assume )

[0286] Calculate micro-batch adjustment contribution:

[0287]

[0288] The total optimization parameter changes are:

[0289]

[0290] Explanation of the example results:

[0291] It represents a small optimization of drilling parameters under the current feedback mechanism, showing a slight improvement in drilling efficiency after fine-tuning.

[0292] See also Figure 6 , the specific steps for obtaining the drilling efficiency results are:

[0293] Start the drilling operation, use the optimized drilling parameters, and monitor the drilling speed and drill wear data in real time, using the formula:

[0294]

[0295] generating adjusted drilling speed data;

[0296] in, is the adjusted drilling speed data, is the initial drilling speed, , Adjust the drilling depth and the influence of drill bit wear respectively, Indicates the current drilling depth. is the drill wear index;

[0297] The speed-wear correlation analysis was performed using the adjusted drilling speed data using the formula:

[0298]

[0299] Calculate the correlation between speed and wear, and generate correlation analysis results between drilling speed and wear data;

[0300] in, is the result of correlation analysis between drilling speed and wear data, is the drilling speed data, is the corresponding drill wear rate;

[0301] According to the correlation analysis results between drilling speed and wear data, the drilling efficiency is calculated using the formula:

[0302]

[0303] Get drilling efficiency results;

[0304] in, Drilling efficiency results.

[0305] Detailed explanation and derivation of the formula for step 1

[0306] formula:

[0307]

[0308] Parameter Description:

[0309] : Adjusted drilling speed, indicating the depth The drilling speed at .

[0310] : Initial drilling speed, assumed to be .

[0311] : Depth attenuation coefficient, assumed to be .

[0312] : Drill wear adjustment factor, assumed to be 0.03.

[0313] : Drilling depth, assumed to be .

[0314] : Drill wear index, assumed to be 0.5.

[0315] Calculation process:

[0316] calculate

[0317] calculate

[0318] calculate

[0319] Exponential part of comprehensive formula

[0320] Calculate drilling speed

[0321] Result description:

[0322] Calculation results It means that at a depth of 100 meters, after considering the influence of drill bit wear, the drilling speed slows down significantly and almost stops, showing extremely high attenuation.

[0323] Detailed explanation and derivation of the formula for step 2

[0324] formula:

[0325]

[0326] Parameter Description:

[0327] : Correlation coefficient between speed and wear.

[0328] : Drilling speed data, assumed to be .

[0329] : Drill wear rate data, assumed to be or 0.5.

[0330] Calculation process:

[0331] calculate

[0332] calculate

[0333] calculate

[0334] calculate

[0335] calculate

[0336] Result description:

[0337] result It indicates that there is a perfect positive correlation between drilling speed and drill bit wear, i.e., the slower the speed, the greater the wear, or vice versa, which in practical applications may indicate that drilling parameters need to be adjusted to optimize performance.

[0338] Detailed explanation and derivation of the formula for step 3

[0339] formula:

[0340]

[0341] Parameter Description:

[0342] : Drilling efficiency.

[0343] : The correlation coefficient calculated in step 2 is known to be 1.

[0344] Calculation process:

[0345] calculate

[0346] use

[0347] calculate

[0348] Result description:

[0349] result This indicates that current drilling efficiency is extremely low, which may be caused by the perfect positive correlation between drilling speed and drill bit wear.

[0350] See also Figure 7 , the specific steps for obtaining parameter efficiency evaluation are:

[0351] Based on the drilling efficiency results, the data of drilling speed and drill bit wear rate are collected from the historical database and real-time monitoring system, and the formula is used:

[0352]

[0353] Calculate and generate real-time drilling parameter data;

[0354] in, For real-time drilling parameter data, represents the weight coefficient of speed, represents the weight coefficient of wear, represents the square root of the real-time drilling speed, Represents the natural logarithmic increment of the real-time drill wear rate;

[0355] The real-time drilling parameter data is compared with the stored historical drilling parameter data using the formula:

[0356]

[0357] Calculate the difference between historical data and real-time data to obtain parameter difference analysis results;

[0358] in, is the result of parameter difference analysis. represents the historical drilling parameters, Represents the number of data points;

[0359] According to the results of parameter difference analysis, the formula is used:

[0360]

[0361] The average difference value was calculated to obtain parameter efficiency evaluation;

[0362] in, is the parameter efficiency evaluation, Represents the total number of data points involved in the calculation.

[0363] Assume the following numbers:

[0364] Drilling speed In real-time monitoring, it is usually between 100 and (revolutions per minute).

[0365] Drill wear rate Usually in arrive between.

[0366] Weight coefficient and , usually depends on the actual equipment and operating conditions, it is assumed here and .

[0367] formula:

[0368]

[0369] Parameter details:

[0370] : The weight coefficient of speed, which determines the influence of speed in the overall formula.

[0371] : The weight coefficient of wear determines the influence of wear rate in the overall formula.

[0372] : Real-time drilling speed.

[0373] Or 0.0005: Real-time drill wear rate.

[0374] Calculation process and examples:

[0375] Impact of calculation speed:

[0376] Calculate the effect of wear rate: (Using natural logarithms)

[0377] Putting these two parts together:

[0378] This result It is the weighted result of the calculated real-time drilling parameters, reflecting the combined effect of drilling speed and wear rate.

[0379] formula:

[0380]

[0381] Parameter details:

[0382] : Historical drilling parameters.

[0383] : Real-time drilling parameters calculated above.

[0384] : The total number of data points.

[0385] Calculation process and examples:

[0386] For a single data point the difference is calculated:

[0387] Assuming that all points have similar differences, then

[0388] This result It is the sum of the differences of all data points, reflecting the total amount of difference between historical data and real-time data.

[0389] formula:

[0390]

[0391] Parameter details:

[0392] : The sum of the differences calculated above.

[0393] : The total number of data points.

[0394] Calculation process and examples:

[0395] Calculate the mean difference:

[0396] This result It is an evaluation of the efficiency of drilling parameters, showing the average degree of difference between historical and real-time data. This helps to evaluate the accuracy and reliability of the real-time monitoring system.

[0397] The above are only preferred embodiments of the present invention and are not intended to limit the present invention in other forms. Any technician familiar with the profession may use the technical contents disclosed above to change or modify them into equivalent embodiments with equivalent changes and apply them to other fields. However, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention without departing from the technical solution of the present invention still falls within the protection scope of the technical solution of the present invention.

Claims

1. A well logging operation automation system in plateau areas, characterized in that: The system comprises: The rock formation characteristic analysis module performs pixel detection and region annotation based on the rock formation texture and fault line image data of the plateau area, and performs feature extraction and node association analysis on each node to establish a structured geological map; The data purification module identifies and filters noise on the structured geological map, performs quality assessment on each data point and removes erroneous data, obtains a clean geological map, performs data enhancement processing on the basis of the clean geological map, adjusts the relationship strength between data points, and creates an enhanced geological map; The parameter optimization module extracts rock hardness and wear rate from the enhanced geological map, performs parameter impact analysis and sensitivity test in combination with current drilling parameters, formulates a preliminary adjustment plan, and dynamically adjusts the learning rate and micro-batch data in the preliminary adjustment plan according to real-time drilling feedback to obtain optimized drilling parameters; The efficiency evaluation module uses the optimized drilling parameters to perform drilling operations, collects drilling speed and drill bit wear data in real time, performs correlation analysis between speed and wear, generates drilling efficiency results, and evaluates the effectiveness of drilling parameters by comparing historical data with real-time data to obtain parameter efficiency evaluation; The steps for obtaining the structured geological map are specifically as follows: Load the rock texture and fault line image data of the plateau area, and use the edge enhancement algorithm to process the image using the formula: Calculate edge strength to extract rock texture and fault lines and obtain characteristic boundary map; in, Indicates location The edge strength, , is a parameter to adjust the sensitivity of edge detection, represents the image intensity, , Respectively represent the images in and Partial derivatives of direction; The feature boundary map is used to perform pixel detection and region annotation, and an image segmentation algorithm based on brightness and texture is used to refine the region division, using the formula: Perform weighted summation on the pixels in each area to obtain a regional map of geological characteristics; in, represents the score of the region, It is The importance weight of each pixel, It is The attribute value of each pixel, is the brightness reference value, is the brightness adjustment factor; Based on the geological feature area map, node association analysis is performed to construct a connection network between nodes using the formula: Evaluate the strength of associations between nodes to obtain a structured geological map; in, Representation Node For nodes The strength of association, Is a node and The relationship weight between is the stability parameter; The steps for obtaining the clean geological map are specifically as follows: The structured geological map is subjected to noise identification, and the image data is analyzed using the frequency domain filtering method, using the formula: in, Indicates location The image intensity after noise attenuation, represents the Fourier transform, is the image at position The original strength, is the attenuation parameter, , The coordinates representing the pixel location; The data point quality assessment was performed using the geological map after the preliminary filtering to exclude erroneous data using the formula: Calculate the consistency of each data point with its neighboring points, exclude abnormal data points by setting thresholds, and generate a clean geological map; in, Indicates The quality score of the data point, , Represents data points and neighboring points The attribute value of Yes The set of neighboring points of It is a small amount, avoiding the denominator being zero; The steps for obtaining the enhanced geological map are specifically as follows: Starting from the clean geological map, the relationship between data points is strengthened using the formula: By adjusting the weights between data points, a data relationship reinforcement graph is generated; in, Represents data points and The strength of the relationship between , is the attribute value of the corresponding data point, is the stability coefficient, avoiding the denominator being zero, is the spatial attenuation coefficient, is the distance between points; Using the data relationship enhancement diagram, spatial continuity analysis is performed using the formula: By referring to the effect of geometric distance on the strength of relationship, a spatial continuity enhancement graph is generated; in, Represents data points The spatial continuity score of is the adjustment parameter of spatial attenuation; Integrate the spatial continuity enhancement map and perform image rendering using the formula: Enhance the readability and detail of geological maps to obtain enhanced geological maps; in, represents an enhanced geological map, It is a data point The original attribute value of is the logarithmic gain parameter; The steps for obtaining the optimized drilling parameters are specifically as follows: From the enhanced geological map, the data of rock hardness and wear rate are extracted, and the rock formation characteristics are analyzed using image processing technology, using the formula: Generate rock hardness and wear data; in, Represents the light intensity value of the rock layer, , Represent the texture and edge strength information of the rock formation, Indicates location The hardness rating, is the weight coefficient of multiple attributes; Based on the rock hardness and wear data, combined with the current drilling parameters, parameter impact analysis and sensitivity testing were performed using the formula: Determine key adjustment points and generate preliminary drilling adjustment plans; in, represents the drilling depth, represents the drilling performance prediction, Represent the depth and hardness parameters respectively, is a parameter , The influence function of , are the adjustment factors for impact and depth attenuation, respectively; The learning rate and micro-batch data are dynamically adjusted according to the preliminary drilling adjustment plan and real-time drilling feedback, using the formula: After adjustment, the optimized drilling parameters are obtained; in, is the learning rate, is a fine-tuning parameter, represents the adjustment amount from drilling feedback, is with Associated drilling performance indicators, Indicates the optimization results of drilling parameters; The steps for obtaining the drilling efficiency result are specifically as follows: Start the drilling operation, use the optimized drilling parameters, monitor the drilling speed and drill bit wear data in real time, and use the formula: generating adjusted drilling speed data; in, is the adjusted drilling speed data, is the initial drilling speed, , Adjust the drilling depth and the influence of drill bit wear respectively, Indicates the current drilling depth. is the drill wear index; The adjusted drilling speed data was used to perform a speed-wear correlation analysis using the formula: Calculate the correlation between speed and wear, and generate correlation analysis results between drilling speed and wear data; in, is the result of correlation analysis between drilling speed and wear data, is the drilling speed data, is the corresponding drill wear rate; According to the correlation analysis results of the drilling speed and wear data, the drilling efficiency is calculated using the formula: Get drilling efficiency results; in, for drilling efficiency results; The steps for obtaining the parameter efficiency evaluation are specifically as follows: Based on the drilling efficiency results, the data of drilling speed and drill bit wear rate were collected from the historical database and real-time monitoring system, and the formula was used: Calculate and generate real-time drilling parameter data; in, For real-time drilling parameter data, represents the weight coefficient of speed, represents the weight coefficient of wear, represents the square root of the real-time drilling speed, Represents the natural logarithmic increment of the real-time drill wear rate; The real-time drilling parameter data is compared with the stored historical drilling parameter data using the formula: Calculate the difference between historical data and real-time data to obtain parameter difference analysis results; in, is the result of parameter difference analysis. represents the historical drilling parameters, Represents the number of data points; According to the results of the parameter difference analysis, the formula is used: The average difference value was calculated to obtain parameter efficiency evaluation; in, is the parameter efficiency evaluation, Represents the total number of data points involved in the calculation.

Images