A calibration method and calibration device for a distributed optical fiber temperature measurement system

Abstract

The application discloses a kind of calibration method and calibration device of distributed optical fiber temperature measurement system.Method includes: the scattering response data along temperature measurement fiber under preset calibration condition is collected, and original calibration response sequence is constructed in combination with reference temperature zone temperature information;Candidate link disturbance position is extracted based on original calibration response sequence, target disturbance position caused by joint or fusion point is identified, and corresponding disturbance type is determined;According to disturbance type, local response distortion variable, inter-segment baseline offset, subsequent segment slope disturbance and calibration weight adjustment amount are determined, and the influence propagation relationship of each joint or fusion point on subsequent segment calibration parameter is constructed;With target disturbance position as section division reference, based on influence propagation relationship, temperature measurement fiber is segmented calibration function reconstruction, and combined with section continuity constraint and reference temperature zone constraint is jointly solved, and the calibrated distributed optical fiber temperature inversion model is obtained.The application improves temperature calibration accuracy under complex disturbance condition.

Description

Technical Field

[0001] This application relates to the field of fiber optic temperature measurement technology, specifically to a calibration method and calibration device for a distributed fiber optic temperature measurement system. Background Technology

[0002] Distributed Temperature Sensing (DTS) systems utilize the Raman or Brillouin scattering effects in optical fibers to achieve continuous temperature monitoring along the entire length of the fiber. They are widely used in scenarios such as power cable monitoring, tunnel fire early warning, and oil and gas pipeline safety monitoring. During system operation, to ensure measurement accuracy, it is typically necessary to calibrate the system by setting a reference temperature zone or standard temperature source, thereby establishing a mapping relationship between the scattered signal and the actual temperature.

[0003] However, in practical engineering applications, temperature-measuring optical fibers typically have numerous structural nodes along their length, such as fusion splices, connectors, and maintenance joints. These nodes introduce local insertion loss, enhanced reflection, and abrupt changes in scattering characteristics, causing abrupt changes or spikes in the scattering response curve along the fiber at corresponding locations. Existing calibration methods often use threshold detection or simple smoothing to correct these anomalies, but these methods only compensate for local data and fail to effectively consider the propagation effects of joint disturbances on the signal attenuation trend, signal-to-noise ratio distribution, and calibration parameter fitting process in subsequent sections.

[0004] Furthermore, when real temperature changes occur near the joint, the joint loss characteristics and temperature change characteristics overlap, which can easily lead to misjudgment, resulting in an overall deviation or even failure of the calibration model.

[0005] Therefore, how to accurately identify and describe the structural impact of connector disturbances on the calibration curve during the calibration process, and on this basis, achieve effective reconstruction of calibration parameters, has become an urgent technical problem to be solved. Summary of the Invention

[0006] To address the technical problems mentioned in the background section, the purpose of this application is to provide a calibration method, apparatus, electronic device, storage medium, and computer program product for a distributed fiber optic temperature measurement system.

[0007] According to a first aspect of this application, a calibration method for a distributed optical fiber temperature measurement system is provided, comprising the following steps: S1: Collect the scattering response data along the temperature measuring fiber of the distributed fiber temperature measurement system under the preset calibration conditions, and construct the original calibration response sequence by combining it with the temperature information of the reference temperature zone. S2, based on the original calibration response sequence, extract candidate link disturbance locations, identify target disturbance locations caused by joints or fusion splices among the candidate link disturbance locations, and determine the corresponding joint or fusion splice disturbance type according to the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics and channel response asymmetry characteristics of each target disturbance location; S3. Based on the disturbance type of each joint or weld point, determine the corresponding local response distortion, inter-segment baseline offset, subsequent segment slope disturbance, and calibration weight adjustment, and construct the influence propagation relationship of each joint or weld point on the calibration parameters of the subsequent segment. S4. Using the target disturbance location as a segmentation reference, the temperature measuring fiber is reconstructed into a segmented calibration function based on the influence propagation relationship. The calibration parameters of each segment are jointly solved by combining the segment continuity constraint and the reference temperature zone constraint to obtain the calibrated distributed fiber temperature inversion model.

[0008] According to a second aspect of this application, a calibration device for a distributed optical fiber temperature measurement system is provided, the device comprising: The original sequence construction module is used to collect the scattering response data along the temperature measuring fiber of the distributed fiber temperature measurement system under the preset calibration conditions, and to construct the original calibration response sequence by combining the temperature information of the reference temperature zone. The disturbance identification module is used to extract candidate link disturbance locations based on the original calibration response sequence, identify target disturbance locations caused by joints or fusion splices among the candidate link disturbance locations, and determine the corresponding joint or fusion splice disturbance type based on the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics and channel response asymmetry characteristics of each target disturbance location. The propagation relationship construction module is used to determine the corresponding local response distortion, inter-segment baseline offset, subsequent segment slope disturbance and calibration weight adjustment based on the disturbance type of each joint or weld point, and to construct the propagation relationship of the influence of each joint or weld point on the calibration parameters of the subsequent segment. The model reconstruction module is used to reconstruct the segmented calibration function of the temperature measuring fiber based on the target disturbance location as the segment division reference and the influence propagation relationship. It also combines the segment continuity constraint and the reference temperature zone constraint to jointly solve the calibration parameters of each segment, thereby obtaining the calibrated distributed fiber temperature inversion model.

[0009] According to a third aspect of this application, an electronic device is provided, including a memory and a processor, wherein the memory stores a computer program that, when executed by the processor, implements the method as described in any of the preceding claims.

[0010] According to a fourth aspect of this application, a storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the method as described in any of the preceding claims.

[0011] According to a fifth aspect of this application, a computer program product is provided, comprising a computer program that, when executed by a processor, implements the method as described in any of the preceding claims.

[0012] Compared with existing technologies, this invention identifies link disturbances caused by joints or fusion splices in a typological manner, and extracts local response distortions, inter-segment baseline offsets, subsequent segment slope disturbances, and calibration weight adjustments for different disturbance types. It further constructs the propagation relationship of the disturbance's influence on subsequent segment calibration parameters. This propagation relationship is introduced into the segmented calibration function reconstruction process and combined with segment continuity constraints and reference temperature zone constraints for joint solution. This effectively characterizes the cumulative effect of multiple disturbance sources superimposed and propagating along the optical fiber, reducing misjudgments of local anomalies by traditional overall calibration or single-point correction methods. This significantly improves the accuracy and stability of distributed optical fiber temperature measurement results under complex link conditions, while also enhancing the model's adaptability to different installation conditions and long-term operational drift. Attached Figure Description

[0013] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0014] Other features, objects, and advantages of this application will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings: Figure 1 A schematic flowchart illustrating a calibration method for a distributed optical fiber temperature measurement system provided in this application embodiment; Figure 2 This is a schematic diagram of the structure of a calibration device for a distributed optical fiber temperature measurement system provided in an embodiment of this application; Figure 3 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Detailed Implementation

[0015] To make the objectives, technical solutions, and advantages of this application clearer, the technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of this application without creative effort are within the scope of protection of this application.

[0016] like Figure 1 As shown in the figure, this application discloses a calibration method for a distributed optical fiber temperature measurement system, which is applied to the distributed optical fiber temperature measurement system. The system can obtain temperature distribution information along the temperature measurement optical fiber based on the principles of Raman scattering or Brillouin scattering.

[0017] The method includes the following steps: S1: Collect the scattering response data along the temperature measuring fiber of the distributed fiber temperature measurement system under the preset calibration conditions, and construct the original calibration response sequence by combining it with the temperature information of the reference temperature zone. Specifically, data is acquired from the distributed fiber optic temperature measurement system under preset calibration conditions. These preset calibration conditions involve a known or controllable temperature field; for example, several reference temperature zones are set along the temperature-measuring fiber, and each zone is kept at a stable and known temperature value. By emitting a probe light signal and receiving the scattered signals generated along the fiber (such as Raman anti-Stokes light and Stokes light signals or Brillouin frequency shift signals), scattering response data corresponding to each sampling position along the fiber's length is obtained.

[0018] After obtaining the scattering response data, it is preprocessed, including mapping the time axis to the spatial location, normalizing the signal intensity, and performing necessary noise reduction, thereby forming a response sequence that corresponds one-to-one with the spatial location of the optical fiber. It should be noted that this response sequence can be represented as a set of discrete sampling points arranged along the length of the optical fiber, with each sampling point corresponding to a response value.

[0019] The above response sequence is aligned with the temperature information of the reference temperature region. Specifically, this can be done by identifying the position interval of the reference temperature region in the response sequence and establishing a mapping relationship between the response values ​​within the corresponding interval and the known temperature values, thereby introducing temperature calibration anchor points into the response sequence. In some embodiments, interpolation or fitting methods can be used to extend the temperature information of the discrete reference temperature region to the entire response sequence to form a continuous temperature reference constraint.

[0020] Through the above processing, the original calibration response sequence is constructed, which not only contains the scattering response information distributed along the optical fiber space, but also implicitly contains the temperature reference relationship with the reference temperature region.

[0021] S2, based on the original calibration response sequence, extract candidate link disturbance locations, identify target disturbance locations caused by joints or fusion splices among the candidate link disturbance locations, and determine the corresponding joint or fusion splice disturbance type according to the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics and channel response asymmetry characteristics of each target disturbance location; In the constructed original calibration response sequence, anomalous locations that structurally affect the temperature inversion model are identified. Further analysis is conducted to determine whether these anomalous locations are caused by joints or fusion splices, and to identify the specific type of disturbance. It is understandable that in the actual deployment of distributed fiber optic temperature measurement systems, multiple joints, fusion splices, or maintenance connection points often exist along the temperature-measuring fiber. These locations not only cause local optical power loss but may also be accompanied by enhanced reflection, changes in subsequent attenuation trends, and inconsistencies in multi-channel responses. If these locations are not identified separately from the overall calibration process, and a unified fitting method is still used for parameter solving, anomalies caused by changes in link structure are easily misjudged as temperature changes, leading to calibration function shifts and reducing the accuracy of the temperature inversion results.

[0022] Therefore, this step extracts various local feature parameters that can reflect the characteristics of link disturbances at each sampling location. Based on these feature parameters, candidate link disturbance locations suspected of having structural disturbances are identified. Furthermore, locations that consistently appear in multiple calibration acquisitions and exhibit typical response characteristics of joints or fusion points are selected as target disturbance locations. Subsequently, the disturbance type is determined based on the local step characteristics, reflection spike characteristics, attenuation gradient abrupt changes, and channel response asymmetry characteristics of the target disturbance location.

[0023] The above method can distinguish between structural disturbances in the link and actual temperature changes during the calibration phase, enabling subsequent processing operations to determine the impact amount based on the disturbance type and construct the impact propagation relationship.

[0024] In some embodiments, candidate link disturbance locations are extracted based on the original calibration response sequence, target disturbance locations caused by joints or fusion splices are identified among the candidate link disturbance locations, and the corresponding joint or fusion splice disturbance type is determined based on the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics, and channel response asymmetry characteristics of each target disturbance location, including: S21. Based on the original calibration response sequence, calculate the mean difference of the response before and after each sampling position, the local peak amplitude, the change in attenuation slope, and the difference in the dual-channel response within a preset window to obtain the disturbance characterization parameters corresponding to each sampling position. Specifically, each sampling position in the original calibration response sequence is taken as the analysis center, and a local window of a preset length is selected before and after the sampling position. The response change characteristics within the window are statistically calculated to form a disturbance characterization parameter that can characterize whether there is a link disturbance at that position.

[0025] The difference in the mean response is used to characterize whether there is a significant step change in the signal baseline before and after the sampling position. For abrupt changes in insertion loss caused by a joint or fusion splice, this is often manifested as a significant difference in the average response level before and after that position. Therefore, by calculating the difference between the mean response within the window before and after the sampling position, it is possible to preliminarily determine whether a step-type disturbance exists.

[0026] The local peak amplitude is used to characterize whether there is a significant reflection spike near the sampling location. Typically, when there is a mismatch at the joint end face, unstable welding quality, or local reflection enhancement, a local narrow peak will appear in the response sequence. Therefore, the intensity of the reflection spike can be characterized by extracting the difference between the local maximum peak near this location and its neighborhood background level.

[0027] The change in attenuation slope is used to reflect whether there is an abrupt change in the attenuation trend of the signal along the fiber length before and after the sampling position. If a joint or fusion splice not only causes local power loss but also changes the effective transmission state of its subsequent segments, then there will usually be a detectable difference in the response attenuation slope before and after that position.

[0028] The dual-channel response difference is used to characterize the change in response consistency between the two detection channels. For example, in a Raman scattering-based temperature measurement system, the response values ​​of the Stokes channel and the anti-Stokes channel can be collected separately, and the relative changes between the two channels before and after the position can be compared. If the difference increases significantly, it indicates that the position may have triggered mode redistribution, polarization change, or channel imbalance.

[0029] The aforementioned disturbance characterization parameters are represented as feature vectors containing multiple components, each component of which corresponds to the aforementioned mean difference in response, local peak amplitude, change in attenuation slope, and difference in dual-channel response. It should be understood that using multiple feature quantities in combination can avoid false detection problems caused by relying solely on a single anomaly indicator. For example, the mere presence of a mean difference in response does not necessarily indicate a joint disturbance; it could also be a response change caused by a local temperature gradient. Similarly, the mere presence of a local peak does not necessarily indicate an abnormal reflection; it could also be caused by transient noise.

[0030] S22, Identify candidate link disturbance locations that meet preset conditions based on the disturbance characterization parameters, and select stable abnormal locations that repeatedly appear in multiple calibration acquisition results from the candidate link disturbance locations as the target disturbance locations; Specifically, based on the disturbance characterization parameters obtained above, preliminary anomaly screening is performed on each sampling position in the original calibration response sequence to identify candidate link disturbance positions that may have link disturbances.

[0031] Specifically, corresponding judgment thresholds or joint judgment conditions can be set for different disturbance characterization parameters. For example, when the mean difference of the response at a certain sampling location exceeds a first preset threshold, it can be determined that the location has a significant step anomaly tendency; when the local peak amplitude exceeds a second preset threshold, it can be determined that the location has a reflection enhancement tendency; when the change in attenuation slope exceeds a third preset threshold, it can be determined that the location may cause changes in the transmission state of subsequent segments; when the difference in the dual-channel response exceeds a fourth preset threshold, it can be determined that the location may be accompanied by changes in mode distribution or channel coupling state. Sampling locations that meet one or more of these conditions can be used as candidate link disturbance locations.

[0032] It should be noted that during actual measurement, the original calibration response sequence may contain spurious disturbance points caused by transient noise, environmental perturbations, detector fluctuations, or local short-term anomalies. These spurious disturbance points typically only appear in individual acquisition cycles and lack spatial repeatability. Therefore, this embodiment further utilizes multiple calibration acquisition results for stability screening. Specifically, candidate link disturbance positions in multiple acquisition cycles can be spatially aligned, and the recurrence frequency or probability of each candidate position in different acquisition cycles can be counted. For abnormal positions that appear continuously in multiple calibration acquisition results and whose spatial position deviation remains within a preset tolerance range, they can be identified as stable abnormal positions and used as target disturbance positions. Correspondingly, abnormal points that only appear sporadically or have large positional drift can be considered as spurious anomalies caused by random noise or transient fluctuations and removed from the candidate set.

[0033] The target disturbance location obtained through the above processing is a structural anomaly location with strong spatial stability and repeatability. It is more consistent with the fixed distribution characteristics of joints or fusion points in the actual link, thus helping to improve the accuracy of subsequent disturbance type identification.

[0034] S23. For each of the target disturbance locations, based on the corresponding local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics, and channel response asymmetry characteristics, determine the corresponding joint or fusion point disturbance type, including at least one of pure attenuation step disturbance, reflection spike coupled attenuation disturbance, and mode redistribution disturbance.

[0035] Specifically, the previously identified target disturbance locations need further typological determination to clarify their specific impact mechanism on the calibration process. It is understandable that the response characteristics caused by different abnormal states of joints or weld points are not entirely the same; therefore, classification based on the combined relationships of multiple characteristic quantities is necessary.

[0036] Specifically, when a target disturbance location is characterized by a significant difference in the mean values ​​of the responses before and after the disturbance, while the local peak amplitude is not significant and the difference in the dual-channel response does not increase significantly, the disturbance location can be identified as a pure attenuation step-type disturbance. This type of disturbance typically corresponds to joints or welds where insertion loss increases but reflection is weak. Its main effects are abrupt changes in the local response baseline and discontinuities in energy levels between segments.

[0037] When a target disturbance location not only exhibits a significant step change in response but also displays a marked local spike, and the attenuation slope of its subsequent segments changes, the disturbance location can be identified as a reflection spike coupled attenuation type disturbance. This type of disturbance typically corresponds to end-face mismatch, enhanced local reflection, or unstable junction conditions. Its impact is not only reflected in the local anomaly at that location but also further affects the signal attenuation trend of subsequent segments.

[0038] When a target perturbation location exhibits a significant increase in the difference between the two-channel responses, while the fluctuation level, signal-to-noise ratio, or inter-channel consistency of its subsequent segments change, and local steps or reflection spikes are not necessarily significant, the perturbation location can be identified as a mode redistribution type perturbation. This type of perturbation is usually related to energy redistribution caused by changes in mode coupling, polarization state changes, or local microbending, and its main impact is reflected in changes in the response stability and calibration reliability of subsequent segments.

[0039] The type of disturbance can be determined using either a rule-based approach or a classification model built based on training samples. If a rule-based approach is used, the determination conditions can be set according to the logical combination relationships between various features. If a classification model is used, local step features, reflection spike features, attenuation gradient abrupt changes, and channel response asymmetry features can be used as model inputs to output corresponding disturbance type labels. This invention does not impose specific limitations on these methods.

[0040] S3. Based on the disturbance type of each joint or weld point, determine the corresponding local response distortion, inter-segment baseline offset, subsequent segment slope disturbance, and calibration weight adjustment, and construct the influence propagation relationship of each joint or weld point on the calibration parameters of the subsequent segment. This step is used to determine how different types of disturbances affect the temperature measurement signal, the extent of their impact, and how the impact is transmitted to subsequent segments. This requires extracting multidimensional influence quantities that can characterize local anomalies, segment offsets, attenuation changes, and reliability changes based on the physical mechanisms of different types of disturbances.

[0041] Specifically, different types of joint or weld point disturbances have significantly different effects on the temperature measurement signal. For pure attenuation step disturbances, the main manifestation is a sudden change in the response level before and after the disturbance location. This type of disturbance usually does not form a significant reflection peak, but it will cause a jump in the local energy baseline. Therefore, it is necessary to focus on extracting the local response distortion variable that can reflect the degree of response anomaly near the location, as well as the inter-segment baseline offset that can characterize the change in the overall response reference level of the segment before and after the disturbance.

[0042] For reflection spike coupling attenuation type disturbances, not only are spike-shaped local anomalies introduced at the disturbance location, but the response attenuation trend of subsequent segments is also changed. Therefore, in addition to local response distortion and inter-segment baseline offset, it is also necessary to extract the slope disturbance of subsequent segments to characterize the continuous influence of the disturbance on the variation law of the response of subsequent segments along the fiber direction.

[0043] For pattern redistribution type perturbations, their effect on local amplitude abrupt changes may not be significant, but they often lead to a decrease in channel consistency, increased volatility, or worse signal stability in subsequent segments. Therefore, this type of perturbation is more suitable for characterization by adjusting the calibration weights to reflect the change in the reliability of the corresponding segment in the subsequent parameter fitting and joint solution process.

[0044] Furthermore, the local response distortion variable is used to describe the instantaneous deviation of the response near the disturbance location from the normal state; the inter-segment baseline offset is used to describe the difference in the overall response level between adjacent segments before and after the disturbance; the subsequent segment slope disturbance is used to describe the change in the decay rate or trend of the response of the subsequent segment due to the disturbance; and the calibration weight adjustment is used to describe how the participation intensity of certain segment data in the subsequent calibration model should be adjusted under the influence of the disturbance.

[0045] In some embodiments, constructing the propagation relationship of the influence of each joint or weld point on the calibration parameters of subsequent sections includes: S31, taking each of the target disturbance positions as the propagation starting point, determine the corresponding subsequent influence section according to the extension direction of the temperature measuring fiber; determine the influence intensity of each target disturbance position on each subsequent influence section based on the disturbance type of the joint or fusion splice and the corresponding local response distortion, inter-segment baseline offset, subsequent section slope disturbance amount and calibration weight adjustment amount. Specifically, each target disturbance location is taken as the propagation starting point, and its corresponding subsequent influence segment is determined according to the extension direction of the temperature-measuring optical fiber (i.e., the direction of optical signal propagation). It is understood that in a distributed optical fiber temperature measurement system, the optical signal propagates along the fiber from the transmitting end and gradually attenuates. Changes in the upstream link structure have a cumulative impact on the downstream signal, while downstream disturbances typically do not have a reverse effect on the upstream region. Therefore, this step only considers the unidirectional influence relationship between the target disturbance location and its subsequent segments.

[0046] The subsequent influence section can be divided into several continuous sub-segments according to a fixed spatial length, or it can be adaptively divided based on the attenuation trend change or signal stability change in the original calibration response sequence, so that the segment division is more in line with the actual signal change characteristics.

[0047] After determining the subsequent affected sections, the impact intensity is further determined based on the disturbance type and corresponding impact amount at the joint or weld point. Specifically, the method for constructing the impact intensity differs for different disturbance types: For pure attenuation step disturbances, the main manifestation is that the signal energy changes abruptly at the disturbance location. Therefore, the influence intensity of this type of disturbance is mainly determined by the inter-segment baseline offset, and is corrected by the local response distortion, which is used to describe the degree of offset of the overall signal baseline of the subsequent segment. For reflection spike coupling attenuation type disturbances, since they simultaneously introduce local reflection enhancement and subsequent transmission loss changes, their impact intensity includes not only the inter-segment baseline offset, but also the subsequent segment slope disturbance amount, in order to characterize the continuous impact of the disturbance on the signal attenuation trend of the subsequent segment. At the same time, the local response distortion is used to correct the initial amplitude of the impact at the disturbance initiation point. For pattern redistribution type perturbations, their main impact is on signal channel consistency and stability. Therefore, the intensity of their impact is more reflected by the adjustment of calibration weights, which is used to characterize the change in the degree of participation of subsequent segments in the parameter solution process, that is, to reduce the weight contribution of low confidence segments.

[0048] To facilitate subsequent propagation and allocation, the influence intensity can be expressed as a function related to the segment location. For example, an initial influence weight function can be constructed with the target disturbance location as a reference point. This function takes its maximum value at the propagation starting point and serves as the basic input for decreasing allocation in subsequent steps.

[0049] S32, combining the relative distance between the target disturbance location and each subsequent affected section, the signal quality of the section, and the superimposed influence of other subsequent target disturbance locations, the influence intensity is distributed in a decreasing manner to obtain the influence propagation relationship of each joint or fusion point on the calibration parameters of the subsequent section.

[0050] It should be understood that the impact of disturbances on subsequent segments typically exhibits the following characteristics: first, it gradually weakens with increasing propagation distance; second, it is affected by the signal quality of the segment, and its effective impact should be appropriately reduced in areas with poor signal quality; third, when multiple disturbance locations exist, the impacts of different disturbance sources on the same segment will be superimposed or even compete. Therefore, this invention constructs a method for allocating the impact of distance attenuation, signal quality modulation, and the superposition of multiple disturbances, as detailed below: In some embodiments, the influence intensity is progressively distributed by combining the relative distance between the target disturbance location and each subsequent affected segment, the segment signal quality, and the superimposed influence of other subsequent target disturbance locations, to obtain the propagation relationship of the influence of each joint or weld point on the calibration parameters of the subsequent segment, including: S321, Based on the relative distance between the target disturbance location and each subsequent affected segment and the signal quality parameters of each subsequent affected segment, determine the basic influence intensity of each target disturbance location on each subsequent affected segment; Specifically, the influence attenuation relationship is first established based on the relative distance between the target disturbance location and each subsequent affected segment. Exponential decay, piecewise linear decay, or other monotonically decreasing function forms can be used, so that the influence intensity is greater when close to the disturbance location and gradually decreases when far away from the disturbance location.

[0051] Simultaneously, the intensity of the fundamental influence is modulated using segment signal quality parameters. These signal quality parameters may include signal-to-noise ratio, response fluctuation, and channel consistency indicators. When the signal quality of a segment is high, it indicates that the response data of that segment is more stable and reliable, and a larger weight can be retained for the impact of disturbances. Conversely, when the signal quality is low, to avoid amplifying the noise effect on the calibration process, the intensity of the fundamental influence can be suppressed, thereby reducing the effective impact of the disturbance on that segment.

[0052] The specific process for determining the basic influence intensity includes: first, using the initial influence intensity corresponding to the previously determined target disturbance position as the benchmark intensity, dividing each subsequent influence segment into a near-end segment, an intermediate transition segment, and a far-end segment according to the relative distance between them and the target disturbance position, and assigning a preset proportional coefficient to each segment to obtain the influence intensity after distance allocation; then, obtaining the signal quality parameters of each subsequent influence segment, and performing a weighted correction on the influence intensity after distance allocation based on the signal quality parameters, wherein the correction coefficient corresponding to the segment with higher signal quality is greater than or equal to the preset benchmark value, and the correction coefficient corresponding to the segment with lower signal quality is less than the preset benchmark value, thereby obtaining the basic influence intensity of each target disturbance position on each subsequent influence segment.

[0053] It should be understood that, in order to avoid abrupt changes in the basic influence intensity between adjacent segments, the modified influence intensity can also be smoothed between segments to make it change continuously along the direction of the temperature measuring fiber.

[0054] S322, combining the spatial distribution relationship of other target disturbance positions and the corresponding basic influence intensity, the basic influence intensity is reduced and superimposed to obtain the influence propagation relationship of each target disturbance position on the calibration parameters of each subsequent influence section.

[0055] Specifically, within the same subsequent influence segment, there may be base influence intensities from multiple target disturbance locations simultaneously. Direct linear superposition could lead to an abnormal amplification of the influence intensity, resulting in overcompensation during subsequent calibration. Therefore, this embodiment employs a decreasing allocation and superposition correction mechanism, which includes the following processing methods: (1) Normalization allocation: When multiple disturbance locations affect the same segment, the basic influence intensity of each disturbance is normalized so that the sum is kept within the preset range, thereby avoiding the overall weight loss caused by the superposition of multiple disturbances.

[0056] (2) Prioritize dominant disturbances: Determine the location of dominant disturbances based on the magnitude of the basic influence intensity of each disturbance location. For disturbance sources with a significantly higher basic influence intensity than other disturbances, assign them a higher weight, while attenuating other weak disturbances, thereby highlighting the impact of the main disturbance on the section.

[0057] (3) Spatial competition and shielding: When multiple disturbances are spatially close and their affected areas highly overlap, a competition function can be introduced to redistribute their effects. For example, when an upstream disturbance has already had a significant impact on a certain area, the additional impact of a downstream disturbance on that area can be partially suppressed, thereby avoiding the duplication of similar effects.

[0058] In addition, a coupling attenuation factor based on the perturbation spacing can be introduced. That is, when the distance between two perturbation locations is small, their influence coupling degree is high, and their superimposed influence is compressed as a whole; when the distance is large, their influence is allowed to superimpose relatively independently.

[0059] Specifically, the superposition correction process of the basic influence intensity includes: for any subsequent influence segment, obtaining the basic influence intensity corresponding to all target disturbance positions that affect the segment; determining the priority of each target disturbance position based on the relative distance between each target disturbance position and the segment and their order in the fiber extension direction; taking the basic influence intensity corresponding to the target disturbance position with the highest priority as the dominant influence intensity, and applying an attenuation coefficient to the basic influence intensities corresponding to the remaining target disturbance positions before superimposing them to obtain the comprehensive influence intensity of the segment.

[0060] When the distance between multiple target disturbance locations is less than a preset spacing threshold, they are considered as a coupled disturbance group. The basic influence intensity of the coupled disturbance group is normalized as a whole before being allocated to the corresponding subsequent influence segment, so as to avoid repeated amplification of disturbance influence in local areas.

[0061] In addition, component-level superposition corrections can be performed separately for different types of disturbances. That is, the above superposition correction process is performed on the local response distortion, inter-segment baseline offset, subsequent segment slope disturbance, and calibration weight adjustment, so that the effects of different types of disturbances on each calibration parameter dimension are independent yet synergistically constrained.

[0062] Through the aforementioned decreasing allocation and superposition correction process, the propagation relationship of the influence of each target disturbance location on the calibration parameters of each subsequent affected segment is finally obtained. This propagation relationship not only reflects the spatial competition and coupling characteristics among multiple disturbance sources, but also avoids the overcompensation problem caused by simple superposition, thus providing a more stable and physically consistent constraint basis for the construction of subsequent piecewise calibration functions.

[0063] S4. Using the target disturbance location as a segmentation reference, the temperature measuring fiber is reconstructed into a segmented calibration function based on the influence propagation relationship. The calibration parameters of each segment are jointly solved by combining the segment continuity constraint and the reference temperature zone constraint to obtain the calibrated distributed fiber temperature inversion model.

[0064] Because disturbances caused by joints or welds exhibit significant spatial segmentation, using a single global calibration function would be insufficient to simultaneously account for the response differences across different segments. Therefore, this embodiment divides the system into segments based on the target disturbance location and constructs a separate calibration function for each segment. A constraint mechanism is also used to ensure reasonable connection between the segments.

[0065] In some embodiments, the target disturbance location is used as a segmentation reference, and the temperature-measuring fiber is reconstructed using a segmented calibration function based on the influence propagation relationship. The calibration parameters for each segment are then jointly solved by combining segment continuity constraints and reference temperature zone constraints to obtain a calibrated distributed fiber temperature inversion model, including: S41, the temperature measuring fiber is divided into multiple sections to be calibrated using the target disturbance positions as boundaries, and a piecewise calibration function containing a baseline correction term, a slope correction term, and a weight correction term is constructed for each section to be calibrated; the influence propagation relationship is introduced into the piecewise calibration function of the corresponding section to be calibrated to make preliminary corrections to the calibration parameters of each section to be calibrated. Specifically, the target disturbance positions are first sorted according to their spatial location on the temperature-sensing fiber, and the temperature-sensing fiber is divided into multiple continuous segments to be calibrated, using the interval between adjacent target disturbance positions as the dividing boundary. The interval between the fiber start end and the first target disturbance position, and the interval from the last target disturbance position to the fiber end, are treated as the beginning and end segments respectively, thus forming a complete segment division structure.

[0066] After segmenting the data, a piecewise calibration function is constructed for each segment to be calibrated. This function describes the correction process for the scattering response within that segment and the mapping relationship between the corrected response and temperature. For example, the first segment... The section calibration response function of each section to be calibrated can be expressed as: in, For the first Section in location The calibrated response value at the location, This is the original scattering response. , and These are the weight correction parameter, baseline correction parameter, and slope correction parameter, respectively. This is the starting position of this section. This is the residual error term.

[0067] It should be understood that the aforementioned piecewise calibration response function does not directly correct the temperature value, but rather first performs a structured correction on the original scattering response within the response domain. Specifically, the item... This represents the adjustment result of the overall amplitude of the original scattering response, used to reflect the influence of the disturbance on the strength relationship of the response in this section; Item This represents the compensation result for the overall response baseline of the section, used to eliminate the overall offset of the section caused by joints or weld points; the term represents the correction result for the response within the section as a function of spatial location, where is used to characterize the rate of response change caused by a unit change in spatial location, thereby reducing the positional offset. This is converted into a corresponding response correction value to characterize the impact of the disturbance on the attenuation characteristics of subsequent segments.

[0068] Furthermore, based on the calibrated response value The calibrated temperature of this section is obtained through the corresponding temperature inversion mapping relationship, i.e.:

[0069] Among them, is the first Section in location The calibrated temperature at point is the temperature inversion mapping function corresponding to the segment . It can be understood that the temperature inversion mapping function can be a pre-calibrated response-temperature mapping relationship, or a segment mapping relationship dynamically updated during subsequent joint solution.

[0070] Furthermore, the initial values ​​of the above parameters are constructed from the various disturbance influence quantities determined in step S3. Among them, the baseline correction parameter... The slope correction parameter is determined by the inter-segment baseline offset. The weighting correction parameters are determined by the slope perturbation of subsequent segments. The calibration weight adjustment amount determines the piecewise calibration function so that it can reflect the initial impact of local disturbances on the signal.

[0071] Based on this, the aforementioned influence propagation relationship is introduced into the piecewise calibration function of each segment. Specifically, for any segment to be calibrated, all target perturbation locations that affect that segment are identified, and the corresponding influence intensities are obtained. These intensities are then used as modulation factors applied to each correction parameter. For example, the influence intensities can be applied to the amplitude of the baseline correction parameter, the adjustment range of the slope correction parameter, and the scaling ratio of the weight correction parameter, thereby achieving the mapping of perturbation influence across different parameter dimensions.

[0072] In some embodiments, the parameter update relationship can be expressed as:

[0073]

[0074]

[0075] in, These are the initial parameters. They represent the first The target disturbance location is related to the first The weighting effect of the segment, the baseline offset, and the slope disturbance. This is the corresponding influence propagation coefficient. Further, the influence propagation coefficient is determined by the relative distance between the disturbance location and the segment, the signal quality of the segment, and the superposition relationship of multiple disturbances. For example, the propagation coefficient can be expressed as:

[0076] in, The distance between the disturbance location and the segment. The attenuation coefficient is... These are the signal quality parameters for the section. This is a correction factor for multiple perturbations. Similarly, and It can also be determined by the propagation relationship of the same or similar structure.

[0077] When the same segment is affected by multiple disturbances, the intensity of each disturbance is normalized or prioritized to ensure it meets a preset range, thereby avoiding excessive amplification caused by the superposition of multiple disturbances. For disturbances that are spatially close to each other, coupling processing can also be performed so that their effects participate in parameter correction as a whole.

[0078] S42, combining the segment continuity constraints at the boundary of adjacent calibrated segments and the temperature constraints of the reference temperature zone, the calibration parameters of each calibrated segment are jointly solved to obtain the calibrated distributed fiber temperature inversion model.

[0079] After constructing the segmented calibration functions for each segment and making preliminary parameter corrections, constraints are introduced to jointly solve the parameters of all segments.

[0080] First, a segment continuity constraint is introduced. At the boundary of adjacent segments, the output values ​​of their calibration functions are required to remain consistent, and their changing trends are restricted to avoid abrupt changes in the temperature inversion results at the segment boundaries. In some embodiments, adjacent segments can be constrained to satisfy both the continuity of the calibration response function and the continuity of the temperature inversion results. Specifically, at the segment boundary... The following conditions can be met:

[0081] Furthermore, the corresponding temperature inversion results are constrained to remain consistent at the boundaries or to vary by no more than a preset threshold. This ensures the continuity of the response correction process and avoids unreasonable jumps in the final temperature result.

[0082] Secondly, a reference temperature range constraint is introduced. For sampling points located within the reference temperature range, their calibrated temperature output should be consistent with the known temperature or meet a preset error range. That is, for locations within the reference temperature range, the following conditions must be met: in, The reference temperature is known. This constraint allows for the overall model to be corrected and anchored, preventing overall shifts or accumulated errors.

[0083] Based on this, the calibration parameters of each section are treated as overall variables to construct a joint solution model that satisfies both data fitting requirements and section continuity and reference temperature range constraints. In some embodiments, the solution can be achieved through the following objectives:

[0084] in, For the first Section, The weighting function is based on signal quality. This refers to the observed temperature value at the corresponding location or the temperature observation value obtained from the reference temperature zone and calibration relationship. It should be noted that here... and All terms are defined within the temperature domain to ensure the consistency of the objective function in terms of dimensions.

[0085] Furthermore, the weights of each segment in the solution process can be adjusted according to the signal quality parameters of each segment, so that the segments with higher signal quality have a stronger constraint effect on the final model, thereby improving the stability and anti-interference ability of the overall inversion results.

[0086] Through the above-mentioned piecewise function construction and multi-constraint joint solution process, a distributed optical fiber temperature inversion model that is continuous and consistent throughout the entire optical fiber range and can effectively compensate for the effects of multiple disturbance propagation is finally obtained.

[0087] Reference Figure 2 As shown in the figure, this application embodiment also provides a calibration device 200 for a distributed optical fiber temperature measurement system, the device comprising: The original sequence construction module 10 is used to collect the scattering response data along the temperature measuring fiber of the distributed fiber temperature measurement system under the preset calibration conditions, and to construct the original calibration response sequence by combining the temperature information of the reference temperature zone. The disturbance identification module 20 is used to extract candidate link disturbance locations based on the original calibration response sequence, identify target disturbance locations caused by joints or fusion splices among the candidate link disturbance locations, and determine the corresponding joint or fusion splice disturbance type based on the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics and channel response asymmetry characteristics of each target disturbance location. The propagation relationship construction module 30 is used to determine the corresponding local response distortion, inter-segment baseline offset, subsequent segment slope disturbance and calibration weight adjustment based on the disturbance type of each joint or weld point, and to construct the propagation relationship of the influence of each joint or weld point on the calibration parameters of the subsequent segment. The model reconstruction module 40 is used to reconstruct the segmented calibration function of the temperature measuring fiber based on the target disturbance location as the segment division reference and the influence propagation relationship, and to jointly solve the calibration parameters of each segment by combining the segment continuity constraint and the reference temperature zone constraint, so as to obtain the calibrated distributed fiber temperature inversion model.

[0088] In some embodiments, the disturbance identification module 20 is configured to: Based on the original calibration response sequence, the mean difference of the response, the local peak amplitude, the change in attenuation slope, and the difference in the dual-channel response within the preset window before and after each sampling position are calculated to obtain the disturbance characterization parameters corresponding to each sampling position. Candidate link disturbance locations that meet preset conditions are identified based on the disturbance characterization parameters, and stable abnormal locations that repeatedly appear in multiple calibration acquisition results are selected from the candidate link disturbance locations as the target disturbance locations. For each target disturbance location, based on the corresponding local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics, and channel response asymmetry characteristics, the corresponding joint or fusion point disturbance type is determined, including at least one of pure attenuation step disturbance, reflection spike coupled attenuation disturbance, and mode redistribution disturbance.

[0089] In some embodiments, the propagation relationship construction module 30 is configured to: Using the location of each target disturbance as the propagation starting point, the corresponding subsequent affected sections are determined according to the extension direction of the temperature measuring fiber; based on the type of disturbance at the joint or fusion splice and the corresponding local response distortion, inter-segment baseline offset, subsequent section slope disturbance, and calibration weight adjustment, the intensity of the influence of each target disturbance location on each subsequent affected section is determined. By combining the relative distance between the target disturbance location and each subsequent affected segment, the segment signal quality, and the superimposed influence of other subsequent target disturbance locations, the influence intensity is distributed in a decreasing manner to obtain the influence propagation relationship of each joint or fusion point on the calibration parameters of the subsequent segment.

[0090] In some embodiments, the propagation relationship construction module 30 is configured specifically for: Based on the relative distance between the target disturbance location and each subsequent affected segment and the signal quality parameters of each subsequent affected segment, the basic influence intensity of each target disturbance location on each subsequent affected segment is determined. By combining the spatial distribution of other target disturbance locations and their corresponding basic influence intensity, the basic influence intensity is progressively distributed and superimposed to obtain the influence propagation relationship of each target disturbance location on the calibration parameters of each subsequent influence segment.

[0091] In some embodiments, the model reconstruction module 40 is configured to: The temperature measuring fiber is divided into multiple calibration segments using the target disturbance locations as boundaries, and a piecewise calibration function containing a baseline correction term, a slope correction term, and a weight correction term is constructed for each calibration segment. The influence propagation relationship is introduced into the piecewise calibration function of the corresponding calibration segment to perform preliminary correction of the calibration parameters of each calibration segment. By combining the continuity constraints of adjacent sections at the boundary of the sections to be calibrated and the temperature constraints of the reference temperature zone, the calibration parameters of each section to be calibrated are solved jointly to obtain the calibrated distributed fiber temperature inversion model.

[0092] like Figure 3 As shown, this application embodiment also provides an electronic device 300, including a memory 302 and a processor 301. The memory 302 stores a computer program, which, when executed by the processor 301, implements the method described in any of the preceding claims.

[0093] Specifically, the processor 301 can be a central processing unit (CPU), digital signal processor (DSP), microprocessor, or application-specific integrated circuit (ASIC) with data processing capabilities, used to execute various instructions and complete data processing tasks. The memory 302 can include volatile memory and / or non-volatile memory, such as random access memory (RAM), read-only memory (ROM), flash memory, or solid-state drive (SSD), used to store program code and intermediate data generated during operation.

[0094] In some embodiments, the electronic device 300 may further include a data acquisition interface module (not shown) for receiving scattering response data from a distributed fiber optic temperature measurement system; a communication module for data interaction with an external system; and a display module for outputting calibration results or temperature distribution information.

[0095] This application also provides a storage medium having a computer program stored thereon, which, when executed by a processor, implements the method as described in any of the preceding claims.

[0096] This application also provides a computer program product, including a computer program that, when executed by a processor, implements the method as described in any of the preceding claims.

[0097] The above description represents the preferred embodiments of the present invention. It should be noted that, for those skilled in the art, various improvements and modifications can be made without departing from the principles of the present invention, and these improvements and modifications are also considered to be within the scope of protection of the present invention.

Claims

1. A calibration method for a distributed optical fiber temperature measurement system, characterized in that, Includes the following steps: S1: Collect the scattering response data along the temperature measuring fiber of the distributed fiber temperature measurement system under the preset calibration conditions, and construct the original calibration response sequence by combining it with the temperature information of the reference temperature zone. S2, based on the original calibration response sequence, extract candidate link disturbance locations, identify target disturbance locations caused by joints or fusion splices among the candidate link disturbance locations, and determine the corresponding joint or fusion splice disturbance type according to the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics and channel response asymmetry characteristics of each target disturbance location; S3. Based on the disturbance type of each joint or weld point, determine the corresponding local response distortion, inter-segment baseline offset, subsequent segment slope disturbance, and calibration weight adjustment, and construct the influence propagation relationship of each joint or weld point on the calibration parameters of the subsequent segment. S4. Using the target disturbance location as a segmentation reference, the temperature measuring fiber is reconstructed into a segmented calibration function based on the influence propagation relationship. The calibration parameters of each segment are jointly solved by combining the segment continuity constraint and the reference temperature zone constraint to obtain the calibrated distributed fiber temperature inversion model.

2. The calibration method for a distributed optical fiber temperature measurement system according to claim 1, characterized in that, Candidate link disturbance locations are extracted based on the original calibration response sequence. Target disturbance locations caused by joints or fusion splices are identified among these candidate locations. The corresponding joint or fusion splice disturbance type is determined based on the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics, and channel response asymmetry characteristics of each target disturbance location, including: S21. Based on the original calibration response sequence, calculate the mean difference of the response before and after each sampling position, the local peak amplitude, the change in attenuation slope, and the difference in the dual-channel response within a preset window to obtain the disturbance characterization parameters corresponding to each sampling position. S22, Identify candidate link disturbance locations that meet preset conditions based on the disturbance characterization parameters, and select stable abnormal locations that repeatedly appear in multiple calibration acquisition results from the candidate link disturbance locations as the target disturbance locations; S23. For each of the target disturbance locations, based on the corresponding local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics, and channel response asymmetry characteristics, determine the corresponding joint or fusion point disturbance type, including at least one of pure attenuation step disturbance, reflection spike coupled attenuation disturbance, and mode redistribution disturbance.

3. The calibration method for a distributed optical fiber temperature measurement system according to claim 1, characterized in that, Construct the propagation relationship of the influence of each joint or weld point on the calibration parameters of subsequent sections, including: S31, taking each of the target disturbance positions as the propagation starting point, determine the corresponding subsequent influence section according to the extension direction of the temperature measuring fiber; determine the influence intensity of each target disturbance position on each subsequent influence section based on the disturbance type of the joint or fusion splice and the corresponding local response distortion, inter-segment baseline offset, subsequent section slope disturbance amount and calibration weight adjustment amount. S32, combining the relative distance between the target disturbance location and each subsequent affected section, the signal quality of the section, and the superimposed influence of other subsequent target disturbance locations, the influence intensity is distributed in a decreasing manner to obtain the influence propagation relationship of each joint or fusion point on the calibration parameters of the subsequent section.

4. The calibration method for a distributed optical fiber temperature measurement system according to claim 3, characterized in that, By combining the relative distance between the target disturbance location and each subsequent affected segment, the signal quality of the segment, and the superimposed influence of other subsequent target disturbance locations, the influence intensity is progressively distributed to obtain the propagation relationship of the influence of each joint or fusion point on the calibration parameters of the subsequent segment, including: S321, Based on the relative distance between the target disturbance location and each subsequent affected segment and the signal quality parameters of each subsequent affected segment, determine the basic influence intensity of each target disturbance location on each subsequent affected segment; S322, combining the spatial distribution relationship of other target disturbance positions and the corresponding basic influence intensity, the basic influence intensity is reduced and superimposed to obtain the influence propagation relationship of each target disturbance position on the calibration parameters of each subsequent influence section.

5. The calibration method for a distributed optical fiber temperature measurement system according to claim 1, characterized in that, Using the target disturbance location as a segmentation reference, the temperature-measuring fiber is reconstructed using a segmented calibration function based on the influence propagation relationship. The calibration parameters for each segment are then jointly solved by combining segment continuity constraints and reference temperature zone constraints to obtain the calibrated distributed fiber temperature inversion model, including: S41, the temperature measuring fiber is divided into multiple sections to be calibrated using the target disturbance positions as boundaries, and a piecewise calibration function containing a baseline correction term, a slope correction term, and a weight correction term is constructed for each section to be calibrated; the influence propagation relationship is introduced into the piecewise calibration function of the corresponding section to be calibrated to make preliminary corrections to the calibration parameters of each section to be calibrated. S42, combining the segment continuity constraints at the boundary of adjacent calibrated segments and the temperature constraints of the reference temperature zone, the calibration parameters of each calibrated segment are jointly solved to obtain the calibrated distributed fiber temperature inversion model.

6. A calibration device for a distributed fiber optic temperature measurement system, characterized in that, The device includes: The original sequence construction module is used to collect the scattering response data along the temperature measuring fiber of the distributed fiber temperature measurement system under the preset calibration conditions, and to construct the original calibration response sequence by combining the temperature information of the reference temperature zone. The disturbance identification module is used to extract candidate link disturbance locations based on the original calibration response sequence, identify target disturbance locations caused by joints or fusion splices among the candidate link disturbance locations, and determine the corresponding joint or fusion splice disturbance type based on the local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics and channel response asymmetry characteristics of each target disturbance location. The propagation relationship construction module is used to determine the corresponding local response distortion, inter-segment baseline offset, subsequent segment slope disturbance and calibration weight adjustment based on the disturbance type of each joint or weld point, and to construct the propagation relationship of the influence of each joint or weld point on the calibration parameters of the subsequent segment. The model reconstruction module is used to reconstruct the segmented calibration function of the temperature measuring fiber based on the target disturbance location as the segment division reference and the influence propagation relationship. It also combines the segment continuity constraint and the reference temperature zone constraint to jointly solve the calibration parameters of each segment, thereby obtaining the calibrated distributed fiber temperature inversion model.

7. The calibration device for a distributed optical fiber temperature measurement system according to claim 6, characterized in that, The disturbance identification module is configured to: Based on the original calibration response sequence, the mean difference of the response, the local peak amplitude, the change in attenuation slope, and the difference in the dual-channel response within the preset window before and after each sampling position are calculated to obtain the disturbance characterization parameters corresponding to each sampling position. Candidate link disturbance locations that meet preset conditions are identified based on the disturbance characterization parameters, and stable abnormal locations that repeatedly appear in multiple calibration acquisition results are selected from the candidate link disturbance locations as the target disturbance locations. For each target disturbance location, based on the corresponding local step characteristics, reflection spike characteristics, attenuation gradient abrupt change characteristics, and channel response asymmetry characteristics, the corresponding joint or fusion point disturbance type is determined, including at least one of pure attenuation step disturbance, reflection spike coupled attenuation disturbance, and mode redistribution disturbance.

8. An electronic device, characterized in that, The method includes a memory and a processor, wherein the memory stores a computer program, and the computer program, when executed by the processor, implements the method as described in any one of claims 1-5.

9. A storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the method as described in any one of claims 1-5.

10. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the method as described in any one of claims 1-5.

Images