Pulse identification method for MWD and device and electronic equipment thereof

Abstract

The present invention discloses a pulse identification method for MWD, an apparatus thereof, and an electronic device thereof, relating to the field of pulse technology. The pulse identification method comprises: receiving pulse signal data, and extracting pulse signal segments of a preset length from the pulse signal data in chronological order based on pulse width; processing each pulse signal segment to obtain a target pulse signal segment; identifying the peak pulse on each target pulse signal segment; and determining the target position of each peak pulse on a pulse signal sequence to obtain an initial pulse sequence; removing erroneous pulses on each target pulse signal segment in the initial pulse sequence based on a preset pulse permutation and combination to obtain a target pulse sequence; and converting the target pulse sequence into a decoded value based on a preset pulse parameter table. The present invention solves the technical problem of low accuracy in pulse identification in the related art.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of pulse technology, in particular to a pulse identification method for MWD and a device thereof and an electronic device. BACKGROUND

[0002] In the process of measurement while drilling (MWD) in oil and gas exploration, the downhole measurement data at the present stage is mostly transmitted in the form of pulse signals relying on mud as a carrier. With the continuous increase of well depth, the pump pressure signal is continuously attenuated, the influence of noise signal is continuously increased, and the decoding accuracy is also decreased. Therefore, in the complex signal, how to improve the signal-to-noise ratio of the pulse signal, quickly identify the effective pulse signal and decode becomes a key problem to be solved in the process of measurement while drilling.

[0003] In the related art, the following methods can be used for pulse identification: (1) the amplitude contrast approximation method is used to obtain the amplitude characteristics of the mud pulse signal, however, this method is sensitive to the noise of the mud pulse signal, if the noise signal in the field is large, it may affect the amplitude of the step-by-step approximation, and may cause distortion of the amplitude, and further cause incorrect identification of the pulse, affecting decoding; (2) a convolutional neural network method is used to realize identification of the mud pulse, however, this method needs to establish a model, and a large amount of model training is performed on the basis of the established model, and the desired effect can be obtained after multiple training, if the quantity of training is not enough or the quality is not good, it may affect the identification and processing of the actual mud pulse signal, affecting decoding; (3) adaptive filtering is used to remove the noise in the mud pulse signal, however, the adaptive filtering is affected by noise and other factors, and the convergence speed will be different, which may deform the actual pulse signal, and further cause deviation in pulse position identification.

[0004] In view of the above problems, no effective solution has been proposed so far. SUMMARY

[0005] The embodiments of the present application provide a pulse identification method for MWD and a device thereof and an electronic device, which at least solve the technical problem of low accuracy in identifying the pulse in the related art.

[0006] According to an aspect of the embodiments of the present application, there is provided a pulse identification method for MWD, comprising: receiving pulse signal data, and extracting pulse signal segments of a preset length from the pulse signal data in time sequence based on pulse width; processing each pulse signal segment to obtain a target pulse signal segment; identifying a peak pulse on each target pulse signal segment and determining a target position of each peak pulse in a pulse signal sequence, to obtain an initial pulse sequence, wherein the pulse signal sequence is a sequence obtained by combining all target pulse signal segments in time sequence; removing false pulses on each target pulse signal segment in the initial pulse sequence based on a preset pulse permutation and combination, to obtain a target pulse sequence, and converting the target pulse sequence into a decoded numerical value based on a preset pulse parameter table.

[0007] Further, the step of processing each pulse signal segment to obtain a target pulse signal segment comprises: performing an open operation on each signal in the pulse signal segment to obtain a first pulse signal segment, and performing a close operation on each signal in the first pulse signal segment to obtain a second pulse signal segment; performing filtering processing on each signal in the second pulse signal segment based on a filter factor array to obtain a third pulse signal segment; sorting signal values of each signal in the third pulse signal segment, and determining a minimum signal value based on the sorting result; taking the minimum signal value as a baseline value, and removing the baseline of each signal in the third pulse signal segment based on the baseline value to obtain a fourth pulse signal segment; performing amplification processing on each signal in the fourth pulse signal segment based on an amplification factor array to obtain the target pulse signal segment.

[0008] Further, the step of performing an open operation on each signal in the pulse signal segment to obtain a first pulse signal segment, and performing a close operation on each signal in the first pulse signal segment to obtain a second pulse signal segment comprises: performing erosion processing on each signal in the pulse signal segment based on a preset circular region to obtain a first erosion pulse signal segment; performing inflation processing on each signal in the first erosion pulse signal segment based on a preset circular region to obtain a first inflation pulse signal segment, and taking the first inflation pulse signal segment as the first pulse signal segment; performing inflation processing on each signal in the first inflation pulse signal segment based on a preset circular region to obtain a second inflation pulse signal segment; performing erosion processing on each signal in the second inflation pulse signal segment based on a preset circular region to obtain a second erosion pulse signal segment, and taking the second erosion pulse signal segment as the second pulse signal segment.

[0009] Further, the signal value of each signal in the first dilated pulse signal segment is a first dilated signal value, and the step of eroding each signal in the pulse signal segment comprises: taking the signal as a current signal, and determining a neighborhood signal set of the current signal, wherein the neighborhood signal set comprises a preset number of signals adjacent to the current signal and located before the current signal, and a preset number of signals adjacent to the current signal and located after the current signal, and the preset number is equal to the radius of the preset circular region; determining a height value of each signal in the neighborhood signal set, and taking a difference value between the current signal value of each signal and the height value as an erosion standard value of each signal; in a case where the current signal value of the signal is greater than or equal to the erosion standard value, taking the erosion standard value as the first dilated signal value of the signal, and in a case where the current signal value of the signal is less than the erosion standard value, taking the current signal value as the first dilated signal value of the signal.

[0010] Further, the signal value of each signal in the first dilated pulse signal segment is a first dilated signal value, and the step of eroding each signal in the pulse signal segment comprises: taking the signal as a current signal, and determining a neighborhood signal set of the current signal, wherein the neighborhood signal set comprises a preset number of signals adjacent to the current signal and located before the current signal, and a preset number of signals adjacent to the current signal and located after the current signal, and the preset number is equal to the radius of the preset circular region; determining a height value of each signal in the neighborhood signal set, and taking a difference value between the current signal value of each signal and the height value as an erosion standard value of each signal; in a case where the current signal value of the signal is greater than or equal to the erosion standard value, taking the erosion standard value as the first dilated signal value of the signal, and in a case where the current signal value of the signal is less than the erosion standard value, taking the current signal value as the first dilated signal value of the signal.

[0011] Further, the step of identifying the peak pulse on each target pulse signal segment comprises: traversing each preset signal on the target pulse signal segment, wherein the preset signal is a signal on the target pulse signal segment other than the beginning signal and the end signal; for each preset signal, determining a previous signal and a next signal adjacent to the preset signal; determining a peak signal value of the preset signal based on the signal value of the preset signal, the signal value of the previous signal adjacent to the preset signal, and the signal value of the next signal adjacent to the preset signal; in a case where the peak signal value is less than 0, determining the preset signal as the peak pulse on the target pulse signal segment, wherein each target pulse signal segment has at most one peak pulse.

[0012] Further, the step of determining the target position of each peak pulse on the pulse signal sequence comprises: identifying a synchronization header on the pulse signal sequence, and taking a pulse signal indicated by a preset synchronization mark on the synchronization header as an initial pulse signal of the pulse signal sequence; determining a length between each peak pulse and the initial pulse signal, and taking the length as an absolute position of the peak pulse; in a case where the absolute position is a floating point value, performing an integer conversion on the floating point value to obtain an integer value of the peak pulse, and determining an error of the peak pulse based on the floating point value and the integer value; determining a prediction error based on the errors of all the peak pulses; and correcting the absolute position of each peak pulse based on the prediction error to obtain the target position of the peak pulse on the pulse signal sequence.

[0013] Further, the preset pulse permutation combination comprises a plurality of preset pulse permutations, and the step of removing the error pulses on each target pulse signal segment in the initial pulse sequence based on the preset pulse permutation combination comprises: determining to-be-decoded pulses on each target pulse signal segment; in a case where the number of the to-be-decoded pulses is greater than a preset decoding number, performing permutation and combination on all the to-be-decoded pulses to obtain a plurality of to-be-decoded pulse permutations; matching each to-be-decoded pulse permutation with the plurality of preset pulse permutations, and determining a to-be-decoded pulse permutation matched with any one preset pulse permutation as a target to-be-decoded pulse permutation; determining, as error pulses, to-be-decoded pulses on the target pulse signal segment except for target to-be-decoded pulses included in the target to-be-decoded pulse permutation, and removing the error pulses on the target pulse signal segment.

[0014] According to another aspect of the embodiment of the present application, there is also provided a pulse identification device for MWD, comprising: an extraction unit configured to receive pulse signal data, and extract pulse signal segments of a preset length from the pulse signal data in chronological order based on pulse widths; a processing unit configured to process each pulse signal segment to obtain target pulse signal segments; an identification unit configured to identify peak pulses on each target pulse signal segment, and determine target positions of each peak pulse on a pulse signal sequence to obtain an initial pulse sequence, wherein the pulse signal sequence is a sequence obtained by combining all the target pulse signal segments in chronological order; and a removal unit configured to remove error pulses on each target pulse signal segment in the initial pulse sequence based on a preset pulse permutation combination to obtain a target pulse sequence, and convert the target pulse sequence into a decoded numerical value based on a preset pulse parameter table.

[0015] Further, the processing unit comprises: a first processing module, configured to perform an opening operation on each signal in the pulse signal segment to obtain a first pulse signal segment, and perform a closing operation on each signal in the first pulse signal segment to obtain a second pulse signal segment; a second processing module, configured to perform filtering processing on each signal in the second pulse signal segment based on the filter factor array to obtain a third pulse signal segment; a first determining module, configured to sort signal values of each signal in the third pulse signal segment, and determine a minimum signal value based on a sorting result; a first removing module, configured to take the minimum signal value as a baseline value, and remove the baseline of each signal in the third pulse signal segment based on the baseline value to obtain a fourth pulse signal segment; and a third processing module, configured to perform amplification processing on each signal in the fourth pulse signal segment based on the amplification factor array to obtain the target pulse signal segment.

[0016] Further, the first processing module comprises: a first erosion module, configured to perform erosion processing on each signal in the pulse signal segment based on the preset circular region to obtain a first erosion pulse signal segment; a first inflation module, configured to perform inflation processing on each signal in the first erosion pulse signal segment based on the preset circular region to obtain a first inflation pulse signal segment, and take the first inflation pulse signal segment as the first pulse signal segment; a second inflation module, configured to perform inflation processing on each signal in the first inflation pulse signal segment based on the preset circular region to obtain a second inflation pulse signal segment; and a second erosion module, configured to perform erosion processing on each signal in the second inflation pulse signal segment based on the preset circular region to obtain a second erosion pulse signal segment, and take the second erosion pulse signal segment as the second pulse signal segment.

[0017] Further, a signal value of each signal in the first inflation pulse signal segment is a first inflation signal value, the first erosion module comprises: a first determining submodule, configured to take the signal as a current signal, and determine a neighbor signal set of the current signal, wherein the neighbor signal set comprises: a preset number of signals adjacent to the current signal and located before the current signal, and a preset number of signals adjacent to the current signal and located after the current signal, and the preset number is equal to a radius of the preset circular region; a second determining submodule, configured to determine a height value of each signal in the neighbor signal set, and take a difference between a current signal value of each signal and the height value as an erosion standard value of each signal; and a first taking submodule, configured to take the erosion standard value as the first inflation signal value of the signal in a case where the current signal value of the signal is greater than or equal to the erosion standard value, and take the current signal value as the first inflation signal value of the signal in a case where the current signal value of the signal is less than the erosion standard value.

[0018] Further, the signal value of each signal in the first erosion pulse signal segment is a first erosion signal value, and the first expansion module comprises: a third determination submodule, configured to take the signal as a current signal, and determine a neighbor signal set of the current signal; a fourth determination submodule, configured to determine a height value of each signal in the neighbor signal set, and take a sum value between the current signal value of each signal and the height value as an expansion standard value of each signal; and a second taking submodule, configured to take the expansion standard value as the first erosion signal value of the signal in a case where the current signal value of the signal is less than or equal to the expansion standard value, and take the current signal value as the first erosion signal value of the signal in a case where the current signal value of the signal is greater than the expansion standard value.

[0019] Further, the recognition unit comprises: a first traversal module, configured to traverse each preset signal on the target pulse signal segment, wherein the preset signal is a signal on the target pulse signal segment except for a beginning signal and an ending signal; a second determination module, configured to determine, for each preset signal, a previous signal and a next signal adjacent to the preset signal; a third determination module, configured to determine, based on a signal value of the preset signal, a signal value of the previous signal adjacent to the preset signal, and a signal value of the next signal adjacent to the preset signal, a peak signal value of the preset signal; and a fourth determination module, configured to determine, in a case where the peak signal value is less than 0, the preset signal as a peak pulse on the target pulse signal segment, wherein each target pulse signal segment has at most one peak pulse.

[0020] Further, the recognition unit further comprises: a first recognition module, configured to recognize a sync head on the pulse signal sequence, and take a pulse signal indicated by a preset sync identifier on the sync head as an initial pulse signal of the pulse signal sequence; a fifth determination module, configured to determine a length between each peak pulse and the initial pulse signal, and take the length as an absolute position of the peak pulse; a sixth determination module, configured to, in a case where the absolute position is a floating-point value, perform an integer conversion on the floating-point value to obtain an integer value of the peak pulse, and determine an error of the peak pulse based on the floating-point value and the integer value; a seventh determination module, configured to determine, based on the errors of all the peak pulses, a prediction error by using a least square method; and a first correction module, configured to correct the absolute position of each peak pulse based on the prediction error to obtain a target position of the peak pulse on the pulse signal sequence.

[0021] Further, the preset pulse arrangement combination includes a plurality of preset pulse arrangements, the removing unit includes an eighth determination module, configured to determine to-be-decoded pulses on each target pulse signal segment, a first sorting module, configured to, in a case where the number of to-be-decoded pulses is greater than the preset decoding number, arrange and combine all to-be-decoded pulses to obtain a plurality of to-be-decoded pulse arrangements, a first matching module, configured to match each to-be-decoded pulse arrangement with the plurality of preset pulse arrangements, and determine a to-be-decoded pulse arrangement matched with any one preset pulse arrangement as a target to-be-decoded pulse arrangement, and a ninth determination module, configured to determine, as error pulses, to-be-decoded pulses on the target pulse signal segment other than target to-be-decoded pulses included in the target to-be-decoded pulse arrangement, and remove the error pulses on the target pulse signal segment.

[0022] According to another aspect of the embodiments of the present application, a computer program product is also provided, including a non-volatile computer readable storage medium, the non-volatile computer readable storage medium storing a computer program, the computer program being executed by a processor to implement the pulse identification method for MWD of any one of the above.

[0023] According to another aspect of the embodiments of the present application, an electronic device is also provided, including one or more processors and a memory, the memory being configured to store one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the pulse identification method for MWD of any one of the above.

[0024] In the present application, pulse signal data is received, and based on pulse width, pulse signal segments of a preset length are extracted from the pulse signal data in time sequence, each pulse signal segment is processed to obtain a target pulse signal segment, peak pulses on each target pulse signal segment are identified, and a target position of each peak pulse on a pulse signal sequence is determined to obtain an initial pulse sequence, based on a preset pulse arrangement combination, error pulses on each target pulse signal segment in the initial pulse sequence are removed to obtain a target pulse sequence, and based on a preset pulse parameter table, the target pulse sequence is converted into a decoding value, thereby solving the technical problem of low accuracy in identifying pulses in the related art.

[0025] In the present application, each pulse signal segment can be processed in time sequence first to remove noise with a large amplitude in the signal and improve the signal-to-noise ratio of the signal, then the accurate position of each peak pulse is identified through the peak value method, and the arrangement combination method is used to exclude error pulses, and then the identified target pulse sequence can be converted into a decoding value, achieving the technical effects of accurately and efficiently identifying pulses and improving pulse decoding accuracy and efficiency. BRIEF DESCRIPTION OF DRAWINGS

[0026] The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiments of the application and together with the description serve to explain the application. In the drawings:

[0027] Figure 1 is a flow chart of an optional pulse identification method for MWD according to an embodiment of the application;

[0028] Figure 2 is a schematic diagram of an optional signal sampling and filtering and identification flow according to an embodiment of the application;

[0029] Figure 3 is a schematic diagram of simulation test results of an optional pulse identification method for MWD according to an embodiment of the application;

[0030] Figure 4 is a schematic diagram of water circulation test results of an optional pulse identification method for MWD according to an embodiment of the application;

[0031] Figure 5 is a schematic diagram of an optional pulse identification device for MWD according to an embodiment of the application;

[0032] Figure 6 is a hardware structure block diagram of an electronic device (or mobile device) for a pulse identification method for MWD according to an embodiment of the application. DETAILED DESCRIPTION

[0033] In order to make the persons skilled in the art better understand the application scheme, the technical scheme in the embodiments of the application will be described clearly and completely below in conjunction with the drawings in the embodiments of the application. Obviously, the described embodiments are only a part of the embodiments of the application, rather than all the embodiments. Based on the embodiments in the application, all other embodiments obtained by those of ordinary skill in the art without creative work should belong to the protection scope of the application.

[0034] It should be noted that the terms "first", "second", and the like in the present application are used to distinguish similar objects, and do not necessarily indicate a specific order or a chronological sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the application described herein can be implemented in an order other than that illustrated or described herein. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but can include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0035] It should be noted that the information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) collected and related to the present application are all information and data authorized by the user or authorized by all parties, and the collection, storage, use, processing, transmission, provision, disclosure and application of related data comply with relevant laws, regulations and standards in the relevant region, necessary security measures are taken, public order and good customs are not violated, and appropriate operation portals are provided for users to choose authorization or refusal. For example, an interface is provided between the system and the related users or agencies, and before obtaining the relevant information, the interface needs to send a request to the aforementioned user or agency, and after receiving the consent information feedback from the aforementioned user or agency, the relevant information is obtained.

[0036] The present application relates to the measurement-while-drilling technology in the field of oil and gas exploration and development, and proposes a pulse identification method. First, the open-close operation method can be used to remove the pulse signals with large amplitude (i.e. pump pressure mutation signals) in the pulse signals, and the low-pass filtering method is used to further process the signals after the open-close operation, so as to retain the signals of the low-frequency part. Then, the sliding median filtering algorithm is used to remove the baseline to retain the change trend of the signals. Then, the long-time slot correlation algorithm is used to improve the signal-to-noise ratio of the pulse signals, so that the effective signals obtain the maximum gain. Then, the peak value method is used to obtain the pulse, and the accumulated error is removed by the error compensation method, so as to ensure that the error of the pulse from the ideal position is within the minimum range. Then, the permutation and combination method is used to exclude the error pulse, and the error pulse signal similar to the normal pulse signal can be excluded. In this way, the effective pulse can be quickly identified, the complexity of the calculation is reduced, and the effective pulse signal can be quickly identified for subsequent pulse decoding.

[0037] The present application will be described in detail below in conjunction with various embodiments.

[0038] Embodiment one

[0039] According to the embodiments of the present application, an embodiment of a pulse identification method for MWD is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0040] Figure 1 An optional flowchart of a pulse identification method for MWD according to an embodiment of the present application is shown in FIG. 1, which includes the following steps: Figure 1

[0041] ​Step S101, receiving pulse signal data, and extracting a preset length of pulse signal segment from the pulse signal data in time sequence based on the pulse width.

[0042] In the embodiment of the present application, in the process of measurement while drilling in oil and gas exploration, the collected downhole measurement data can be processed by the ground processing instrument to obtain analog signals, and the analog signals are converted into digital signals (i.e. pulse signal data), and then transmitted to the host computer through the bus. The host computer avoids missing data through the cache mechanism. Then, according to the actual pulse width, a pulse signal segment of a preset length (consistent with the pulse width, so as to facilitate subsequent signal processing and decoding) is extracted from the pulse signal data in time sequence.

[0043] In some optional embodiments, the data (pulse signal data) received by the host computer is stored in the buffer area. The size of the buffer area is theoretically unlimited and can be unlimitedly stored. The data can be taken out of the buffer area (i.e. extracting data according to the actual slot width (determined according to the pulse width)) in each time slot, and the time of storing data in the buffer area is substantially the same as the time of taking out data. In this way, the data storage and data extraction are strictly separated, which ensures that the buffer area will not increase unlimitedly and there will be no missing of transmission data.

[0044] Exemplarily, the data storage structure is: list={X0, X1, X2, …, Xn-1, Xn}, wherein n represents the length of the collected data, and Xn represents the pulse signal at the nth moment.

[0045] Here, the pulse width refers to the length of time during which the pulse signal lasts. The slot width refers to the length of the signal time window for analysis and decoding, also known as time slot. The setting of the slot width is determined based on the frequency and duration of the pulse signal. For example, if the frequency of the pulse signal is 10 times per second, the slot width can be set to 0.1 second, and a segment of data can be extracted from the signal every 0.1 second for processing to identify and decode the pulse signal.

[0046] Step S102, processing each pulse signal segment to obtain a target pulse signal segment.

[0047] In the embodiment of the present application, processing each pulse signal segment includes open-close operation processing, low-pass filtering processing, sliding median baseline removal processing, and signal-to-noise ratio improvement processing, etc. to obtain a target pulse signal segment. In this way, after filtering out large burr signals and low-pass filtering processing, the low-frequency part of the signal is retained and the high-frequency noise signal is removed. After using the median algorithm to remove the baseline, the signal-to-noise ratio of the signal is improved through the correlation algorithm, which can ensure that the trend of signal change is not destroyed, provide more effective pulse signals for the subsequent, and ensure the effectiveness and accuracy of decoding calculation.

[0048] In the embodiment of the present application, the open-close operation processing refers to removing the spike noise in the pulse signal segment through the open operation, and then restoring the continuity and smoothness of the signal through the close operation, so as to ensure the quality of the signal. The low-pass filtering processing can use a finite impulse response filter (such as a FIR (Finite Impulse Response) low-pass filter) to remove the high-frequency noise component in the signal and retain the low-frequency signal, thereby improving the signal-to-noise ratio. The sliding median baseline removal processing eliminates the baseline drift by finding the median value in a sliding window, and is suitable for processing signals with random noise.

[0049] In step S103, the peak pulse on each target pulse signal segment is identified, and the target position of each peak pulse on the pulse signal sequence is determined, to obtain an initial pulse sequence, wherein the pulse signal sequence is a sequence obtained by combining all target pulse signal segments in time sequence.

[0050] In the embodiment of the present application, the peak pulse on each target pulse signal segment can be identified first, and then the position of each peak pulse on the pulse signal sequence (i.e., a sequence obtained by combining all target pulse signal segments in time sequence) is determined, and the difference between the cumulative error and the predicted pulse is counted, the linear coefficient is calculated using the least square method, and the actual position of the peak pulse is corrected to obtain the target position of each peak pulse on the pulse signal sequence, thereby obtaining the initial pulse sequence (i.e., a sequence containing the accurate position of the peak pulse).

[0051] In step S104, based on the preset pulse permutation combination, the error pulse on each target pulse signal segment in the initial pulse sequence is removed to obtain a target pulse sequence, and the target pulse sequence is converted into a decoded numerical value based on the preset pulse parameter table.

[0052] In the embodiment of the present application, the corrected pulse signal (i.e., the initial pulse sequence) is converted into a signal permutation that can be used for decoding, and the valid pulses within the time slot length are arranged and combined according to the preset pulse permutation combination (i.e., the legal pulse permutation set according to past experience or theory, which is used to exclude error pulses), and the error pulses are excluded to obtain the target pulse sequence. Then, the corresponding pulse data can be converted into the actual attitude parameter (i.e., the decoded numerical value) using the table lookup method according to the preset pulse parameter table (i.e., a table containing the correspondence between the pulse sequence and the specific decoding information).

[0053] In some optional embodiments, the possible pulse combinations of various lengths can be written in a dictionary in advance, and then the arrangement result (i.e., the target pulse sequence) is used to query the dictionary, the calculation parameters required for the current decoding content are queried according to the query result, and then the actual decoding numerical value is converted through proportional calculation and displayed.

[0054] In summary, each pulse signal segment can be processed in time sequence to remove noise with large amplitude in the signal and improve the signal-to-noise ratio of the signal, then the accurate position of each peak pulse is identified by calculating the peak value method, and the wrong pulses are excluded by using the permutation and combination method, and then the identified target pulse sequence can be decoded and converted into a numerical value, so that the technical effects of accurately and efficiently identifying pulses and improving the accuracy and efficiency of pulse decoding are achieved, thereby solving the technical problem of low accuracy of identifying pulses in the related art.

[0055] To improve the accuracy of determining the target pulse signal segment, in the pulse identification method for MWD provided in Embodiment One of the present application, each signal in the pulse signal segment is subjected to open operation processing to obtain a first pulse signal segment, and each signal in the first pulse signal segment is subjected to closed operation processing to obtain a second pulse signal segment; each signal in the second pulse signal segment is subjected to filtering processing based on a filter factor array to obtain a third pulse signal segment; the signal values of each signal in the third pulse signal segment are sorted, and based on the sorting result, a minimum signal value is determined; the minimum signal value is taken as a baseline value, and based on the baseline value, the baseline of each signal in the third pulse signal segment is removed to obtain a fourth pulse signal segment; each signal in the fourth pulse signal segment is subjected to amplification processing based on an amplification factor array to obtain a target pulse signal segment.

[0056] In the embodiment of the present application, according to the actual pulse width, data of a fixed length (i.e., a pulse signal segment) is extracted, the extracted data can be amplified according to the requirements of open and closed operations, then the processed data is subjected to open operation to remove large peaks and burrs in the data, and then the data after open operation processing is shrunk by using closed operation to ensure the smoothness of the data, and then the data is subjected to reduction processing to obtain the data after open and closed operation processing.

[0057] Specifically, the extracted pulse signal segment of the preset length is subjected to open operation processing to remove noise or peak interference with large amplitude in the signal and ensure the continuity and smoothness of the signal. Subsequently, the first pulse signal segment after open operation processing is subjected to closed operation processing to further eliminate small-range signal fluctuations, make the signal shape more smooth, and reduce the misjudgment in subsequent processing.

[0058] In the embodiment of the present application, a finite impulse response (FIR) low-pass filter can be used to filter the second pulse signal segment after open and closed operation processing to obtain a third pulse signal segment. The filter factor array defines the frequency response characteristics of the filter to further eliminate high-frequency noise and retain the low-frequency components in the signal, i.e., the main information of the pulse signal. The processing of low-pass filtering is represented as follows:

[0059] ;

[0060] wherein the input signal is , the unit impulse response of the FIR filter is (i.e. filter factor), the length is k=0, 1, …, n-1, n represents the number of signals, and the output signal is y(k). After processing, the noise of the high-frequency part in the signal can be processed, only the signal and noise of the low-frequency part are retained, and the effectiveness and smoothness of the signal are further ensured.

[0061] In the embodiment of the present application, the signal value y(k) is sorted and the minimum value is found as the baseline value Xbase of the current signal segment in the filtered third pulse signal segment, and then, based on the baseline value, the baseline of each signal in the third pulse signal segment is removed to obtain the fourth pulse signal segment, wherein the formula for removing the baseline is as follows:

[0062] ;

[0063] By removing the baseline, the baseline drift is eliminated, the data offset after the fluctuation of the field pump pressure data is avoided, and a unified data standard is provided for pulse decoding.

[0064] In the embodiment of the present application, after the data (i.e. the fourth pulse signal segment) of is obtained, because the signal still has part of the low-frequency noise signal, the curve does not look very perfect and needs to be further processed to improve the signal-to-noise ratio. Therefore, an amplification processing algorithm is adopted, a preset amplification factor array (which defines the amplification ratio of the signal at different positions or different frequencies, is used to enhance the effective signal and improve the signal-to-noise ratio) is used to further remove the noise of the signal to obtain the target pulse signal segment, wherein the amplification processing of the signal is represented as follows:

[0065] ;

[0066] wherein, represents the amplification factor, and the input data , k=0, 1, …, n-1 are operated to amplify the signal , and improve the usability of the signal.

[0067] In the embodiment, the signal-to-noise ratio of the pulse signal in the measurement while drilling is effectively improved, and the continuity and stability of the signal are ensured. The open-close operation removes large burrs and small fluctuations, the low-pass filter eliminates high-frequency noise, the baseline removal makes the signal fluctuation relatively stable, and the amplification processing further highlights the effective pulse signal. This series of processing makes the pulse signal clearer and more reliable, provides a high-quality signal basis for subsequent pulse recognition and decoding, and improves the accuracy and efficiency of the measurement while drilling.

[0068] In order to remove the pump pressure mutation in the signal, the open-close operation processing is performed on the signal, in the pulse identification method for MWD provided in the embodiment one of the present application, based on the preset circular region, the erosion processing is performed on each signal in the pulse signal segment, to obtain the first erosion pulse signal segment; based on the preset circular region, the inflation processing is performed on each signal in the first erosion pulse signal segment, to obtain the first inflation pulse signal segment, and the first inflation pulse signal segment is taken as the first pulse signal segment; based on the preset circular region, the inflation processing is performed on each signal in the first inflation pulse signal segment, to obtain the second inflation pulse signal segment; based on the preset circular region, the erosion processing is performed on each signal in the second inflation pulse signal segment, to obtain the second erosion pulse signal segment, and the second erosion pulse signal segment is taken as the second pulse signal segment.

[0069] In the embodiment of the present application, a circular data region (i.e. the preset circular region) with a radius of r can be set first to realize the erosion and inflation processing of the signal through the circular region.

[0070] In the embodiment of the present application, the open operation is first performed on the signal, assuming that the original data is A and the preset circular region is B, the erosion processing can be first performed on A by B to obtain data C (i.e. the first erosion pulse signal segment), so as to remove the burr and peak data greater than or equal to r in the data, and then the data C is inflated with a radius of r, so as to restore the processed data to the original position and size, to obtain data D (i.e. the first inflation pulse signal segment). In the processing process, the peak and burr are removed, and then the over-processed data is restored, so as to ensure the change trend of the data.

[0071] Exemplarily, the open operation can be represented as: Wherein, A represents the original data, and B is the circular data region with a radius of r.

[0072] In the embodiment of the present application, the closed operation processing is performed on the data D processed by the open operation, B can be used to perform the inflation processing on D first to obtain data E (i.e. the second inflation pulse signal segment), so as to amplify the smaller fluctuation in the data, and then the erosion processing is performed on the data E, so as to remove the amplified burr or peak data, to ensure that the data after processing is effective and available, and finally data F (i.e. the second erosion pulse signal segment) is obtained.

[0073] Exemplarily, the closed operation can be represented as: Wherein, D represents the original data, and B is the circular data region with a radius of r.

[0074] In this way, the data after the opening and closing operation processing can process the large peaks or burrs in the data, and process the small fluctuations, so that the data becomes smoother, that is, the change trend of the data is ensured, and the fluctuations in the data are processed.

[0075] In the embodiment, the noise and large burrs in the signal are removed by the corrosion processing, then the eroded part is recovered by the expansion processing, the pulse feature is strengthened again by the expansion, and finally the excess peaks introduced in the secondary processing are removed by the corrosion again, so that the signal is clean and the pulse is clear, the signal noise ratio is improved, the pulse signal can still be accurately captured and recognized in a complex environment, and the data quality and decoding accuracy in the MWD process are improved.

[0076] Optionally, the signal value of each signal in the first expansion pulse signal segment is a first expansion signal value, in order to perform accurate corrosion processing, in the pulse identification method for MWD provided in Embodiment One of the application, the signal is taken as a current signal, and a neighborhood signal set of the current signal is determined, wherein the neighborhood signal set includes a preset number of signals adjacent to the current signal and located before the current signal, and a preset number of signals adjacent to the current signal and located after the current signal, and the preset number is equal to the radius of the preset circular region; the height value of each signal in the neighborhood signal set is determined, and the difference between the current signal value of each signal and the height value is taken as a corrosion standard value of each signal; in the case that the current signal value of the signal is greater than or equal to the corrosion standard value, the corrosion standard value is taken as the first expansion signal value of the signal, and in the case that the current signal value of the signal is less than the corrosion standard value, the current signal value is taken as the first expansion signal value of the signal.

[0077] In the embodiment of the application, when processing each signal in the pulse signal segment, a certain signal is first selected as a current signal, and then a neighborhood signal set is determined around the signal. The set not only includes a preset number of signals adjacent to the current signal and located before the current signal (that is, the number of adjacent signals before the current signal is equal to the preset number), but also includes a preset number of signals adjacent to the current signal and located after the current signal (that is, the number of adjacent signals after the current signal is equal to the preset number), and the preset number is equal to the radius value of the preset circular region.

[0078] In some optional embodiments, if the number of signals adjacent to the current signal and located before the current signal is less than the preset number, the previous pulse signal segment adjacent to the pulse signal segment is spliced with the pulse signal segment to obtain the preset number of signals adjacent to the current signal and located before the current signal; if the number of signals adjacent to the current signal and located after the current signal is less than the preset number, the next pulse signal segment adjacent to the pulse signal segment is spliced with the pulse signal segment to obtain the preset number of signals adjacent to the current signal and located after the current signal. If the pulse signal segment is a starting pulse signal segment or an ending pulse signal segment of the collected pulse signal data, the rth signal on the starting pulse signal segment is taken as the current signal, and the last r signals on the ending pulse signal segment are skipped.

[0079] In the embodiment of the application, the height value of each signal in the neighborhood signal set can be calculated. Specifically, each signal in the neighborhood signal set can be numbered in order, ranging from -r to +r, and the height value of the signal can be calculated based on the value corresponding to the number of the signal and the radius. For example, each number i corresponds to a value, and the height value of the signal can be calculated according to the formula The height value of the signal can be obtained, where |i| represents the value corresponding to the number i. Then, the difference between the current signal value and the height value of each signal is taken as the corrosion standard value of each signal. If the current signal value of the signal is greater than or equal to the corrosion standard value, the corrosion standard value is taken as the first dilatation signal value of the signal. This is a signal weakening process, which is equivalent to reducing the amplitude of the signal to the corrosion standard value, thereby removing the spikes or noise in the signal. If the current signal value of the signal is less than the corrosion standard value, the signal value remains unchanged (i.e., the current signal value is taken as the first dilatation signal value of the signal).

[0080] In the embodiment, the corrosion processing process is a dynamic and local signal optimization strategy, which determines the processing method of the signal point by evaluating the environment around the signal point, helps to remove local abnormal high points, that is, noise or spikes, and ensures that the signal is smoother and more continuous. By setting the neighborhood signal set and the corrosion standard value, the signal can be accurately controlled, avoiding unnecessary damage to the effective signal while removing noise. In addition, it can also adapt to different signal environments, because the corrosion standard value is dynamically calculated according to the local signal characteristics, rather than a fixed threshold.

[0081] Optionally, the signal value of each signal in the first erosion pulse signal segment is a first erosion signal value, in order to perform accurate dilation processing, in the pulse identification method for MWD provided by Embodiment One of the present application, the signal is taken as a current signal, and a neighborhood signal set of the current signal is determined; the height value of each signal in the neighborhood signal set is determined, and the sum value between the current signal value of each signal and the height value is taken as a dilation standard value of each signal; in the case where the current signal value of the signal is less than or equal to the dilation standard value, the dilation standard value is taken as the first erosion signal value of the signal, and in the case where the current signal value of the signal is greater than the dilation standard value, the current signal value is taken as the first erosion signal value of the signal.

[0082] In the embodiment of the present application, the signal point to be processed can be selected as a “current signal” first, and a “neighborhood signal set” is determined around this point. This set not only includes the signal points immediately adjacent to and located before and after the current signal (referred to as neighborhood signals), but also the number of these neighborhood signals is equal to the radius of the preset circular region. For example, if the radius of the preset circular region is 3, the neighborhood signal set will include the first 3 and last 3 signal points before and after the current signal.

[0083] Then, the “height value” (i.e. signal intensity or amplitude) of each signal point in the neighborhood signal set is calculated, and the current signal value of each signal is added to the height value of each signal to obtain a sum value as the “dilation standard value” of the signal. This standard value is used for subsequent dilation processing judgment. A high dilation standard value means that the signal point is more susceptible to the influence of dilation processing, and vice versa.

[0084] After that, according to the relationship between the value of the current signal point and the dilation standard value, it is decided whether to perform dilation processing. Specifically, if the value of the current signal point is less than or equal to the dilation standard value, the dilation standard value is taken as the “first erosion signal value”, i.e. dilation processing is performed; if the value of the current signal point is greater than the dilation standard value, the current signal value remains unchanged, i.e. dilation processing is not performed.

[0085] In the embodiment, local peaks (i.e. noise) in the signal are removed by erosion processing, and then the areas in the signal that may have been excessively eroded are restored by dilation processing, maintaining the continuity of the signal profile, so that the interference of noise on signal recognition can be effectively reduced, while the true characteristics of the signal are preserved as much as possible, ensuring the accuracy of subsequent pulse identification.

[0086] In order to improve the accuracy of identifying the peak pulse, in the pulse identification method for MWD provided in Embodiment One of the present application, each preset signal on the target pulse signal segment is traversed, wherein the preset signal is a signal on the target pulse signal segment except the beginning signal and the ending signal; for each preset signal, a previous signal and a next signal adjacent to the preset signal are determined; based on the signal value of the preset signal, the signal value of the previous signal adjacent to the preset signal and the signal value of the next signal, the peak signal value of the preset signal is determined; in the case that the peak signal value is less than 0, the preset signal is determined as the peak pulse on the target pulse signal segment, wherein each target pulse signal segment has at most one peak pulse.

[0087] In the embodiment of the present application, all the signal points except the beginning and ending of the signal segment can be checked one by one. These internal signal points are referred to as “preset signals”. For each preset signal point being checked, the signal points adjacent to it are determined. According to the signal value (i.e. amplitude value) of the current preset signal point, and the signal values of the previous signal and the next signal, the peak signal value of the preset signal is calculated, the peak signal value The calculation formula is as follows:

[0088] ;

[0089] Wherein, represents the signal value of the current preset signal, represents the signal value of the previous signal of the current preset signal, represents the signal value of the next signal of the current preset signal.

[0090] In the embodiment of the present application, if the peak signal value is less than 0, the preset signal is determined as the peak pulse on the target pulse signal segment. If there is no preset signal with a peak signal value less than 0 on the target pulse signal segment, it indicates that the target pulse signal segment has no peak pulse.

[0091] It should be noted that at most one peak pulse can be identified in each target pulse signal segment, so as to ensure the uniqueness and accuracy of pulse identification.

[0092] In the embodiment, by comparing the signal values of the signal points from the preprocessed signal segment, the possible peak pulse is accurately located, which not only can effectively identify the key turning point in the pulse signal, but also can ensure the accuracy of the identification process and the uniqueness of the pulse.

[0093] In order to improve the accuracy of determining the target position of each peak pulse on the pulse signal sequence, in the pulse identification method for MWD provided by the embodiment one of the present application, the preamble is identified on the pulse signal sequence, and a pulse signal indicated by a preset synchronization mark on the preamble is taken as the initial pulse signal of the pulse signal sequence; the length between each peak pulse and the initial pulse signal is determined, and the length is taken as the absolute position of the peak pulse; in the case that the absolute position is a floating point value, the floating point value is converted by rounding, to obtain an integer value of the peak pulse, and based on the floating point value and the integer value, the error of the peak pulse is determined; based on the errors of all peak pulses, the least square method is used to determine the prediction error; based on the prediction error, the absolute position of each peak pulse is corrected, to obtain the target position of the peak pulse on the pulse signal sequence.

[0094] In the embodiment of the present application, the beginning of the received pulse signal data has a preamble, therefore, by identifying the preamble in the pulse signal sequence, the starting position of the pulse signal to be decoded can be determined. The preamble is the starting mark of the signal sequence, has a unique pulse feature or sequence, and is easy to distinguish and locate. After the preamble is found, according to the preset synchronization mark, a specific pulse signal indicated in the preamble is determined as the initial pulse signal. In this way, the reference point for subsequent pulse positioning is accurate, and the decoding error caused by the positioning error of the starting point is avoided.

[0095] Then, for each identified peak pulse, the length, that is, the time interval or the distance of the signal segment, between the peak pulse and the initial pulse signal is calculated. The length value represents the absolute position of the peak pulse relative to the initial pulse signal, and is an important parameter for pulse signal decoding.

[0096] In the embodiment of the present application, the current time slot length is used to make a preliminary pulse value calculation, to calculate the floating point pulse value corresponding to the absolute position of each peak pulse, and then the floating point pulse value is converted by pulse rounding, to obtain the standard ideal position value of the peak pulse (that is, the integer value obtained by conversion, for example, if the floating point pulse value is 9.9, the standard ideal position value is 10). Then, the error between the integer value and the floating point pulse value is calculated. By calculating the difference between the actual position and the predicted position of each pulse, the linear conversion parameters are obtained by using the least square method, to calculate the prediction error of the current pulse. Then, the time slot length used to calculate the pulse value is updated based on the prediction error, to supplement (that is, correct) the absolute position of each peak pulse, and the pulse data supplemented (that is, the target position of each peak pulse on the pulse signal sequence) is used as the actual pulse for decoding.

[0097] In the embodiment, the position deviation commonly existing in the pulse signal sequence is automatically corrected by error calculation and least square method, and the reliability and accuracy of pulse signal decoding are improved.

[0098] Optionally, the preset pulse arrangement combination includes a plurality of preset pulse arrangements. In the pulse identification method for MWD provided in Embodiment One of the present application, the to-be-decoded pulses on each target pulse signal segment are determined. In the case that the number of to-be-decoded pulses is greater than the preset decoding number, the to-be-decoded pulses are arranged and combined to obtain a plurality of to-be-decoded pulse arrangements. Each to-be-decoded pulse arrangement is matched with the plurality of preset pulse arrangements, and the to-be-decoded pulse arrangement matched with any one of the preset pulse arrangements is determined as a target to-be-decoded pulse arrangement. The to-be-decoded pulses other than the target to-be-decoded pulses included in the target to-be-decoded pulse arrangement on the target pulse signal segment are determined as error pulses, and the error pulses on the target pulse signal segment are removed.

[0099] In the embodiment, the to-be-decoded pulses (judged as pulse signals possibly carrying valid information) on each target pulse signal segment can be determined. If the total number of detected to-be-decoded pulses exceeds the preset decoding number, the to-be-decoded pulses are arranged and combined according to the actual number of pulses to be used for decoding, to generate a plurality of possible to-be-decoded pulse arrangements. Here, the preset decoding number is usually set based on the encoding rule and actual application requirements, and the pulses exceeding this number may contain error or redundant information.

[0100] In the embodiment, each to-be-decoded pulse arrangement is matched with the plurality of preset pulse arrangements in the preset pulse arrangement combination. Here, the preset pulse arrangement is a legal pulse sequence template defined in advance according to the encoding protocol. If a to-be-decoded pulse arrangement completely matches any template in the preset pulse arrangement, the to-be-decoded pulse arrangement is taken as a target to-be-decoded pulse arrangement, i.e., a valid sequence conforming to the encoding rule. Then, all pulses not belonging to the target to-be-decoded pulse arrangement are marked as error pulses, and these error pulses are excluded from the target pulse signal segment. In this way, it is ensured that the decoded pulse sequence completely conforms to the encoding protocol, thereby avoiding decoding errors and improving the accuracy and efficiency of decoding.

[0101] In the embodiment, through arrangement and combination and matching verification, not only the pulse sequence conforming to the encoding rule can be identified, but also the error pulses can be effectively excluded, thereby ensuring the accuracy and reliability of the decoding result. In this way, the accuracy of pulse signal decoding is improved, and especially in a complex signal and high noise environment, the valid pulse information can be more efficiently screened out, and the information transmission and decoding capability of the MWD system is enhanced.

[0102] The following will be described in detail in combination with another alternative specific embodiment.

[0103] In the embodiment of the present application, a signal sampling filtering and identification method is provided, Figure 2 is a schematic diagram of an alternative signal sampling filtering and identification process according to the embodiment of the present application, as shown in the figure, comprising the following processes: Figure 2

[0104] (1) Data buffering and extraction: the upper computer collects the pump pressure pulse signal and saves it in the buffer area, and extracts the data according to the actual slot width, and the data saving and data extraction are strictly separated;

[0105] (2) Opening and closing operation and FIR low-pass processing: the extracted data is processed by opening and closing operation to remove the pump pressure mutation, and a finite impulse response (FIR) low-pass filter is used to filter out high-frequency signals and retain low-frequency data;

[0106] (3) Sliding median baseline removal and related processing: the sliding median algorithm is used to remove the baseline of the signal, and related operations are performed to amplify the effective signal to ensure that the signal change can be accurately identified;

[0107] (4) Peak value acquisition, cumulative error calculation and pulse correction: the pulse peak value is calculated, the difference between the cumulative error and the predicted pulse is calculated, the linear coefficient is calculated using the least square method, and the actual pulse signal is corrected;

[0108] (5) Arrangement and combination to exclude error pulses: the corrected pulse signal is converted into a signal arrangement that can be used for decoding, and the pulses in the effective time slot length are arranged and combined to exclude error pulses;

[0109] (6) Combined type lookup table and data conversion: the corresponding pulse data is converted into actual attitude parameters using the lookup table method.

[0110] In the embodiment, the change trend of the signal can be retained while the noise data with large amplitude is removed, so that in the process of low-pass filtering, the pulse signal waveform or the relative position of the pulse after low-pass filtering will not be affected by the noise with large amplitude. Then, the error compensation algorithm is used to compensate the cumulative error of the actual collected pulse, so that the pulse position used for actual decoding is the most reasonable.

[0111] In the embodiment of the present application, the effectiveness of the filtering effect can be verified by using simulation and water circulation test methods.

[0112] (1) Simulation.

[0113] ​The pump pressure data was collected from the actual drilling site, and the original signal was processed using the above-mentioned filtering method in a 1:1 restoration method. Finally, after processing, the effective signal can be identified from the pulse signal with peak changes and various noises. Figure 3 FIG. 1 is a schematic diagram of simulation test results of an optional pulse identification method for MWD according to an embodiment of the present invention, as shown in FIG. Figure 3 As shown, there are respectively displayed the curve graph drawn based on the original data, the curve graph drawn based on the data after opening and closing operation processing, the curve graph drawn based on the data after low-pass processing, the curve graph drawn based on the data after baseline removal processing, and the curve graph drawn based on the data after correlation processing.

[0114] from Figure 3 It can be seen that in the original data, there is a pump pressure signal with a large amplitude change at the peak. After the opening and closing operation, the signal with a large amplitude can be successfully processed, and the signal-to-noise ratio of the signal is further improved after low-pass processing. After baseline removal and correlation processing, the pulse signal has been clearly displayed and has the same form as the pulse in the original data, which is conducive to subsequent pulse decoding.

[0115] (2) Water circulation test.

[0116] The pulse recognition method can be integrated into the host computer decoding software to conduct actual water circulation tests. On-site data collection, real-time filtering and recognition are performed. The water circulation test results show that the pulse recognition method can meet the application scenario. Using this method to filter and recognize pulse signals in real time, a clear and effective pulse signal is obtained. Figure 4 FIG. 1 is a schematic diagram of water circulation test results of an optional pulse identification method for MWD according to an embodiment of the present invention, as shown in FIG. Figure 4 As shown in the figure, it shows the curve diagram after filtering and identification through water circulation test.

[0117] In this embodiment of the present invention, large-amplitude noise in the signal is first removed through opening and closing operations, and high-frequency signals are filtered out using a low-pass filter. A median filter removes the baseline while a correlation algorithm is used to improve the signal-to-noise ratio. The changing trend of the pulse signal after autocorrelation is then used to identify peak values ​​as valid pulses. Invalid pulses are filtered out using a threshold filter. Finite accumulation of pulse positions is used to obtain precise pulse positions, providing a valid signal for decoding and improving decoding efficiency.

[0118] The following describes it in detail with reference to another embodiment.

[0119] Example 2

[0120] The pulse recognition device for MWD provided in the embodiment comprises a plurality of implementation units, each of which corresponds to the implementation steps in Embodiment 1.

[0121] Figure 5 is a schematic diagram of an optional pulse recognition device for MWD according to an embodiment of the present application, as shown in the figure, the pulse recognition device can comprise an extraction unit 50, a processing unit 51, an identification unit 52, and a removal unit 53. Figure 5

[0122] The extraction unit 50 is configured to receive pulse signal data, and extract pulse signal segments of a preset length from the pulse signal data in time sequence based on pulse width.

[0123] The processing unit 51 is configured to process each pulse signal segment to obtain a target pulse signal segment.

[0124] The identification unit 52 is configured to identify a peak pulse on each target pulse signal segment, and determine a target position of each peak pulse in a pulse signal sequence to obtain an initial pulse sequence, wherein the pulse signal sequence is a sequence obtained by combining all target pulse signal segments in time sequence.

[0125] The removal unit 53 is configured to remove an error pulse on each target pulse signal segment in the initial pulse sequence based on a preset pulse permutation and combination to obtain a target pulse sequence, and convert the target pulse sequence into a decoded value based on a preset pulse parameter table.

[0126] The pulse recognition device can first process each pulse signal segment in time sequence to remove noise with a large amplitude in the signal and improve the signal-to-noise ratio of the signal, then identify the accurate position of each peak pulse by calculating the peak value method, and exclude error pulses by using the permutation and combination method, and finally convert the identified target pulse sequence into a decoded value, thereby achieving the technical effects of accurately and efficiently recognizing pulses and improving the accuracy and efficiency of pulse decoding, and further solving the technical problem of low accuracy of recognizing pulses in the related art.

[0127] ​Optionally, the processing unit comprises: a first processing module, configured to perform an opening operation on each signal in the pulse signal segment to obtain a first pulse signal segment, and perform a closing operation on each signal in the first pulse signal segment to obtain a second pulse signal segment; a second processing module, configured to perform filtering processing on each signal in the second pulse signal segment based on the filter factor array to obtain a third pulse signal segment; a first determining module, configured to sort signal values of each signal in the third pulse signal segment, and determine a minimum signal value based on a sorting result; a first removing module, configured to take the minimum signal value as a baseline value, and remove the baseline of each signal in the third pulse signal segment based on the baseline value to obtain a fourth pulse signal segment; and a third processing module, configured to perform amplification processing on each signal in the fourth pulse signal segment based on the amplification factor array to obtain the target pulse signal segment.

[0128] Optionally, the first processing module comprises: a first erosion module, configured to perform erosion processing on each signal in the pulse signal segment based on the preset circular region to obtain a first erosion pulse signal segment; a first inflation module, configured to perform inflation processing on each signal in the first erosion pulse signal segment based on the preset circular region to obtain a first inflation pulse signal segment, and take the first inflation pulse signal segment as the first pulse signal segment; a second inflation module, configured to perform inflation processing on each signal in the first inflation pulse signal segment based on the preset circular region to obtain a second inflation pulse signal segment; and a second erosion module, configured to perform erosion processing on each signal in the second inflation pulse signal segment based on the preset circular region to obtain a second erosion pulse signal segment, and take the second erosion pulse signal segment as the second pulse signal segment.

[0129] Optionally, a signal value of each signal in the first inflation pulse signal segment is a first inflation signal value, the first erosion module comprises: a first determining submodule, configured to take the signal as a current signal, and determine a neighborhood signal set of the current signal, wherein the neighborhood signal set comprises: a preset number of signals adjacent to the current signal and located before the current signal, and a preset number of signals adjacent to the current signal and located after the current signal, and the preset number is equal to the radius of the preset circular region; a second determining submodule, configured to determine a height value of each signal in the neighborhood signal set, and take a difference between a current signal value of each signal and the height value as an erosion standard value of each signal; and a first as submodule, configured to, in a case where the current signal value of the signal is greater than or equal to the erosion standard value, take the erosion standard value as the first inflation signal value of the signal, and in a case where the current signal value of the signal is less than the erosion standard value, take the current signal value as the first inflation signal value of the signal.

[0130] Optionally, the signal value of each signal in the first erosion pulse signal segment is a first erosion signal value, and the first expansion module comprises: a third determination submodule, configured to take the signal as a current signal and determine a neighbor signal set of the current signal; a fourth determination submodule, configured to determine a height value of each signal in the neighbor signal set and take a sum value between the current signal value of each signal and the height value as an expansion standard value of each signal; and a second taking submodule, configured to take the expansion standard value as the first erosion signal value of the signal in a case where the current signal value of the signal is less than or equal to the expansion standard value, and take the current signal value as the first erosion signal value of the signal in a case where the current signal value of the signal is greater than the expansion standard value.

[0131] Optionally, the identification unit comprises: a first traversal module, configured to traverse each preset signal on the target pulse signal segment, wherein the preset signal is a signal on the target pulse signal segment other than the beginning signal and the ending signal; a second determination module, configured to determine, for each preset signal, a previous signal and a next signal adjacent to the preset signal; a third determination module, configured to determine, based on the signal value of the preset signal, the signal value of the previous signal adjacent to the preset signal and the signal value of the next signal, a peak signal value of the preset signal; and a fourth determination module, configured to determine the preset signal as a peak pulse on the target pulse signal segment in a case where the peak signal value is less than 0, wherein each target pulse signal segment has at most one peak pulse.

[0132] Optionally, the identification unit further comprises: a first identification module, configured to identify a sync head on the pulse signal sequence and take a pulse signal indicated by a preset sync identifier on the sync head as an initial pulse signal of the pulse signal sequence; a fifth determination module, configured to determine a length between each peak pulse and the initial pulse signal and take the length as an absolute position of the peak pulse; a sixth determination module, configured to, in a case where the absolute position is a floating-point value, perform an integer conversion on the floating-point value to obtain an integer value of the peak pulse and determine an error of the peak pulse based on the floating-point value and the integer value; a seventh determination module, configured to determine a prediction error by using a least square method based on the errors of all the peak pulses; and a first correction module, configured to correct the absolute position of each peak pulse based on the prediction error to obtain a target position of the peak pulse on the pulse signal sequence.

[0133] Optionally, the preset pulse arrangement combination includes a plurality of preset pulse arrangements, and the removing unit includes an eighth determining module, configured to determine the to-be-decoded pulses on each target pulse signal segment; a first sorting module, configured to, in a case where the number of the to-be-decoded pulses is greater than the preset decoding number, arrange and combine all the to-be-decoded pulses to obtain a plurality of to-be-decoded pulse arrangements; a first matching module, configured to match each to-be-decoded pulse arrangement with the plurality of preset pulse arrangements, and determine the to-be-decoded pulse arrangement matched with any one preset pulse arrangement as a target to-be-decoded pulse arrangement; and a ninth determining module, configured to determine the to-be-decoded pulses other than the target to-be-decoded pulses included in the target to-be-decoded pulse arrangement on the target pulse signal segment as error pulses, and remove the error pulses on the target pulse signal segment.

[0134] The pulse recognition device described above can further include a processor and a memory, and the extraction unit 50, the processing unit 51, the recognition unit 52, and the removing unit 53 are all stored in the memory as program units, and the corresponding functions are realized by the processor executing the program units stored in the memory.

[0135] The processor described above includes a kernel, and the kernel calls the corresponding program units from the memory. The kernel can be one or more, and the error pulses on each target pulse signal segment in the initial pulse sequence are removed based on the preset pulse arrangement combination by adjusting the kernel parameters to obtain a target pulse sequence, and the target pulse sequence is converted into a decoding value based on the preset pulse parameter table.

[0136] The memory described above can include a non-permanent memory in a computer readable medium, a random access memory (RAM), and / or a non-volatile memory such as a read-only memory (ROM) or a flash memory (flash RAM), and the memory includes at least one memory chip.

[0137] The application further provides a computer program product adapted to execute the program of the following method steps when executed on a data processing device: receiving pulse signal data, and extracting pulse signal segments of a preset length from the pulse signal data in chronological order based on pulse width, processing each pulse signal segment to obtain a target pulse signal segment, identifying peak pulses on each target pulse signal segment and determining the target position of each peak pulse in a pulse signal sequence to obtain an initial pulse sequence, removing error pulses on each target pulse signal segment in the initial pulse sequence based on a preset pulse arrangement combination to obtain a target pulse sequence, and converting the target pulse sequence into a decoding value based on a preset pulse parameter table.

[0138] According to another aspect of the embodiments of the present application, there is also provided a computer program product comprising a non-transitory computer readable medium storing a computer program which, when executed by a processor, implements any of the above pulse identification methods for MWD.

[0139] According to another aspect of the embodiments of the present application, there is also provided an electronic device comprising one or more processors and a memory for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to implement any of the above pulse identification methods for MWD.

[0140] Figure 6 is a hardware structural block diagram of an electronic device (or a mobile device) for a pulse identification method for MWD according to an embodiment of the present application. As shown in Figure 6 , the electronic device can include one or more processors (for example, processors 602a, 602b, …, 602n, etc. in the processor 602 of Figure 6 , which can include but are not limited to processing devices such as microprocessors MCU or programmable logic devices FPGA, etc.), a memory 604 for storing data. In addition, it can also include a display, an input / output interface (I / O interface), a universal serial bus (USB) port (which can be included as one of the ports of the I / O interface), a network interface, a keyboard, a power supply and / or a camera. Those skilled in the art can understand that Figure 6 , the structure shown is only schematic, which does not limit the structure of the above-mentioned electronic device. For example, the electronic device can include more or less components than those shown in Figure 6 , or have a different configuration from Figure 6 .

[0141] The above-mentioned serial numbers of the embodiments of the present application are only for description, and do not represent the advantages and disadvantages of the embodiments.

[0142] The embodiments or examples of the present disclosure are not exhaustive, and are only a part of the embodiments or examples, and are not specific limitations on the protection scope of the present disclosure. Each step in a certain embodiment or example can be implemented as an independent example, and the steps can be combined arbitrarily, for example, a scheme after removing some steps in a certain embodiment or example can also be implemented as an independent example, and the order of the steps in a certain embodiment or example can be exchanged arbitrarily, in addition, the optional ways or optional examples in a certain embodiment or example can be combined arbitrarily; in addition, the embodiments or examples can be combined arbitrarily, for example, the steps of different embodiments or examples can be combined arbitrarily, a certain embodiment or example can be combined with the optional ways or optional examples of other embodiments or examples.

[0143] In the above-described embodiments of the present application, the description of each embodiment has its own focus, and the parts not described in detail in a certain embodiment can be referred to the related description of other embodiments.

[0144] In the several embodiments of the present application, it should be understood that the disclosed technology can be implemented in other ways. Of course, the embodiment described above is only a schematic, for example, the division of the units can be a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some interface, unit or module, and can be electrical or other forms.

[0145] The units described as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, that is, they can be located in one place, or can be distributed on multiple units. Part or all of the units can be selected according to actual needs to achieve the purpose of the present embodiment scheme.

[0146] In addition, each functional unit in each embodiment of the present application can be integrated in one processing unit, or each unit can exist physically, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of a software functional unit.

[0147] The integrated unit, if implemented in the form of a software function unit and sold or used as an independent product, can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application, essentially or in other words, the part that contributes to the prior art or the whole or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, including a number of instructions to make a computer device (which can be a personal computer, a server or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present application. The aforementioned storage medium includes: a U disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a mobile hard disk, a magnetic disk or an optical disk, and various media that can store program codes.

[0148] The above is only the preferred embodiment of the present application, it should be pointed out that, for those skilled in the art, without departing from the principles of the present application, can make a number of improvements and refinements, these improvements and refinements should also be considered as the protection scope of the present application.

Claims

1. A pulse identification method for MWD, characterized in that: include: receiving pulse signal data, and extracting pulse signal segments of a preset length from the pulse signal data in a time sequence based on the pulse width; Processing each of the pulse signal segments to obtain a target pulse signal segment; Identify the peak pulse on each target pulse signal segment, and determine the target position of each peak pulse on the pulse signal sequence to obtain an initial pulse sequence, wherein the pulse signal sequence is a sequence of all the target pulse signal segments combined in time order; determine the target position of each peak pulse on the pulse signal sequence, including: identifying the synchronization head on the pulse signal sequence, and using the pulse signal indicated by the preset synchronization mark on the synchronization head as the initial pulse signal of the pulse signal sequence; determine the length between each peak pulse and the initial pulse signal, and use the length as the absolute position of the peak pulse; when the absolute position is a floating-point value, round the floating-point value to obtain the integer value of the peak pulse, and determine the error of the peak pulse based on the floating-point value and the integer value; determine the prediction error based on the errors of all the peak pulses using the least squares method; correct the absolute position of each peak pulse based on the prediction error to obtain the target position of the peak pulse on the pulse signal sequence; Based on the preset pulse arrangement and combination, the erroneous pulses on each target pulse signal segment in the initial pulse sequence are removed to obtain a target pulse sequence, and based on the preset pulse parameter table, the target pulse sequence is converted into a decoding value; wherein, removing the erroneous pulses on each target pulse signal segment in the initial pulse sequence includes: determining the pulses to be decoded on each target pulse signal segment; when the number of the pulses to be decoded is greater than the preset decoding number, arranging and combining all the pulses to be decoded to obtain a plurality of pulse arrangements to be decoded; matching each pulse arrangement to be decoded with a plurality of the preset pulse arrangements, and determining the pulse arrangement to be decoded that matches any one of the preset pulse arrangements as the target pulse arrangement to be decoded; determining the other pulses to be decoded on the target pulse signal segment except the target pulse to be decoded included in the target pulse arrangement to be decoded as the erroneous pulses, and removing the erroneous pulses on the target pulse signal segment.

2. The pulse identification method according to claim 1, characterized in that The step of processing each of the pulse signal segments to obtain a target pulse signal segment comprises: Performing an opening operation on each signal in the pulse signal segment to obtain a first pulse signal segment, and performing a closing operation on each signal in the first pulse signal segment to obtain a second pulse signal segment; Based on the filter factor array, filtering each signal in the second pulse signal segment to obtain a third pulse signal segment; sorting the signal value of each signal in the third pulse signal segment, and determining a minimum signal value based on the sorting result; Taking the minimum signal value as a baseline value, and based on the baseline value, removing the baseline of each signal in the third pulse signal segment to obtain a fourth pulse signal segment; Based on the amplification factor array, each signal in the fourth pulse signal segment is amplified to obtain the target pulse signal segment.

3. The pulse identification method according to claim 2, characterized in that The steps of performing an opening operation on each signal in the pulse signal segment to obtain a first pulse signal segment, and performing a closing operation on each signal in the first pulse signal segment to obtain a second pulse signal segment include: Based on a preset circular area, performing corrosion processing on each signal in the pulse signal segment to obtain a first corrosion pulse signal segment; Based on the preset circular area, each signal in the first corrosion pulse signal segment is expanded to obtain a first expanded pulse signal segment, and the first expanded pulse signal segment is used as the first pulse signal segment; Based on the preset circular area, performing expansion processing on each signal in the first expansion pulse signal segment to obtain a second expansion pulse signal segment; Based on the preset circular area, each signal in the second expansion pulse signal segment is corroded to obtain a second corroded pulse signal segment, and the second corroded pulse signal segment is used as the second pulse signal segment.

4. The pulse identification method according to claim 3, characterized in that: The signal value of each signal in the first expansion pulse signal segment is the first expansion signal value, and the step of performing corrosion processing on each signal in the pulse signal segment includes: The signal is used as a current signal, and a neighborhood signal set of the current signal is determined, wherein the neighborhood signal set includes: a preset number of signals adjacent to the current signal and located before the current signal, and a preset number of signals adjacent to the current signal and located after the current signal, wherein the preset number is equal to the radius of the preset circular area; Determine the height value of each signal in the neighborhood signal set, and use the difference between the current signal value of each signal and the height value as the corrosion standard value of each signal; wherein the height value is calculated by sequentially numbering each signal in the neighborhood signal set, and calculating the height value of each signal based on the value corresponding to the signal number and the radius; the signal value refers to the amplitude value of the signal; When the current signal value of the signal is greater than or equal to the corrosion standard value, the corrosion standard value is used as the first expansion signal value of the signal, and when the current signal value of the signal is less than the corrosion standard value, the current signal value is used as the first expansion signal value of the signal.

5. The pulse identification method according to claim 3, characterized in that The signal value of each signal in the first corrosion pulse signal segment is a first corrosion signal value, and the step of performing expansion processing on each signal in the first corrosion pulse signal segment includes: Taking the signal as a current signal, and determining a neighborhood signal set of the current signal; Determine a height value of each signal in the neighborhood signal set, and use the sum of the current signal value of each signal and the height value as the expansion standard value of each signal; When the current signal value of the signal is less than or equal to the expansion standard value, the expansion standard value is used as the first corrosion signal value of the signal, and when the current signal value of the signal is greater than the expansion standard value, the current signal value is used as the first corrosion signal value of the signal.

6. The pulse identification method according to claim 1, characterized in that The step of identifying the peak pulse on each target pulse signal segment comprises: Traversing each preset signal on the target pulse signal segment, wherein the preset signal is a signal on the target pulse signal segment excluding the beginning signal and the ending signal; For each of the preset signals, determining a previous signal and a next signal adjacent to the preset signal; Determine the peak signal value of the preset signal based on the signal value of the preset signal, the signal value of the previous signal adjacent to the preset signal, and the signal value of the next signal; wherein the calculation formula of the peak signal value F(k) is as follows: F(k)=(f(k-1)-f(k))*(f(k)-f(k+1)), f(k) represents the signal value of the current preset signal, f(k-1) represents the signal value of the previous signal of the current preset signal, and f(k+1) represents the signal value of the next signal of the current preset signal; When the peak signal value is less than 0, the preset signal is determined as the peak pulse on the target pulse signal segment, wherein each target pulse signal segment has at most one peak pulse.

7. A pulse identification device for MWD, characterized in that: include: an extraction unit, configured to receive pulse signal data and extract pulse signal segments of a preset length from the pulse signal data in a time sequence based on the pulse width; a processing unit, configured to process each of the pulse signal segments to obtain a target pulse signal segment; an identification unit, configured to identify a peak pulse on each target pulse signal segment, and determine a target position of each peak pulse on a pulse signal sequence to obtain an initial pulse sequence, wherein the pulse signal sequence is a sequence of all target pulse signal segments combined in time order; determining the target position of each peak pulse on a pulse signal sequence comprises: identifying a synchronization header on the pulse signal sequence, and using a pulse signal indicated by a preset synchronization identifier on the synchronization header as an initial pulse signal of the pulse signal sequence; determining a length between each peak pulse and the initial pulse signal, and using the length as an absolute position of the peak pulse; when the absolute position is a floating-point value, rounding the floating-point value to obtain an integer value of the peak pulse, and determining an error of the peak pulse based on the floating-point value and the integer value; determining a prediction error based on the errors of all the peak pulses using a least squares method; and correcting the absolute position of each peak pulse based on the prediction error to obtain the target position of the peak pulse on the pulse signal sequence; A removal unit is used to remove the erroneous pulses on each target pulse signal segment in the initial pulse sequence based on preset pulse permutations and combinations to obtain a target pulse sequence, and convert the target pulse sequence into a decoding value based on a preset pulse parameter table; wherein, removing the erroneous pulses on each target pulse signal segment in the initial pulse sequence includes: determining the pulses to be decoded on each target pulse signal segment; when the number of the pulses to be decoded is greater than the preset decoding number, arranging and combining all the pulses to be decoded to obtain multiple pulse arrangements to be decoded; matching each pulse arrangement to be decoded with multiple preset pulse arrangements, and determining the pulse arrangement to be decoded that matches any one of the preset pulse arrangements as a target pulse arrangement to be decoded; determining the other pulses to be decoded on the target pulse signal segment except the target pulse to be decoded included in the target pulse arrangement to be decoded as the erroneous pulses, and removing the erroneous pulses on the target pulse signal segment.

8. An electronic device, characterized in that: The device comprises one or more processors and a memory, wherein the memory is used to store one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the pulse identification method for MWD according to any one of claims 1 to 6.

Images