Deep well anode air resistance identification method

Abstract

The invention discloses a deep well anode air resistance identification method. The method comprises steps S10 to S110. In the embodiment of the invention, in a period of time before the constant potential rectifier is abnormal, according to the loop resistors at the head end and the tail end of the selected characteristic part, the length is calculated by using a grounding resistance formula, and then a reference array formed by the loop resistors arranged from small to large is obtained through back-stepping; and obtaining an actual measurement array formed by the loop resistors which are directly calculated according to the time sequence, and calculating a correlation coefficient between the two arrays, so as to determine whether the operation abnormality of the potentiostat is caused by the air resistance problem or not, so that manual judgment according to experience is not needed, the air resistance problem is quickly and accurately identified, and the accuracy of the air resistance problem is improved. Therefore, whether the system abnormity is caused by the air resistance problem or not is rapidly discharged, and the system abnormity can be rapidly and accurately identified.

Description

Technical Field

[0001] The invention belongs to the technical field of cathode protection and anti-corrosion, and in particular relates to a method for identifying deep well anode gas resistance. Background Art

[0002] In smart pipelines, intelligent cathodic protection is an important part of ensuring pipeline anti-corrosion performance and operational safety. The intelligence of cathodic protection relies on the intelligent management system to correctly identify various abnormal conditions in the system and pipelines. Only by accurately identifying these abnormal conditions can targeted treatment measures be formulated to improve the safety and reliability of pipeline operation.

[0003] The anode bed (also known as auxiliary anode) in the impressed current cathodic protection system is one of the important equipment for achieving cathodic protection. The anode bed introduces the protection current from the power supply into the soil through the conductor, among which the deep well anode is one of the common auxiliary anode forms. The deep well anode adopts a structural design in which the grounding electrodes are arranged vertically at intervals, so as to ensure that there is sufficient spacing between the grounding electrodes while occupying a small area, thereby providing a wider distribution area. This design can effectively ensure that the protection current enters the soil evenly and achieve good protection for the pipeline.

[0004] However, when the current is sent into the soil by the deep well anode, the ions in the soil will produce gas due to electrolysis. When the gas covers the surface of the grounding electrode, the grounding resistance will increase. Even when the constant potential instrument increases the output voltage to the maximum, it still cannot provide the current that meets the cathodic protection requirements, thus affecting the anti-corrosion effect of the protected pipeline. In order to solve this problem, deep well anodes are usually equipped with exhaust pipes to discharge the gas generated by electrolysis to avoid gas blockage. However, due to construction problems or equipment aging, the exhaust pipe may be blocked or the exhaust may not be smooth, resulting in gas blockage. The gas blockage phenomenon will directly affect the normal operation of the constant potential instrument, making the system unable to provide sufficient protection current, which will lead to the failure of the cathodic protection of the pipeline.

[0005] In addition to the gas blockage problem, there are some factors that can cause the constant potential instrument to increase its output voltage to the maximum, causing abnormal output of the constant potential instrument. Therefore, identifying whether a gas blockage problem occurs is an important part of efficiently discovering abnormalities in the cathodic protection system. Traditional gas blockage identification relies entirely on manual experience for judgment, which has obvious deficiencies in accuracy and timeliness. Especially under the demand for intelligent cathodic protection, manual judgment is difficult to meet the requirements for rapid and accurate identification of abnormalities. Therefore, the existing technology urgently needs a method that can intelligently and timely identify deep well anode gas blockage problems, so as to improve the operating efficiency and reliability of the cathodic protection system, fully meet the needs of smart pipeline construction, and quickly and correctly identify abnormal conditions in the system. Summary of the invention

[0006] In view of the shortcomings existing in the related art, the present invention provides a deep well anode gas resistance identification method to solve the problem that the current gas resistance identification relies on manual experience judgment, resulting in difficulty in meeting abnormal fast and accurate identification requirements.

[0007] The present invention provides a deep well anode gas resistance identification method, which is characterized by comprising the following steps:

[0008] S10: Determine the current output voltage V of the constant potential instrument out Is it greater than or equal to its rated output voltage U en If the judgment result is yes, the current moment is regarded as an abnormal moment and the process goes to step S20;

[0009] S20: extracting the identification period, the end of the identification period is the abnormal moment, and dividing the identification period into a starting section, a middle section and a terminal section in order from the beginning to the end of time;

[0010] S30: Get the initial loop resistance R b_1 , V out_b and I out_b are the output voltage and output current of the potentiostat in the initial stage;

[0011] S40: Get the final loop resistance R e_1 , V out_e and I out_e are the output voltage and output current of the potentiostat in the final stage;

[0012] S50: According to the grounding resistance formula, using the initial loop resistance R b_1 , obtain the equivalent soil resistivity ρ;

[0013] S60: Get the final maximum loop resistance R e , e is the base of natural logarithm, D a is the deep well anode diameter;

[0014] S70: Obtain the initial starting loop resistance Rb, Rb = min (α × Rb _1 , Rb _1 +β×(Re-Rb _1 )), the first calculation coefficient α is greater than 1, and the second calculation coefficient β is greater than 0 and less than 1;

[0015] S80: According to the grounding resistance formula, the maximum loop resistance R at the end is used. e and initial starting circuit resistance R b , respectively obtain the minimum equivalent length L a_e and the maximum equivalent length L a_b ;

[0016] S90: L a_e and L a_b Divide it into X equal parts, and get L a_0_0 ~L a_0_X , L a_0_0 =L a_e , L a_0_X =L a_b According to the grounding resistance formula, using L a_0_0 ~L a_0_X , and get R 0_0 ~R 0_X , and get R 0_0 ~R 0_X An array of references arranged from small to large;

[0017] S100: Extract the measured time period, the end of which is the abnormal time. Divide the measured time period into X parts to obtain the measured times T0~T X , obtain the measured time T0~T X The output voltage and output current of the potentiostat are used to calculate the measured loop resistance R 1-0 ~R 1-X , and get R 1-0 ~R 1-X An array of measurements arranged in chronological order;

[0018] S110: Calculate the correlation coefficient R between the reference array and the measured array according to the Pearson correlation coefficient formula, and determine whether the correlation coefficient R is greater than or equal to a set threshold value γ. If the determination result is yes, it is determined that an air blockage problem occurs; if the determination result is no, it is determined that no air blockage problem occurs;

[0019] The grounding resistance formula is: R d is the grounding resistance, L a is the deep well anode length.

[0020] In some embodiments, the deep well anode has a plurality of coaxially arranged ground electrodes, each of which has a diameter D a , L a is the sum of the lengths of all ground electrodes.

[0021] In some embodiments, in step S20, if the continuous operation time of the constant potential instrument is greater than 1 year, the duration of the identification period is 1 year; if the continuous operation time of the constant potential instrument is less than 1 year, the beginning of the identification segment is the start-up time of the constant potential instrument.

[0022] In some embodiments, the duration of the initial phase is 7 to 15 days, and the duration of the terminal phase is 1 to 7 days.

[0023] In some embodiments, 1<α<3, 0<β<0.5.

[0024] In some embodiments, in step S80, when hour,

[0025] In some of these embodiments, X≥20.

[0026] In some embodiments, in step S100, when extracting the measured time period, searching forward from the abnormal time, when the quotient of the output voltage and output current of the constant potential instrument at a certain time is less than or equal to the initial starting loop resistance R b , then this moment is taken as the beginning of the measured period.

[0027] In some of the embodiments, in step S100, if the output voltage and output current of the constant potentiostat are not collected at a certain measuring moment, the output voltage and output current at the measuring moment are obtained by linear interpolation method according to the output voltages and output currents at the adjacent two collecting moments.

[0028] In some of these embodiments, 0.5≤γ≤1.

[0029] Compared with the prior art, the beneficial effect of the present application lies in that: in the embodiment of the present invention, in the period of time before the abnormality of the constant potential instrument occurs, the loop resistance at the first and second ends of the characteristic part are selected, and the length is calculated using the grounding resistance formula, and then the reference array composed of the loop resistance arranged from small to large is obtained by reverse calculation, and the measured array composed of the loop resistance directly calculated according to the time sorting is obtained, and the correlation coefficient between the two arrays is calculated, so as to determine whether the operation abnormality of the constant potential instrument is caused by an air resistance problem, without the need for manual judgment based on experience, and to achieve rapid and accurate identification of the air resistance problem, so as to quickly rule out whether the system abnormality is caused by the air resistance problem, so that the system abnormality can be quickly and accurately identified, solving the problem that the current air resistance identification relies on manual experience judgment, making it difficult to meet the requirements of rapid and accurate identification of abnormalities. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of this application. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0031] Figure 1 The present invention is a process flow chart of the deep well anode gas resistance identification method. DETAILED DESCRIPTION

[0032] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0033] In the description of the present invention, it should be understood that the terms "center", "lateral", "longitudinal", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as a limitation on the present invention.

[0034] The terms "first", "second", and "third" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the technical features indicated. Thus, a feature defined as "first", "second", and "third" may explicitly or implicitly include one or more of the features.

[0035] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a direct connection, or an indirect connection through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0036] like Figure 1 As shown, in an exemplary embodiment of the deep well anode gas resistance identification method of the present invention, the deep well anode gas resistance identification method includes the following steps:

[0037] S10: Determine the current output voltage V of the constant potential instrument out Is it greater than or equal to its rated output voltage U en If the result of the judgment is yes, the current moment is regarded as an abnormal moment, and the process proceeds to step S20. More specifically, the potentiostat collects the output voltage and output current during its operation at a certain time frequency. Each time an output voltage is collected, it is compared with the rated output voltage. When the current output voltage is greater than or equal to the rated output voltage, the output voltage is deemed to be out of standard, and the time when the output voltage that exceeds the standard is collected is marked as an abnormal moment, and the process proceeds to the subsequent steps.

[0038] S20: extracting the identification period, the end of the identification period is the abnormal moment, that is, the time from the abnormal moment forward is taken as the identification period, and the identification period is divided into a starting section, a middle section and a terminal section in order from the beginning to the end according to time.

[0039] S30: Get the initial loop resistance R b_1 , V out_b and I out_b is the output voltage and output current of the constant potential instrument in the initial stage. More specifically, in the initial stage, the constant potential instrument collects a set of output voltage and output current at each collection moment, and divides the output voltage and output current at that moment to obtain the loop resistance. The average value of all the loop resistances obtained in the initial stage is taken as the initial loop resistance R b_1 .

[0040] S40: Get the final loop resistance R e_1 , V out_e and I out_e is the output voltage and output current of the potentiostat in the final stage. More specifically, in the final stage, the potentiostat collects a set of output voltage and output current at each collection moment, and divides the output voltage and output current at that moment to obtain the loop resistance. The average value of all loop resistances obtained in the final stage is taken as the final loop resistance R e_1 .

[0041] S50: According to the grounding resistance formula, using the initial loop resistance R b_1 , obtain the equivalent soil resistivity ρ. More specifically, the ground resistance formula is: Set the initial loop resistance R b_1 As the grounding resistor R d Substitute into the formula, the deep well anode length L in the formula a and deep well anode diameter D a All are known deep well anode parameters, thus obtaining the equivalent soil resistivity ρ.

[0042] S60: Get the final maximum loop resistance R e , That is, the final maximum loop resistance R e Select the final loop resistance R e_1 and through The minimum value among the values ​​found, where e is the base of the natural logarithm.

[0043] S70: Obtain the initial starting loop resistance Rb, Rb = min (α × Rb _1 , Rb _1 +β×(Re-Rb _1)), that is, the initial starting circuit resistance Rb is selected as α×Rb _1 The value obtained and Rb _1 +β×(Re-Rb _1 ) is the minimum value among the numerical values ​​obtained, the first calculation coefficient α is greater than 1, and the first calculation coefficient β is greater than 0 and less than 1.

[0044] S80: According to the grounding resistance formula, the maximum loop resistance R at the end is used. e and initial starting circuit resistance R b , respectively obtain the minimum equivalent length L a_e and the maximum equivalent length L a_b More specifically, the deep well anode diameter D in the grounding resistance formula a and the equivalent soil resistivity ρ are known, the final maximum loop resistance R e As the grounding resistor R d Substitute into the formula and we get the minimum equivalent length L a_e In the same way, the initial starting circuit resistance R b As the grounding resistor R d Substitute into the formula and we get the maximum equivalent length L a_b .

[0045] S90: L a_e and L a_b Divide it into X equal parts, and get L a_0_0 ~L a_0_X , L a_0_0 =L a_e , L a_0_X =L a_b According to the grounding resistance formula, using L a_0_0 ~L a_0_X , and get R 0_0 ~R 0_X , obtained by L a_0_0 ~L a_0_X and R 0_0 ~R 0_X Composition and according to R 0_0 ~R 0_X A reference array arranged from small to large. More specifically, L a_0_0 , L a_0_1 , L a_0_2 ...L a_0_X Arrange from small to large, and the difference between adjacent lengths is L a_e With L a_b The difference is divided by X. The diameter of the deep well anode D in the grounding resistance formula a and the equivalent soil resistivity ρ are known, L a_0_0 , L a_0_1 , L a_0_2 ...La_0_X As the deep well anode length L a Substitute the grounding resistance formula into the grounding resistance, and the grounding resistance is R 0_0 , R 0_1 , R 0_2 ……R 0_X . R 0_0 ~R 0_X Arrange them in ascending order and use them as a reference array.

[0046] S100: extract the measured time period, the end of which is the abnormal time, that is, the time from the abnormal time forward is taken as the measured time period, and the measured time period is divided into X parts to obtain the measured times T0, T1, T2...T X , that is, the time of all the separation points of the measured time period and the time of the first two ends, to obtain the measured time T0~T X The output voltage and output current of the constant potential instrument are obtained by dividing the output voltage at each measured moment by the output current, that is, the loop resistance at each measured moment is obtained, and thus the measured loop resistance R is calculated. 1-0 ~R 1-X . The measured loop resistance R 1-0 ~R 1-X According to the corresponding measured time T0~T X Arrange them in order from earliest to latest time and use them as the measured array.

[0047] S110: According to the Pearson correlation coefficient formula, calculate the correlation coefficient R between the reference array and the measured array, that is, the R 0_0 ~R 0_X The reference array and the R array which will be arranged from the beginning to the end in time 1-0 ~R 1-X The measured array is brought into the Pearson correlation coefficient formula for calculation to obtain the correlation coefficient R. It is determined whether the correlation coefficient R is greater than or equal to the set threshold γ. If the determination result is yes, it is determined that an air lock problem occurs; if the determination result is no, it is determined that no air lock problem occurs.

[0048] In step S10, the output voltage exceeding the standard of the constant potential instrument is used as a condition for starting the subsequent steps, so as to ensure that the air blockage problem is identified when the operation state of the constant potential instrument is abnormal.

[0049] In steps S30 and S40, the average loop resistance of the head section and the tail end is selected as the initial loop resistance R in a period of time before the constant potential instrument becomes abnormal. b_1 and the final loop resistance R e_1 , eliminate the fluctuation influence of the output voltage and output current of the constant potential instrument, and ensure the accuracy of the acquired data.

[0050] Since soil resistivity changes with time, the soil resistivity is obtained by obtaining the loop resistance in the initial stage before the constant potential instrument in step S50 as the equivalent soil resistivity ρ, thereby ensuring the accuracy of the resistance and length calculated using the ground resistance formula.

[0051] Due to the relationship curve between the grounding resistance and the electrode length reflected in the grounding resistance formula, the grounding resistance will first increase as the electrode length decreases, and then begin to decrease after reaching an extreme value. In step S60, the smaller one of the final loop resistance and the loop resistance extreme value is selected, that is, the final maximum loop resistance R e The selected position is reasonably forward, so as to ensure the maximum loop resistance R obtained at the end e The characteristic part of the relationship curve where the electrode length and the grounding resistance are positively correlated. In step S70, the initial loop resistance R is calculated using the first calculation coefficient α. b_1 amplify, or use the initial loop resistance R b_1 and the final loop resistance R e_1 The difference between the initial loop resistance R and the second calculation coefficient β b_1 Amplify the initial starting loop resistance R b The selection position is reasonably far back, so as to ensure that the initial starting circuit resistance R b The characteristic part of the relationship curve where the electrode length and ground resistance are positively correlated.

[0052] In step S80, the initial starting loop resistance R of the first and last ends of the characteristic part in the relationship curve is characterized. b and the final maximum loop resistance R e , the electrode length variation range that best reflects the characteristics of the relationship curve is inferred from the ground resistance formula, that is, the minimum equivalent length L a_e To the maximum equivalent length L a_b . Further, the minimum equivalent length L a_e To the maximum equivalent length L a_b Divide it into multiple parts, calculate the loop resistance at the dividing point and both ends, and arrange the loop resistance from small to large to obtain a reference array.

[0053] When the gas blockage problem occurs, as the gas generated by electrolysis increases, the surface area of ​​the electrode covered by the gas gradually increases, the effective length of the electrode gradually decreases, and the grounding resistance of the electrode gradually increases, which in turn causes the loop resistance to gradually increase. Therefore, when the gas blockage problem occurs, the increase in grounding resistance is positively correlated with the decrease in the effective length of the electrode, that is, the increase in loop resistance is positively correlated with the decrease in the effective length of the electrode.

[0054] In step S100, the loop resistance obtained by dividing the output voltage and output current directly from the output voltage and output current at a certain time before the abnormality occurs is arranged in chronological order as a measured array. In step S110, the reference array and the measured array are subjected to correlation coefficient calculation. When the correlation coefficient exceeds the threshold value, it is proved that there is a high correlation between the loop resistance calculated by directly taking the quotient of the output voltage and the output current according to the time sequence before the abnormality occurs and the loop resistance calculated by the grounding resistance formula according to the soil resistivity, the diameter and length of the deep well anode, that is, before the abnormality occurs, the increase of the loop resistance is positively correlated with the decrease of the electrode length, that is, the increase of the loop resistance is caused by the decrease of the effective electrode length due to the air resistance problem, and thus it is determined that the air resistance problem has occurred.

[0055] The output voltage of the constant potentiostat exceeded the standard, and it was determined that the cathodic protection system was abnormal. Then it was judged whether it was a gas blockage problem, thus effectively ruling out whether the system abnormality was caused by a gas blockage problem. The real cause of the abnormality of the cathodic protection system was discovered efficiently and quickly, and targeted treatment measures could be formulated as soon as possible to improve the safety and reliability of pipeline operation.

[0056] In the above-mentioned exemplary embodiment, the deep well anode gas resistance identification method selects a time before the abnormality occurs, selects the two ends of the most characteristic section of the loop resistance, and calculates the effective length of the electrode using the grounding resistance formula, and then reversely infers the reference array composed of loop resistances arranged from small to large based on the obtained effective length of the electrode, and calculates the correlation coefficient with the measured array composed of loop resistances directly calculated at one end of the time directly selected in the same time period, so as to determine whether the increase in loop resistance is positively correlated with the decrease in effective length, and then determine whether such change is caused by gas resistance problems, thereby determining whether the operation abnormality of the constant potential instrument is caused by gas resistance problems, without the need for manual judgment based on experience, to achieve rapid and accurate identification of gas resistance problems, thereby quickly eliminating whether the system abnormality is caused by gas resistance problems, so that the system abnormality can be quickly and accurately identified, solving the problem that the current gas resistance identification relies on manual experience judgment and is difficult to meet the requirements of rapid and accurate identification of abnormalities.

[0057] In some embodiments, the deep well anode has a plurality of coaxially arranged ground electrodes, each of which has a diameter D a , L a is the sum of the lengths of each grounding electrode, thus ensuring that when the grounding resistance formula is used for calculation, the deep well anode diameter D a And deep well anode length is accurate.

[0058] In some embodiments, in step S20, if the continuous operation time of the constant potential instrument is greater than 1 year, the duration of the identification period is 1 year; if the continuous operation time of the constant potential instrument is less than 1 year, the beginning of the identification segment is the start-up time of the constant potential instrument.

[0059] Since soil resistivity may change over time, the identification period is selected based on the length of time the potentiostat has been running, and a period of one year before the abnormality occurs is selected to ensure that the equivalent soil resistivity ρ calculated from the data in the selected time period can accurately reflect the actual resistivity of the current soil, thereby ensuring the accuracy of the electrode length and resistance calculated subsequently using the ground resistance formula.

[0060] In some embodiments, the duration of the initial phase is 7 to 15 days, and the duration of the final phase is 1 to 7 days. b_1 When obtaining, the average value of the loop resistance over a long period of time is used to ensure that the initial loop resistance R b_1 The final loop resistance R is the loop resistance feedback during normal operation of the potentiostat. e_1 When obtaining, the average value of the loop resistance in a short period of time is used to ensure that the final loop resistance R e_1 It provides loop resistance feedback when the potentiostat approaches abnormal operation, thereby eliminating the loop resistance fluctuation caused by output current and voltage fluctuation at both ends to the greatest extent, and improving the accuracy of data acquisition.

[0061] In some embodiments, 1<α<3, 0<β<0.5, so as to obtain the initial starting loop resistance R b When selecting the loop resistance, avoid selecting the position too far back to prevent the data feature segment from being shorter and the inability to fully evaluate the correlation, thereby ensuring the accuracy of the correlation evaluation.

[0062] In some embodiments, between ρ and D a Under the condition of being determined unchanged, the grounding resistance formula can reflect the relationship between the grounding resistance and the effective length of the deep well anode. However, the formula for obtaining the grounding resistance cannot directly inverse the function to obtain the functional relationship of the deep well anode length from the known grounding resistance. Therefore, in step S80, when When the resistance is at the extreme value of the relationship curve presented by the grounding resistance formula, Therefore, the anode length can be directly obtained by using the extreme value of the function, which improves the convenience of data acquisition. It is necessary to obtain the anode length through numerical methods such as binary search.

[0063] In some embodiments, in steps S90 and S100, the time period is divided into the same X parts, and X≥20, so as to ensure that the number of selected length points and time points is sufficient, so that the obtained reference array and measured array can accurately reflect the relationship between the loop resistance and the change of the effective length of the deep well anode.

[0064] In some embodiments, in step S100, when extracting the measured time period, searching forward from the abnormal time, when the quotient of the output voltage and output current of the constant potential instrument at a certain time is less than or equal to the initial starting loop resistance R b The acquisition of the actual measurement period is to search for the loop resistance less than or equal to the initial starting loop resistance R b The interference of loop resistance fluctuation on the selection of the first end of the actual measurement period is eliminated, ensuring that the actual measurement period is consistent with R 0_0 ~R 0_X The correlation of the effective length range of the deep well anode is shown.

[0065] In some embodiments, in step S100, if the output voltage and output current of the constant potentiostat are not collected at a certain measuring moment, the output voltage and output current at the measuring moment are obtained by linear interpolation method according to the output voltages and output currents at the adjacent two collecting moments.

[0066] Since the output voltage and output current of the constant potentiostat are not collected continuously, when the output voltage and output current are not collected at a certain moment in the actual measurement, the data of the adjacent collection points before and after are used to fit the output voltage and output current, ensuring the feasibility and representativeness of the actual loop resistance calculation.

[0067] In some embodiments, when 0.5≤γ≤1, that is, when the correlation between the reference array and the measured array is not less than 50%, it is determined that an air blockage problem has occurred, ensuring that the correlation can accurately identify whether an air blockage problem has occurred.

[0068] Finally, it should be noted that: the various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0069] The above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, ordinary technicians in the field should understand that the specific implementation methods of the present invention can still be modified or some technical features can be replaced by equivalents without departing from the spirit of the technical solution of the present invention, which should be included in the scope of the technical solution for protection of the present invention.

Claims

1. A deep well anode gas resistance identification method, characterized in that: The steps include: S10: Determine the current output voltage V of the constant potential instrument out Is it greater than or equal to its rated output voltage U en If the judgment result is yes, the current moment is regarded as an abnormal moment and the process goes to step S20; S20: extracting the identification period, the end of the identification period is the abnormal moment, and dividing the identification period into a starting section, a middle section and a terminal section in order from the beginning to the end of time; S30: Get the initial loop resistance R b_1 , V out_b and I out_b are the output voltage and output current of the potentiostat in the initial stage; S40: Get the final loop resistance R e_1 , V out_e and I out_e are the output voltage and output current of the potentiostat in the final stage; S50: According to the grounding resistance formula, using the initial loop resistance R b_1 , obtain the equivalent soil resistivity ρ; S60: Get the final maximum loop resistance R e , e is the base of natural logarithm, D a is the deep well anode diameter; S70: Obtain the initial starting loop resistance Rb, Rb = min (α × Rb _1 , Rb _1 +β×(Re-Rb _1 )), the first calculation coefficient α is greater than 1, and the second calculation coefficient β is greater than 0 and less than 1; S80: According to the grounding resistance formula, the maximum loop resistance R at the end is used. e and initial starting circuit resistance R b , respectively obtain the minimum equivalent length L a_e and the maximum equivalent length L a_b ; S90: L a_e and L a_b Divide it into X equal parts, and get L a_0_0 ~L a_0_X , L a_0_0 =L a_e , L a_0_X =L a_b According to the grounding resistance formula, using L a_0_0 ~L a_0_X , and get R 0_0 ~R 0_X , and get R 0_0 ~R 0_X An array of references arranged from small to large; S100: Extract the measured time period, the end of which is the abnormal time, and divide the measured time period into X parts to obtain the measured times T0~T X , obtain the measured time T0~T X The output voltage and output current of the potentiostat are used to calculate the measured loop resistance R 1-0 ~R 1-X , and get R 1-0 ~R 1-X An array of measurements arranged in chronological order; S110: Calculate the correlation coefficient R between the reference array and the measured array according to the Pearson correlation coefficient formula, and determine whether the correlation coefficient R is greater than or equal to a set threshold value γ. If the determination result is yes, it is determined that an air blockage problem occurs; if the determination result is no, it is determined that no air blockage problem occurs; The grounding resistance formula is: R d is the grounding resistance, L a is the deep well anode length.

2. The deep well anode gas resistance identification method according to claim 1, characterized in that: The deep well anode has multiple coaxially arranged grounding electrodes, each with a diameter of D a , L a is the sum of the lengths of all ground electrodes.

3. The deep well anode gas resistance identification method according to claim 1, characterized in that: In step S20, if the continuous operation time of the constant potential instrument is greater than 1 year, the duration of the identification period is 1 year; if the continuous operation time of the constant potential instrument is less than 1 year, the beginning of the identification section is the start time of the constant potential instrument.

4. The deep well anode gas resistance identification method according to claim 1, characterized in that: The duration of the initial stage is 7 to 15 days, and the duration of the terminal stage is 1 to 7 days.

5. The deep well anode gas resistance identification method according to claim 1, characterized in that: 1＜α＜3，0＜β＜0.5。 6. The deep well anode gas resistance identification method according to claim 1, characterized in that: In step S80, when hour, 7. The deep well anode gas resistance identification method according to claim 1, characterized in that: X≥20。 8. The deep well anode gas resistance identification method according to claim 1, characterized in that: In step S100, when extracting the measured time period, search forward from the abnormal time, when the quotient of the output voltage and output current of the constant potential instrument at a certain time is less than or equal to the initial starting loop resistance R b , then this moment is taken as the beginning of the measured period.

9. The deep well anode gas resistance identification method according to claim 1, characterized in that: In step S100, if the output voltage and output current of the constant potential instrument are not collected at a certain measuring moment, the output voltage and output current at the measuring moment are obtained by linear interpolation method according to the output voltages and output currents at the adjacent two collecting moments.

10. The deep well anode gas resistance identification method according to claim 1, characterized in that: 0.5≤γ≤1.

Images