Vibration signal acquisition method, device and equipment and computer readable storage medium

Abstract

The application discloses an acquisition method and device of a vibration signal, equipment and a computer readable storage medium, and belongs to the technical field of computers. The method comprises the following steps: acquiring the weight of each sampling point of an optical fiber; determining the first channel signal and the second channel signal of each operation unit of the optical fiber at each sampling time in a first time period according to the weight of each sampling point; determining the first phase angle of each operation unit at each sampling time in the first time period according to the first channel signal and the second channel signal of each operation unit at each sampling time in the first time period; performing phase error elimination processing on the first phase angle of each operation unit at each sampling time in the first time period to obtain the second phase angle of each operation unit at each sampling time in the first time period; and determining the vibration signal of each operation unit in the first time period according to the second phase angle of each operation unit at each sampling time in the first time period. The accuracy of the determined vibration signal is improved.

Description

Technical Field

[0001] This application relates to the field of computer technology, and in particular to a method, apparatus, device, and computer-readable storage medium for acquiring vibration signals. Background Technology

[0002] Distributed fiber optic vibration sensing technology utilizes optical fibers laid in the same trench as buried pipelines as vibration sensors to acquire vibration signals from the buried pipelines. Distributed fiber optic vibration sensing systems are divided into intensity demodulation type and phase demodulation type. Among them, the phase demodulation type is widely used.

[0003] However, the vibration signal of the optical fiber acquired by the phase demodulation type distributed optical fiber vibration sensing system has a relatively serious phase fading, resulting in poor accuracy of the acquired vibration signal and thus low accuracy of the vibration signal of the buried pipeline. Summary of the Invention

[0004] This application provides a method, apparatus, device, and computer-readable storage medium for acquiring vibration signals, which can be used to solve problems in related technologies. The technical solution is as follows:

[0005] On one hand, embodiments of this application provide a method for acquiring vibration signals, the method comprising:

[0006] Obtain the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline;

[0007] Based on the weights of each sampling point, the first channel signal and the second channel signal of each operating unit of the optical fiber are determined for each sampling in the first time period. Each operating unit includes at least one sampling point. The first channel signal of the operating unit at any sampling in the first time period is used to indicate the vibration of the operating unit at any sampling. The second channel signal of the operating unit at any sampling in the first time period is used to indicate the vibration of the operating unit after the phase angle changes to the target angle at any sampling.

[0008] Based on the first channel signal and the second channel signal of each operation unit during each sampling in the first time period, the first phase angle of each operation unit during each sampling in the first time period is determined;

[0009] Phase error elimination processing is performed on the first phase angle of each operation unit during each sampling in the first time period to obtain the second phase angle of each operation unit during each sampling in the first time period;

[0010] Based on the second phase angle of each sampling of each operating unit in the first time period, the vibration signal of each operating unit in the first time period is determined. The vibration signal of any operating unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as the operating unit in the first time period.

[0011] In one possible implementation, obtaining the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline includes:

[0012] For any sampling point, acquire the first channel signal and the second channel signal of that sampling point at each sampling time in the second time period, where the second time period is earlier than the first time period. The first channel signal of that sampling point at any sampling time in the second time period is used to indicate the vibration status of that sampling point at that sampling time, and the second channel signal of that sampling point at any sampling time in the second time period is used to indicate the vibration status of that sampling point after the phase angle changes to the target angle at that sampling time. Based on the first channel signal and the second channel signal of that sampling point at each sampling time in the second time period, determine the complex number of that sampling point at each sampling time in the second time period. The real part of the complex number of a sample point at any sampling time in the second time period is the first channel signal of that sample point at that sampling time, and the imaginary part of the complex number of that sample point at any sampling time in the second time period is the second channel signal of that sample point at that sampling time. Based on the complex number of that sample point at each sampling time in the second time period, a target complex number is determined. The real part of the target complex number is used to indicate the first channel signal of that sample point in the second time period, and the imaginary part of the target complex number is used to indicate the second channel signal of that sample point in the second time period. The angle between the point indicated by the target complex number and the horizontal number axis is determined. Based on the angle, the weight of that sample point is determined.

[0013] In one possible implementation, determining the target complex number based on the complex number at each sampling point in the second time period includes:

[0014] Determine the modulus of the complex number at each sampling point in the second time period; based on the modulus of the complex number at each sampling point in the second time period, determine a reference complex number among the complex numbers at each sampling point in the second time period, wherein the modulus of the reference complex number satisfies the target requirement; take the average of the real parts of the reference complex number as the real part of the target complex number, and take the average of the imaginary parts of the reference complex number as the imaginary part of the target complex number, thereby obtaining the target complex number.

[0015] In one possible implementation, determining a reference complex number from the complex number at each sampling time of the second time period based on the modulus of the complex number at each sampling time of the any sampling point in the second time period includes:

[0016] Arrange the moduli of the complex numbers at each sampling point in the second time period according to the target order; and take the complex number corresponding to the moduli that are placed before the target position in the arranged moduli as the reference complex number.

[0017] In one possible implementation, determining the first channel signal and the second channel signal of each operating unit of the optical fiber at each sampling time in the first time period based on the weights of the respective sampling points includes:

[0018] For any sampling of any of the operation units in the first time period, acquire the first channel signal and the second channel signal of the sampling point included in the operation unit at any sampling of the first time period; determine the first channel weighted value and the second channel weighted value of the sampling point included in the operation unit at any sampling of the first time period based on the weight of the sampling point included in the operation unit and the first channel signal and the second channel signal of the sampling point included in the operation unit at any sampling of the first time period; determine the first channel signal and the second channel signal of the operation unit at any sampling of the first time period based on the first channel weighted value and the second channel weighted value of the sampling point included in the operation unit at any sampling of the first time period.

[0019] In one possible implementation, determining the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time within the first time period, based on the weights of the sampling points included in any operation unit and the first channel signal and the second channel signal of the sampling points included in any operation unit at any sampling time within the first time period, includes:

[0020] For any sampling point included in any operation unit, the product of the weight of any sampling point and the first channel signal of any sampling point in any sampling of the first time period is used as the first channel weight value of any sampling point in any sampling of the first time period; the product of the weight of any sampling point and the second channel signal of any sampling point in any sampling of the first time period is used as the second channel weight value of any sampling point in any sampling of the first time period.

[0021] In one possible implementation, determining the first channel signal and the second channel signal of any operation unit at any sampling time within the first time period based on the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time within the first time period includes:

[0022] The weighted values ​​of the first channel of the sampling points included in any one of the operation units during any one sampling time period are summed to obtain the first channel signal of the operation unit during any one sampling time period; the weighted values ​​of the second channel of the sampling points included in any one of the operation units during any one sampling time period are summed to obtain the second channel signal of the operation unit during any one sampling time period.

[0023] In one possible implementation, determining the first phase angle of each operation unit during each sampling of the first time period based on the first channel signal and the second channel signal during each sampling of the first time period includes:

[0024] For any sampling of any of the operation units in the first time period, the first phase angle of any operation unit at any sampling time in the first time period is obtained by using the arctangent function based on the first channel signal and the second channel signal at any sampling time in the first time period.

[0025] In one possible implementation, the step of performing phase error elimination processing on the first phase angle of each operation unit at each sampling time in the first time period to obtain the second phase angle of each operation unit at each sampling time in the first time period includes:

[0026] Obtain a first phase array including the first phase angle of each operation unit at any sampling time in the first time period; perform linear fitting on the first phase array to obtain a fitted straight line; take the difference between the first phase angle of any operation unit at any sampling time in the first time period and the reference phase angle as the second phase angle of any operation unit at any sampling time in the first time period, and the reference phase angle is the phase angle on the fitted straight line that corresponds to the first phase angle of any operation unit at any sampling time in the first time period.

[0027] In one possible implementation, determining the vibration signal of each operating unit during the first time period based on the second phase angle of each sampling of the respective operating unit during the first time period includes:

[0028] For any of the operation units, a difference sequence corresponding to the operation unit is obtained based on the second phase angle of each sampling of the operation unit in the first time period. The difference sequence includes multiple values, where any value is the difference between the second phase angles of two adjacent samplings. The difference sequence is processed to obtain the vibration signal of the operation unit in the first time period.

[0029] In one possible implementation, processing the differential sequence to obtain the vibration signal of any operating unit in the first time period includes:

[0030] For a second value in the difference sequence that is greater than a first value, a third value is determined based on the second value and the first value, wherein the third value is less than the first value; for a fifth value in the difference sequence that is less than a fourth value, a sixth value is determined based on the fifth value and the fourth value, wherein the sixth value is greater than the fourth value; the third value replaces the second value in the difference sequence, and the sixth value replaces the fourth value in the difference sequence, to obtain the vibration signal of any operating unit in the first time period.

[0031] On the other hand, embodiments of this application provide a vibration signal acquisition device, the device comprising:

[0032] The acquisition module is used to acquire the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline;

[0033] The determining module is used to determine the first channel signal and the second channel signal of each operation unit of the optical fiber at each sampling time in the first time period according to the weight of each sampling point. Each operation unit includes at least one sampling point. The first channel signal of the operation unit at any sampling time in the first time period is used to indicate the vibration of the operation unit at any sampling time. The second channel signal of the operation unit at any sampling time in the first time period is used to indicate the vibration of the operation unit at any sampling time after the phase angle changes to the target angle.

[0034] The determining module is further configured to determine the first phase angle of each operation unit during each sampling in the first time period based on the first channel signal and the second channel signal during each sampling in the first time period.

[0035] The determining module is further configured to perform phase error elimination processing on the first phase angle of each operation unit during each sampling in the first time period to obtain the second phase angle of each operation unit during each sampling in the first time period.

[0036] The determining module is further configured to determine the vibration signal of each operating unit in the first time period based on the second phase angle of each sampling of each operating unit in the first time period, wherein the vibration signal of any operating unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as the operating unit in the first time period.

[0037] In one possible implementation, the acquisition module is configured to, for any sampling point, acquire a first channel signal and a second channel signal of that sampling point at each sampling time in a second time period, wherein the second time period is earlier than the first time period; the first channel signal of the sampling point at any sampling time in the second time period is used to indicate the vibration status of the sampling point at that sampling time; and the second channel signal of the sampling point at any sampling time in the second time period is used to indicate the vibration status of the sampling point after the phase angle changes by a target angle at that sampling time; and, based on the first channel signal and the second channel signal of the sampling point at each sampling time in the second time period, determine the sampling time of the sampling point in the second time period. The complex number is defined as follows: the real part of the complex number of any sampling point during any sampling in the second time period is the first channel signal of the sampling point during that sampling, and the imaginary part of the complex number of any sampling point during any sampling in the second time period is the second channel signal of the sampling point during that sampling. Based on the complex number of any sampling point during each sampling in the second time period, a target complex number is determined. The real part of the target complex number indicates the first channel signal of the sampling point in the second time period, and the imaginary part of the target complex number indicates the second channel signal of the sampling point in the second time period. The angle between the point indicated by the target complex number and the horizontal number axis is determined. Based on the angle, the weight of the sampling point is determined.

[0038] In one possible implementation, the acquisition module is configured to determine the modulus of the complex number at each sampling time of any sampling point in the second time period; determine a reference complex number from the complex numbers at each sampling time of any sampling point in the second time period based on the modulus of the complex number at each sampling time of any sampling point in the second time period, wherein the modulus of the reference complex number satisfies the target requirement; and obtain the target complex number by taking the average of the real parts of the reference complex number as the real part of the target complex number and taking the average of the imaginary parts of the reference complex number as the imaginary part of the target complex number.

[0039] In one possible implementation, the acquisition module is used to arrange the moduli of the complex numbers of any sampling point in each sampling time period according to the target order; and to take the complex number corresponding to the moduli that are placed before the target position in the arranged moduli as the reference complex number.

[0040] In one possible implementation, the determining module is configured to, for any sampling of any of the operation units in a first time period, acquire the first channel signal and the second channel signal of the sampling points included in the any operation unit at any sampling time in the first time period; determine the first channel weighted value and the second channel weighted value of the sampling points included in the any operation unit at any sampling time in the first time period based on the weight of the sampling points included in the any operation unit and the first channel signal and the second channel signal of the sampling points included in the any operation unit at any sampling time in the first time period; and determine the first channel signal and the second channel signal of the any operation unit at any sampling time in the first time period based on the first channel weighted value and the second channel weighted value of the sampling points included in the any operation unit at any sampling time in the first time period.

[0041] In one possible implementation, the determining module is configured to, for any sampling point included in any operation unit, take the product of the weight of the sampling point and the first channel signal of the sampling point at any sampling time in the first time period as the first channel weighted value of the sampling point at any sampling time in the first time period; and take the product of the weight of the sampling point and the second channel signal of the sampling point at any sampling time in the first time period as the second channel weighted value of the sampling point at any sampling time in the first time period.

[0042] In one possible implementation, the determining module is configured to add the weighted values ​​of the first channel of the sampling points included in any one operation unit during any one sampling in the first time period to obtain the first channel signal of the any one operation unit during any one sampling in the first time period; and add the weighted values ​​of the second channel of the sampling points included in any one operation unit during any one sampling in the first time period to obtain the second channel signal of the any one operation unit during any one sampling in the first time period.

[0043] In one possible implementation, the determining module is configured to, for any sampling of any of the operation units in the first time period, obtain the first phase angle of any operation unit at any sampling time in the first time period based on the first channel signal and the second channel signal at any sampling time in the first time period using an arctangent function.

[0044] In one possible implementation, the determining module is configured to acquire a first phase array including the first phase angle of each operation unit at any sampling time in the first time period; perform linear fitting on the first phase array to obtain a fitted straight line; and take the difference between the first phase angle of any operation unit at any sampling time in the first time period and a reference phase angle as the second phase angle of any operation unit at any sampling time in the first time period, wherein the reference phase angle is the phase angle on the fitted straight line corresponding to the first phase angle of any operation unit at any sampling time in the first time period.

[0045] In one possible implementation, the determining module is configured to, for any one of the operation units, obtain a difference sequence corresponding to any one operation unit based on the second phase angle of each sampling of the operation unit in the first time period, wherein the difference sequence includes multiple values, and any one value is the difference between the second phase angles of two adjacent samplings; and process the difference sequence to obtain the vibration signal of the operation unit in the first time period.

[0046] In one possible implementation, the determining module is configured to: for a second value in the difference sequence that is greater than a first value, determine a third value based on the second value and the first value, wherein the third value is less than the first value; for a fifth value in the difference sequence that is less than a fourth value, determine a sixth value based on the fifth value and the fourth value, wherein the sixth value is greater than the fourth value; replace the second value in the difference sequence with the third value, and replace the fourth value in the difference sequence with the sixth value, to obtain the vibration signal of any operating unit in the first time period.

[0047] On the other hand, embodiments of this application provide a computer device, the computer device including a processor and a memory, the memory storing at least one piece of program code, the at least one piece of program code being loaded and executed by the processor, so that the computer device implements any of the vibration signal acquisition methods described above.

[0048] On the other hand, a computer-readable storage medium is also provided, wherein at least one piece of program code is stored in the computer-readable storage medium, the at least one piece of program code being loaded and executed by a processor to enable a computer to implement any of the above-described methods for acquiring vibration signals.

[0049] On the other hand, a computer program or computer program product is also provided, wherein the computer program or computer program product stores at least one computer instruction, which is loaded and executed by a processor to enable the computer to implement any of the above-mentioned methods for acquiring vibration signals.

[0050] The technical solution provided in this application has at least the following beneficial effects:

[0051] The technical solution provided in this application obtains the weights of each sampling point of the optical fiber. When determining the first channel signal and the second channel signal of each operating unit of the optical fiber at each sampling time in the first time period, the weights of each sampling point included in the operating unit are taken into account, making the accuracy of the vibration signal of the operating unit in the first time period more accurate based on the first channel signal and the second channel signal of the operating unit at each sampling time in the first time period. Since the vibration signal of the operating unit in the first time period is the same as the vibration signal of the buried pipeline at the same location as the operating unit in the first time period, the accuracy of the determined vibration signal of the buried pipeline in the first time period is also higher. The vibration signal of the buried pipeline is used for geological exploration of the location of the buried pipeline, thereby making the geological exploration of the location of the buried pipeline more accurate. Attached Figure Description

[0052] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0053] Figure 1 This is a schematic diagram of the implementation environment of a vibration signal acquisition method provided in an embodiment of this application;

[0054] Figure 2 This is a flowchart of a method for acquiring vibration signals provided in an embodiment of this application;

[0055] Figure 3 This is a flowchart illustrating how to determine the weight of any sampling point, as provided in an embodiment of this application.

[0056] Figure 4 This is a schematic diagram showing the angle between a point indicating a target complex number and a horizontal number line, as provided in an embodiment of this application.

[0057] Figure 5 This is a schematic diagram showing a fitted straight line provided in an embodiment of this application;

[0058] Figure 6 This is a schematic diagram of a vibration signal acquisition system for buried pipelines provided in an embodiment of this application;

[0059] Figure 7 This is a schematic diagram of the structure of a vibration signal acquisition device provided in an embodiment of this application;

[0060] Figure 8 This is a schematic diagram of the structure of a terminal device provided in an embodiment of this application;

[0061] Figure 9 This is a schematic diagram of the structure of a server provided in an embodiment of this application. Detailed Implementation

[0062] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0063] It should be noted that the terms "first," "second," etc., used in this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such terms can be used interchangeably where appropriate so that the embodiments of this application described herein can be implemented in orders other than those illustrated or described herein. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0064] Figure 1 This is a schematic diagram illustrating the implementation environment of a vibration signal acquisition method provided in this application embodiment, as shown below. Figure 1 As shown, the implementation environment includes a computer device 101, which can be a terminal device or a server; this embodiment does not limit the specific type of device. The computer device 101 is used to execute the vibration signal acquisition method provided in this embodiment.

[0065] Optionally, computer device 101 is a terminal device. A terminal device can be any electronic device that allows human-computer interaction with a user through one or more methods such as a keyboard, touchpad, remote control, voice interaction, or handwriting device. Examples include PCs (Personal Computers), mobile phones, smartphones, PDAs (Personal Digital Assistants), wearable devices, PPCs (Pocket PCs), tablets, smart car systems, smart TVs, smart speakers, and smartwatches.

[0066] A terminal device can refer to one of multiple terminal devices; this embodiment uses only one terminal device as an example. Those skilled in the art will understand that the number of terminal devices can be more or less. For example, there may be only one terminal device, or there may be dozens or hundreds, or even more. This application embodiment does not limit the number or type of terminal devices.

[0067] When computer device 101 is a server, the server can be a single server, a server cluster consisting of multiple servers, or any of the following: a cloud computing platform or a virtualization center. This application embodiment does not limit this. The server and terminal devices communicate via a wired or wireless network. The server has data receiving, data processing, and data sending functions. Of course, the server may also have other functions, which this application embodiment does not limit.

[0068] Those skilled in the art should understand that the above-described terminal devices and servers are merely illustrative examples. Other existing or future terminal devices or servers that are applicable to this application should also be included within the scope of protection of this application, and are hereby incorporated by reference.

[0069] This application provides a method for acquiring vibration signals, which can be applied to the above-mentioned... Figure 1 The implementation environment shown is as follows: Figure 2 The flowchart shown in this embodiment of the present application illustrates a method for acquiring vibration signals. This method can be implemented by... Figure 1 The computer device 101 in the middle performs the operation. For example... Figure 2 As shown, the method includes the following steps 201 to 205.

[0070] In step 201, the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline are obtained.

[0071] In one possible implementation, the optical fiber laid in the same trench as the buried pipeline includes multiple sampling points. The distance between each sampling point is a first distance, which is set based on experience or adjusted according to the implementation environment; this application embodiment does not limit this. For example, the first distance is 1 meter. The number of sampling points included in the optical fiber is the sum of the quotient of the optical fiber length and the first distance, and 1. For example, if the optical fiber length is 100 meters and the first distance is 1 meter, then the optical fiber includes 101 sampling points.

[0072] Each sampling point corresponds to a weight, and the process for determining the weight of each sampling point is similar. This application embodiment does not limit the process for determining the weight of each sampling point. Optionally, for any sampling point among the sampling points of the optical fiber, the process for determining the weight of any sampling point is as follows: Figure 3 As shown, it includes the following steps 301 to 305.

[0073] Step 301: Obtain the first channel signal and the second channel signal of any sampling point at each sampling time in the second time period. The first channel signal of any sampling point at any sampling time in the second time period is used to indicate the vibration of any sampling point at any sampling time. The second channel signal of any sampling point at any sampling time in the second time period is used to indicate the vibration of any sampling point after the phase angle changes to the target angle at any sampling time.

[0074] The target angle is set based on experience or adjusted according to the implementation environment; this embodiment does not limit this. For example, the target angle is 90 degrees. The second time period can be any time period, and its duration can be any duration; this embodiment does not limit this either. The number of samples in the second time period can be any number; this embodiment does not limit this either. For example, the duration of the second time period is 1 second, and the number of samples in the second time period is 10.

[0075] Step 302: Based on the first channel signal and the second channel signal of any sampling point in each sampling of the second time period, determine the complex number of any sampling point in each sampling of the second time period. The real part of the complex number of any sampling point in the second time period is the first channel signal of any sampling point in any sampling, and the imaginary part of the complex number of any sampling point in the second time period is the second channel signal of any sampling point in any sampling.

[0076] In one possible implementation, when sampling any sampling point at any time during the second time period, the first channel signal and the second channel signal of that sampling point at that time can be obtained. For example, if the first channel signal of any sampling point at any time during the second time period is A and the second channel signal is B, then the complex number of any sampling point at any time during the second time period is A+Bi, where i is the unit of the imaginary part of the complex number.

[0077] Step 303: Determine the target complex number based on the complex number of any sampling point during each sampling in the second time period.

[0078] Optionally, the process of determining the target complex number based on the complex number of any sampling point at each sampling time in the second time period includes: determining the modulus of the complex number of any sampling point at each sampling time in the second time period; determining a reference complex number from the complex numbers of any sampling point at each sampling time in the second time period based on the modulus of the complex number of any sampling point at each sampling time in the second time period, wherein the modulus of the reference complex number satisfies the target requirement; taking the average of the real parts of the reference complex number as the real part of the target complex number, and taking the average of the imaginary parts of the reference complex number as the imaginary part of the target complex number, thus obtaining the target complex number. The real part of the target complex number is used to indicate the first channel signal of any sampling point in the second time period, and the imaginary part of the target complex number is used to indicate the second channel signal of any sampling point in the second time period.

[0079] Optionally, the modulus L of the complex number at any sampling point in the second time period is determined according to the following formula (1).

[0080]

[0081] Where A is the real part of the complex number of any sampling point in any sampling time period, that is, the first channel signal of any sampling point in any sampling time period; B is the imaginary part of the complex number of any sampling point in any sampling time period, that is, the second channel signal of any sampling point in any sampling time period.

[0082] In one possible implementation, the process of determining a reference complex number from the complex numbers of any sampling point in each sampling time period based on the modulus of the complex number of any sampling point in each sampling time period includes: arranging the modulus of the complex numbers of any sampling point in each sampling time period according to the target order; and taking the complex number corresponding to the modulus that ranks before the target position in the arranged modulus as the reference complex number.

[0083] The target order can be from largest to smallest or from smallest to largest; this embodiment does not limit this. The target position can be any position; this embodiment does not limit this either. Optionally, the target position can be determined based on the number of samples in the second time period. For example, the target position is the product of the number of samples in the second time period and the target value. The target value is set based on experience or adjusted according to the implementation environment; this embodiment does not limit this. For example, if the target value is 20% and the number of samples in the second time period is 10, then the target position is 10 * 20% = 2, which means that the complex numbers corresponding to the first two moduli in the arranged moduli are used as reference complex numbers.

[0084] For example, if the target complex numbers are A+Bi and C+Di, then the target complex number is

[0085] Step 304: Determine the angle between the point indicated by the target complex number and the horizontal number line.

[0086] The process of determining the angle between the point indicated by the target complex number and the horizontal number axis includes: using the real part of the target complex number as the abscissa and the imaginary part as the ordinate to obtain the point indicated by the target complex number; determining the point indicated by the target complex number in a two-dimensional coordinate system; and using the angle between the line connecting the point indicated by the target complex number and the origin and the horizontal number axis as the angle between the point indicated by the target complex number and the horizontal number axis. The horizontal number axis is the x-axis in the two-dimensional coordinate system. Figure 4 This is a schematic diagram illustrating the angle between a point indicating a target complex number and a horizontal number line, as provided in an embodiment of this application. Figure 4 In the diagram, 401 is the point indicated by the target complex number, the horizontal number line is the x-axis, and 402 is the angle between the point indicated by the target complex number and the horizontal number line.

[0087] Step 305: Determine the weight of any sampling point based on the included angle.

[0088] Based on the included angle, the weight of any sampling point is determined as follows: Where i is the unit of the imaginary part of the complex number. The included angle.

[0089] In step 202, based on the weight of each sampling point, the first channel signal and the second channel signal of each operation unit of the optical fiber are obtained during each sampling in the first time period. Each operation unit includes at least one sampling point.

[0090] In this embodiment, the first time period is later than the second time period. The number of samples in the first time period and the number of samples in the second time period can be the same or different, and this embodiment does not limit this. The first channel signal of any operating unit at any sampling time in the first time period is used to indicate the vibration of any operating unit at any sampling time, and the second channel signal of any operating unit at any sampling time in the first time period is used to indicate the vibration of any operating unit after the phase angle changes to the target angle at any sampling time.

[0091] The optical fiber includes multiple operating units, each including multiple sampling points. The distance between any two operating units is greater than the distance between any two sampling points, and is an integer multiple of the distance between any two sampling points. The distance between any two operating units is related to the optical pulse width; a pulse width of 200 ns corresponds to a distance of 20 meters. The process of determining the first and second channel signals for each sampling time period of each operating unit is similar. This embodiment only uses the first and second channel signals of any operating unit in any sampling time period as an example for illustration.

[0092] Optionally, for any sampling of any operation unit in the first time period, the first channel signal and the second channel signal of the sampling points included in any operation unit at any sampling of the first time period are acquired; based on the weight of the sampling points included in any operation unit and the first channel signal and the second channel signal of the sampling points included in any operation unit at any sampling of the first time period, the first channel weight value and the second channel weight value of the sampling points included in any operation unit at any sampling of the first time period are determined; based on the first channel weight value and the second channel weight value of the sampling points included in any operation unit at any sampling of the first time period, the first channel signal and the second channel signal of any operation unit at any sampling of the first time period are determined.

[0093] The process of determining the first channel weight value and the second channel weight value of the sampling point included in any operation unit during any sampling time period, based on the weight value of the sampling point included in any operation unit and the first channel signal and the second channel signal of the sampling point included in any operation unit during any sampling time period, includes: for any sampling point included in any operation unit, the product between the weight value of any sampling point and the first channel signal of any sampling point during any sampling time period is taken as the first channel weight value of any sampling point during any sampling time period; the product between the weight value of any sampling point and the second channel signal of any sampling point during any sampling time period is taken as the second channel signal of any sampling point during any sampling time period.

[0094] In one possible implementation, the process of determining the first channel signal and the second channel signal of any operation unit at any sampling time in the first time period, based on the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time in the first time period, includes: adding the first channel weighted values ​​of the sampling points included in any operation unit at any sampling time in the first time period to obtain the first channel signal of any operation unit at any sampling time in the first time period; and adding the second channel weighted values ​​of the sampling points included in any operation unit at any sampling time in the first time period to obtain the second channel signal of any operation unit at any sampling time in the first time period.

[0095] For example, if the weight of any sampling point in any operation unit is Q, and the first channel signal of any sampling point in any sampling of the first time period is M and the second channel signal is N, then the weighted value of the first channel of any sampling point in any sampling of the first time period is Q*M, and the weighted value of the second channel is Q*N. If any operation unit includes 5 sampling points, where the weighted value of the first channel of the first sampling point in any sampling of the first time period is P1 and the weighted value of the second channel is W1; the weighted value of the first channel of the second sampling point in any sampling of the first time period is P2 and the weighted value of the second channel is W2; the weighted value of the first channel of the third sampling point in any sampling of the first time period is P3 and the weighted value of the second channel is W3; the weighted value of the first channel of the fourth sampling point in any sampling of the first time period is P4 and the weighted value of the second channel is W4; and the weighted value of the first channel of the fifth sampling point in any sampling of the first time period is P5 and the weighted value of the second channel is W5. Then, the first channel signal of any operation unit during any sampling in the first time period is P1+P2+P3+P4+P5, and the second channel signal is W1+W2+W3+W4+W5.

[0096] In step 203, the first phase angle of each operation unit during each sampling in the first time period is determined based on the first channel signal and the second channel signal during each sampling in the first time period.

[0097] In one possible implementation, the process of determining the first phase angle for each sampling of each operation unit in the first time period is similar. This embodiment only illustrates the process of determining the first phase angle for any operation unit in any sampling of the first time period. Optionally, for any sampling of any operation unit in the first time period, the first phase angle for any sampling of any operation unit in the first time period is obtained by using an arctangent function based on the first channel signal and the second channel signal of any operation unit in the first time period.

[0098] Optionally, based on the first channel signal and the second channel signal of any operation unit at any sampling time in the first time period, the first phase angle α of any operation unit at any sampling time in the first time period is obtained by the arctangent function according to the following formula (2).

[0099]

[0100] Where U is the first phase angle of any operation unit during any sampling in the first time period, and V is the second channel signal of any operation unit during any sampling in the first time period. arctan is the arctangent function.

[0101] It should be noted that the first phase angle of any operating unit during any sampling in the first time period can also be determined by other means, and the embodiments of this application do not limit this.

[0102] In step 204, the first phase angle of each operation unit during each sampling in the first time period is processed to eliminate phase error, so as to obtain the second phase angle of each operation unit during each sampling in the first time period.

[0103] As the laser propagates along the optical fiber, its phase changes cumulatively, which is not conducive to detecting phase changes. Therefore, it needs to be eliminated. In this embodiment, a fitting algorithm is used to eliminate phase errors.

[0104] In one possible implementation, the process of eliminating phase error by processing the first phase angle of each operation unit at each sampling time in the first time period to obtain the second phase angle of each operation unit at each sampling time in the first time period includes: obtaining a first phase array including the first phase angle of each operation unit at any sampling time in the first time period; performing linear fitting on the first phase array to obtain a fitted straight line; taking the difference between the first phase angle of any operation unit at any sampling time in the first time period and a reference phase angle as the second phase angle of any operation unit at any sampling time in the first time period, where the reference phase angle is the phase angle on the fitted straight line corresponding to the first phase angle of any operation unit at any sampling time in the first time period.

[0105] For example, if the optical fiber includes 10 operating units, then the first phase array includes the first phase angle of the 10 operating units at any sampling time in the first time period, and the first phase array is (α1, α2, α3, α4, α5, α6, α7, α8, α9, α10).

[0106] The process of fitting a straight line to the first phase array to obtain the fitted straight line includes: determining the coordinate points corresponding to each operation unit based on each operation unit and its first phase angle at any sampling time in the first time period. The x-coordinate of the coordinate point corresponding to any operation unit is the position of that operation unit on the optical fiber, and the y-coordinate is the first phase angle of that operation unit at any sampling time in the first time period. The coordinate points corresponding to each operation unit are determined in a two-dimensional coordinate system; the coordinate points corresponding to each operation unit are fitted to obtain the fitted straight line; the y-coordinate of the point on the fitted straight line whose x-coordinate is the position of any operation unit is used as the reference phase angle of that operation unit at any sampling time in the first time period; and the difference between the first phase angle of that operation unit at any sampling time in the first time period and the reference phase angle is used as the second phase angle of that operation unit at any sampling time in the first time period.

[0107] like Figure 5 This is a schematic diagram showing a fitted straight line provided in an embodiment of this application. Figure 5 The white dots in the diagram represent the coordinates of each operation unit. 501 is the fitted straight line. The x-coordinate of the fitted straight line represents the position of any operation unit (the black dots), and the y-coordinate represents the reference phase angle of any operation unit during any sampling in the first time period.

[0108] For example, the first phase angle of operation unit 1 during any sampling in the first time period is α1, and the ordinate of the point of operation unit 1 on the fitted straight line is β. That is, the reference phase angle of operation unit 1 during any sampling in the first time period is β, and the second phase angle of operation unit 1 during any sampling in the first time period is α1-β.

[0109] In step 205, the vibration signal of each operating unit in the first time period is determined based on the second phase angle of each sampling in the first time period. The vibration signal of any operating unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as any operating unit in the first time period.

[0110] In one possible implementation, the process of determining the vibration signal of each operating unit in the first time period is similar. This embodiment only illustrates the process of determining the vibration signal of any one operating unit in the first time period as an example. Optionally, an unwinding algorithm is used to determine the vibration signal of each operating unit in the first time period.

[0111] For any one of the operation units, the process of determining the vibration signal of any one operation unit in the first time period based on the second phase angle of each sampling in the first time period includes: obtaining the difference sequence corresponding to any one operation unit based on the second phase angle of each sampling in the first time period, wherein the difference sequence includes multiple values, and any value is the difference between the second phase angles of two adjacent samplings; and processing the difference sequence to obtain the vibration signal of any one operation unit in the first time period.

[0112] The process of obtaining the difference sequence corresponding to any operation unit based on the second phase angle of each sampling in the first time period includes: obtaining a second phase array including the second phase angle of any operation unit in each sampling in the first time period; determining the difference between two adjacent second phase angles in the second phase array; obtaining the difference sequence corresponding to any operation unit based on the difference between two adjacent second phase angles, wherein the difference sequence corresponding to any operation unit includes the difference between two adjacent second phase angles.

[0113] For example, the first time period was sampled 10 times. The second phase angle of any operation unit during the first sampling of the first time period is γ1, the second phase angle of any operation unit during the second sampling of the first time period is γ2, the second phase angle of any operation unit during the third sampling of the first time period is γ3, the second phase angle of any operation unit during the fourth sampling of the first time period is γ4, the second phase angle of any operation unit during the fifth sampling of the first time period is γ5, the second phase angle of any operation unit during the sixth sampling of the first time period is γ6, the second phase angle of any operation unit during the seventh sampling of the first time period is γ7, the second phase angle of any operation unit during the eighth sampling of the first time period is γ8, the second phase angle of any operation unit during the ninth sampling of the first time period is γ9, and the second phase angle of any operation unit during the tenth sampling of the first time period is γ10. Then the difference sequence corresponding to any operation unit is (γ2-γ1, γ3-γ2, γ4-γ3, γ5-γ4, γ6-γ5, γ7-γ6, γ8-γ7, γ9-γ8, γ10-γ9).

[0114] After obtaining the differential sequence corresponding to any operating unit, the process of processing the differential sequence to obtain the vibration signal of any operating unit in the first time period includes: for the second value in the differential sequence that is greater than the first value, a third value is determined based on the second value and the first value, and the third value is less than the first value; for the fifth value in the differential sequence that is less than the fourth value, a sixth value is determined based on the fifth value and the fourth value, and the sixth value is greater than the fourth value; the third value replaces the second value in the differential sequence, and the sixth value replaces the fourth value in the differential sequence, to obtain the vibration signal of any operating unit in the first time period.

[0115] The first phase angle of each sampling point at any sampling time in the first time period is obtained through the arctangent function. Since the output range of the arctangent function is (-π / 2, π / 2), the first phase angle of each sampling point at any sampling time in the first time period is within the output range of the arctangent function. The first value is the maximum value of the output range of the arctangent function, which is π / 2, and the fourth value is the minimum value of the output range of the arctangent function, which is -π / 2.

[0116] Optionally, based on the second and first values, the third value is determined as X - 2 * Y, where X is the second value and Y is the first value. Based on the fifth and fourth values, the sixth value is determined as Z - 2 * C, where Z is the fifth value and C is the fourth value. Through the above operations, the discontinuities caused by phase periodicity can be gradually eliminated, restoring a continuous phase diagram that more closely approximates the actual physical phenomenon.

[0117] In one possible implementation, the process of replacing the second value in the difference sequence with the third value and the fourth value in the difference sequence to obtain the vibration signal of any operating unit in the first time period includes: replacing the second value in the difference sequence with the third value and the fourth value in the difference sequence with the sixth value to obtain an adjusted difference sequence, and using the adjusted difference sequence as the vibration signal of any operating unit in the first time period. Alternatively, the number of fading events in the adjusted difference sequence is determined; based on the fact that the number of fading events in the adjusted difference sequence is less than a threshold, the adjusted difference sequence is used as the vibration signal of any operating unit in the first time period; based on the fact that the number of fading events in the adjusted difference sequence is not less than the threshold, steps 201 to 205 are re-executed to obtain a difference sequence with a number of fading events less than the threshold, and the difference sequence with a number of fading events less than the threshold is used as the vibration signal of any operating unit in the first time period.

[0118] The number of times threshold is set based on experience or adjusted according to the implementation environment, and this application embodiment does not limit it in this way.

[0119] Optionally, the process of determining the number of fading values ​​in the adjusted difference sequence includes: determining the number of values ​​in the adjusted difference sequence that are greater than the fading threshold, and using the number of values ​​in the adjusted difference sequence that are greater than the fading threshold as the number of fading values ​​in the adjusted difference sequence.

[0120] The fading threshold is set based on experience or adjusted according to the implementation environment; this application does not limit this.

[0121] In one possible implementation, after obtaining the adjusted difference sequence, a binary classification algorithm can be used to process the adjusted difference sequence to obtain the number of fading events. The binary classification algorithm can be a Support Vector Machine (SVM).

[0122] Optionally, the vibration signal of any determined operating unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as the operating unit in the first time period, and the vibration signal of the buried pipeline in the first time period is used to conduct geological exploration around the buried pipeline.

[0123] In one possible implementation, after determining the vibration signal of any operating unit in the first time period, the vibration signal of any operating unit in other time periods can also be determined in the same way as the vibration signal of any operating unit in the first time period, so as to realize the real-time detection of the vibration signal of the buried pipeline.

[0124] The above method obtains the weights of each sampling point in the optical fiber. When determining the first and second channel signals of each operating unit in the first time period, it considers the weights of each sampling point included in the operating unit, thus improving the accuracy of the vibration signal of the operating unit in the first time period determined based on the first and second channel signals of each sampling in the first time period. Since the vibration signal of the operating unit in the first time period is the same as the vibration signal of the buried pipeline at the same location as the operating unit in the first time period, the accuracy of the determined vibration signal of the buried pipeline in the first time period is also improved. The vibration signal of the buried pipeline is used for geological exploration of the location of the buried pipeline, thereby making the geological exploration of the location of the buried pipeline more accurate.

[0125] Figure 6 This is a schematic diagram of a vibration signal acquisition system for buried pipelines provided in an embodiment of this application. Figure 6 In this diagram, 1 represents a narrow-linewidth laser, used to provide stable, ultra-narrow linewidth laser light. 2 represents a pulsed laser modulator, typically an acousto-optic modulator or an electro-optic modulator, which converts continuous light into pulsed light. 3 represents a circulator. 4 represents a photodetector module, which can be categorized into heterodyne detection and homodyne detection based on its structure. 5 represents a data acquisition module, where a data acquisition card performs the analog-to-digital signal conversion. 6 represents the first channel signal. 7 represents the second channel signal. The phase difference between 6 and 7 is 90 degrees. 8 represents an optical fiber, and 10 represents an operating unit on the optical fiber. 11, 12, 13, 14, and 15 represent sampling points on an operating unit.

[0126] Figure 7 The diagram shown is a structural schematic of a vibration signal acquisition device provided in an embodiment of this application. Figure 7 As shown, the device includes:

[0127] The acquisition module 701 is used to acquire the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline.

[0128] The determining module 702 is used to determine the first channel signal and the second channel signal of each operation unit of the optical fiber at each sampling time in the first time period according to the weight of each sampling point. Each operation unit includes at least one sampling point. The first channel signal of any operation unit at any sampling time in the first time period is used to indicate the vibration of any operation unit at any sampling time. The second channel signal of any operation unit at any sampling time in the first time period is used to indicate the vibration of any operation unit after the phase angle changes to the target angle at any sampling time.

[0129] The determining module 702 is further configured to determine the first phase angle of each operation unit during each sampling in the first time period based on the first channel signal and the second channel signal during each sampling in the first time period.

[0130] The determining module 702 is also used to perform phase error elimination processing on the first phase angle of each operation unit during each sampling in the first time period to obtain the second phase angle of each operation unit during each sampling in the first time period.

[0131] The determining module 702 is further configured to determine the vibration signal of each operating unit in the first time period based on the second phase angle of each sampling in the first time period. The vibration signal of any operating unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as any operating unit in the first time period.

[0132] In one possible implementation, the acquisition module 701 is configured to, for any given sampling point, acquire the first channel signal and the second channel signal of that sampling point at each sampling time within a second time period, wherein the second time period is earlier than the first time period. The first channel signal of any sampling point at any sampling time within the second time period is used to indicate the vibration status of that sampling point at any sampling time, and the second channel signal of any sampling point at any sampling time within the second time period is used to indicate the vibration status of that sampling point after the phase angle changes to the target angle at any sampling time. Based on the first channel signal and the second channel signal of any sampling point at each sampling time within the second time period, the module determines the vibration status of that sampling point at the second time period. The complex number is determined for each sampling of a segment. The real part of the complex number of any sampling point in any sampling of the second time segment is the first channel signal of any sampling point in any sampling of the second time segment, and the imaginary part of the complex number of any sampling point in any sampling of the second time segment is the second channel signal of any sampling point in any sampling of the second time segment. Based on the complex number of any sampling point in each sampling of the second time segment, a target complex number is determined. The real part of the target complex number is used to indicate the first channel signal of any sampling point in the second time segment, and the imaginary part of the target complex number is used to indicate the second channel signal of any sampling point in the second time segment. The angle between the point indicated by the target complex number and the horizontal number axis is determined. Based on the angle, the weight of any sampling point is determined.

[0133] In one possible implementation, the acquisition module 701 is used to determine the modulus of the complex number at each sampling time of any sampling point in the second time period; based on the modulus of the complex number at each sampling time of any sampling point in the second time period, a reference complex number is determined from the complex numbers at each sampling time of any sampling point in the second time period, and the modulus of the reference complex number satisfies the target requirement; the average value of the real part of the reference complex number is used as the real part of the target complex number, and the average value of the imaginary part of the reference complex number is used as the imaginary part of the target complex number to obtain the target complex number.

[0134] In one possible implementation, the acquisition module 701 is used to arrange the modulus of the complex number at each sampling point in the second time period according to the target order; and to take the complex number corresponding to the modulus that is placed before the target position in the arranged modulus as the reference complex number.

[0135] In one possible implementation, the determining module 702 is configured to, for any sampling of any operation unit in the first time period, acquire the first channel signal and the second channel signal of the sampling points included in any operation unit at any sampling time in the first time period; determine the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time in the first time period based on the weight of the sampling points included in any operation unit and the first channel signal and the second channel signal of the sampling points included in any operation unit at any sampling time in the first time period; and determine the first channel signal and the second channel signal of any operation unit at any sampling time in the first time period based on the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time in the first time period.

[0136] In one possible implementation, the determining module 702 is configured to, for any sampling point included in any operation unit, take the product between the weight of any sampling point and the first channel signal of any sampling point during any sampling in the first time period as the first channel weight value of any sampling point during any sampling in the first time period; and take the product between the weight of any sampling point and the second channel signal of any sampling point during any sampling in the first time period as the second channel weight value of any sampling point during any sampling in the first time period.

[0137] In one possible implementation, the determining module 702 is used to add the first channel weights of the sampling points included in any operation unit at any sampling time in the first time period to obtain the first channel signal of any operation unit at any sampling time in the first time period; and to add the second channel weights of the sampling points included in any operation unit at any sampling time in the first time period to obtain the second channel signal of any operation unit at any sampling time in the first time period.

[0138] In one possible implementation, the determining module 702 is used to obtain the first phase angle of any operation unit at any sampling time in the first time period based on the first channel signal and the second channel signal at any sampling time in the first time period using an arctangent function.

[0139] In one possible implementation, the determining module 702 is used to obtain a first phase array including the first phase angle of each operation unit at any sampling time in the first time period; to perform linear fitting on the first phase array to obtain a fitted straight line; and to take the difference between the first phase angle of any operation unit at any sampling time in the first time period and a reference phase angle as the second phase angle of any operation unit at any sampling time in the first time period, wherein the reference phase angle is the phase angle on the fitted straight line corresponding to the first phase angle of any operation unit at any sampling time in the first time period.

[0140] In one possible implementation, the determining module 702 is used to, for any one of the operation units, obtain the difference sequence corresponding to any operation unit based on the second phase angle of each sampling of any operation unit in the first time period. The difference sequence includes multiple values, where any value is the difference between the second phase angles of two adjacent samplings. The difference sequence is then processed to obtain the vibration signal of any operation unit in the first time period.

[0141] In one possible implementation, the determining module 702 is used to determine a third value, which is less than the first value, based on the second value and the first value in the difference sequence; to determine a sixth value, which is greater than the fourth value, based on the fifth value and the fourth value in the difference sequence; to replace the second value in the difference sequence with the third value and to replace the fourth value in the difference sequence with the sixth value, thereby obtaining the vibration signal of any operating unit in the first time period.

[0142] The aforementioned device acquires the weights of each sampling point of the optical fiber. When determining the first and second channel signals of each operating unit of the optical fiber during each sampling in the first time period, it considers the weights of each sampling point included in the operating unit, thus improving the accuracy of the vibration signal of the operating unit in the first time period determined based on the first and second channel signals of each sampling in the first time period. Since the vibration signal of the operating unit in the first time period is the same as the vibration signal of the buried pipeline at the same location as the operating unit in the first time period, the accuracy of the determined vibration signal of the buried pipeline in the first time period is also improved. The vibration signal of the buried pipeline is used for geological exploration of the location of the buried pipeline, thereby making the geological exploration of the location of the buried pipeline more accurate.

[0143] It should be understood that the above-described apparatus is only illustrated by the division of the functional modules described above when implementing its functions. In practical applications, the functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. In addition, the apparatus and method embodiments provided in the above embodiments belong to the same concept, and their specific implementation process can be found in the method embodiments, which will not be repeated here.

[0144] Figure 8 This illustration shows a structural block diagram of a terminal device 800 provided in an exemplary embodiment of this application. The terminal device 800 can be any electronic device product capable of human-computer interaction with a user through one or more methods such as a keyboard, touchpad, remote control, voice interaction, or handwriting device. Examples include PCs (Personal Computers), mobile phones, smartphones, PDAs (Personal Digital Assistants), wearable devices, PPCs (Pocket PCs), tablet computers, smart car systems, smart TVs, smart speakers, and smartwatches.

[0145] Typically, terminal device 800 includes a processor 801 and a memory 802.

[0146] Processor 801 may include one or more processing cores, such as a quad-core processor, an octa-core processor, etc. Processor 801 may be implemented using at least one hardware form selected from DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). Processor 801 may also include a main processor and a coprocessor. The main processor, also known as a CPU (Central Processing Unit), is used to process data in the wake-up state; the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, processor 801 may integrate a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed on the screen. In some embodiments, processor 801 may also include an AI (Artificial Intelligence) processor, which is used to handle computational operations related to machine learning.

[0147] The memory 802 may include one or more computer-readable storage media, which may be non-transitory. The memory 802 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices or flash memory devices. In some embodiments, the non-transitory computer-readable storage media in the memory 802 are used to store at least one instruction, which is executed by the processor 801 to implement the vibration signal acquisition method provided in the method embodiments of this application.

[0148] In some embodiments, the terminal device 800 may also optionally include a peripheral device interface 803 and at least one peripheral device. The processor 801, memory 802, and peripheral device interface 803 can be connected via a bus or signal line. Each peripheral device can be connected to the peripheral device interface 803 via a bus, signal line, or circuit board. Specifically, the peripheral device includes at least one of the following: a radio frequency circuit 804, a display screen 805, a camera assembly 806, an audio circuit 807, and a power supply 809.

[0149] Peripheral device interface 803 can be used to connect at least one I / O (Input / Output) related peripheral device to processor 801 and memory 802. In some embodiments, processor 801, memory 802 and peripheral device interface 803 are integrated on the same chip or circuit board; in some other embodiments, any one or two of processor 801, memory 802 and peripheral device interface 803 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.

[0150] The radio frequency (RF) circuit 804 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 804 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 804 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals back into electrical signals. Optionally, the RF circuit 804 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, etc. The RF circuit 804 can communicate with other terminal devices through at least one wireless communication protocol. This wireless communication protocol includes, but is not limited to: the World Wide Web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 804 may also include circuitry related to NFC (Near Field Communication), which is not limited in this application.

[0151] Display screen 805 is used to display a UI (User Interface). This UI may include graphics, text, icons, videos, and any combination thereof. When display screen 805 is a touch display screen, it also has the ability to collect touch signals on or above its surface. These touch signals can be input as control signals to processor 801 for processing. In this case, display screen 805 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, there may be one display screen 805, disposed on the front panel of terminal device 800; in other embodiments, there may be at least two display screens, disposed on different surfaces of terminal device 800 or in a folded design; in still other embodiments, display screen 805 may be a flexible display screen, disposed on a curved or folded surface of terminal device 800. Furthermore, display screen 805 may be configured as a non-rectangular irregular shape, i.e., a non-rectangular screen. Display screen 805 may be made of materials such as LCD (Liquid Crystal Display) or OLED (Organic Light-Emitting Diode).

[0152] The camera assembly 806 is used to acquire images or videos. Optionally, the camera assembly 806 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is located on the front panel of the terminal device 800, and the rear-facing camera is located on the back of the terminal device 800. In some embodiments, there are at least two rear-facing cameras, which are any one of a main camera, a depth-sensing camera, a wide-angle camera, and a telephoto camera, to achieve background blurring by fusion of the main camera and the depth-sensing camera, panoramic shooting by fusion of the main camera and the wide-angle camera, VR (Virtual Reality) shooting, or other fusion shooting functions. In some embodiments, the camera assembly 806 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash refers to a combination of a warm-light flash and a cool-light flash, which can be used for light compensation at different color temperatures.

[0153] The audio circuit 807 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, converting the sound waves into electrical signals that are input to the processor 801 for processing, or input to the radio frequency circuit 804 for voice communication. For stereo sound acquisition or noise reduction purposes, multiple microphones may be used, each located at a different part of the terminal device 800. The microphone may also be an array microphone or an omnidirectional microphone. The speaker is used to convert the electrical signals from the processor 801 or the radio frequency circuit 804 into sound waves. The speaker may be a conventional diaphragm speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can convert electrical signals not only into audible sound waves but also into inaudible sound waves for purposes such as distance measurement. In some embodiments, the audio circuit 807 may also include a headphone jack.

[0154] Power supply 809 is used to supply power to the various components in terminal device 800. Power supply 809 can be AC ​​power, DC power, a disposable battery, or a rechargeable battery. When power supply 809 includes a rechargeable battery, the rechargeable battery can be a wired rechargeable battery or a wireless rechargeable battery. A wired rechargeable battery is a battery that is charged via a wired line, and a wireless rechargeable battery is a battery that is charged via a wireless coil. The rechargeable battery can also be used to support fast charging technology.

[0155] In some embodiments, the terminal device 800 further includes one or more sensors 810. The one or more sensors 810 include, but are not limited to, an accelerometer 811, a gyroscope 812, a pressure sensor 813, an optical sensor 815, and a proximity sensor 816.

[0156] Accelerometer 811 can detect the magnitude of acceleration along the three coordinate axes of a coordinate system established by terminal device 800. For example, accelerometer 811 can be used to detect the components of gravitational acceleration along the three coordinate axes. Processor 801 can control display screen 805 to display the user interface in either a landscape or portrait view based on the gravitational acceleration signal acquired by accelerometer 811. Accelerometer 811 can also be used for games or for acquiring user motion data.

[0157] The gyroscope sensor 812 can detect the orientation and rotation angle of the terminal device 800. The gyroscope sensor 812, in conjunction with the accelerometer sensor 811, can collect 3D motion data from the user on the terminal device 800. Based on the data collected by the gyroscope sensor 812, the processor 801 can perform the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.

[0158] The pressure sensor 813 can be disposed on the side bezel of the terminal device 800 and / or on the lower layer of the display screen 805. When the pressure sensor 813 is disposed on the side bezel of the terminal device 800, it can detect the user's grip signal on the terminal device 800, and the processor 801 can perform left / right hand recognition or quick operation based on the grip signal collected by the pressure sensor 813. When the pressure sensor 813 is disposed on the lower layer of the display screen 805, the processor 801 can control the operable controls on the UI interface based on the user's pressure operation on the display screen 805. The operable controls include at least one of button controls, scroll bar controls, icon controls, and menu controls.

[0159] An optical sensor 815 is used to collect ambient light intensity. In one embodiment, the processor 801 can control the display brightness of the display screen 805 based on the ambient light intensity collected by the optical sensor 815. Specifically, when the ambient light intensity is high, the display brightness of the display screen 805 is increased; when the ambient light intensity is low, the display brightness of the display screen 805 is decreased. In another embodiment, the processor 801 can also dynamically adjust the shooting parameters of the camera assembly 806 based on the ambient light intensity collected by the optical sensor 815.

[0160] The proximity sensor 816, also known as a distance sensor, is typically located on the front panel of the terminal device 800. The proximity sensor 816 is used to detect the distance between the user and the front of the terminal device 800. In one embodiment, when the proximity sensor 816 detects that the distance between the user and the front of the terminal device 800 is gradually decreasing, the processor 801 controls the display screen 805 to switch from a screen-on state to a screen-off state; when the proximity sensor 816 detects that the distance between the user and the front of the terminal device 800 is gradually increasing, the processor 801 controls the display screen 805 to switch from a screen-off state to a screen-on state.

[0161] Those skilled in the art will understand that Figure 8 The structure shown does not constitute a limitation on the terminal device 800, and may include more or fewer components than shown, or combine certain components, or use different component arrangements.

[0162] Figure 9This is a schematic diagram of the server structure provided in the embodiments of this application. The server 900 can vary considerably due to different configurations or performance. It may include one or more Central Processing Units (CPUs) 901 and one or more memories 902. The one or more memories 902 store at least one line of program code, which is loaded and executed by the one or more processors 901 to implement the vibration signal acquisition method provided in the various method embodiments described above. Of course, the server 900 may also have wired or wireless network interfaces, a keyboard, and input / output interfaces for input and output. The server 900 may also include other components for implementing device functions, which will not be elaborated here.

[0163] In an exemplary embodiment, a computer-readable storage medium is also provided, which stores at least one piece of program code that is loaded and executed by a processor to enable a computer to implement any of the above-described methods for acquiring vibration signals.

[0164] Optionally, the aforementioned computer-readable storage medium may be a read-only memory (ROM), a random access memory (RAM), a compact disc read-only memory (CD-ROM), magnetic tape, floppy disk, and optical data storage device, etc.

[0165] In an exemplary embodiment, a computer program or computer program product is also provided, which stores at least one computer instruction, which is loaded and executed by a processor to enable the computer to implement any of the above-described methods for acquiring vibration signals.

[0166] It should be noted that all information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data used for analysis, stored data, displayed data, etc.), and signals involved in this application are authorized by the user or fully authorized by all parties, and the collection, use, and processing of related data must comply with the relevant laws, regulations, and standards of the relevant countries and regions. For example, the first channel signal and the second channel signal of any sampling point in this application at any time period during any sampling are obtained under fully authorized conditions.

[0167] It should be understood that "multiple" as used in this article refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, or B alone. The character " / " generally indicates that the preceding and following related objects have an "or" relationship.

[0168] The sequence numbers of the embodiments in this application are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.

[0169] The above description is merely an exemplary embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the principles of this application should be included within the protection scope of this application.

Claims

1. A method for acquiring vibration signals, characterized in that, The method includes: Obtain the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline; Based on the weights of each sampling point, the first channel signal and the second channel signal of each operating unit of the optical fiber are determined for each sampling in the first time period. Each operating unit includes at least one sampling point. The first channel signal of the operating unit at any sampling in the first time period is used to indicate the vibration of the operating unit at any sampling. The second channel signal of the operating unit at any sampling in the first time period is used to indicate the vibration of the operating unit after the phase angle changes to the target angle at any sampling. Based on the first channel signal and the second channel signal of each operation unit during each sampling in the first time period, the first phase angle of each operation unit during each sampling in the first time period is determined; Phase error elimination processing is performed on the first phase angle of each operation unit during each sampling in the first time period to obtain the second phase angle of each operation unit during each sampling in the first time period; Based on the second phase angle of each sampling of each operation unit in the first time period, the vibration signal of each operation unit in the first time period is determined. The vibration signal of any operation unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as the operation unit in the first time period. The process of obtaining the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline includes: For any sampling point, acquire the first channel signal and the second channel signal of the sampling point at each sampling time in the second time period, wherein the second time period is earlier than the first time period; Based on the first channel signal and the second channel signal of any sampling point in the second time period, the complex number of any sampling point in the second time period is determined. The real part of the complex number of any sampling point in the second time period is the first channel signal of any sampling point in the second time period, and the imaginary part of the complex number of any sampling point in the second time period is the second channel signal of any sampling point in the second time period. A target complex number is determined based on the complex number of any sampling point at each sampling time in the second time period. The real part of the target complex number is used to indicate the first channel signal of any sampling point in the second time period, and the imaginary part of the target complex number is used to indicate the second channel signal of any sampling point in the second time period. Determine the angle between the point indicated by the target complex number and the horizontal number line. ; Based on the included angle, the weight of any sampling point is determined to be... The i is the unit of the imaginary part of the complex number.

2. The method according to claim 1, characterized in that, The first channel signal of any sampling point during any sampling in the second time period is used to indicate the vibration of any sampling point during any sampling. The second channel signal of any sampling point during any sampling in the second time period is used to indicate the vibration of any sampling point after the phase angle changes to the target angle during any sampling.

3. The method according to claim 1, characterized in that, Determining the target complex number based on the complex number of any sampling point at each sampling time in the second time period includes: Determine the modulus of the complex number for each sampling point in the second time period; Based on the modulus of the complex number at each sampling point in the second time period, a reference complex number is determined from the complex numbers at each sampling point in the second time period, and the modulus of the reference complex number satisfies the target requirement. The target complex number is obtained by taking the average of the real parts of the reference complex number as the real part and taking the average of the imaginary parts of the reference complex number as the imaginary part.

4. The method according to claim 3, characterized in that, The step of determining a reference complex number from the complex numbers at each sampling point in the second time period based on the modulus of the complex number at each sampling point in the second time period includes: Arrange the modulus of the complex number of any sampling point in the second time period according to the target order; The complex number corresponding to the modulus that precedes the target position in the arranged modulus is taken as the reference complex number.

5. The method according to any one of claims 1 to 4, characterized in that, The step of determining the first channel signal and the second channel signal of each operating unit of the optical fiber at each sampling time in the first time period based on the weight of each sampling point includes: For any sampling of any of the operation units in the first time period, the first channel signal and the second channel signal of the sampling point included in the operation unit at any sampling of the first time period are obtained. Based on the weight of the sampling point included in any operation unit, and the first channel signal and second channel signal of the sampling point included in any operation unit during any sampling in the first time period, determine the first channel weighted value and second channel weighted value of the sampling point included in any operation unit during any sampling in the first time period. Based on the first channel weighted value and the second channel weighted value of the sampling points included in any one of the operation units during any one sampling in the first time period, the first channel signal and the second channel signal of any one of the operation units during any one sampling in the first time period are determined.

6. The method according to claim 5, characterized in that, The step of determining the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time within the first time period, based on the weights of the sampling points included in any operation unit and the first channel signal and the second channel signal of the sampling points included in any operation unit at any sampling time within the first time period, includes: For any sampling point included in any operation unit, the product of the weight of any sampling point and the first channel signal of any sampling point in any sampling time of the first time period is used as the first channel weight value of any sampling point in any sampling time of the first time period. The product of the weight of any sampling point and the second channel signal of any sampling point in any sampling time of the first time period is used as the second channel weight value of any sampling point in any sampling time of the first time period.

7. The method according to claim 5, characterized in that, The step of determining the first channel signal and the second channel signal of any operation unit at any sampling time within the first time period based on the first channel weighted value and the second channel weighted value of the sampling points included in any operation unit at any sampling time within the first time period includes: The weighted values ​​of the first channel of the sampling points included in any one operation unit during any sampling in the first time period are added together to obtain the first channel signal of any one operation unit during any sampling in the first time period. The weighted values ​​of the second channel of the sampling points included in any one of the operation units during any sampling in the first time period are added together to obtain the second channel signal of any one of the operation units during any sampling in the first time period.

8. The method according to any one of claims 1 to 4, characterized in that, The step of determining the first phase angle of each operation unit during each sampling in the first time period based on the first channel signal and the second channel signal during each sampling in the first time period includes: For any sampling of any operation unit in the first time period, the first phase angle of any operation unit at any sampling time in the first time period is obtained by using the arctangent function based on the first channel signal and the second channel signal at any sampling time in the first time period.

9. The method according to any one of claims 1 to 4, characterized in that, The step of eliminating phase error by processing the first phase angle of each operation unit during each sampling in the first time period to obtain the second phase angle of each operation unit during each sampling in the first time period includes: Obtain a first phase array including the first phase angle of each operation unit at any sampling time in the first time period; The first phase array is fitted with a straight line to obtain the fitted straight line; The difference between the first phase angle and the reference phase angle of any operation unit at any sampling time in the first time period is taken as the second phase angle of the operation unit at any sampling time in the first time period. The reference phase angle is the phase angle on the fitted straight line that corresponds to the first phase angle of the operation unit at any sampling time in the first time period.

10. The method according to any one of claims 1 to 4, characterized in that, The step of determining the vibration signal of each operating unit in the first time period based on the second phase angle of each sampling of each operating unit in the first time period includes: For any of the operation units, a difference sequence corresponding to the operation unit is obtained based on the second phase angle of the operation unit at each sampling time in the first time period. The difference sequence includes multiple values, and any value is the difference between the second phase angles at two adjacent sampling times. The differential sequence is processed to obtain the vibration signal of any of the operating units in the first time period.

11. The method according to claim 10, characterized in that, The process of processing the differential sequence to obtain the vibration signal of any operating unit in the first time period includes: For a second value in the difference sequence that is greater than a first value, a third value is determined based on the second value and the first value, wherein the third value is less than the first value; For the fifth value in the difference sequence that is less than the fourth value, a sixth value is determined based on the fifth value and the fourth value, wherein the sixth value is greater than the fourth value; The third value replaces the second value in the difference sequence, and the sixth value replaces the fourth value in the difference sequence, to obtain the vibration signal of any operating unit in the first time period.

12. A device for acquiring vibration signals, characterized in that, The device includes: The acquisition module is used to acquire the weights of each sampling point of the optical fiber laid in the same trench as the buried pipeline; The determining module is used to determine the first channel signal and the second channel signal of each operation unit of the optical fiber at each sampling time in the first time period according to the weight of each sampling point. Each operation unit includes at least one sampling point. The first channel signal of the operation unit at any sampling time in the first time period is used to indicate the vibration of the operation unit at any sampling time. The second channel signal of the operation unit at any sampling time in the first time period is used to indicate the vibration of the operation unit at any sampling time after the phase angle changes to the target angle. The determining module is further configured to determine the first phase angle of each operation unit during each sampling in the first time period based on the first channel signal and the second channel signal during each sampling in the first time period. The determining module is further configured to perform phase error elimination processing on the first phase angle of each operation unit during each sampling in the first time period to obtain the second phase angle of each operation unit during each sampling in the first time period. The determining module is further configured to determine the vibration signal of each operating unit in the first time period based on the second phase angle of each sampling of each operating unit in the first time period. The vibration signal of any operating unit in the first time period is used to indicate the vibration signal of the buried pipeline at the same location as the operating unit in the first time period. The acquisition module is configured to, for any sampling point, acquire the first channel signal and the second channel signal of that sampling point at each sampling time in a second time period, wherein the second time period is earlier than the first time period; determine a complex number of that sampling point at each sampling time in the second time period based on the first channel signal and the second channel signal of that sampling point at each sampling time in the second time period, wherein the real part of the complex number of that sampling point at any sampling time in the second time period is the first channel signal of that sampling point at that sampling time, and the imaginary part of the complex number of that sampling point at any sampling time in the second time period is the second channel signal of that sampling point at that sampling time; determine a target complex number based on the complex number of that sampling point at each sampling time in the second time period, wherein the real part of the target complex number is used to indicate the first channel signal of that sampling point in the second time period, and the imaginary part of the target complex number is used to indicate the second channel signal of that sampling point in the second time period; and determine the angle between the point indicated by the target complex number and the horizontal number axis. Based on the included angle, the weight of any sampling point is determined to be... The i is the unit of the imaginary part of the complex number.

13. A computer device, characterized in that, The computer device includes a processor and a memory, the memory storing at least one piece of program code, which is loaded and executed by the processor to enable the computer device to implement the vibration signal acquisition method as described in any one of claims 1 to 11.

14. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one piece of program code, which is loaded and executed by a processor to enable the computer to implement the vibration signal acquisition method as described in any one of claims 1 to 11.

15. A computer program product, characterized in that, The computer program product stores at least one computer instruction, which is loaded and executed by a processor to enable the computer to implement the vibration signal acquisition method as described in any one of claims 1 to 11.

Images