A memory type time-depth logging method

Abstract

The application provides a storage type time-depth logging method, and relates to the technical field of logging control.The method comprises the following steps: obtaining encoder, tension and acceleration data, drawing an encoder curve, a tension curve and an acceleration curve; receiving an effective motion time interval of a single complete drill rod operation according to curve characteristics; keeping the total change amount and total length of the encoder value unchanged for each interval original encoder curve segment, re-matching the encoder value and time point according to the layered scale relationship between the encoder data and the actual displacement of the hook, and obtaining a corrected target encoder curve; converting all target encoder curves into time-depth curves according to the depth interval corresponding to each drill rod, and splicing the time-depth curves according to the operation sequence of the drill rod operation event interval to generate a time-depth file.The method solves the technical problem of how to improve the depth measurement precision of the storage type time-depth logging and the reliability of the system as a whole.

Description

Technical Field

[0001] This application relates to the field of well logging control technology, specifically to a storage-based time-depth well logging method. Background Technology

[0002] In oil and gas exploration, for complex well conditions such as highly deviated wells and horizontal wells, pushing a storage-type time-depth logging instrument to the target formation via a rigidly connected drill pipe is an effective operational method. During this process, a reliable depth measurement scheme is needed to accurately correlate the formation parameters stored in the instrument with the depth.

[0003] To achieve the above objectives, one existing technical solution is based on depth information provided by a surface logging system. This solution directly uses the drill pipe depth recorded by the logging system as the depth of the logging instrument. Another improved solution is to add a dedicated magnetic positioning sub to the downhole instrument string, which assists in depth correction and segmentation by identifying magnetic marks pre-installed on the drill pipe coupling.

[0004] However, the aforementioned existing technical solutions have inherent limitations in practical applications. Directly using the depth provided by the logging system results in a non-negligible systematic error between the provided depth and the actual depth of the logging instrument, easily causing local distortion of the logging curve. While using a single magnetic positioning instrument for auxiliary correction improves accuracy, the reliability of the entire depth system depends entirely on this dedicated positioning section. In harsh downhole environments, if this positioning instrument malfunctions or its signal is interfered with, depth information can be completely lost or produce significant errors, leading to a high risk of single-point failure and impacting the success rate of logging operations. Summary of the Invention

[0005] This application provides a storage-based time-depth logging method to address the technical problem of how to improve the accuracy of depth measurement and the overall reliability of the system in existing technologies.

[0006] Firstly, this application provides a storage-based time-depth logging method, including:

[0007] During the process of raising or lowering the drill pipe, encoder data from the encoder installed on the winch cable drum, tension data from the tension sensor connected to the hook cable, and acceleration data from the accelerometer built into the downhole logging instrument are acquired. The winch cable drum is used to raise or lower the drill pipe.

[0008] The encoder curve is plotted based on the encoder data, the tension curve is plotted based on the tension data, and the acceleration curve is plotted based on the acceleration data.

[0009] The user receives a selection of several effective motion time intervals for a single complete drill pipe operation based on the graphic features of the encoding curve, the tension curve, and the acceleration curve. The encoding curve, the tension curve, and the acceleration curve include multiple drill pipe operation event intervals, and each drill pipe operation event interval corresponds to the lifting or lowering process of one drill pipe.

[0010] For any original encoded curve segment within any effective motion time interval, keep the total change and total duration of the encoded value of the original encoded curve segment unchanged, and re-match the encoded value and time point on the original encoded curve segment according to the layered scale relationship between the encoder data and the actual displacement of the hook, to obtain the corrected target encoded curve.

[0011] All the target encoded curves are converted into multiple time-depth curves according to the depth range corresponding to each drill pipe. All the time-depth curves are then spliced ​​together according to the operation sequence of the corresponding drill pipe operation event range to generate a time-depth file. The time-depth curves represent the correspondence between time and depth.

[0012] Optionally, the step of re-matching the encoding values ​​and time points on the original encoding curve segment based on the pre-acquired layered scale relationship between the encoder data and the actual displacement of the hook to obtain the corrected target encoding curve further includes, before:

[0013] The winch is operated to continuously raise the hook from the starting height point to the ending height point, and the encoder count value sequence and the real-time hook height sequence are continuously recorded.

[0014] Align the count value sequence and the real-time height sequence of the large hook in time to generate an original scatter dataset with the encoder count value as the independent variable and the real-time height of the large hook as the dependent variable;

[0015] Piecewise linear regression or nonlinear curve fitting is performed on the original scatter dataset to construct a piecewise function model for quantifying the nonlinear correspondence between the encoder count value and the actual displacement of the hook. The piecewise function model is used to represent the hierarchical scale relationship.

[0016] Optionally, the step of re-matching the encoding values ​​and time points on the original encoding curve segment based on the pre-acquired layered scale relationship between the encoder data and the actual displacement of the hook to obtain the corrected target encoding curve specifically includes:

[0017] Based on the aforementioned layered scale relationship, calculate the theoretical total displacement of the large hook corresponding to the original encoded curve segment;

[0018] The original encoded curve segment is divided into several consecutive encoded increments, and the physical displacement corresponding to each encoded increment is calculated based on the hierarchical scale relationship.

[0019] Based on the proportion of each physical displacement in the total theoretical displacement, the ideal time allocation proportion corresponding to each coding increment is calculated;

[0020] Based on the ideal time allocation ratio corresponding to each coding increment, the coding value on the original coding curve segment is calculated and matched with the corrected time point to obtain the target coding curve.

[0021] Optionally, after the step of generating the time-depth file, the method further includes:

[0022] The time-depth file is matched with the formation data measured by the downhole logging instrument based on the same data timestamp to generate a logging curve with depth on the horizontal axis and formation data on the vertical axis.

[0023] Optionally, the step of converting all the target encoded curves into multiple time-depth curves according to the depth interval corresponding to each drill pipe specifically includes:

[0024] For any of the target encoding curves, the time-depth curve is converted by the following steps:

[0025] Construct a mapping relationship between each encoded value of the target encoded curve and the depth range of the drill pipe corresponding to the target encoded curve;

[0026] Through the mapping relationship, each encoded value of the target encoding curve is converted into an actual depth value, thereby generating a time-depth curve corresponding to the target encoding curve.

[0027] Optionally, the step of plotting the tension curve based on the tension data specifically includes:

[0028] Real-time acquisition of displacement data of the hook's lateral swing;

[0029] Calculate the oscillation period and oscillation amplitude based on the displacement data;

[0030] The additional force of the tension data is calculated based on the oscillation period and the oscillation amplitude.

[0031] When the additional force exceeds the correction threshold, the tension data is compensated and corrected according to the additional force to obtain the target tension data;

[0032] The tension curve is plotted using the target tension data.

[0033] Optionally, the step of plotting the acceleration curve based on the acceleration data specifically includes:

[0034] Acquire attitude data of the downhole logging instrument during the lifting or lowering process;

[0035] The tilt angle and azimuth angle of the downhole logging instrument are calculated based on the attitude data.

[0036] Perform three-dimensional spatial coordinate transformation on the acceleration data based on the tilt angle and the azimuth angle;

[0037] Gravity component compensation is performed on the converted acceleration data to obtain the target acceleration data;

[0038] The acceleration curve is plotted using the target acceleration data.

[0039] In a second aspect, embodiments of this application provide a storage-type time-depth logging device, which includes: one or more processors and a memory; the memory is coupled to the one or more processors, and the memory is used to store computer program code, which includes computer instructions, and the one or more processors call the computer instructions to cause the storage-type time-depth logging device to perform the method described in the first aspect and any possible implementation thereof.

[0040] Thirdly, embodiments of this application provide a computer program product containing instructions that, when the computer program product is run on a stored time-depth logging device, cause the stored time-depth logging device to perform the method described in the first aspect and any possible implementation thereof.

[0041] Fourthly, embodiments of this application provide a computer-readable storage medium including instructions that, when executed on a stored time-depth logging device, cause the stored time-depth logging device to perform the method described in the first aspect and any possible implementation thereof.

[0042] In summary, one or more technical solutions provided in the embodiments of this application have at least the following technical effects or advantages:

[0043] 1. By adopting the above technical solution, this method integrates three different sources and physical dimensions of information: encoder data from the surface, tension data, and acceleration data from the downhole instrument itself. Acceleration data accurately reflects the actual motion state of the downhole instrument, while tension data characterizes the stress state of the drill pipe during lifting or lowering. By combining these two signals with encoder data, which is the primary measure of displacement, the actual motion state of the downhole logging instrument can be effectively identified. Furthermore, internal time rematching of the encoding curve based on the layered scale relationship further eliminates the error between the encoder's encoded value and the actual displacement of the hook caused by the layered winding of the drum. This collaborative working and correction mechanism of multiple information sources not only significantly improves the accuracy of depth measurement under complex well conditions, but also, due to the diversity of information sources, allows for inference based on other data even if a sensor temporarily fails, thus greatly enhancing the reliability and stability of depth measurement.

[0044] 2. By adopting the above technical solution, this method obtains a large amount of raw data by synchronously acquiring continuous encoder count sequences and real-time hook height sequences as the true values. Compared with the traditional method of manually recording only a few layer-changing points, this method based on piecewise linear regression or nonlinear curve fitting of a large amount of scattered data can more accurately capture the error caused by the different real-time hook height displacements corresponding to one turn of the cable in different layers due to the layered winding of the drum, while also reducing human observation errors. Therefore, the layered scale relationship model constructed by this method has higher accuracy.

[0045] 3. By adopting the above technical solution, this method acknowledges and bases on the physical fact revealed by the layered scale relationship that "equal coding increments are not equal to equal physical displacements." By calculating the proportion of each microscopic physical displacement in the theoretical total displacement and using this as a weight to redistribute the total duration, this method reconstructs a primitive, non-uniform process in terms of physical displacement into an ideal process that is non-linear in time but physically more consistent with a uniform motion model. This reconstruction of the internal time axis accurately restores the instrument's actual movement rhythm downhole, resulting in a target coding curve with extremely high time fidelity.

[0046] 4. By adopting the above technical solution, the process of converting the target encoded curve into a time-depth curve is specified by a conversion method based on mapping relationships. This method constructs and applies a direct mapping from the encoded value of the target encoded curve to its corresponding drill pipe depth range, converting the abstract encoded value into an actual depth value with clear physical meaning. This direct mapping conversion method ensures that the high-precision time relationship obtained in the internal time rematching step can be completely inherited and transferred to the final time-depth curve, thereby guaranteeing the accuracy of the time-depth file. Attached Figure Description

[0047] Figure 1 This is a flowchart illustrating a storage-based time-depth logging method in an embodiment of this application;

[0048] Figure 2 This is a schematic diagram of the stored time-depth logging method in the embodiments of this application;

[0049] Figure 3 This is a schematic diagram of the physical structure of a storage-type time-depth logging device in the embodiments of this application. Detailed Implementation

[0050] To enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments.

[0051] In the description of the embodiments of this application, the words "for example" or "for instance" are used to indicate examples, illustrations, or explanations. Any embodiment or design that is described as "for example" or "for instance" in the embodiments of this application should not be construed as being more preferred or advantageous than other embodiments or design options. Rather, the use of the words "for example" or "for instance" is intended to present the relevant concepts in a specific manner.

[0052] In the description of the embodiments of this application, the term "multiple" means two or more. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the indicated technical features. Thus, a feature defined as "first," "second," or "third" may explicitly or implicitly include one or more of that feature. The terms "comprising," "including," "having," and variations thereof all mean "including but not limited to," unless otherwise specifically emphasized.

[0053] This application provides a storage-based time-depth logging method, referencing... Figure 1 , Figure 1 This is a flowchart of a storage-based time-depth logging method provided in an embodiment of this application. The method includes:

[0054] Step S101: During the process of raising or lowering the drill pipe, acquire encoder data from the encoder installed on the winch cable drum, tension data from the tension sensor connected to the hook cable, and acceleration data from the accelerometer built into the downhole logging instrument.

[0055] Drill pipe typically refers to a hollow steel pipe of standard length used in drilling or logging operations to connect surface equipment to downhole tools. A winch cable drum is a power device used to raise or lower the entire drill pipe by winding up and unwinding a steel cable. An encoder is a sensor mounted on the rotating shaft of the winch cable drum to measure the drum's angular displacement; its output is a series of pulse counts or digital values ​​proportional to the drum's rotation angle. A hook cable is a steel cable connecting the winch and the hook. A tension sensor is a sensor installed on the load-bearing path of the hook cable to measure the tensile force on the cable; it can be, for example, a strain gauge tension meter or a piezomagnetic sensor, and its output tension data represents the real-time load value of the cable. A downhole logging instrument is a device that does not rely on cables, has its own power supply and data storage module, and is used to measure formation physical parameters. An accelerometer is an inertial sensor integrated within the downhole logging instrument to measure the instrument's acceleration in at least one axis; its output acceleration data represents changes in the instrument's motion state.

[0056] Specifically, this step is executed throughout the entire process of raising and lowering the drilling rig in storage-based depth logging. Before the operation begins, the encoder, tension sensor, and surface data acquisition system are connected, and the time of the downhole logging instrument is synchronized with the time of the surface system. During the operation, the surface data acquisition system continuously acquires the count values ​​output by the encoder and the voltage or current signals output by the tension sensor at a preset sampling period (e.g., 100ms), and records these signals along with timestamps. Simultaneously, the downhole logging instrument, at its internally set sampling period (e.g., 100ms), independently acquires the acceleration values ​​of its built-in accelerometer in multiple dimensions, such as the axial and radial directions, and stores these acceleration values ​​along with timestamps provided by the instrument's internal clock in its built-in storage medium. These three sets of data streams with independent timestamps together constitute the basic raw data required for subsequent depth calculation and correction.

[0057] Step S102: Draw an encoding curve based on the encoder data, draw a tension curve based on the tension data, and draw an acceleration curve based on the acceleration data.

[0058] The encoding curve refers to a two-dimensional curve with the operation time as the horizontal axis and the encoder's encoding value as the vertical axis. The slope of this curve directly reflects the rotational speed of the winch drum. The tension curve refers to a two-dimensional curve with the operation time as the horizontal axis and the tension value of the tension sensor as the vertical axis. The value and shape changes of this curve reflect the stress state of the entire drill pipe, such as stationary, lifting, lowering, or encountering resistance. The acceleration curve refers to a two-dimensional curve with the operation time as the horizontal axis and the acceleration value of the accelerometer as the vertical axis. The vibration amplitude and frequency characteristics of this curve directly characterize the actual motion state of the downhole logging instrument, such as stationary, uniform motion, or variable speed motion. The plotting can be an actual graphical display or a digital construction during data processing.

[0059] Specifically, this step is the core of data preprocessing and visualization. After acquiring the raw data, the system first needs to perform necessary processing, such as converting the raw electrical signal from the tension sensor into a tension value in tons or kilonewtons using calibration coefficients, and converting the raw output value from the accelerometer using calibration coefficients and performing gravity component compensation to obtain a net acceleration value in m / s². Then, the system digitally presents the three sets of processed data, based on their respective timestamps, in one or more associated time coordinate systems. These three curves can be simultaneously displayed on a graphical user interface for operators to observe and interact with. For example, in a typical drill pipe hoisting operation, the user will typically see the tension curve climb from a lower, stable value and enter a higher plateau, while the coded curve begins to rise steadily, followed by a significant group of vibrational peaks in the acceleration curve.

[0060] Step S103: Receive the effective motion time interval of a selected number of single complete drill pipe operations based on the graphical features of the encoding curve, the tension curve, and the acceleration curve.

[0061] The graphic features refer to the typical morphological combinations on the three curves that visually reflect the start and end of a complete drill pipe operation. For example, the initial graphic features may include the tension curve starting to deviate from the stationary baseline and climbing upwards, followed by severe vibrations in the acceleration curve; the ending graphic features may include the tension curve significantly declining from a high plateau, and the vibrations of the acceleration curve returning to background noise levels. The effective motion time interval refers to a continuous period of time selected by the user on the timeline through human-computer interaction. The start and end points of this interval define, in time, an independent drill pipe lifting or lowering operation used for logging data acquisition. The drill pipe operation event interval refers to a data segment in the continuous data stream that contains the aforementioned graphic features, corresponding to each physical drill pipe operation (such as lifting one drill pipe or lowering one drill pipe).

[0062] For details, please refer to Figure 2 , Figure 2 This step defines a core human-computer interaction process, illustrating the encoding curve, tension curve, and acceleration curve. On the visualization interface generated in step S102, the following is displayed: Figure 2 The system displays three curves: the encoding curve, the tension curve, and the acceleration curve. Interactive tools are also provided, such as two vertical marker lines that users can drag using a mouse or touch. Users identify a complete operation by observing the interaction between these three curves. Typically, after the cable is taut, the tension curve reaches a high, stable plateau. Only then will the hook forcefully pull the drill bit, causing the instrument to move. Acceleration then occurs, and the acceleration curve begins to fluctuate significantly. Simultaneously, the encoding curve begins to change. Based on this logic, for example, a user might drag the starting marker line to a point where they determine the tension is about to reach a high, stable plateau, while observing whether the acceleration curve begins to fluctuate significantly and the encoding curve begins to change. When the tension is about to reach a high, stable plateau, and both the acceleration curve and the encoding curve begin to change, this can be used to determine the starting point for the cutoff. Users can also select the cutoff endpoint based on characteristics such as the tension curve starting to decline and the acceleration curve becoming flatter at the endpoint. Once the user confirms the selection, the system records the start and end times defined by the two marked lines, forming a valid movement time interval, and associates it with the corresponding drill pipe number and length for subsequent processing.

[0063] Step S104: For any original encoded curve segment within any effective motion time interval, keep the total change in the encoded value and the total duration of the original encoded curve segment unchanged. Based on the layered scale relationship between the encoder data and the actual displacement of the hook, re-match the encoded value and time point on the original encoded curve segment to obtain the corrected target encoded curve.

[0064] The original encoded curve segment refers to the time-encoded value correspondence extracted from the effective motion time interval determined in step S103, without any correction. The layered calibration relationship refers to a mathematical model, pre-established through field calibration, used to describe the nonlinear mapping between encoder count values ​​and the actual physical displacement of the hook in the vertical direction. This model quantifies the measurement error caused by the different effective diameters of the winch cable drum at different winding layers. Re-matching refers to an algorithmic process that nonlinearly reconstructs the timestamps corresponding to each encoded value point within the interval, while keeping the total interval duration and total encoding change constant. The corrected target encoded curve refers to the new time-encoded value correspondence whose internal shape, after the above re-matching, more realistically reflects the uniform physical displacement of the downhole instrument.

[0065] Specifically, this step aims to eliminate depth measurement errors caused by layered winding of the rollers. Its core is to correct the linear or irregular distribution of the original coding curve on the time axis to a distribution that better reflects physical reality. After receiving a valid motion time interval determined in step S103, the system invokes the layered scale relationship. This layered scale relationship reveals that the physical displacement represented by each coding increment is unequal at different stages of coding value change. The purpose of this rematching process is to ensure that, on the corrected target coding curve, within any equal time interval, the corresponding physical displacement (converted through the layered scale relationship) tends to be consistent or conforms to a preset ideal motion model. For example, if the layered scale relationship shows that the displacement corresponding to a certain coding increment is large, then after rematching, the time occupied by that coding increment on the time axis will be correspondingly lengthened, and vice versa, thus obtaining the corrected target coding curve.

[0066] Step S105: Convert all the target encoded curves into multiple time-depth curves according to the depth interval corresponding to each drill pipe, and stitch all the time-depth curves together according to the operation sequence of the corresponding drill pipe operation event interval to generate a time-depth file.

[0067] The depth range refers to a depth interval corresponding to a specific drill pipe operation, with a clearly defined starting and ending depth, and its length is equal to the actual length of one or more drill pipes involved in that operation. The time-depth curve is a two-dimensional curve with operation time as the horizontal axis and the absolute depth value of the downhole instrument as the vertical axis. The splicing refers to connecting multiple time-depth curves representing different drill pipe operation intervals, which may be discontinuous in time, according to the chronological order of their physical operations, into a continuous time-depth relationship that runs through the entire logging section. The time-depth file is the final generated data file recording the downhole depth corresponding to each point in time during the entire logging operation.

[0068] Specifically, this step aims to convert the internally time-corrected target coded curves, with the coded value as the vertical axis, into a final result with absolute depth as the vertical axis. For each target coded curve, the system first determines its corresponding depth interval. The starting depth of this interval is the cumulative length of all previously operated drill pipe, and the ending depth is the starting depth plus the actual length of the drill pipe operated in this operation. Then, the system establishes and applies a mapping relationship from the range of coded values ​​of the target coded curve to its corresponding depth interval, converting the coded value corresponding to each time point on the curve into an absolute depth value within that depth interval. By performing this conversion operation on all target coded curves, multiple independent time-depth curves can be obtained. Finally, the system stitches these time-depth curves together according to their actual chronological order in the downhole operation, thereby forming a unique and continuous time-depth file covering the entire logging operation timeline.

[0069] The following is a more detailed description of the process of the method provided in this implementation.

[0070] Optionally, steps S106-S108 are the steps for obtaining the layer scale relationship in this scheme.

[0071] Step S106: Operate the winch to continuously raise the hook from the starting height point to the ending height point, and continuously record the encoder count value sequence and the real-time hook height sequence.

[0072] The starting and ending height points refer to the physical limit positions of the hook within the derrick guide rails, or a sufficiently large, manually set travel range covering all possible winding layers of the winch drum during logging operations. The count value sequence refers to a series of continuously changing count values ​​output by the encoder during the hoisting process. The real-time hook height sequence refers to the real-time vertical position data of the hook, obtained as the true value, through an independent, high-precision reference measurement system during the hoisting process. This reference measurement system can be a laser rangefinder, a wire displacement sensor, or a machine vision-based measurement system.

[0073] Specifically, this step aims to acquire all the raw data used to calculate the stratified scale relationships through a standardized physical operation. Before conducting formal logging operations, field operators perform a hook lifting operation under no-load or light-load conditions. During the operation, the surface data acquisition system records two independent data channels simultaneously and in parallel at a high sampling frequency: one is the sequence of count values ​​generated by the encoder connected to the winch drum; the other is the real-time hook height sequence measured by a high-precision reference measurement system. Both data sequences are associated with a unified timestamp, ensuring the feasibility of subsequent time alignment. This lifting process typically requires as uniform a speed as possible to minimize the interference of dynamic effects on the calibration results.

[0074] Step S107: Align the count value sequence and the real-time height sequence of the large hook in time to generate an original scatter dataset with the encoder count value as the independent variable and the real-time height of the large hook as the dependent variable.

[0075] The temporal alignment refers to pairing each encoder count value with the real-time hook height value measured at the same moment by matching the same timestamps in two data sequences. The original scatter dataset is a two-dimensional dataset consisting of a large number of data points (E, D), where E represents the encoder count value at a certain moment, and D represents the actual hook height at the same moment.

[0076] Specifically, this step is the core of data preparation and preprocessing. After receiving the two time-series data collected in step S106, the system executes a data alignment algorithm. This algorithm iterates through each timestamp of one sequence (e.g., the encoder count sequence) and searches for data points with the same or closest timestamps in the other sequence (the real-time hook height sequence). If the sampling frequencies of the two sequences are different, the system will also use linear interpolation or higher-order interpolation methods to calculate the precise hook height corresponding to each encoder sampling time. By finding the true height of each encoder count at the same time, the system generates an original scatter dataset that reflects the intrinsic functional relationship between the two.

[0077] Step S108: Perform piecewise linear regression or nonlinear curve fitting on the original scatter dataset to construct a piecewise function model for quantifying the nonlinear correspondence between the encoder count value and the actual displacement of the hook. The piecewise function model is used to represent the hierarchical scale relationship.

[0078] Piecewise linear regression refers to a statistical modeling method that divides a dataset into multiple segments and performs linear regression within each segment. Nonlinear curve fitting refers to using a nonlinear function (such as a polynomial function, exponential function, etc.) to approximate the distribution trend of the entire scatter dataset. The piecewise function model refers to a mathematical model composed of multiple sub-functions that are effective in different encoding value intervals; this model is the final output hierarchical scaling relationship.

[0079] Specifically, this step is the core modeling process of extracting physical laws from data using mathematical tools. The system takes the original scatter plot dataset generated in step S107 as input. A typical implementation is that the system first identifies "inflection points" where the slope changes significantly in the scatter plot; these inflection points physically correspond to the layer-changing positions of the winch drum. Then, using these inflection points as boundaries, the dataset is divided into multiple subsets. The system applies least squares linear regression to each subset to obtain the slope and intercept within that segment, thus constructing an expression for a sub-function. Combining the expressions of all sub-functions and their effective coded value ranges forms a complete piecewise function model. This model can accurately quantify the nonlinear relationship of "how many meters the hook actually moves for each additional count when the coded value is within a certain range."

[0080] Optionally, steps S10401-S10404 are more specific steps than step S104.

[0081] Step S10401: Based on the layered scale relationship, calculate the theoretical total displacement of the large hook corresponding to the original coded curve segment;

[0082] Step S10402: Divide the encoding value of the original encoding curve segment into several consecutive encoding increments, and calculate the physical displacement corresponding to each encoding increment based on the hierarchical scale relationship.

[0083] Step S10403: Calculate the ideal time allocation ratio corresponding to each coding increment based on the proportion of each physical displacement in the total theoretical displacement;

[0084] Step S10404: Based on the ideal time allocation ratio corresponding to each coding increment, calculate and match the corrected time points for the coding values ​​on the original coding curve segment to obtain the target coding curve;

[0085] Specifically, the following example will be used to explain steps S10401-S10404. First, the system executes step S10401 to calculate the theoretical total displacement.

[0086] Suppose a valid motion time interval is captured, with a total duration of 10 seconds and a starting time of 120 seconds. The starting code value of the original coded curve segment corresponding to this interval is E_start = 100000, and the ending code value is E_end = 100003. The system calls the hierarchical scale relationship D = f(E) to query or calculate: f(100000) = 10.000 meters; f(100003) = 10.334 meters. Therefore, the theoretical total displacement of this interval is: 10.334 - 10.000 = 0.334 meters.

[0087] The system then executes step S10402 to calculate the physical displacement corresponding to each encoding increment.

[0088] The encoded value range [100000, 100003] is divided into three consecutive encoded increments of size 1. Based on the hierarchical scale relationship model D=f(E), the system calculates the microscopic physical displacement corresponding to each increment. Since the slope increases on the outer layer of the roller, the displacement corresponding to each increment is not equal.

[0089] The physical displacement corresponding to the first coding increment [100000-100001] is: f(100001)-f(100000)=10.110-10.000=0.110 meters.

[0090] The physical displacement corresponding to the second coding increment [100001-100002] is: f(100002)-f(100001)=10.222-10.110=0.112 meters.

[0091] The physical displacement corresponding to the third encoding increment [100002-100003] is: f(100003) - f(100002) = 10.334 - 10.222 = 0.112 meters. (For simplicity, assume the slope here is the same as in the previous step.)

[0092] Next, step S10403 is executed to calculate the ideal time allocation percentage.

[0093] Based on the calculation results from the previous step, the system determines the proportion of time that should be allocated to each coding increment:

[0094] The time percentage of the first increment is approximately 32.9% (0.110 / 0.334).

[0095] The second increment's time percentage is approximately 33.5% (0.112 / 0.334).

[0096] The third increment's time percentage is approximately 33.5% (0.112 / 0.334).

[0097] Finally, step S10404 is executed to calculate and match the corrected time points.

[0098] The system redistributes the total duration of 10 seconds according to the calculated proportions and calculates a new timestamp for each coded value.

[0099] The time point for the encoded value 100000 is the start time: 120.00 seconds.

[0100] The new time point for the encoded value 100001 is calculated as follows: The time that the first increment should take is 10 seconds × 32.9% = 3.29 seconds. The matching time point is 120.00 seconds + 3.29 seconds = 123.29 seconds.

[0101] The new time point for the encoded value 100002 is calculated as follows: the time that the second increment should take is 10 seconds × 33.5% = 3.35 seconds. The matching time point is 123.29 seconds + 3.35 seconds = 126.64 seconds.

[0102] The new time point for the encoded value 100003 is calculated as follows: the time that the third increment should take is 10 seconds × 33.5% = 3.35 seconds. The matching time point is 126.64 seconds + 3.35 seconds ≈ 130.00 seconds, which is consistent with the total duration.

[0103] Through this refinement process, we calculated and matched a new, corrected time point for each coding value from 100000 to 100003 on the original coding curve, thereby obtaining the corrected target coding curve.

[0104] Of course, the actual effective motion time interval is not limited to three encoded values. The actual effective motion time interval can be divided into encoded increments based on the actual data sampling period.

[0105] Optionally, after step S105, this method may also perform step S109.

[0106] Step S109: Match the time-depth file with the formation data measured by the downhole logging instrument based on the same data timestamp to generate a logging curve with depth on the horizontal axis and formation data on the vertical axis.

[0107] Formation data refers to the sequence of raw physical parameters collected and stored by various measuring sensors (such as natural gamma, resistivity, and acoustic sensors) on downhole logging instruments, synchronized with the instrument's internal clock. Essentially, it is a multidimensional dataset of "time-formation parameters." A logging curve is a two-dimensional graph or data set that visually displays the variation of one or more formation physical parameters with depth as the reference axis. It forms the basis for geological evaluation and reservoir interpretation. Data timestamps refer to the time points corresponding to each data point in the raw physical parameter sequence and the time points corresponding to each data point in the time-depth file.

[0108] Specifically, this step involves data fusion and final result generation. The system first loads the time-depth file generated in step S105 and the formation data file exported from the downhole logging instrument simultaneously. The system can perform point-by-point matching based on the unique correspondence between data points with the same timestamp. For any record in the formation data file (e.g., time: 14:35:10.010, gamma value: 85 API, resistivity: 25 ohm.m), the system will search for the exact same timestamp 14:35:10.010 in the time-depth file and extract its corresponding depth value (e.g., 1500.00 meters). By merging this depth value with the formation parameters in that record, the system generates a new, complete depth-based logging data record (depth: 1500.00 meters, gamma value: 85 API, resistivity: 25 ohm.m). The system repeats this matching process for all formation data records, ultimately generating a final dataset indexed by depth and containing all formation parameters. Based on this final dataset, standard well logging curves for geological analysis can be generated.

[0109] Optionally, steps S10501-S10502 are more specific steps than step S105.

[0110] Step S10501: Construct a mapping relationship between each encoded value of the target encoded curve and the depth range of the drill pipe corresponding to the target encoded curve;

[0111] The depth range refers to a continuous absolute depth range from the starting depth to the ending depth that the corresponding drill pipe actually occupies downhole, such as 950 meters to 980 meters. The depth range is determined by the actual length of the drill pipe and its position in the wellbore. The mapping relationship refers to a mathematical function or table used to associate any coded value on the target coded curve with a specific depth value within the absolute depth range of the corresponding drill pipe, thereby providing a rule for converting the coded value into an actual absolute depth value.

[0112] Specifically, this step aims to establish a quantification relationship so that each coded value on the target coded curve can be accurately assigned an absolute depth value in the wellbore. For example, suppose a drill pipe operation corresponds to an effective movement time interval, which generates a target coded curve segment, and it is known that the drill pipe movement starts at 950 meters downhole and ends at 980 meters (either raised or lowered). Then, the absolute depth interval corresponding to this drill pipe operation is 950 meters to 980 meters. First, the system determines the starting absolute depth (e.g., 950 meters) and ending absolute depth (e.g., 980 meters) of the drill pipe operation, thus deriving a total depth change of 30 meters. Second, the system calculates the total change between the starting and ending coded values ​​of the target coded curve within the effective movement time interval. For example, if the starting coded value is 200000 and the ending coded value is 200005, then the total change in coded values ​​is 5 coded units. Next, the system calculates a depth-to-code-value ratio based on the total depth change of 30 meters and the total change in code values ​​over 5 code units. For example, the ratio is calculated as 30 meters / 5 code units = 6 meters / code unit. This ratio, along with the initial absolute depth, constitutes the mapping relationship between the code values ​​and the absolute depth range within that specific effective operating time interval. This mapping relationship is a localized proportional relationship established for a single drill pipe operating interval. It linearly distributes the changes in code values ​​within that interval to the known absolute depth range, thus providing a specific calculation basis for subsequently converting each code value into an actual absolute depth value.

[0113] Step S10502: Through the mapping relationship, each coding value of the target coding curve is converted into an actual depth value to generate the time depth curve corresponding to the target coding curve;

[0114] The actual depth value refers to the numerical value representing the true absolute vertical position of the instrument in the wellbore, which is converted from the encoded value based on the mapping relationship. It is usually expressed in meters, such as 953.8 meters. The time-depth curve refers to a two-dimensional dataset, with time on the horizontal axis and the actual depth value corresponding to that time on the vertical axis. It is used to describe the trajectory of the logging instrument's downhole position changing with time within a certain effective movement time interval, and is a correspondence between time and depth.

[0115] Specifically, this step aims to convert each coded value on the target coded curve point by point into its corresponding actual absolute depth value based on the established mapping relationship. For example, continuing the example from step S10501, we have established a mapping relationship for the target coded curve (coded value range [200000, 200005]) within a certain effective movement time interval, where the depth-coded value ratio coefficient is 6 meters / coded unit, and the initial absolute depth of the drill pipe operation is 950 meters. For the coded value 200002 appearing on the target coded curve at a certain time t1, the system will apply the mapping relationship for conversion: First, calculate the change in the coded value relative to the initial coded value, i.e., 200002 - 200000 = 2 coded units; then multiply this change in coded value by the previously determined depth-coded value ratio coefficient, i.e., 2 coded units * 6 meters / coded unit = 12 meters. Finally, add this calculated depth increment to the initial absolute depth of the drill pipe operation, i.e., 950 meters + 12 meters = 962 meters. Therefore, at time t1, the actual absolute depth of the instrument is 962 meters. For the coded value 200003 appearing at time t2, its corresponding actual absolute depth will be 950 meters + (200003 - 200000) * 6 meters / coded unit = 950 meters + 3 * 6 meters = 950 meters + 18 meters = 968 meters. By performing this operation on all time points and their corresponding coded values ​​on the target coded curve, converting the coded value at each time point into an actual absolute depth value, the system generates a new dataset. This dataset contains the actual absolute depth value corresponding to each time point (i.e., the time point on the target coded curve) within the effective movement time interval, thus forming a local time-depth curve for this specific drill pipe operation. This local time-depth curve accurately reflects the trajectory of the absolute depth change of the logging instrument during the corresponding drill pipe operation, providing basic data for subsequently stitching together the local time-depth curves of all drill pipe operations into a full-well section time-depth curve according to the drill pipe lowering or raising operation sequence.

[0116] Optionally, steps S10101-S10105 are more specific steps for drawing a tension curve based on the tension data;

[0117] Step S10101: Real-time acquisition of displacement data of the lateral swing of the hook;

[0118] The lateral oscillation of the hook refers to the pendulum-like, non-axial reciprocating motion of the hook suspending the entire drill pipe during the lifting or lowering of the drill string, caused by factors such as wind load, equipment vibration, or unstable operation. The displacement data refers to the physical quantity used to describe this lateral oscillation state. For example, it can be obtained by installing laser displacement sensors, ultrasonic rangefinders, or machine vision systems on the drilling rig or derrick to acquire a real-time sequence of displacement changes of the hook in one or two dimensions of the horizontal plane.

[0119] Specifically, this step aims to acquire raw input data for quantifying the hook's oscillation state. Simultaneously with the logging operation, one or more non-contact displacement sensors are deployed at fixed locations capable of effectively monitoring changes in the hook's horizontal position. These sensors continuously measure the hook's lateral offset relative to its static equilibrium point at a preset sampling frequency and transmit this displacement data, along with a timestamp, to the surface data acquisition system. This sequence of displacement data provides the necessary foundation for subsequent calculations of the additional dynamic forces generated by the oscillation.

[0120] Step S10102: Calculate the oscillation period and oscillation amplitude based on the displacement data;

[0121] The oscillation period refers to the time required for the hook to complete one full reciprocating oscillation, and is a physical quantity describing the oscillation frequency. The oscillation amplitude refers to the maximum distance the hook's horizontal displacement deviates from the equilibrium point during the oscillation process, and is a physical quantity describing the intensity of the oscillation.

[0122] Specifically, this step aims to extract key swing characteristic parameters from the raw displacement data. After receiving the displacement data time series acquired in step S10101, the system performs signal processing. For example, the system can identify the dominant frequency in the signal by performing a Fast Fourier Transform (FFT) or autocorrelation analysis on the displacement series, and thereby calculate the swing period. Simultaneously, the system determines the maximum swing amplitude within a sliding time window by finding local maxima and minima of the displacement data. These two parameters—swing period and swing amplitude—together define the main kinematic characteristics of the hook's swing at the current moment.

[0123] Step S10103: Calculate the additional force of the tension data based on the oscillation period and the oscillation amplitude;

[0124] The additional force refers to the dynamic force component, excluding the static load of the drill pipe, generated in the vertical direction on the tension sensor due to the swing of the hook. This additional force is mainly contributed by the centripetal force generated by the swing and the projection of the tangential acceleration in the vertical direction.

[0125] Specifically, this step is the core of establishing a quantitative relationship between oscillation and tension error. The system is based on a simplified physical model of a simple or compound pendulum, using the oscillation period and amplitude calculated in step S10102 as input parameters. According to the principles of oscillation dynamics, the system can calculate the vertical component of the centripetal force generated by the maximum velocity at the lowest point of the oscillation trajectory, and the additional force generated by the tangential component of gravity and angular acceleration at the highest point of the oscillation trajectory. By modeling the entire oscillation process, the system can estimate a time-varying sequence of additional forces caused by the oscillation.

[0126] Step S10104: When the additional force exceeds the correction threshold, the tension data is compensated and corrected according to the additional force to obtain the target tension data;

[0127] The correction threshold is a pre-set critical value used to determine whether the hook oscillation is severe enough to significantly affect the tension measurement. When the calculated additional force is less than this threshold, its influence is considered negligible. The compensation correction refers to the process of subtracting the additional force calculated in step S10103 from the original tension sensor reading to eliminate oscillation interference. The target tension data refers to higher-quality tension data that, after compensation correction, more purely reflects the axial stress state of the drill pipe.

[0128] Specifically, this step performs actual error correction. The system compares the additional force sequence calculated in step S10103 with the correction threshold in real time. Once the instantaneous value of the additional force exceeds the threshold, the system activates the compensation correction module. This module subtracts the corresponding calculated additional force value from the original tension data at the same timestamp. For example, if the original tension reading at a certain moment is 50.5 tons, and the calculated upward additional force due to violent swaying is 0.2 tons, then the corrected target tension value is 50.5 - 0.2 = 50.3 tons. This correction process is dynamic and is triggered only when the swaying is violent, thus accurately eliminating abnormal interference without affecting normal data.

[0129] Step S10105: Plot the tension curve using the target tension data;

[0130] Specifically, this step is the final output of this optional process. The system no longer uses the raw tension data, which may contain oscillation noise, but instead uses the target tension data processed in step S10104 to perform the tension curve plotting operation described in step S102. The resulting tension curve has a smoother shape and can more accurately reflect the actual axial force changes of the drill pipe. For example, its high-value plateau period will be more stable, and the inflection points of ascent and descent will be clearer, providing higher-quality input for subsequent user interception or automatic identification steps, thereby improving the accuracy and reliability of the entire depth determination method.

[0131] Optionally, steps S10106-S10110 are more specific steps for plotting acceleration curves based on acceleration data;

[0132] Step S10106: Obtain the attitude data of the downhole logging instrument during the lifting or lowering process;

[0133] The attitude data refers to a series of parameters used to describe the orientation and tilt state of the downhole logging instrument in three-dimensional space. This data is typically acquired by attitude sensors (such as a triaxial magnetometer and a triaxial gyroscope) integrated into the instrument along with the accelerometer. The attitude data can be expressed as Euler angles (such as roll, pitch, and yaw) or quaternions.

[0134] Specifically, this step aims to acquire the fundamental data used for coordinate system transformation and gravity compensation. While the downhole logging instrument is acquiring data, its internal attitude sensor synchronously measures the instrument's real-time attitude relative to a geographic coordinate system (e.g., N-E-Ground) or the direction of the gravity field at a preset sampling frequency. This attitude data, along with acceleration data, formation parameter data, etc., shares the same internal clock and is recorded together in the instrument's storage medium, ensuring strict temporal synchronization of all downhole data.

[0135] Step S10107: Calculate the tilt angle and azimuth angle of the downhole logging instrument based on the attitude data;

[0136] The tilt angle refers to the angle between the axis of the downhole logging instrument and its projection in the vertical direction, used to describe the degree of instrument tilt, and its value varies between 0° (vertical) and 90° (horizontal). The azimuth angle refers to the angle between the projection of the axis of the downhole logging instrument on the horizontal plane and the geographic north direction, used to describe the orientation of the instrument.

[0137] Specifically, this step aims to calculate the key angular parameters required for gravity compensation from the raw attitude data. After reading the downhole stored data, the system processes the attitude data acquired in step S10106. For example, if the attitude data is provided by an accelerometer and a magnetometer, the system can calculate the tilt and roll angles relative to the gravitational field using the triaxial readings of the accelerometer, and then, after magnetic declination correction, calculate the azimuth angle relative to geographic north using the triaxial readings of the magnetometer. These two angles—tilt and azimuth—together define the precise orientation of the instrument axis in the geographic coordinate system.

[0138] Step S10108: Perform three-dimensional spatial coordinate transformation on the acceleration data based on the tilt angle and the azimuth angle;

[0139] The three-dimensional spatial coordinate transformation refers to a mathematical process of converting a vector from one coordinate system to another. In this step, the process specifically refers to converting the original triaxial acceleration data measured in the instrument's own coordinate system to a fixed geodetic coordinate system (e.g., the North-East-Earth coordinate system).

[0140] Specifically, the purpose of this step is to unify all acceleration measurements to a fixed reference frame for gravity compensation. Based on the real-time tilt and azimuth angles calculated in step S10107, the system constructs a rotation matrix from the instrument coordinate system to the geodetic coordinate system. Then, for each time point, the three-axis acceleration data vector [ax, ay, az] is transformed into a new vector [a_north, a_east, a_down] in the geodetic coordinate system by left-multiplying it by the inverse (or transpose) of this rotation matrix. After this step, the original acceleration data is freed from its own coordinate system, which was affected by the instrument's rotation and rotation, and its components acquire a clear geographical orientation.

[0141] Step S10109: Perform gravity component compensation on the converted acceleration data to obtain the target acceleration data;

[0142] The gravity component compensation refers to subtracting the projection of the gravity acceleration vector g in the coordinate system from the three-axis acceleration data after coordinate transformation, in order to eliminate the static influence of gravity on acceleration measurement. The target acceleration data refers to the actual motion acceleration of the downhole instrument caused only by external forces (such as the pushing and pulling force of the drill pipe, frictional resistance, etc.) after compensation.

[0143] Specifically, this step is the core of eliminating gravity interference. In the geodetic coordinate system, the direction of the gravitational acceleration vector g is constant, usually represented as [0, 0, +g] (assuming the "ground" direction is positive). The system subtracts this gravity vector from the triaxial acceleration data [a_north, a_east, a_down] converted in step S10108. For example, the compensated Z-axis (vertical direction) acceleration is a_down - g. After this calculation, gravity is accurately eliminated regardless of the inclination or azimuth angle of the downhole instrument. The final target acceleration data can more purely and sensitively reflect the true motion state changes of the instrument.

[0144] Step S10110: Use the target acceleration data to plot the acceleration curve;

[0145] The term "use" refers to using the target acceleration data sequence obtained in step S10109, after attitude correction and gravity compensation, as the final data source for plotting the acceleration curve. Typically, the component most relevant to the instrument axis (e.g., the vertical component in vertical to highly inclined wells, and the horizontal component in horizontal wells) is selected to plot the single-axis acceleration curve.

[0146] Specifically, this step is the final output of this optional process. The system uses the target acceleration data processed in step S10109 to perform the acceleration curve plotting operation described in step S102. The resulting acceleration curve will have a baseline value that is stable near zero and will no longer drift with changes in the instrument's tilt angle. Its vibration peaks can more realistically and clearly reflect the abrupt changes in the instrument's motion state caused by changes in force, resulting in a higher signal-to-noise ratio. This provides a more reliable and clearer basis for subsequent user-selected or automatically identified steps, thereby improving the accuracy of the entire depth determination method.

[0147] The storage-type time-depth logging device in the embodiments of this invention is described below from the perspective of hardware processing. Please refer to [link / reference]. Figure 3 This is a schematic diagram of the physical structure of a storage-type time-depth logging device in the embodiments of this application.

[0148] It should be noted that, Figure 3 The structure of the storage-type time-depth logging device shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of the present invention.

[0149] like Figure 3As shown, the storage-type time-depth logging device includes a CPU 301, which can perform various appropriate actions and processes according to a program stored in the read-only memory ROM 302 or a program loaded from the storage section 308 into the random access memory RAM 303, such as performing the methods described in the above embodiments. The RAM 303 also stores various programs and data required for system operation. The CPU 301, ROM 302, and RAM 303 are interconnected via a bus 304. An I / O interface 305 is also connected to the bus 304.

[0150] The following components are connected to I / O interface 305: input section 306 including audio input devices, push-button switches, etc.; output section 307 including a liquid crystal display (LCD) and audio output devices, indicator lights, etc.; storage section 308 including a hard disk, etc.; and communication section 309 including a network interface card such as a LAN (Local Area Network) card, modem, etc. Communication section 309 performs communication processing via a network such as the Internet. Drive 310 is also connected to I / O interface 305 as needed. Removable media 311, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., are installed on drive 310 as needed so that computer programs read from them can be installed into storage section 308 as needed.

[0151] In particular, according to embodiments of the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a computer-readable medium, the computer program containing computer programs for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication section 309, and / or installed from removable medium 311. When the computer program is executed by CPU 301, it performs the various functions defined in the present invention.

[0152] It should be noted that specific examples of computer-readable storage media may include, but are not limited to: electrical connections having one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), flash memory, optical fiber, portable compact disc read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this invention, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0153] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. Each block in a flowchart or block diagram may represent a module, program segment, or portion of code, which contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those shown in the drawings.

[0154] Specifically, the storage-type time-depth logging device of this embodiment includes a processor and a memory. The memory stores a computer program. When the computer program is executed by the processor, it implements the storage-type time-depth logging method provided in the above embodiment.

[0155] In another aspect, the present invention also provides a computer-readable storage medium, which may be included in the storage-based time-depth logging device described in the above embodiments; or it may exist independently and not assembled into the storage-based time-depth logging device. The storage medium carries one or more computer programs, which, when executed by a processor of the storage-based time-depth logging device, cause the storage-based time-depth logging device to implement the storage-based time-depth logging method provided in the above embodiments.

[0156] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. This program can be stored in a computer-readable storage medium, and when executed, it can include the processes described in the above method embodiments. The aforementioned storage medium includes various media capable of storing program code, such as ROM or random access memory (RAM), magnetic disks, or optical disks.

Claims

1. A storage-based time-depth logging method, characterized in that, Includes the following steps: During the process of raising or lowering the drill pipe, encoder data from the encoder installed on the winch cable drum, tension data from the tension sensor connected to the hook cable, and acceleration data from the accelerometer built into the downhole logging instrument are acquired. The winch cable drum is used to raise or lower the drill pipe. The encoder curve is plotted based on the encoder data, the tension curve is plotted based on the tension data, and the acceleration curve is plotted based on the acceleration data. The user receives a selection of several effective motion time intervals for a single complete drill pipe operation based on the graphic features of the encoding curve, the tension curve, and the acceleration curve. The encoding curve, the tension curve, and the acceleration curve include multiple drill pipe operation event intervals, and each drill pipe operation event interval corresponds to the lifting or lowering process of one drill pipe. For any original encoded curve segment within any effective motion time interval, keep the total change and total duration of the encoded value of the original encoded curve segment unchanged, and re-match the encoded value and time point on the original encoded curve segment according to the layered scale relationship between the encoder data and the actual displacement of the hook to obtain the corrected target encoded curve. All the target encoded curves are converted into multiple time-depth curves according to the depth range corresponding to each drill pipe. All the time-depth curves are then spliced ​​together according to the operation sequence of the corresponding drill pipe operation event range to generate a time-depth file. The time-depth curves represent the correspondence between time and depth. The step of re-matching the encoding values ​​and time points on the original encoding curve segment based on the layered scale relationship between the pre-acquired encoder data and the actual displacement of the hook, to obtain the corrected target encoding curve, specifically includes: Based on the aforementioned layered scale relationship, calculate the theoretical total displacement of the large hook corresponding to the original encoded curve segment; The original encoded curve segment is divided into several consecutive encoded increments, and the physical displacement corresponding to each encoded increment is calculated based on the hierarchical scale relationship. Based on the proportion of each physical displacement in the total theoretical displacement, the ideal time allocation proportion corresponding to each coding increment is calculated; Based on the ideal time allocation ratio corresponding to each coding increment, the coding value on the original coding curve segment is calculated and matched with the corrected time point to obtain the target coding curve. The step of plotting the acceleration curve based on the acceleration data specifically includes: Acquire attitude data of the downhole logging instrument during the lifting or lowering process; The tilt angle and azimuth angle of the downhole logging instrument are calculated based on the attitude data. Perform three-dimensional spatial coordinate transformation on the acceleration data based on the tilt angle and the azimuth angle; Gravity component compensation is performed on the converted acceleration data to obtain the target acceleration data; The acceleration curve is plotted based on the target acceleration data.

2. The method according to claim 1, characterized in that, The step of re-matching the encoding values ​​and time points on the original encoding curve segment based on the layered scale relationship between the pre-acquired encoder data and the actual displacement of the hook to obtain the corrected target encoding curve also includes the following steps before: The winch is operated to continuously raise the hook from the starting height point to the ending height point, and the encoder count value sequence and the real-time hook height sequence are continuously recorded. Align the count value sequence and the real-time height sequence of the large hook in time to generate an original scatter dataset with the encoder count value as the independent variable and the real-time height of the large hook as the dependent variable; Piecewise linear regression or nonlinear curve fitting is performed on the original scatter dataset to construct a piecewise function model for quantifying the nonlinear correspondence between the encoder count value and the actual displacement of the hook. The piecewise function model is used to represent the hierarchical scale relationship.

3. The method according to claim 1, characterized in that, After the step of generating the time-depth file, the method further includes: The time-depth file is matched with the formation data measured by the downhole logging instrument based on the same data timestamp to generate a logging curve with depth on the horizontal axis and formation data on the vertical axis.

4. The method according to claim 1, characterized in that, The step of converting all the target encoded curves into multiple time-depth curves according to the depth range corresponding to each drill pipe specifically includes: For any of the target encoding curves, the time-depth curve is converted by the following steps: Construct a mapping relationship between each encoded value of the target encoded curve and the depth range of the drill pipe corresponding to the target encoded curve; Through the mapping relationship, each encoded value of the target encoding curve is converted into an actual depth value, thereby generating a time-depth curve corresponding to the target encoding curve.

5. The method according to claim 1, characterized in that, The step of plotting the tension curve based on the tension data specifically includes: Real-time acquisition of displacement data of the hook's lateral swing; Calculate the oscillation period and oscillation amplitude based on the displacement data; The additional force of the tension data is calculated based on the oscillation period and the oscillation amplitude. When the additional force exceeds the correction threshold, the tension data is compensated and corrected according to the additional force to obtain the target tension data; The tension curve is plotted based on the target tension data.

Images