Pipeline siltation risk assessment method and system

Abstract

The invention relates to the technical field of pipeline siltation risk assessment, and particularly discloses a pipeline siltation risk assessment method and system, and the method comprises the steps: setting each transmitting angle, operating a sound wave transmitting device of each sub-pipeline section to transmit corresponding sound waves into each sub-pipeline section according to each transmitting angle, measuring a time difference value from sound wave emission to sound wave receiving in each sub-pipeline section by the sound wave emission device according to each emission angle, thereby determining a preliminary sedimentation evaluation coefficient corresponding to each emission angle in each sub-pipeline section; meanwhile, whether deposits exist in all the sub-pipeline sections or not is determined; according to the method, the sub-pipeline sections with the deposits are marked as target pipeline sections, then the deposition risk assessment coefficients of all the target pipeline sections are comprehensively analyzed and fed back to the pipeline background terminal, detection can be completed under the condition that normal operation of the pipeline is not interfered, and the real-time performance and safety of detection are ensured; and accurate positioning of a problem area is facilitated, and unnecessary maintenance of the whole pipeline is avoided, so that the detection efficiency and the resource utilization rate are improved.

Description

Technical Field

[0001] The invention belongs to the technical field of pipeline siltation risk assessment, and relates to a pipeline siltation risk assessment method and system. Background Art

[0002] Pipeline systems have been widely used in industries such as industry and agriculture. However, the siltation problem in pipelines often leads to poor fluid flow and even causes pipeline blockage, which seriously affects production efficiency and the stability of system operation. Traditional pipeline siltation detection methods mainly rely on manual inspection or flow monitoring, but these methods have many technical defects: first, manual inspection usually requires shutdown operation and cannot achieve online monitoring, resulting in greater interference of the detection process on pipeline operation; second, manual inspection relies on experience judgment, and the accuracy and consistency of the detection results are difficult to guarantee, especially in complex pipeline environments, which are prone to missed detection or misjudgment; and although the flow monitoring-based method can provide a certain degree of siltation prompts, its detection depends on the change of the overall flow rate, and it is impossible to accurately locate the location and degree of siltation. At the same time, in the case of slight siltation or local siltation, the flow change may not be obvious, resulting in insufficient detection sensitivity. In addition, these methods generally lack real-time performance, making it difficult to continuously monitor and early warn the dynamic changes of siltation in the pipeline, and cannot meet the needs of modern industry and agriculture for efficient, accurate and intelligent management of pipeline systems. Summary of the invention

[0003] In view of the above problems existing in the prior art, the present invention provides a pipeline sedimentation risk assessment method and system for solving the above technical problems.

[0004] In order to achieve the above purpose and other purposes, the technical solution adopted by the present invention is as follows:

[0005] In one aspect, the present invention provides a pipeline siltation risk assessment method, the method comprising the following steps:

[0006] Step 1: Obtain the total number of sections of the target pipeline, and redivide the target pipeline into sub-pipeline sections based on the set number of sections, respectively arrange a sound wave transmitting device and a sound wave receiving device at the starting position of each sub-pipeline section, set each emission angle corresponding to each sub-pipeline section, operate the sound wave transmitting device of each sub-pipeline section to transmit corresponding sound waves into each sub-pipeline section according to each emission angle, and place a sound wave reflecting device at the end position of each sub-pipeline section;

[0007] Step 2: measuring the time difference from the emission of the sound wave to the reception of the sound wave in each sub-pipeline section according to each emission angle, thereby determining the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section; and determining whether there is sediment in each sub-pipeline section;

[0008] Step 3: Record the sub-pipeline section with sediment as the target pipeline section, and then comprehensively analyze the sedimentation risk assessment coefficient of each target pipeline section and feed it back to the pipeline backend terminal.

[0009] Set each sub-pipeline segment to correspond to each emission angle. The specific setting process includes the following steps:

[0010] The test launch angle is recorded as θ i , where i is the number of each test launch, and θ i ∈[θ min ,θ max ],θ min ,θ max They are the set angle lower threshold and angle upper threshold respectively;

[0011] Select any one of the sub-pipeline segments as the target sub-pipeline segment;

[0012] The digital twin technology is used to generate a model of the target sub-pipeline segment, and then the target sub-pipeline segment model is simulated at each test launch angle. The comprehensive evaluation value β(θ i )=ω1 i β1(θ i )-ω2 i β2(θ i )-ω3 i β3(θ i )-ω4 i β4(θ i )+τ, where ω1 i ω2 i ,ω3 i ,ω4 i Respectively represent the preset weight factors for each test emission angle, which are used to balance the influence of different characteristics and satisfy ω1 i +ω2 i +ω3 i +ω4 i =1, τ is the set correction coefficient;

[0013] In the above formula, β1(θ i ) is the total energy receiving intensity corresponding to the i-th test emission angle in the target sub-pipeline segment, and the specific calculation process is as follows:

[0014] Import the actual length L' of the target sub-pipeline segment into the strength calculation model The total energy receiving intensity β1(θ i ), R1(x,θ i) represents the receiving intensity of the sound wave at the first reflection at position x in the target sub-pipeline segment under the i-th test transmission angle; R n (x,θ i ) represents the receiving intensity of the sound wave at the position x in the target sub-pipeline segment at the nth reflection under the i-th test transmission angle; n is the number of reflections of the sound wave, e is a natural constant, d is used to represent the continuous change of x in the integral, φ is the reflection attenuation factor, which is used to describe the attenuation of energy after each reflection, N max is the maximum number of reflections of the sound wave;

[0015] The total attenuation factor β2(θ i ), ψ1 and ψ2 decibels represent the set medium absorption coefficient and angle attenuation correction coefficient, d(x,θ i ) is the propagation distance of the sound wave in the target sub-pipeline section at the i-th test emission angle;

[0016] The propagation path complexity factor β3(θ i ), is the path complexity weight factor;

[0017] Then the propagation time corresponding to each test launch angle in the target sub-pipeline section is calculated. and They respectively represent the propagation speed of the sound wave in the current medium and the correction coefficient of the path angle to the propagation time;

[0018] The comprehensive evaluation value corresponding to each test emission angle in the target sub-pipeline segment is compared with the preset evaluation threshold. If the comprehensive evaluation value corresponding to a test emission angle in the target sub-pipeline segment is greater than or equal to the preset evaluation threshold, the corresponding test emission angle in the target sub-pipeline segment is used as the emission angle of the corresponding sub-pipeline segment, otherwise the angle is discarded;

[0019] The emission angles corresponding to the sub-pipeline segments are set in the same manner as the target sub-pipeline segments.

[0020] The calculation formula for the preset weight factor for each test launch angle is as follows:

[0021] ω1 i =ω10·(1+k1·cos(θ i ));

[0022] ω2 i =ω20·(1+k2·sin(θ i ));

[0023] ω3 i =ω30·(1+k3·tan(θ i ));

[0024] ω4i=ω40·(1-k4·cos 2 (θi));

[0025] In the above formula, ω10, ω20, ω30, and ω40 are the initial values ​​of the weight factors at the reference angle, and k1, k2, k3, and k4 are the coefficients for controlling the range of weight change.

[0026] Determine the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section. The specific determination process includes:

[0027] Combined with the average temperature and pressure of each sub-pipeline segment during the corresponding time period from the emission of the sound wave to the reception of the sound wave, the sound wave propagation speed in each sub-pipeline segment is corrected.

[0028]

[0029] In the above formula, j is the number of each sub-pipeline section, υ0 represents the speed of sound under standard conditions, α1 and α2 represent the correction coefficients of temperature and pressure respectively, T and P represent the temperature and pressure under standard conditions respectively, are the mean temperature and mean pressure of the j-th sub-pipeline segment in the corresponding time period from the emission of the sound wave to the reception of the sound wave;

[0030] Then calculate the propagation distance of the sound wave in each sub-pipeline section from each emission angle Δt jp is the time difference from the emission of the sound wave to the reception of the sound wave in the jth sub-pipeline section from the pth emission angle, and p is the number of each emission angle;

[0031] The preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section is calculated from this Among them, judgment condition A means that the difference between the propagation distance of the sound wave in each sub-pipeline segment at each emission angle and the actual distance of the pipeline length of each sub-pipeline segment is less than the set allowable threshold, ZF j is the reference amplitude in the jth sub-pipeline section when there is no siltation, ZF j ' p represents the amplitude of the sound wave received in the jth sub-pipeline section from the pth emission angle; Indicates the symbol NOT.

[0032] The logic for determining whether there is sediment in each sub-pipeline section is:

[0033] The preliminary sedimentation assessment coefficient in each sub-pipeline segment is compared with the set sedimentation assessment threshold coefficient. If the preliminary sedimentation assessment coefficient in a sub-pipeline segment is greater than the set sedimentation assessment threshold coefficient, it is determined that there is sediment in the sub-pipeline segment; otherwise, it is determined that there is no sediment in the sub-pipeline segment.

[0034] Comprehensively analyze the sedimentation risk assessment coefficient of each target pipeline section, including:

[0035] Select any one of the target pipeline segments as the target calculation pipeline segment;

[0036] Based on this target, the sedimentation risk assessment coefficient of the pipeline section is calculated L is the actual distance of the pipeline length of the target calculation pipeline segment, x is the coordinate of any position in the target calculation pipeline segment, Φ(x) is the estimated risk assessment function, and the specific definition formula is as follows:

[0037] Φ(x)=ζ1×ΔS1(x)+ζ1×ΔS2(x)+ζ1×ΔS3(x)+ζ1×ΔS4(x)+τ1;

[0038] Wherein, ζ1, ζ2, ζ3 and ζ4 represent the set weighted indexes respectively, and satisfy ζ1+ζ2+ζ3+ζ4=1; ΔS1(x), ΔS2(x), ΔS3(x), ΔS4(x) are the change of the acoustic wave reflection intensity, the change of the acoustic wave reflection time, the change of the acoustic wave scattering intensity and the distribution value of the inner wall reflection coefficient of the target calculated pipeline section respectively, and τ1 is the set environmental correction coefficient;

[0039] The siltation risk assessment coefficient of each target pipeline section is calculated in the same way as the calculation method of the siltation risk assessment coefficient of the target pipeline section.

[0040] The goal is to calculate the change in the sound wave reflection intensity, the change in the sound wave reflection time, the change in the sound wave scattering intensity and the distribution of the inner wall reflection coefficient of the pipeline section. The specific calculation formula is:

[0041] Objective Calculate the change in the sound wave reflection intensity of the pipeline section Where b is the number of each emission angle, b = 1, 2, ... N, N is the total number of emission angles, S1 b (x) represents the reflection intensity measured at position x inside the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S1(x) is the reflection intensity baseline at position x inside the target calculation pipeline segment under normal conditions;

[0042] Objective Calculate the change in the acoustic reflection time of the pipeline section S2 b(x) represents the sound wave reflection time measured at the internal position x of the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S2(x) is the sound wave reflection time baseline at the internal position x of the target calculation pipeline segment under normal conditions;

[0043] Objective Calculate the change in the sound wave scattering intensity of the pipeline section S3 i (x) represents the sound wave scattering intensity measured at position x inside the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S3(x) is the sound wave scattering intensity baseline at position x inside the target calculation pipeline segment under normal conditions;

[0044] Objective Calculate the inner wall reflection coefficient distribution value of the pipeline segment S4(x) is the maximum reflection intensity at position x inside the target calculated pipeline segment under normal conditions. Calculate the average reflection intensity at different angles inside the pipe segment for the target.

[0045] Another aspect of the present invention provides a pipeline siltation risk assessment system, which includes a launch angle setting module, a siltation preliminary judgment module and a siltation risk assessment module, wherein the above modules are connected by wired and / or wireless connection to achieve data transmission between the modules;

[0046] Emission angle setting module: obtain the total number of sections of the target pipeline, and re-divide the target pipeline into sub-pipeline sections based on the set number of sections, and set the emission angles corresponding to each sub-pipeline section, and operate the sound wave emission device of each sub-pipeline section to emit corresponding sound waves into each sub-pipeline section according to each emission angle;

[0047] Preliminary sedimentation judgment module: measures the time difference from the emission of the sound wave to the reception of the sound wave in each sub-pipeline section at each emission angle, thereby determining the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section; and at the same time determines whether there is sediment in each sub-pipeline section;

[0048] Sedimentation risk assessment module: The sub-pipeline segment with sediment is recorded as the target pipeline segment, and then the sedimentation risk assessment coefficient of each target pipeline segment is comprehensively analyzed and fed back to the pipeline background terminal.

[0049] As described above, the pipeline siltation risk assessment method and system provided by the present invention have at least the following beneficial effects:

[0050] The present invention provides a pipeline siltation risk assessment method and system, which transmits sound waves into each sub-pipeline segment by setting each emission angle and operating the sound wave transmitting device, measures the time difference between the sound wave emission and reception to determine the preliminary siltation assessment coefficient, further identifies the target pipeline segment, and then comprehensively analyzes the siltation risk assessment coefficients of these target pipeline segments and feeds them back to the pipeline background terminal. This process has important technical necessity and significant technical advantages: first, through multi-angle sound wave detection, the pipeline cross section can be fully covered, and the recognition accuracy of the distribution and thickness of the silt can be improved; secondly, using the time difference to evaluate the siltation situation is a non-contact detection method, which can complete the detection without interfering with the normal operation of the pipeline, ensuring the real-time and safety of the detection; thirdly, the segmented evaluation of the siltation situation of the sub-pipeline segment is helpful to accurately locate the problem area and avoid unnecessary maintenance of the entire pipeline, thereby improving the detection efficiency and resource utilization; it can realize dynamic monitoring and data management of the pipeline status, provide a scientific basis for subsequent maintenance decisions, further support the intelligent management and predictive maintenance of the pipeline, and reduce the risk of pipeline blockage or reduced operating efficiency due to siltation. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for describing the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0052] Figure 1 It is a schematic diagram of the connection of each step of the method of the present invention.

[0053] Figure 2 It is a schematic diagram of the connection of various modules of the system of the present invention. DETAILED DESCRIPTION

[0054] The above contents in combination with the implementation of the present invention are merely examples and explanations of the concept of the present invention. The technicians in this technical field may make various modifications or additions to the specific embodiments described or replace them in a similar manner. As long as they do not deviate from the concept of the invention or exceed the scope defined by the claims, they shall all fall within the protection scope of the present invention.

[0055] Example 1

[0056] See also Figure 1 As shown, a pipeline siltation risk assessment method comprises the following steps:

[0057] Step 1: Obtain the total number of sections of the target pipeline, and redivide the target pipeline into sub-pipeline sections based on the set number of sections, respectively arrange a sound wave transmitting device and a sound wave receiving device at the starting position of each sub-pipeline section, set each emission angle corresponding to each sub-pipeline section, operate the sound wave transmitting device of each sub-pipeline section to transmit corresponding sound waves into each sub-pipeline section according to each emission angle, and place a sound wave reflecting device at the end position of each sub-pipeline section;

[0058] Step 2: measuring the time difference from the emission of the sound wave to the reception of the sound wave in each sub-pipeline section according to each emission angle, thereby determining the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section; and determining whether there is sediment in each sub-pipeline section;

[0059] Step 3: Record the sub-pipeline section with sediment as the target pipeline section, and then comprehensively analyze the sedimentation risk assessment coefficient of each target pipeline section and feed it back to the pipeline backend terminal.

[0060] Exemplarily, each sub-pipeline segment is set to correspond to each emission angle, and the specific setting process includes the following steps:

[0061] The test launch angle is recorded as θ i , where i is the number of each test launch, and θ i ∈[θ min ,θ max ],θ min ,θ max are the set angle lower threshold and angle upper threshold respectively; θ min =0,θ max =360;

[0062] Select any one of the sub-pipeline segments as the target sub-pipeline segment;

[0063] The digital twin technology is used to generate a model of the target sub-pipeline segment, and then the target sub-pipeline segment model is simulated at each test launch angle. The comprehensive evaluation value β(θ i )=ω1 i β1(θ i )-ω2 i β2(θ i )-ω3 i β3(θ i )-ω4 i β4(θ i )+τ, where ω1 i ω2 i ,ω3 i ,ω4 iRespectively represent the preset weight factors for each test emission angle, which are used to balance the influence of different characteristics and satisfy ω1 i +ω2 i +ω3 i +ω4 i =1, τ is the set correction coefficient;

[0064] In the above formula, β1(θ i ) is the total energy receiving intensity corresponding to the i-th test emission angle in the target sub-pipeline segment, and the specific calculation process is as follows:

[0065] Import the actual length L' of the target sub-pipeline segment into the strength calculation model The total energy receiving intensity β1(θ i ), R1(x,θ i ) represents the receiving intensity of the sound wave at the first reflection at position x in the target sub-pipeline segment under the i-th test transmission angle; R n (x,θ i ) represents the receiving intensity of the sound wave at the position x in the target sub-pipeline segment at the nth reflection under the i-th test transmission angle; n is the number of reflections of the sound wave, e is a natural constant, d is used to represent the continuous change of x in the integral, φ is the reflection attenuation factor, which is used to describe the attenuation of energy after each reflection, N max is the maximum number of reflections of the sound wave;

[0066] In the formula, n represents the number of reflections of the sound wave. In actual applications, when the sound wave propagates in the pipeline, it may be reflected multiple times on the pipeline wall. Each reflection will cause the attenuation of the sound wave energy, so when calculating the sound wave reception intensity, it is necessary to consider all possible reflection times. Specifically, n = 1 represents the first reflection, the intensity of the sound wave emitted from the transmitting device, directly reaching the reflecting device and returning; n = 2 represents the second reflection, the intensity of the sound wave reflected once, reflected again and returned; n = 3 represents the third reflection, the intensity of the sound wave reflected again after two reflections.

[0067] x can be understood as a position coordinate within the pipeline segment, usually a variable along the length of the pipeline. This variable is used to describe the specific position of the sound wave when it propagates in the pipeline. In practical applications, the pipeline is a one-dimensional space, and the sound wave propagates from the transmitter to different positions along the pipeline. x is the position on this propagation path; assuming the pipeline length is L, the value range of x is usually [0, L]. When calculating, we will integrate the sound wave reflection intensity at different positions within this range to obtain the total receiving intensity of the entire pipeline segment.

[0068] The total attenuation factor β2(θ i ), ψ1 and ψ2 decibels represent the set medium absorption coefficient and angle attenuation correction coefficient, d(x,θ i ) is the propagation distance of the sound wave in the target sub-pipeline section at the i-th test emission angle;

[0069] The propagation path complexity factor β3(θ i ), is the path complexity weight factor;

[0070] Then the propagation time corresponding to each test launch angle in the target sub-pipeline section is calculated. and They respectively represent the propagation speed of the sound wave in the current medium and the correction coefficient of the path angle to the propagation time;

[0071] The comprehensive evaluation value corresponding to each test emission angle in the target sub-pipeline segment is compared with the preset evaluation threshold. If the comprehensive evaluation value corresponding to a test emission angle in the target sub-pipeline segment is greater than or equal to the preset evaluation threshold, the corresponding test emission angle in the target sub-pipeline segment is used as the emission angle of the corresponding sub-pipeline segment, otherwise the angle is discarded;

[0072] The emission angles corresponding to the sub-pipeline segments are set in the same manner as the target sub-pipeline segments.

[0073] The design of the comprehensive evaluation value is based on the actual physical characteristics of sound wave propagation in pipes and application requirements. The construction of this formula takes into account several key factors to ensure that the effectiveness and applicability of each launch angle can be fully evaluated:

[0074] The total energy reception intensity is the core indicator for evaluating the reception effect of sound waves at a specific angle. It reflects the amount of energy that the sound wave finally reaches the receiving device after propagation and reflection. A higher energy reception intensity means that the angle is more effective in practical applications;

[0075] The total attenuation factor is the energy attenuation caused by the absorption and obstruction of the medium during the propagation of sound waves. Therefore, the attenuation factor is used to quantify the energy lost by sound waves during propagation. A higher attenuation factor will reduce the reception intensity, so this value needs to be subtracted in the comprehensive evaluation to reflect the actual reception effect;

[0076] The propagation path complexity factor is that the path of sound waves in the pipeline may not be a straight line, especially when there are bends or changes in diameter. The complex propagation path will cause reflection and scattering of sound waves, affecting the final reception effect. The complexity factor is used to evaluate the nonlinearity of the path and the impact of multiple reflections, so this value also needs to be subtracted in the evaluation.

[0077] The propagation time reflects the time required for the sound wave to be transmitted and received. A longer propagation time may mean signal delay, affecting real-time performance, especially in dynamic monitoring. Therefore, incorporating the propagation time into the evaluation and subtracting it will help select a more efficient transmission angle.

[0078] Exemplarily, the calculation formula of the preset weight factor for each test emission angle is as follows:

[0079] ω1 i =ω10·(1+k1·cos(θ i ));

[0080] ω2 i =ω20·(1+k2·sin(θ i ));

[0081] ω3 i =ω30·(1+k3·tan(θ i ));

[0082] ω4i=ω40·(1-k4·cos 2 (θi));

[0083] In the above formula, ω10, ω20, ω30, and ω40 are the initial values ​​of the weight factors at the reference angle, and k1, k2, k3, and k4 are the coefficients for controlling the range of weight change.

[0084] Making the weight factor change dynamically with the emission angle has important benefits and necessity in the calculation of the comprehensive evaluation value. This dynamic weight adjustment can more accurately reflect the propagation characteristics and environmental influences of sound waves at different angles, thereby improving the accuracy and reliability of the evaluation. For example, when the emission angle changes, the attenuation, reflection and scattering characteristics of the sound wave will also change. By dynamically adjusting the weight factor, the impact of these changes on the sound wave propagation effect can be better captured. In addition, this method can also effectively consider the interference of environmental noise, ensure the detectability and signal-to-noise ratio of the signal under complex conditions, and thus improve the overall performance and stability of the system. Therefore, the design of dynamic weight factors not only enhances the adaptability of the evaluation model, but also provides a more scientific basis for optimizing sound wave detection and analysis.

[0085] Exemplarily, the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section is determined, and the specific determination process includes:

[0086] Combined with the average temperature and pressure of each sub-pipeline segment during the corresponding time period from the emission of the sound wave to the reception of the sound wave, the sound wave propagation speed in each sub-pipeline segment is corrected.

[0087]

[0088] In the above formula, j is the number of each sub-pipeline section, υ0 represents the speed of sound under standard conditions, α1 and α2 represent the correction coefficients of temperature and pressure respectively, T and P represent the temperature and pressure under standard conditions respectively, are the mean temperature and mean pressure of the j-th sub-pipeline segment in the corresponding time period from the emission of the sound wave to the reception of the sound wave;

[0089] Then calculate the propagation distance of the sound wave in each sub-pipeline section from each emission angle Δt jp is the time difference from the emission of the sound wave to the reception of the sound wave in the jth sub-pipeline section from the pth emission angle, and p is the number of each emission angle;

[0090] The preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section is calculated from this Among them, judgment condition A means that the difference between the propagation distance of the sound wave in each sub-pipeline segment at each emission angle and the actual distance of the pipeline length of each sub-pipeline segment is less than the set allowable threshold, ZF j is the reference amplitude in the jth sub-pipeline section when there is no siltation, ZF j ' p represents the amplitude of the sound wave received in the jth sub-pipeline section from the pth emission angle; Indicates the symbol NOT.

[0091] Exemplarily, the logic for determining whether there is sediment in each sub-pipeline segment is:

[0092] The preliminary sedimentation assessment coefficient in each sub-pipeline segment is compared with the set sedimentation assessment threshold coefficient. If the preliminary sedimentation assessment coefficient in a sub-pipeline segment is greater than the set sedimentation assessment threshold coefficient, it is determined that there is sediment in the sub-pipeline segment; otherwise, it is determined that there is no sediment in the sub-pipeline segment.

[0093] Exemplarily, the siltation risk assessment coefficient of each target pipeline section is comprehensively analyzed, including:

[0094] Select any one of the target pipeline segments as the target calculation pipeline segment;

[0095] Based on this target, the sedimentation risk assessment coefficient of the pipeline section is calculated L is the actual distance of the target calculation pipeline segment, x is the coordinate of any position in the target calculation pipeline segment, Φ(x) is the estimated risk assessment function, and the specific definition formula is as follows:

[0096] Φ(x)=ζ1×ΔS1(x)+ζ1×ΔS2(x)+ζ1×ΔS3(x)+ζ1×ΔS4(x)+τ1;

[0097] Wherein, ζ1, ζ2, ζ3 and ζ4 represent the set weighted indexes respectively, and satisfy ζ1+ζ2+ζ3+ζ4=1; ΔS1(x), ΔS2(x), ΔS3(x), ΔS4(x) are the change of the acoustic wave reflection intensity, the change of the acoustic wave reflection time, the change of the acoustic wave scattering intensity and the distribution value of the inner wall reflection coefficient of the target calculated pipeline section respectively, and τ1 is the set environmental correction coefficient;

[0098] The siltation risk assessment coefficient of each target pipeline section is calculated in the same way as the calculation method of the siltation risk assessment coefficient of the target pipeline section.

[0099] Exemplarily, the target is to calculate the change in the sound wave reflection intensity, the change in the sound wave reflection time, the change in the sound wave scattering intensity and the distribution of the inner wall reflection coefficient of the pipeline section. The specific calculation formula is:

[0100] Objective Calculate the change in the sound wave reflection intensity of the pipeline section Where b is the number of each emission angle, b = 1, 2, ... N, N is the total number of emission angles, S1 b (x) represents the reflection intensity measured at position x inside the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S1(x) is the reflection intensity baseline at position x inside the target calculation pipeline segment under normal conditions;

[0101] The physical meaning of the above calculation formula is that when sound waves propagate in a pipe, they will encounter the pipe wall or sediment, and part of the sound waves will be reflected back. Reflection intensity S1 b (x) is the energy or amplitude of the reflected sound wave, which is usually related to the material, shape, and sediment distribution of the inner wall of the pipe. Different locations in the pipe may have different reflection intensities. For example, in places with sediment, the reflection intensity may be stronger or weaker than in places without sediment.

[0102] Objective Calculate the change in the acoustic reflection time of the pipeline section S2 b (x) represents the sound wave reflection time measured at the internal position x of the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S2(x) is the sound wave reflection time baseline at the internal position x of the target calculation pipeline segment under normal conditions;

[0103] When sound waves propagate in the pipe, they encounter the pipe wall or other obstacles (such as silt) and are reflected back. b (x) is the time it takes for the sound wave to travel from the transmitting point to position x and back to the receiver; the propagation time of the sound wave is related to the reflection position x and the speed of the sound wave (usually the speed of sound in the pipe). If there is sediment or other abnormalities in the pipe, the reflection position may change, thus affecting S2 b (x); Different emission angles will cause the sound waves to propagate in different paths in the pipe, so S2 b (x) will be affected by the angle θ b and change;

[0104] If there is sediment in the pipe, the sound waves may reflect on the surface of the sediment instead of on the pipe wall, which will shorten the propagation distance and make the reflection time S2 b (x) decreases. By comparing the actual measured S2 b (x) Compared with the normal sound wave reflection time baseline, it can detect whether there is siltation or other abnormalities in the pipeline.

[0105] Objective Calculate the change in the sound wave scattering intensity of the pipeline section S3 i (x) represents the sound wave scattering intensity measured at position x inside the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S3(x) is the sound wave scattering intensity baseline at position x inside the target calculation pipeline segment under normal conditions;

[0106] Objective Calculate the inner wall reflection coefficient distribution value of the pipeline segment S4(x) is the maximum reflection intensity at position x inside the target calculated pipeline segment under normal conditions. Calculate the average reflection intensity at different angles inside the pipe segment for the target.

[0107] Example 2

[0108] like Figure 2 As shown, a pipeline siltation risk assessment system includes a launch angle setting module, a siltation preliminary judgment module and a siltation risk assessment module, wherein the above modules are connected by wired and / or wireless connection to achieve data transmission between the modules;

[0109] Emission angle setting module: obtain the total number of sections of the target pipeline, and re-divide the target pipeline into sub-pipeline sections based on the set number of sections, and set the emission angles corresponding to each sub-pipeline section, and operate the sound wave emission device of each sub-pipeline section to emit corresponding sound waves into each sub-pipeline section according to each emission angle;

[0110] Preliminary sedimentation judgment module: measures the time difference from the emission of the sound wave to the reception of the sound wave in each sub-pipeline section at each emission angle, thereby determining the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section; and at the same time determines whether there is sediment in each sub-pipeline section;

[0111] Sedimentation risk assessment module: The sub-pipeline segment with sediment is recorded as the target pipeline segment, and then the sedimentation risk assessment coefficient of each target pipeline segment is comprehensively analyzed and fed back to the pipeline background terminal.

[0112] It should be understood that in the various embodiments of the present application, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0113] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art who is familiar with the present technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

[0114] Finally: The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A pipeline siltation risk assessment method, characterized in that: The steps include: Step 1: Obtain the total number of sections of the target pipeline, and re-divide the target pipeline into sub-pipeline sections based on the set number of sections, and set the emission angles corresponding to the sub-pipeline sections, and operate the sound wave emission devices of the sub-pipeline sections to emit corresponding sound waves into the sub-pipeline sections according to the emission angles; Step 2: measuring the time difference from the emission of the sound wave to the reception of the sound wave in each sub-pipeline section according to each emission angle, thereby determining the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section; and determining whether there is sediment in each sub-pipeline section; Step 3: Record the sub-pipeline section with sediment as the target pipeline section, and then comprehensively analyze the sedimentation risk assessment coefficient of each target pipeline section and feed it back to the pipeline backend terminal.

2. A pipeline siltation risk assessment method according to claim 1, characterized in that: Set each sub-pipeline segment to correspond to each emission angle. The specific setting process includes the following steps: The test launch angle is recorded as θ i , where i is the number of each test launch, and θ i ∈[θ min ,θ max ],θ min ,θ max They are the set angle lower threshold and angle upper threshold respectively; Select any one of the sub-pipeline segments as the target sub-pipeline segment; The digital twin technology is used to generate a model of the target sub-pipeline segment, and then the target sub-pipeline segment model is simulated at each test launch angle. The comprehensive evaluation value β(θ i )=ω1 i β1(θ i )-ω2 i β2(θ i )-ω3 i β3(θ i )-ω4 i β4(θ i )+τ, where ω1 i ω2 i ,ω3 i ,ω4 i Respectively represent the preset weight factors for each test emission angle, which are used to balance the influence of different characteristics and satisfy ω1 i +ω2 i +ω3 i +ω4 i =1, τ is the set correction coefficient; In the above formula, β1(θ i ) is the total energy receiving intensity corresponding to the i-th test emission angle in the target sub-pipeline segment, and the specific calculation process is as follows: Import the actual length L' of the target sub-pipeline segment into the strength calculation model The total energy receiving intensity β1(θ i ), R1(x,θ i ) represents the receiving intensity of the sound wave at the first reflection at position x in the target sub-pipeline segment under the i-th test transmission angle; R n (x,θ i ) represents the receiving intensity of the sound wave at the position x in the target sub-pipeline segment at the nth reflection under the i-th test transmission angle; n is the number of reflections of the sound wave, e is a natural constant, d is used to represent the continuous change of x in the integral, φ is the reflection attenuation factor, which is used to describe the attenuation of energy after each reflection, N max is the maximum number of reflections of the sound wave; The total attenuation factor β2(θ i ), ψ1 and ψ2 decibels represent the set medium absorption coefficient and angle attenuation correction coefficient, d(x,θ i ) is the propagation distance of the sound wave in the target sub-pipeline section at the i-th test emission angle; The propagation path complexity factor β3(θ i ), is the path complexity weight factor; Then the propagation time corresponding to each test launch angle in the target sub-pipeline section is calculated. and They respectively represent the propagation speed of the sound wave in the current medium and the correction coefficient of the path angle to the propagation time; The comprehensive evaluation value corresponding to each test emission angle in the target sub-pipeline segment is compared with the preset evaluation threshold. If the comprehensive evaluation value corresponding to a test emission angle in the target sub-pipeline segment is greater than or equal to the preset evaluation threshold, the corresponding test emission angle in the target sub-pipeline segment is used as the emission angle of the corresponding sub-pipeline segment, otherwise the angle is discarded; The emission angles corresponding to the sub-pipeline segments are set in the same manner as the target sub-pipeline segments.

3. A pipeline siltation risk assessment method according to claim 2, characterized in that: The calculation formula for the preset weight factor for each test launch angle is as follows: ω1 i =ω10·(1+k1·cos(θ i )); ω2 i =ω20·(1+k2·sin(θ i )); ω3 i =ω30·(1+k3·tan(θ i )); ω4i=ω40·(1-k4·cos 2 (θi)); In the above formula, ω10, ω20, ω30, and ω40 are the initial values ​​of the weight factors at the reference angle, and k1, k2, k3, and k4 are the coefficients for controlling the range of weight change.

4. A pipeline siltation risk assessment method according to claim 1, characterized in that: Determine the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section. The specific determination process includes: Combined with the average temperature and pressure of each sub-pipeline segment during the corresponding time period from the emission of the sound wave to the reception of the sound wave, the sound wave propagation speed in each sub-pipeline segment is corrected. In the above formula, j is the number of each sub-pipeline section, υ0 represents the speed of sound under standard conditions, α1 and α2 represent the correction coefficients of temperature and pressure respectively, T and P represent the temperature and pressure under standard conditions respectively, are the mean temperature and mean pressure of the j-th sub-pipeline segment in the corresponding time period from the emission of the sound wave to the reception of the sound wave; Then calculate the propagation distance of the sound wave in each sub-pipeline section from each emission angle Δt jp is the time difference from the sound wave emission to the sound wave reception in the jth sub-pipeline section from the pth emission angle, and p is the number of each emission angle; The preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section is calculated from this Among them, judgment condition A means that the difference between the propagation distance of the sound wave in each sub-pipeline segment at each emission angle and the actual distance of the pipeline length of each sub-pipeline segment is less than the set allowable threshold, ZF j is the reference amplitude in the jth sub-pipeline section when there is no siltation, ZF j ' p represents the amplitude of the sound wave received in the jth sub-pipeline section from the pth emission angle; Indicates the symbol NOT.

5. A pipeline siltation risk assessment method according to claim 1, characterized in that: The logic for determining whether there is sediment in each sub-pipeline section is: The preliminary sedimentation assessment coefficient in each sub-pipeline segment is compared with the set sedimentation assessment threshold coefficient. If the preliminary sedimentation assessment coefficient in a sub-pipeline segment is greater than the set sedimentation assessment threshold coefficient, it is determined that there is sediment in the sub-pipeline segment; otherwise, it is determined that there is no sediment in the sub-pipeline segment.

6. A pipeline siltation risk assessment method according to claim 1, characterized in that: Comprehensively analyze the sedimentation risk assessment coefficient of each target pipeline section, including: Select any one of the target pipeline segments as the target calculation pipeline segment; Based on this target, the sedimentation risk assessment coefficient of the pipeline section is calculated L is the actual distance of the target calculation pipeline segment, x is the coordinate of any position in the target calculation pipeline segment, Φ(x) is the estimated risk assessment function, and the specific definition formula is as follows: Φ(x)=ζ1×ΔS1(x)+ζ1×ΔS2(x)+ζ1×ΔS3(x)+ζ1×ΔS4(x)+τ1; Wherein, ζ1, ζ2, ζ3 and ζ4 represent the set weighted indexes respectively, and satisfy ζ1+ζ2+ζ3+ζ4=1; ΔS1(x), ΔS2(x), ΔS3(x), ΔS4(x) are the change of the acoustic wave reflection intensity, the change of the acoustic wave reflection time, the change of the acoustic wave scattering intensity and the distribution value of the inner wall reflection coefficient of the target calculated pipeline section respectively, and τ1 is the set environmental correction coefficient; The siltation risk assessment coefficient of each target pipeline section is calculated in the same way as the calculation method of the siltation risk assessment coefficient of the target pipeline section.

7. A pipeline siltation risk assessment method according to claim 6, characterized in that: The goal is to calculate the change in the sound wave reflection intensity, the change in the sound wave reflection time, the change in the sound wave scattering intensity and the distribution of the inner wall reflection coefficient of the pipeline section. The specific calculation formula is: Objective Calculate the change in the sound wave reflection intensity of the pipeline section Where b is the number of each emission angle, b = 1, 2, ... N, N is the total number of emission angles, S1 b (x) represents the reflection intensity measured at position x inside the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S1(x) is the reflection intensity baseline at position x inside the target calculation pipeline segment under normal conditions; Objective Calculate the change in the acoustic reflection time of the pipeline section S2 b (x) represents the sound wave reflection time measured at the internal position x of the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S2(x) is the sound wave reflection time baseline at the internal position x of the target calculation pipeline segment under normal conditions; Objective Calculate the change in the sound wave scattering intensity of the pipeline section S3 i (x) represents the sound wave scattering intensity measured at position x inside the target calculation pipeline segment when the sound wave is emitted at the bth emission angle; S3(x) is the sound wave scattering intensity baseline at position x inside the target calculation pipeline segment under normal conditions; Objective Calculate the inner wall reflection coefficient distribution value of the pipeline segment S4(x) is the maximum reflection intensity at position x inside the target calculated pipeline segment under normal conditions. Calculate the average reflection intensity at different angles inside the pipe segment for the target.

8. A pipeline siltation risk assessment system, characterized in that: The method is implemented based on a pipeline siltation risk assessment method according to any one of claims 1 to 7, and comprises: Emission angle setting module: obtain the total number of sections of the target pipeline, and re-divide the target pipeline into sub-pipeline sections based on the set number of sections, and set the emission angles corresponding to each sub-pipeline section, and operate the sound wave emission device of each sub-pipeline section to emit corresponding sound waves into each sub-pipeline section according to each emission angle; Preliminary sedimentation judgment module: measures the time difference from the emission of the sound wave to the reception of the sound wave in each sub-pipeline section at each emission angle, thereby determining the preliminary sedimentation assessment coefficient corresponding to each emission angle in each sub-pipeline section; and at the same time determines whether there is sediment in each sub-pipeline section; Sedimentation risk assessment module: The sub-pipeline segment with sediment is recorded as the target pipeline segment, and then the sedimentation risk assessment coefficient of each target pipeline segment is comprehensively analyzed and fed back to the pipeline background terminal.

Images