CN108020713A - A kind of non-contact measurement apparatus and method of buried pipeline geomagnetic induction current GIC - Google Patents

Abstract

The present invention relates to a kind of non-contact measurement apparatus and method of buried pipeline geomagnetic induction current GIC, the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC includes compound flux gate probe, power module, Signal-regulated kinase, microprocessor, memory module, local communication module and remote communication module；The characteristics of for spiral-selded pipe road geomagnetic induction current GIC, the magnetic induction intensity produced using measuring device to spiral line current in spiral weld pipeline is acquired and the calculating of electric current, to realize the accurate measurement of spiral-selded pipe road geomagnetic induction current GIC.The shortcomings that present invention is for present in prior art, give the current model and GIC measuring methods of spiral weld pipeline, this method is according to the manufacture craft of spiral weld pipeline, the electric current flowed in pipeline is modeled with the track of helix, improve the accuracy of model, make the GIC that measurement is drawn and actual current value closer, more really reflect the GIC sizes in pipeline.

Description

A kind of non-contact measurement apparatus and method of buried pipeline geomagnetic induction current GIC

Technical field

The present invention relates to buried pipeline Defense natural calamities technical field, and in particular to a kind of buried pipeline earth induction electricity Flow the non-contact measurement apparatus and method of GIC.

Background technology

From the law of electromagnetic induction, when magnetic storm occurs, in point that buried metal pipeline and its coating are formed with the earth Sensing electric current, i.e. geomagnetic induction current (GIC) are produced in cloth closure conductor circuit.Since GIC is harmful to pipeline, interference pipe Current potential (PSP) measures, and makes cathodic protection (CP) device out of control or even has the opposite effect, causes corrosive pipeline to aggravate, so as to reduce pipe Road service life etc., so, after finding GIC in the duct from 1873, pipeline is influenced GIC and protection question is always by pass Note.Up to the present, almost in the pipeline of each country it has been found that cross GIC, in the Space Weather Program of various countries, equal handles GIC influences pipeline to be classified as important research content.

Currently more for detection buried pipeline and buried depth measuring study both at home and abroad, in the market also has more detection to set It is standby.But it is less for the measuring apparatus of geomagnetic induction current (GIC), following method is currently mainly used：

(1) pipeline is established into DSTL models, Geomagnetic Observational Data and pipeline relevant parameter near input channel, final meter Calculate pipe to soil potential and pipe current.Chinese invention patent application 201510579331.X propose a kind of buried oil-gas pipeline by The GIC and the computational methods of PSP that geomagnetic storm influences, it is proposed that a kind of new pipeline model LZS-DSTL, by the ground of Near Pipelines Magnetic storm data, relevant environmental parameter and pipe parameter input, and then calculate GIC and PSP.The sensing electric current that this method calculates Accuracy dependent on its mode input parameter.

(2) standard pipe to soil potential (PSP) test, close spacing potential test technology (CIPS), DC potential gradient method (DCVG), Pearson came (Person) method, frequency conversion frequency selection method, DC current-potentiometry etc..Chinese invention patent 200810119164.0 propose a kind of measuring device and method of buried pipeline geomagnetic induction current and pipe to soil potential, in pipeline On test pile on install GPS synchronous recording devices, by the pipe of each measurement point transient potential recorded and calculate ground Table electric field strength and pipe parameter., must be under cutting status by Current Mutual Inductance if wishing to install measurement point in pipeline optional position Device entangles pipeline, and current transformer otherwise can only be installed at the pipeline jumper at station end.

(3) the sensing electric current of cordless detection buried pipeline, i.e. magnetic survey are used under not cutting status.Britain The Detector for Stray Currents SCM of Lei Di companies production, mainly for detection of the hot spot of pipeline, flowing in and out for electric current is a little fixed Position, can only be roughly calculated out pipe current, precision is very poor.Chinese invention patent 201110097516.9 proposes a kind of underground pipe Road buried depth and the measuring method of stray electrical current size direction and magnetic azinuth, by linearly arrange in ground level 4 A fluxgate sensor, monitors the magnetic field intensity above buried pipeline in real time, and calculate the buried depth of pipeline, current direction and Size.This method has that measuring method is simple, extra signal source, real-time measuring data is not required and can arbitrarily install monitoring The advantages that point, but the precision of its measurement data is relatively low, and error is larger.

In conclusion the ideal method of detection buried pipeline sensing electric current GIC should be contactless magnetic method and survey Amount.But bigger to the measurement error of buried pipeline GIC using the device of non-contact measuring technology both at home and abroad at present, its reason It is that the sensing electric current of buried pipeline is equivalent to an endless constant current for being located at the pipeline center of circle by it, passes through Biot-Sa Cut down your law and draw relation between electric current and earth's surface magnetic field intensity, to obtain the unknown quantitys such as electric current.This model establishes bar Part is：The distribution of resistance of metallic conduit is uniform, in pipeline everywhere GIC along duct orientation streamlined flow.And actual metal pipeline Resistance and non-uniform Distribution, while influenced be subject to pipeline-weld, cause the flow direction of its interior electric current not along duct orientation straight line Flowing.

But the manufacture craft of pipeline has many kinds at present, is broadly divided into seamless steel pipe, steel pipes with straight and spiral weld steel Pipe, as shown in Figure 1.The pipeline material distribution of resistance of seamless steel pipe and steel pipes with straight is relatively uniform, and this kind of pipeline GIC can use nothing Limit for length's constant current model equivalent conduit electric current.And the weld seam of spiral welded pipe is the helix wound along pipeline, due to weld material Inconsistent with pipeline material, the direction of current flow for causing such pipeline is no longer straight line, but along spiral pipeline flow forward.

The buried long oil and gas pipeline in China largely uses spiral weld pipeline at present, uses the prior art and device measuring The sensing electric current of spiral weld pipeline will necessarily produce large error, can not accurately measure pipeline GIC.

Therefore, during the buried pipeline GIC produced by magnetic survey by different modes, electricity should be calculated using different models Flow valuve.

The content of the invention

In order to solve the above technical problem, the present invention provides a kind of the non-contact of buried pipeline geomagnetic induction current GIC Formula measuring device and method.For the problem of spiral-selded pipe road geomagnetic induction current GIC errors are big is measured in the prior art, originally Invention provides a kind of non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC, and gives spiral weld pipeline Current model and GIC measuring methods, the measuring method is according to the manufacture craft of spiral weld pipeline, the electric current that will be flowed in pipeline It is modeled with the track of helix, improves the accuracy of model, makes the GIC that measurement is drawn and actual current value closer, More really reflect the GIC sizes in pipeline.

To achieve these goals, the present invention adopts the following technical scheme that.

First, the invention discloses a kind of non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC, for adopting The magnetic induction intensity data that spiral line current produces in collection spiral weld pipeline include with the calculating for realizing electric current, the measuring device Compound flux gate probe, power module, Signal-regulated kinase, microprocessor, memory module, local communication module and telecommunication Module；

The Signal-regulated kinase, microprocessor, memory module, local communication module, remote communication module form the survey Measure the signal processing module of device；

The compound flux gate probe connects Signal-regulated kinase, and the Signal-regulated kinase connects microprocessor；It is described Signal-regulated kinase is used to be nursed one's health the signal that compound flux gate probe gathers, and is converted to digital signal and is sent to micro- place Manage device；

The microprocessor connects memory module, local communication module, remote communication module respectively；The microprocessor from The digital signal of magnetic induction intensity is obtained in Signal-regulated kinase, and is further handled；The memory module is used to store micro- place The magnetic induction intensity data that reason device is handled well；The local communication module be used for local peripheral device communication, by the number of storage According to being sent to local external equipment；The remote communication module is used to communicate with remote server, and the data of storage are sent to Remote server is further handled；

The power module connects compound flux gate probe, Signal-regulated kinase, microprocessor, local communication mould respectively Block, remote communication module, memory module, stabilized power source is provided for whole measuring device.

Preferably, the compound flux gate probe includes at least one bar shaped bar, at least two simple component fluxgate senses Device and at least two level meter.

In any of the above-described technical solution preferably, it is described by a bar shaped bar, 2 simple component fluxgate sensors and The compound flux gate probe of 2 level instrument devices composition, a length of 2 (L+d3) of bar shaped bar, it is wide and it is high be D, central point O Point, it is each at distance O points both sides L to fix simple component fluxgate sensor M1, a M2；Simple component fluxgate sensor M1, M2's Length and width is respectively d1, d2；The central point of simple component fluxgate sensor M1 is d3 apart from the left hand edge of bar shaped bar, the center of M2 Point is d3 apart from the right hand edge of bar shaped bar, its measurement direction is x coordinate direction；Surface fixes two at bar shaped bar left hand edge E A horizon rule S1, S2, according to xyz coordinate directions, horizon rule S1 is located at identical yoz with simple component fluxgate sensor M1, M2 Face, horizon rule S2 are located at proximal surface xoz faces, and horizon rule length is D, takes L=b pi/2s, wherein 2 π b are defined as the spiral shell of helix current line Away from.

In any of the above-described technical solution preferably, the compound flux gate probe connects Signal-regulated kinase, described Signal-regulated kinase connect microprocessor, the Signal-regulated kinase by the analog signal of 2 simple component fluxgate sensors into Row conditioning, and be converted to digital signal and be sent to microprocessor.

The invention also discloses a kind of contactless measurement of buried pipeline geomagnetic induction current GIC.For spiral The characteristics of welded tube road geomagnetic induction current GIC, using the non-of as above any one of them buried pipeline geomagnetic induction current GIC Contact measuring apparatus, the magnetic induction intensity produced to spiral line current in spiral weld pipeline is acquired and the meter of electric current Calculate, to realize the accurate measurement of spiral-selded pipe road geomagnetic induction current GIC, which includes：

(1) mathematical model of spiral-selded pipe road geomagnetic induction current GIC is established；

(2) orientation of buried spiral weld pipeline is obtained；

(3) influence of the interference field to GIC computational solution precisions is filtered out；

(4) GIC is calculated using numerical integrating.

In any of the above-described technical solution preferably, in (1) described above item, for spiral-selded pipe road earth magnetism The characteristics of sensing electric current GIC, establishes the mathematical model of spiral-selded pipe road geomagnetic induction current GIC, according to spiral weld pipeline The mathematical model of geomagnetic induction current GIC calculates electric current I in spiral weld pipeline_GICWith a point P (x above pipeline₀,y₀,z₀) Tri- axis magnetic induction density B of x, y, z_x,B_y,B_zRelation, electric current I is calculated_GIC, spiral-selded pipe road geomagnetic induction current GIC The modeling process of mathematical model specifically include：

Electric current in spiral weld pipeline is threadingly advanced around pipe circumference, it is assumed that pipeline is endless, then its single electricity Streamline equation meets that helix parametric equation is

Wherein a is the radius of spiral weld pipeline,It is vectorial in xoy planes for what is a little formed on helix with o points Projection and the angle of x-axis, 2 π b are the pitch of helix, and z directions are moved towards for spiral weld pipeline, and x directions are perpendicular to ground level Upwardly direction, y directions are moved towards perpendicular to spiral weld pipeline and parallel to ground levels；

Assuming that the spiral line current GIC flowed in spiral weld pipeline is the spiral line current of single endless, its electric current I_GICAdvance along the z-axis direction, then first according to Biot -- Sa cuts down your theorem and lists a current elements on helixIn pipeline Point P (the x of top one₀,y₀,z₀) caused by magnetic induction intensity be

Wherein μ₀For the magnetic conductivity in vacuum,For current elementsThe vector of P points is directed toward, R is P points to current elementsMinimum distance；

IfRespectively x, y, the unit vector of z-axis, then point P and current elementsRespectively can table with vector form It is shown as

Point P is relative to current elementsPosition vector be represented by

AndComponent form be represented by

Then

WillFormula (2) is substituted into, current elements is obtained and divides what point P was produced along x, the magnetic induction intensity in tri- directions of y, z Amount is respectively

Wherein

Magnetic induction intensity of the above current elements in all directions integrates endless helix and obtains whole spiral shell Weld seam pipe current is revolved in induction level caused by P points, i.e.,

Wherein L represents the spiral line current of whole spiral weld pipeline, all directions that whole spiral line current is produced in point P Magnetic induction intensity value is respectively

In any of the above-described technical solution preferably, in (2) described above item, the tool of spiral weld pipeline is determined The GIC optimum measuring points of body position, trend and spiral weld pipeline, specifically comprise the following steps：

Step 1, the definite spiral weld pipeline in scene is substantially moved towards as z ' directions, by buried pipeline geomagnetic induction current GIC Non-contact measurement yoz faces as mobile face, make horizon rule S2 above, and adjust the axial direction of S2 and bar shaped bar, protect The measurement direction for demonstrate,proving fluxgate sensor M1 is z ' directions；

Step 2, a path M is selected, its direction is moved simple component fluxgate sensor M1 along M perpendicular to z ' It is dynamic, ensure that movement track is fully passed over directly over buried spiral weld pipeline, and record Component of magnetic flux door sensor M1 in real time Surveyed B_z′Curve；According to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′Curve, judges whether path M passes through Symmetric points (x₀, 0,2nb π) or (x₀, 0, (2n+1) b π) and record corresponding symmetric points position；

Step 3, path M ' parallel and distance l on the right side of it with selected path in above-mentioned steps 2 is selected, by single point Amount fluxgate sensor M1 is moved along M ', ensures that movement track is fully passed over directly over buried spiral weld pipeline, and in fact When record Component of magnetic flux door sensor M1 survey B_z′Curve；According to buried spiral weld pipeline B_z′Curve peak valley positioning mode point Analyse B_z′Curve, it is first determined the direction of GIC in spiral weld pipeline, then judge whether path M ' passes through symmetric points (x₀,0,2nb π) or (x₀, 0, (2n+1) b π) and record corresponding symmetric points position；

Step 4, path M " parallel and distance l on the right side of it with 3 path of above-mentioned steps is selected, by simple component fluxgate Sensor M1 is moved along M ", ensures that movement track is fully passed over directly over buried spiral weld pipeline, and real-time record Component of magnetic flux door sensor M1 surveys B_z′Curve；According to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′Curve, First determine whether the direction of GIC in spiral weld pipeline overturns, if upset, rebound step 2 remeasures, if not overturning, Then judge whether path M " passes through symmetric points (x₀, 0,2nb π) or point (x₀, 0, (2n+1) b π) and record corresponding symmetric points institute in place Put；

Step 5,4 are repeated the above steps, according to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′Curve, Until finding by symmetric points (x₀, 0,2nb π) and point (x₀, 0, (2n+1) b π) neighbouring two paths, stop repeat step 4, and determine two points substantially position P₁, P₂, 2 lines are new z ' directions；

Step 6, the distance between two paths l is reduced, 2 is repeated the above steps to step 5, concentrates analysis site P₁, P₂Near The B of path_z′Curve, until finding and B_z′Point (the x that the wave crest or trough of curve are completely superposed₀, 0,2nb π) and point (x₀,0,(2n + 1) b π) untill, this 2 line, that is, pipelines truly move towards, i.e. z directions；It can use n=0 at this time, then point (x₀, 0,0) surveyed to be optimal Amount point.

In any of the above-described technical solution preferably, (3) described above item filters out interference field to GIC result of calculations essence The influence of degree, after described (2) item determines orientation and the optimum measuring point of spiral weld pipeline, by buried pipeline earth induction The central point O of the non-contact measurement apparatus of electric current GIC is placed in optimum measuring point (x₀, 0,0) and place, wherein x₀For measurement point to pipe The sum of the distance in road and pipe radius, are to obtain at the scene, it is believed that it is known constant, using buried pipeline geomagnetic induction current The non-contact measurement apparatus of GIC, 2 simple component fluxgate sensors are placed in above buried spiral weld pipeline, adjust bar The z-axis direction of shape bar, makes it truly move towards to overlap with pipeline, adjusts horizon rule S2, the measurement direction for making fluxgate sensor is x Direction of principal axis, fixed compound flux gate probe, then pop one's head in upper two fluxgate sensor points P_mLAnd P_mRCoordinate value be respectivelyWithThe measurement data B of two sensors of measurement in real time_xmL,B_xmR, and measurement data is substituted into public affairs Formula (56), i.e.,Interference field data is filtered out by microprocessor, calculates point P_mLSecondary field data B_xL, Preserved in real time using memory module, and remote server is uploaded to by remote communication module.

It is preferably in any of the above-described technical solution, (4) described above item calculates GIC using numerical integrating, Specifically include：

Step 1, pitch constant b, the space permeability μ of pipe radius a, helix current line are determined₀And point P_mLCoordinate ValueDeng field pipes data, and as known parameters；

Step 2, the non-contact measurement dress of live buried pipeline geomagnetic induction current GIC remote server received Put the magnetic induction intensity data B sent_xLAfterwards, the parameter and B step 1 determined_xLFormula (57-3) is substituted into, i.e.,Final spiral-selded pipe road geomagnetic induction current GIC is drawn using numerical integration algorithm.

Compared with prior art, above-mentioned technical proposal of the invention has the characteristics that following several：

1st, for the production technology characteristic of spiral weld pipeline, the stream of the geomagnetic induction current GIC in such pipeline is analyzed Dynamic direction is not straight line, but is threadingly advanced.This physical model with pipe current to be equivalent to unlimited long straight conductor at present Grave fault, the GIC resultant errors that spiral weld pipeline is measured by the prior art are very big.

2nd, the mathematical model of the geomagnetic induction current GIC of spiral weld pipeline is established, and is derived on GIC electric currents and pipeline Side surveys clear and definite mathematical relationship between magnetic data.

3rd, the physical model of the geomagnetic induction current GIC of spiral weld pipeline and endless constant current model difference are very big, The method of traditional elimination background magnetic field is not suitable for spiral weld pipeline.The contactless survey that technical scheme proposes Amount method analyzes the symmetry that helix current line produces magnetic field above pipeline, and the magnetic induction by measuring symmetric points is strong Degree, can obtain geomagnetic noise magnetic field, be eliminated in calculating current, substantially reduce the influence of ambient field, and improve calculating Accuracy.

4th, the non-contact measurement apparatus that technical scheme proposes, uses 2 single shaft fluxgate sensors, composition 2 sensors are placed in probe both sides by 1 bar shaped probe, and center probe is placed at the symmetric points of spiral line current, make spy Head is axial axially in parallel with pipeline, is totally different from most of at present that probe vertical is axially placed in pipeline.

5th, the pipeline location and trend obtained at the scene by pipeline detection device is not accurate enough, this also causes magnetic to sense Measuring result error caused by device riding position is unreasonable.Technical scheme proposes one kind and is directed to spiral-selded pipe The measuring device of GIC features in road, can pass through " buried spiral weld pipeline B proposed by the present invention_z′Curve peak valley positioning mode " essence Really obtain the position of pipeline and true trend.

Brief description of the drawings

In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing There is attached drawing needed in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, may be used also To obtain other attached drawings according to these attached drawings.

Fig. 1 is the different welding procedure schematic diagrames of steel pipe in the prior art；

Fig. 2 is a preferred reality of the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC according to the invention Apply the functional block diagram of example；

Fig. 3 is a preferred reality of the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC according to the invention Apply the compound flux gate probe structure diagram of example；

Fig. 4 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the instrumentation plan that the acquisition pipeline of example accurately moves towards；

Fig. 5 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the final disposing way schematic diagram of measuring device in the measurement process of example；

Fig. 6 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the single helix schematic diagram of endless of example；

Fig. 7 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention The current elements for applying example produce the schematic diagram in magnetic field in P points；

Fig. 8 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the schematic diagram for producing magnetic field in P points with symmetrical 2 current elements of x-axis of example；

Fig. 9 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the schematic diagram for producing magnetic field in P points with symmetrical 2 current elements of example；

Figure 10 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the schematic diagram for producing magnetic field in P1 and P2 points with symmetrical 2 current elements of x-axis of example；

Figure 11 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention The x-axis magnetic induction intensity of example is applied directly over pipeline along the schematic diagram of z-axis change curve；

Figure 12 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the z directly over pipeline of example₀Two side draw measurement point P of=2nb π_L,P_RSchematic diagram；

Figure 13 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention The scene for applying example determines that pipeline substantially moves towards and actually moves towards schematic diagram；

Figure 14 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention That applies example measures path schematic diagram in spiral line current diverse location；

Figure 15 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention The schematic diagram of 9 measurement path magnetic induction intensity curves when applying λ=10 ° of example；

Figure 16 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention The schematic diagram of 9 measurement path magnetic induction intensity curves when applying λ=20 ° of example；

Figure 17 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention The schematic diagram of 9 measurement path magnetic induction intensity curves when applying λ=30 ° of example；

Figure 18 is a preferred reality of the contactless measurement of buried pipeline geomagnetic induction current GIC according to the invention Apply the schematic diagram of the optimum measurement mode of the spiral weld pipeline GIC of example；

Embodiment

Below in conjunction with the attached drawing in the embodiment of the present invention, the technical solution in the embodiment of the present invention is carried out clear, complete Site preparation describes, it is clear that described embodiment is only part of the embodiment of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, those of ordinary skill in the art are obtained every other without making creative work Embodiment, belongs to the scope of protection of the invention.

For the problem of spiral-selded pipe road geomagnetic induction current GIC errors are big is measured in the prior art, the present invention is implemented Example provides a kind of measuring device and method of spiral-selded pipe road geomagnetic induction current GIC, for present in prior art Shortcoming, gives the current model and GIC measuring methods of spiral weld pipeline, and this method is according to the making work of spiral weld pipeline Skill, the electric current flowed in pipeline is modeled with the track of helix, the accuracy of model is improved, makes what measurement was drawn GIC is closer with actual current value, more really reflects the GIC sizes in pipeline.

When magnetic storm occurs, earth's magnetic field is interfered, and is then formed in buried metal pipeline and its coating with the earth Geomagnetic induction current (GIC) is produced in distribution closure conductor circuit.The GIC flowed in spiral weld pipeline receives weld seam Influence, form the current path of zigzag shape, and non-directional route.The parametric equation of its current model formula as described below (1) shown in.As shown in Figure 6, coordinate origin is located at the central axes of pipeline its physical model, and x-axis direction is perpendicular to the ground Upwards, y-axis direction is parallel to ground but perpendicular to x-axis, and z-axis direction is the trend of pipeline.

The specific implementation process of the embodiment of the present invention described further below.

In order to realize accurate measurement spiral weld pipeline GIC, the embodiment of the present invention provides a kind of buried pipeline earth induction The non-contact measurement of electric current GIC, as shown in Fig. 2, the measuring device includes compound flux gate probe, power module, signal tune Manage module, microprocessor, memory module, local communication module and remote communication module module；Signal-regulated kinase, microprocessor Device, memory module, local communication module, the signal processing module of the remote communication module module composition measuring device；Compound magnetic Open gate probe connection Signal-regulated kinase, Signal-regulated kinase connection microprocessor；Signal-regulated kinase is used for compound flux The signal of gate probe collection is nursed one's health, and is converted to digital signal and is sent to microprocessor；Microprocessor connects storage respectively Module, local communication module, remote communication module module；Microprocessor obtains the number of magnetic induction intensity from Signal-regulated kinase Word signal, and further handle；Memory module is used for the magnetic induction intensity data that storage microprocessor is handled well；Local communication mould Block be used for the peripheral device communication such as local PC, the data of storage are sent to the external equipments such as local PC；Remote communication module For being communicated with remote server, the data of storage are sent to remote server and are further handled；Power module connects respectively Compound flux gate probe, Signal-regulated kinase, microprocessor, local communication module, remote communication module, memory module, are whole A measuring device provides stabilized power source.

The non-contact measurement of buried pipeline geomagnetic induction current GIC described in the present embodiment, the compound flux door are visited Head includes a bar shaped bar, 2 simple component fluxgate sensors and 2 level meters.By a bar shaped bar, 2 simple component magnetic fluxs The compound flux gate probe of door sensor and 2 level instrument devices composition, a length of 2 (L+d3) of bar shaped bar, it is wide and it is high be D, Central point is O points, each at distance O points both sides L to fix simple component fluxgate sensor M1, a M2；Simple component fluxgate senses The length and width of device M1, M2 are respectively d1, d2；The central point of simple component fluxgate sensor M1 is d3 apart from the left hand edge of bar shaped bar, The central point of M2 is d3 apart from the right hand edge of bar shaped bar, its measurement direction is x coordinate direction；The table at bar shaped bar left hand edge E Two horizon rules S1, S2 are fixed in face, and according to xyz coordinate directions, horizon rule S1 is located at phase with simple component fluxgate sensor M1, M2 Same yoz faces, horizon rule S2 are located at proximal surface xoz faces, and horizon rule length is D, takes L=b pi/2s, wherein 2 π b are defined as helix current The pitch of line.

With reference to Fig. 2 and Fig. 3, compound flux gate probe connection Signal-regulated kinase, Signal-regulated kinase connects microprocessor, Signal-regulated kinase is nursed one's health the analog signal of 2 simple component fluxgate sensors, and is converted to digital signal and is sent to Microprocessor.

The characteristics of for spiral-selded pipe road geomagnetic induction current GIC, can be according to the mathematical model of spiral weld pipeline GIC Formula (1) be calculated current elements as shown in Figure 7 P points produce magnetic field, calculate spiral weld pipeline in electric current I_GICWith pipe A point P (x above road₀,y₀,z₀) x, tri- axis magnetic induction density B of y, z_x,B_y,B_zRelation, such as the number of spiral weld pipeline GIC Shown in the formula (13) (14) (15) for learning model, as long as measuring the x of certain point above spiral weld pipeline, the one of side of y, z To magnetic induction intensity component, substitute into one of formula, you can obtain electric current I_GIC。

But at the scene in the practical operation of environment, it is thus necessary to determine that specific location, trend and the spiral shell of spiral weld pipeline The GIC optimum measuring points of rotation weld seam pipeline can just obtain accurate electric current I_GIC.The x of its midpoint P coordinates₀Value is point P distance pipes The vertical range in road axle center, is considered as constant.

Assuming that spiral weld pipeline determined by scene is substantially moved towards as z ' directions, true trend is z directions, such as Fig. 4 institutes Show.Using spiral-selded pipe road geomagnetic induction current GIC measuring devices, according to buried spiral weld pipeline B_z′Curve peak valley positions Method, obtains the accurate trend and GIC optimum measuring points of buried spiral weld pipeline, finally obtains accurate I_GIC。

Analyzed more than, the embodiment of the present invention is filled using the non-contact measurement of buried pipeline geomagnetic induction current GIC Put, the magnetic induction intensity produced to spiral line current in spiral weld pipeline is acquired and the calculating of electric current, finally realizes spiral shell Revolve the accurate measurement of welded tube road geomagnetic induction current GIC.Below to the contactless survey of buried pipeline geomagnetic induction current GIC Amount method is described in detail.

(1) mathematical model of spiral weld pipeline GIC is established：

For spiral weld pipeline due to the effect of weld seam, the electric current in pipeline is no longer along duct orientation streamlined flow, but It is threadingly advanced around pipe circumference.Assuming that pipeline is endless, then its single current line equation meets helix parametric equation

Wherein a is the radius of spiral weld pipeline,It is vectorial in xoy planes for what is a little formed on helix with o points Projection and the angle of x-axis, 2 π b are the pitch of helix, and z directions are moved towards for spiral weld pipeline, and x directions are perpendicular to ground level Upwardly direction, y directions are moved towards perpendicular to spiral weld pipeline and parallel to ground levels, as shown in Figure 6.

For the ease of research spiral weld pipeline inner vortex electric current in Distribution of Magnetic Field caused by pipeline, we are false If the spiral line current GIC flowed in spiral weld pipeline is the spiral line current of single endless, its electric current I_GICAlong z-axis Advance in direction.

First according to Biot -- Sa cuts down your theorem and lists a current elements on helixA point P (x above pipeline₀, y₀,z₀) caused by magnetic induction intensity, as shown in Figure 7.

Current elementsIt is in the magnetic induction intensity that P points produce

Wherein μ₀For the magnetic conductivity in vacuum,For current elementsThe vector of P points is directed toward, R is P points to current elementsMinimum distance.

IfRespectively x, y, the unit vector of z-axis, then point P and current elementsRespectively can table with vector form It is shown as

Point P is relative to current elementsPosition vector be represented by

AndComponent form be represented by

Then

WillFormula (2) is substituted into, current elements can be obtained and divided what point P was produced along x, the magnetic induction intensity in tri- directions of y, z Amount is respectively

Wherein

Magnetic induction intensity of the above current elements in all directions integrates endless helix can obtain it is whole Spiral weld pipeline electric current is in induction level caused by P points, i.e.,

Wherein L represents the spiral line current of whole spiral weld pipeline.The all directions that whole spiral line current is produced in point P Magnetic induction intensity value is respectively

Specific analytical expression can not be obtained with upper integral, final integral result can be tried to achieve using numerical integrating.

(2) orientation of buried spiral weld pipeline is obtained：

Step 1, as shown in figure 4, spiral weld pipeline determined by scene is substantially moved towards as z ' directions, by the present embodiment institute The yoz faces for the GIC measuring devices stated are used as mobile face, make horizon rule S2 above, and adjust the axial direction of S2 and bar shaped bar, protect The measurement direction for demonstrate,proving simple component fluxgate sensor M1 is z ' directions.

Step 2, a certain bar path M is selected, its direction is carried out simple component fluxgate sensor M1 along M perpendicular to z ' It is mobile, ensure that movement track is fully passed over directly over buried pipeline, and record M1 surveys B in real time_z′Curve.According to " buried spiral Weld seam pipeline B_z′Curve peak valley positioning mode " analyzes B_z′Curve, judges whether path M passes through symmetric points (x₀, 0,2nb π) or (x₀, 0, (2n+1) b π) and record corresponding symmetric points position.

Step 3, path M ' parallel and distance l on the right side of it with selected path in abovementioned steps 2 is selected, by single point Amount fluxgate sensor M1 is moved along M ', ensures that movement track is fully passed over directly over buried pipeline, and record M1 in real time Surveyed B_z′Curve.According to " buried spiral weld pipeline B_z′Curve peak valley positioning mode " analyzes B_z′Curve, it is first determined spiral weld The direction of GIC in pipeline, then judge whether path M ' passes through symmetric points (x₀, 0,2nb π) or (x₀, 0, (2n+1) b π) and record phase Answer symmetric points position.

Step 4, path M " parallel and distance l on the right side of it with 3 path of abovementioned steps is selected, by simple component fluxgate Sensor M1 is moved along M ", ensures that movement track is fully passed over directly over buried pipeline, and record M1 surveys B in real time_z′It is bent Line.According to " buried spiral weld pipeline B_z′Curve peak valley positioning mode " analyzes B_z′Curve, first determines whether the direction of GIC in pipeline Whether overturn.If upset, rebound step 2 remeasure.If not overturning, judge whether path M " passes through symmetric points (x₀,0, 2nb π) or point (x₀, 0, (2n+1) b π) and record corresponding symmetric points position.

Step 5, abovementioned steps 4 are repeated, according to " buried spiral weld pipeline B_z′Curve peak valley positioning mode " analyzes B_z′It is bent Line, until finding by symmetric points (x₀, 0,2nb π) and point (x₀, 0, (2n+1) b π) neighbouring two paths, stop repeating to walk Rapid 4, and determine two points substantially position P₁, P₂, 2 lines are new z ' directions.

Step 6, the distance between two paths l is reduced, abovementioned steps 2 is repeated to step 5, concentrates analysis site P₁, P₂Near The B of path_z′Curve, until finding and B_z′Point (the x that the wave crest (trough) of curve is completely superposed₀, 0,2nb π) and point (x₀,0,(2n + 1) b π) untill, this 2 line, that is, pipelines truly move towards, i.e. z directions.It can use n=0 at this time, then point (x₀, 0,0) surveyed to be optimal Amount point.

(3) influence of the interference field to GIC computational solution precisions is filtered out：

Behind the orientation and optimum measuring point of spiral weld pipeline to be determined, by the GIC measuring devices described in the present embodiment Central point O is placed in optimum measuring point (x₀, 0,0) and place, wherein x₀, can be for the sum of the distance of measurement point to pipeline and pipe radius Scene obtains, it is believed that it is known constant.2 simple component fluxgate sensors are placed in top, adjust the z-axis side of bar shaped bar To making it truly move towards to overlap with pipeline.Horizon rule S2 is adjusted, the measurement direction for making simple component fluxgate sensor is x-axis side To.Fixed compound flux gate probe, as shown in figure 5, the upper two fluxgate sensor points P that then pops one's head in_mLAnd P_mRCoordinate value RespectivelyWithThe measurement data B of two sensors of measurement in real time_xmL,B_xmR, and by measurement data Formula described below (56) is substituted into, interference field data is filtered out by microprocessor, calculates point P_mLSecondary field data B_xL, Preserved in real time using memory module, and remote server is uploaded to by remote communication module module.

In specific implementation, the ground surface magnetic field above spiral weld pipeline is primary field (earth's magnetic field) and secondary field (pipe The magnetic field that road GIC is produced) composition Composite Field, and primary field is interference field.If only the magnetic measured by Magnetic Sensor Induction data substitute into formula (13) (14) (15), and the current value error calculated is very big, so needing to analyze magnetic induction The feature of intensity all directions component, obtains symmetric points, to filter out the error that interference field is brought to result.And due to formula (13) (14) (15) expression formula is complicated, it is difficult to obtains the analytical expression of its integration, can not intuitively observe magnetic field intensity Distribution character.And x the and y function cycles of helix parametric equation are 2 π, and each function is respectively provided with symmetry characteristic.In order to more The magnetic induction intensity expression formula that whole current line is produced in certain point is easily obtained, can first take one for x or y in helix A cycle portions, analysis signal period current line produce the magnetic field expression formula in magnetic field on its symmetric points, so that simplification integrated Journey.But, then can be respectively to a cycle part of x and y functions since the signal period part of sine and cosine functions is misaligned Studied.The filtering method of field data is disturbed to be specifically described in the magnetic field produced below to spiral weld pipeline GIC.

(1) magnetic induction intensity symmetrical analysis is produced for the spiral line current of symmetry axis with x-axis (y=0, z=0)

Take sectionOne section of current line, using x-axis as symmetry axis, take point P being located at pipeline just in x-axis Top, the coordinate for the P that sets up an office is P (x₀, 0,0), a current elements are respectively taken in symmetry axis both sidesWithCalculate two The magnetic induction intensity that current elements are produced in point P, as shown in Figure 8.

Current elementsIn P (x₀, 0,0) produce magnetic field each axis induction weight expression be

Wherein

Current elementsIn P (x₀, 0,0) produce magnetic field each axis induction weight expression be

Wherein

Two current elements positions are symmetrical with x-axis, i.e.,The increment of two current elementsAccording to public affairs Formula (16)~(23), can obtain

According to result above, using x-axis as symmetrical monocycle (single turn) helix current line, two symmetrical current members are right Claim magnetic induction intensity caused by the point P in axis x-axis (y=0, z=0) that there are following characteristics：

1)dB_px=dB_1x+dB_2x=0, i.e. the magnetic induction intensity in x-axis direction is 0；

2)dB_py=2dB_1y=2dB_2y, dB_pz=2dB_1z=2dB_2z, i.e. the magnetic induction intensity in y-axis and z-axis direction is symmetrical 2 times of the magnetic induction intensity that axis single sided current member produces.

Above property can be fromExtend toThe helix current line of whole pipeline can be obtained in x-axis (y =0, z=2nb π) magnetic induction intensity summation caused by upper point P is following result：

WhereinL ' is the helix of z-axis positive axis.

(2) magnetic induction intensity symmetrical analysis is produced with the spiral line current that (y=0, z=b π) is symmetry axis

Take sectionOne section of current line, symmetry axis is used as with the straight line (y=0, z=b π) parallel to x-axis, Point P is taken to be located at the surface of pipeline on this straight line, the coordinate for the P that sets up an office is P (x₀, 0, b π), respectively take one in symmetry axis both sides Current elementsWithThe magnetic induction intensity that two current elements are produced in point P is calculated, as shown in Figure 9.

As shown in Figure 9, it is assumed that current elements'sCurrent elementsThen WithIn P (x₀, 0, b π) produce magnetic field each axis induction weight expression For

Wherein

Wherein

According to result above, R₁=R₂, with straight line (y=0, z=b π) for symmetrical monocycle (single turn) current line, Magnetic induction intensity caused by a point P of the two symmetrical current members on symmetry axis has following characteristics：

1)dB_px=dB_1x+dB_2x=0, i.e. the magnetic induction intensity in x-axis direction is 0；

2)dB_py=2dB_1y=2dB_2y, dB_pz=2dB_1z=2dB_2z, i.e. the magnetic induction intensity in y-axis and z-axis direction is symmetrical 2 times of the magnetic induction intensity that axis single sided current member produces.

Above property can be fromExtend toThe helix current line of whole pipeline can be obtained symmetrical Magnetic induction intensity summation caused by the upper point P of axis (y=0, z=nb π) is following result：

WhereinHelix.

(3) spiral line current point P between (y=0, z=2nb π) and (y=0, z=nb π) two symmetry axis produces magnetic Induction is analyzed

In x-axis (y=0, z=0) two side draws, two symmetric pointsPositioned at spiral-selded pipe The surface in road, and meetTwo current elements are calculated respectivelyIn point P₁, P₂The magnetic induction intensity of generation, As shown in Figure 10.

Two current elements are calculated in point P₁, P₂Caused magnetic induction intensity, can obtain according to formula (2)

It can be obtained according to formula (6)

Assuming thatIt can be obtained according to Fig. 6

R₁₂=R₂₁ (46)

R₁₁=R₂₂ (47)

Result above is substituted into formula (7) (8) can obtain the magnetic induction intensity of P1 and P2 point x-axis

Understand

dB_P1x=-dB_P2x (52)

Then symmetrical current member all on whole helix current line x-axis component magnetic induction caused by the point P1, P2 is strong Degree, i.e., quadrature whole current line, can obtain

Wherein L is whole helix current line.

Formula (53) shows, is pointed to directly over pipeline, and be symmetrical with 2 points of straight line (y=0, z=nb π)The x-axis component induction level phase of spiral weld pipeline GIC generations at this 2 points Deng direction is opposite.

In conclusion spiral weld pipeline induced electricity stream I_GICIt is x with pipeline axle center distance directly over pipeline₀Produce in place Raw x-axis magnetic induction density B x along duct orientation be in cyclically-varying, the symmetric position in straight line (y=0, z=nb π) both sides Place, it may appear that the forward and reverse maximum of Bx, its change curve are as shown in figure 11.

(4) filtering method of field data is disturbed

Conclusion is analyzed more than, directly over pipeline at z=2nb π, i.e. B_xTwo side draws, the two point P of symmetrical shaft position_L,P_R, Make 2 points away from symmetry axis apart from equal and directly over the pipeline, its coordinate is As shown in figure 12.

Assuming that measure P using sensor_L,P_R2 points of x-axis direction magnetic induction intensity is B_xmL,B_xmR, its numerical value is earth's magnetic field The superposition for the secondary field x-axis direction magnetic induction intensity that x-axis direction magnetic induction intensity and spiral line current produce, it is assumed that helix electricity Stream is in P_L,P_RThe x-axis direction magnetic induction intensity of 2 points of generations is respectively B_xL,B_xR, since 2 points closer to the distance, its earth induction is strong Number of degrees value approximately equal, then it is believed that 2 points of earth induction intensity is B_xG, then

B_xmL=B_xG+B_xL (54)

B_xmR=B_xG+B_xR (55)

According to formula (53) (54) (55) the secondary field induction of helix current line generation at 2 points can be derived with surveying The relation of value is

Gained B_xLAs filter out the secondary field magnetic induction intensity after interference field data.By way midpoint P_LLocate the B of sensor_xLGeneration Enter formula (13), finally obtain accurate spiral-selded pipe road geomagnetic induction current GIC, i.e.,

Wherein

In above formula, a is pipe radius, and b is the pitch constant (2 π b are pitch) of helix current line, μ₀For Vacuum Magnetic Conductance, x₀,y₀,z₀For Figure 12 midpoints P_LCoordinate value.Bring the pipeline data at scene and definite measurement point coordinate into above formula, Unknown quantity only hasThe integration in above formula is calculated by the numerical integration method of computer, concrete numerical value can be obtained

I can then be obtained_GICWith B_xLRelation be

(4) GIC is calculated using numerical integrating：

The first step, determines pitch constant b, the space permeability μ of spiral weld pipeline radius a, helix current line₀And point P_mLCoordinate valueDeng field pipes data, and as known parameters；

Second step, the magnetic induction intensity data B that the live GIC measuring devices that remote server is received are sent_xLAfterwards, The parameter and B that the first step is determined_xLFormula (57-3) is substituted into, show that final spiral-selded pipe is genuine using numerical integration algorithm Induced field current GIC.

According to the analysis of the above-mentioned contactless measurement to buried pipeline geomagnetic induction current GIC, if it is desired to essence The true GIC for calculating spiral weld pipeline, it is necessary to first find B_xDirectly over symmetrical shaft position, i.e. pipeline at z=2nb π, then Determine that two optimum measuring points carry out the measurement of magnetic induction intensity and GIC is calculated.But in the operation of actual field, underground pipe Although there are the mark for representing pipeline location, but the exact position of pipeline center, the position of z=2nb π and pipeline in the top in road Trend is all difficult to obtain, though by the equipment of specialty measurement buried metal, also can only obtain the general paving location of pipeline with And rough move towards information.So need to study the Magnetic Induction Density Distribution that buried spiral weld pipeline GIC is produced on the ground Characteristic, to draw a kind of acquisition pipeline accurately trend and the method for optimum measuring point.

Spiral weld pipeline GIC optimum measurement independent positioning methods：

To determine buried spiral weld pipeline exact position, symmetric points (x=x₀, y=0, z=2nb π) and pipeline walk To, it is necessary to when the position for studying a point P on pipeline changes in yoz planes, the regularity of distribution of its magnetic induction intensity.Due to spiral Weld seam pipeline actual position and move towards unknown, in this case, measurement direction and the duct orientation (z-axis of fluxgate sensor Direction) it is inconsistent, for the ease of problem analysis and calculating, now do hypothesis below：Assuming that identified pipeline is substantially walked at the scene To being λ with the actual angle moved towards of buried pipeline, its relation is as shown in figure 13.

Z ' directions are substantially moved towards for buried pipeline determined by scene, and z directions are the actual trend of buried pipeline, Magnetic Sensor Measurement direction be z ' directions magnetic induction intensity, M paths direction is the direction vertical with z ' directions, its path is that detection is buried Measurement track during ground pipeline specific location, its equation of locus are

WhereinThe as intersection point P of M paths and z-axis₀Coordinate.

According to similar triangles theorem, it is known that the magnetic induction intensity value in Magnetic Sensor measurement gained z ' directions should be

B_z′=B_z cosλ-B_y sinλ (59)

If measurement plane be yoz planes above pipeline, study above spiral line current it is mutually parallel, be spaced equal one B on group M paths_z′Numeric distribution.Spiral line current where different M paths is upper to put as shown in figure 14.

9 measurement locus equations are (58) in upper figure, M1 to M9 paths in a z-direction at intervals of π/4, in its equationValue be respectively-π, -3 π/4,-pi/2,-π/4,0, π/4, pi/2,3 π/4, π.

Assuming that pipe radius a=0.5m, pitch parameter b=0.4, yoz plan range pipeline center distance x₀=1m, pipeline Interior electric current flows from left to right, with B_z′To analyze target, ifAccording to formula (59), Mathematica is utilized Numerical integrating in software, studies above the spiral line current B on the M paths of diverse location_z′/B₀Numeric distribution.

When spiral weld pipeline included angle X=10 ° that truly trend is substantially moved towards with scene gained, surveyed since y=-4 Amount, it is z ' directions to make Magnetic Sensor measurement direction, moves to y=-4, the B of 9 paths of the above_z′/B₀Distribution curve such as Figure 15 institutes Show, its abscissa is y values, ordinate B_z′/B₀。

When spiral weld pipeline included angle X=20 ° that truly trend is substantially moved towards with scene gained, surveyed since y=-4 Amount, it is z ' directions to make Magnetic Sensor measurement direction, moves to y=-4, the B of 9 paths of the above_z′/B₀Distribution curve such as Figure 16 institutes Show, its abscissa is y values, ordinate B_z′/B₀。

When spiral weld pipeline included angle X=30 ° that truly trend is substantially moved towards with scene gained, surveyed since y=-4 Amount, it is z ' directions to make Magnetic Sensor measurement direction, moves to y=-4, the B of 9 paths of the above_z′/B₀Distribution curve such as Figure 17 institutes Show, its abscissa is y values, ordinate B_z′/B₀。

According to upper curve, M5 is to measure path by symmetry axis z=2nb π.At the scene implement pipeline detection and In GIC measurement operations, the angle that true pipeline rout and scene gained pipeline substantially move towards will not be excessive, substantially within 30 °. Experiment curv more than under three groups of difference angles, can obtain as drawn a conclusion：

1) in the case where 30 ° of λ ＜, by the magnetic induction density B measured by the path of symmetry axis z=2nb π_z′Curve is Symmetrical " the lowest point type " curve, i.e. Figure 15, the M5 paths in 16,17, the magnetic induction intensity of the lowest point are for other paths Minimum value.And the lowest point position is at z=2nb π, and directly over pipeline, its coordinate points is (x₀,0,2nbπ)；

2) in the case where 30 ° of λ ＜, passing point (x₀, 0, (2n+1) b π) path on B_z′Curve is symmetrical " peak Top-type " curve, i.e. Figure 15, M1 the and M9 paths in 16,17, the magnetic induction intensity of summit is maximum for other paths.And Summit position is at z=(2n+1) b π, and directly over pipeline, its coordinate points is (x₀,0,(2n+1)bπ)；

3) in the case where 30 ° of λ ＜, passing point (x₀, 0, (2n-1) b π) and point (x₀, 0,2nb π) between path on B_z′ Curve is asymmetric " peak valley is staggered " curve, i.e. Figure 15, M2, M3, M4 path in 16,17, its B_z′The lowest point position of curve On the left of curve, summit is located on the right side of curve；

4) in the case where 30 ° of λ ＜, passing point (x₀, 0,2nb π) and point (x₀, 0, (2n+1) b π) between path on B_z′ Curve is asymmetric " peak valley is staggered " curve, i.e. Figure 15, M6, M7, M8 path in 16,17, its B_z′The lowest point position of curve On the right side of curve, summit is located on the left of curve；

5) in the case where 30 ° of λ ＜, according to the measuring sequence of the path from M1 to M9, its B_z′The lowest point variation tendency of curve For：No the lowest point (M1) → peak valley interlocks, the lowest point left avertence, and numerical value is gradually reduced and moves to right the disappearance of (M2, M3, M4) → summit, the lowest point Placed in the middle and all paths less than before the lowest point value (M5) → peak valleys interlock, the lowest point right avertence, and numerical value gradually increases and continues to move to right (M6, M7, M8) → the lowest point disappears, and summit is placed in the middle and more than the summit value (M1) of all paths before.

6) according to conclusions, in the case that GIC senses of current are consistent with Figure 10 and do not overturn all the time, From a certain position of pipeline along z ' directions, the B on multigroup path is measured with fixed intervals_z′Curve, the change of its waveform peak valley become Gesture should be：

Maximum wave crest, without trough → there is left avertence trough and gradually move to right → maximum trough, without wave crest → and wave crest, ripple occur Paddy right avertence and gradually move to right → maximum wave crest, without trough.

When in the experiment curv of path without trough, the position that wave crest peak occurs is the point directly over pipeline (x₀, 0, (2n+1) b π) and place.When in the experiment curv of path without wave crest, the position that trough minimum point occurs is being located at pipeline just Point (the x of top₀, 0,2nb π) and place, this point is the optimum measuring point of spiral-selded pipe road geomagnetic induction current GIC.This 2 points Line parallel to buried pipeline axis, be the true trend of buried pipeline.

7) in the case that GIC senses of current are opposite with Figure 14 and do not overturn all the time, from a certain position of pipeline Along z ' directions, the B on multigroup path is measured with fixed intervals_z′Curve, the variation tendency of its waveform peak valley should be：

Maximum trough, without wave crest → there is left avertence wave crest and gradually move to right → maximum wave crest, without trough → and trough, ripple occur Peak right avertence and gradually move to right → maximum trough, without wave crest.

When in the experiment curv of path without wave crest, the position that trough minimum point occurs is the point directly over pipeline (x₀, 0, (2n+1) b π) and place.

When in the experiment curv of path without trough, the position that wave crest peak occurs is the point directly over pipeline (x₀, 0,2nb π) and place, this point is the optimum measuring point of spiral-selded pipe road geomagnetic induction current GIC.This 2 points line is put down Row is that the true of buried pipeline moves towards in buried pipeline axis.

8) when in measurement process, the B of a certain path is being measured_z′GIC senses of current change during curve, its B_z′It is bent Line will be overturn, its trend must be unsatisfactory for above analysis acquired results, then need to re-start measurement.

9) due to that be able to might not have just past symmetric points (x in selected multigroup path in measurement process₀,0,2nb π) point (the x of sum₀, 0, (2n+1) b π) path, then need to take multiple measurements, continuous contract measurement scope, and reduce path Interval, finally to determine accurate symmetric points position.

The main B by being measured on path of above method_z′The change of the wave crest and trough of curve judges symmetric points and pipe Road moves towards, therefore above measuring method is defined as " buried spiral weld pipeline B_z′Curve peak valley positioning mode ".

The calculating of spiral weld pipeline GIC：

According to upper one section to it is different measurement geometric locuses research it is concluded that, we can smoothly determine buried spiral shell Revolve the optimum measuring point (x of welded tube road geomagnetic induction current GIC₀, 0,2nb π) and buried pipeline true trend.Once really Determine the optimum measuring point (x of spiral weld pipeline GIC₀, 0,2nb π) after, we can be along the axial direction of pipeline, in pipeline just Top optimum measuring point (x₀, 0,2nb π) place both sides symmetric position on respectively place 1 Magnetic Sensor, its measurement direction is adjusted To x directions and fixation, as shown in figure 18, it is assumed that two sensor point P_mL,P_mRCoordinate be

Two sensors measure it after fixing in real time, obtain data B_xmL,B_xmR, formula (56) is substituted into, i.e., Interference field can be filtered out.Result and two point coordinates values are substituted into formula (57) again, you can calculate accurate I_GIC。

Compared with prior art, buried pipeline geomagnetic induction current GIC's provided in an embodiment of the present invention is contactless Measuring device and method, specifically there is following feature：

1st, for the production technology characteristic of spiral weld pipeline, the stream of the geomagnetic induction current GIC in such pipeline is analyzed Dynamic direction is not straight line, but is threadingly advanced.This physical model with pipe current to be equivalent to unlimited long straight conductor at present Grave fault, the GIC resultant errors that spiral weld pipeline is measured by the prior art are very big；

2nd, the mathematical model of the geomagnetic induction current GIC of spiral weld pipeline is established, and is derived on GIC electric currents and pipeline Side surveys clear and definite mathematical relationship between magnetic data；

3rd, the physical model of the geomagnetic induction current GIC of spiral weld pipeline and endless constant current model difference are very big, The method of traditional elimination background magnetic field is not suitable for spiral weld pipeline.This method analyzes helix current line above pipeline The symmetry in magnetic field is produced, and by measuring the magnetic induction intensity of symmetric points, geomagnetic noise magnetic field can be obtained, in calculating current Eliminated, substantially reduce the influence of ambient field, and improve the accuracy of calculating；

4th, using 2 single shaft fluxgate sensors, 1 bar shaped probe is formed, 2 sensors are placed in probe both sides, and Center probe is placed at the symmetric points of spiral line current, makes probe axial axially in parallel with pipeline, be totally different from big at present Part is axially placed in pipeline by probe vertical；

5th, the pipeline location and trend obtained at the scene by pipeline detection device is not accurate enough, this also causes magnetic to sense Measuring result error caused by device riding position is unreasonable.The present invention proposes one kind and is directed to GIC spies in spiral weld pipeline The measuring device of point, can accurately obtain pipeline by " buried spiral weld pipeline curve peak paddy positioning mode " proposed by the present invention Position and true trend.

The above is only that the preferred embodiment of the present invention is described, and is not that the scope of the present invention is limited It is fixed；The foregoing is merely the embodiment of the present invention, the protection domain being not intended to limit the present invention；Do not departing from On the premise of design spirit of the present invention, what this area ordinary skill technical staff made technical scheme any repaiies Change, equivalent substitution, improvement etc., should all fall into the protection domain that claims of the present invention determines.

Patent of the present invention obtains state key research and development planning item (No.2016YFC0800103) fund assistance.

Claims (9)

1. A kind of 1. non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC, for gathering spiral shell in spiral weld pipeline The magnetic induction intensity data that spin line electric current produces and the calculating for realizing electric current, it is characterised in that：

   The measuring device includes compound flux gate probe, power module, Signal-regulated kinase, microprocessor, memory module, local Communication module and remote communication module；

   The Signal-regulated kinase, microprocessor, memory module, local communication module, remote communication module form measurement dress The signal processing module put；

   The compound flux gate probe connects Signal-regulated kinase, and the Signal-regulated kinase connects microprocessor；The signal Conditioning module is used to be nursed one's health the signal that compound flux gate probe gathers, and is converted to digital signal and is sent to microprocessor Device；

   The microprocessor connects memory module, local communication module, remote communication module respectively；The microprocessor is from signal The digital signal of magnetic induction intensity is obtained in conditioning module, and is further handled；The memory module is used for storage microprocessor The magnetic induction intensity data handled well；The local communication module be used for local peripheral device communication, the data of storage are passed Send to local external equipment；The remote communication module is used to communicate with remote server, and the data of storage are sent to remotely Server is further handled；

   The power module connects compound flux gate probe, Signal-regulated kinase, microprocessor, local communication module, remote respectively Journey communication module, memory module, stabilized power source is provided for whole measuring device.
2. 2. the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC as claimed in claim 1, it is characterised in that： It is horizontal that the compound flux gate probe includes at least one bar shaped bar, at least two simple component fluxgate sensor and at least two Instrument.
3. 3. the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC as claimed in claim 2, it is characterised in that： The compound flux gate probe being made of a bar shaped bar, 2 simple component fluxgate sensors and 2 level instrument devices, bar A length of 2 (L+d3) of shape bar, it is wide and it is high be D, central point is O points, each at distance O points both sides L to fix a simple component magnetic flux Door sensor M1, M2；The length and width of simple component fluxgate sensor M1, M2 are respectively d1, d2；Simple component fluxgate sensor M1 Central point apart from the left hand edge of bar shaped bar be d3, the central point of M2 is d3 apart from the right hand edge of bar shaped bar, its measurement direction is x Coordinate direction；Two horizon rules S1, S2 are fixed in surface at bar shaped bar left hand edge E, according to xyz coordinate directions, horizon rule S1 It is located at identical yoz faces with simple component fluxgate sensor M1, M2, horizon rule S2 is located at proximal surface xoz faces, and horizon rule length is D, L=b pi/2s are taken, wherein 2 π b are defined as the pitch of helix current line.
4. 4. the non-contact measurement apparatus of the buried pipeline geomagnetic induction current GIC as described in claim 1 or 3, its feature exist In：The compound flux gate probe connects Signal-regulated kinase, and the Signal-regulated kinase connects microprocessor, the signal tune Reason module is nursed one's health the analog signal of 2 simple component fluxgate sensors, and is converted to digital signal and is sent to microprocessor Device.
5. A kind of 5. contactless measurement of buried pipeline geomagnetic induction current GIC, for spiral-selded pipe road earth induction The characteristics of electric current GIC, using the non-contact of buried pipeline geomagnetic induction current GIC according to any one of claims 1 to 4 Formula measuring device, the magnetic induction intensity produced to spiral line current in spiral weld pipeline is acquired and the calculating of electric current, with Realize the accurate measurement of spiral-selded pipe road geomagnetic induction current GIC, it is characterised in that：The measuring method includes：

   (1) mathematical model of spiral-selded pipe road geomagnetic induction current GIC is established；

   (2) orientation of buried spiral weld pipeline is obtained；

   (3) influence of the interference field to GIC computational solution precisions is filtered out；

   (4) GIC is calculated using numerical integrating.
6. 6. the contactless measurement of buried pipeline geomagnetic induction current GIC as claimed in claim 5, it is characterised in that： In described (1) item, the characteristics of for spiral-selded pipe road geomagnetic induction current GIC, the genuine magnetic strength of spiral-selded pipe is established The mathematical model of induced current GIC, spiral weld pipeline is calculated according to the mathematical model of spiral-selded pipe road geomagnetic induction current GIC Interior electric current I_GICWith a point P (x above pipeline₀,y₀,z₀) x, tri- axis magnetic induction density B of y, z_x,B_y,B_zRelation, be calculated Electric current I_GIC, the modeling process of the mathematical model of spiral-selded pipe road geomagnetic induction current GIC specifically includes：

   Electric current in spiral weld pipeline is threadingly advanced around pipe circumference, it is assumed that pipeline is endless, then its single current line Equation meets that helix parametric equation is

   Wherein a is the radius of spiral weld pipeline,For a little vectorial projection in xoy planes with o points composition on helix With the angle of x-axis, 2 π b be helix pitch, z directions for spiral weld pipeline move towards, x directions for perpendicular to Horizon upwardly Direction, y directions perpendicular to spiral weld pipeline move towards and parallel to ground level；

   Assuming that the spiral line current GIC flowed in spiral weld pipeline is the spiral line current of single endless, its electric current I_GICEdge The advance of z-axis direction, then first according to Biot -- Sa cuts down your theorem and lists a current elements on helixOne above pipeline Point P (x₀,y₀,z₀) caused by magnetic induction intensity be

   <mrow> <mi>d</mi> <mover> <mi>B</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <mfrac> <msub> <mi>&amp;mu;</mi> <mn>0</mn> </msub> <mrow> <mn>4</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msub> <mi>I</mi> <mrow> <mi>G</mi> <mi>I</mi> <mi>C</mi> </mrow> </msub> <mi>d</mi> <mover> <mi>l</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;times;</mo> <mover> <mi>R</mi> <mo>&amp;RightArrow;</mo> </mover> </mrow> <msup> <mi>R</mi> <mn>3</mn> </msup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

   Wherein μ₀For the magnetic conductivity in vacuum,For current elementsThe vector of P points is directed toward, R is P points to current elements's Minimum distance；

   IfRespectively x, y, the unit vector of z-axis, then point P and current elementsIt is represented by respectively with vector form

   <mrow> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mover> <mi>i</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>+</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> <mover> <mi>j</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>+</mo> <msub> <mi>z</mi> <mn>0</mn> </msub> <mover> <mi>k</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

   Point P is relative to current elementsPosition vector be represented by

   AndComponent form be represented by

   Then

   WillFormula (2) is substituted into, current elements is obtained and divides what point P was produced along x, the magnetic induction intensity component in tri- directions of y, z It is not

   Wherein

   Magnetic induction intensity of the above current elements in all directions integrates endless helix and obtains whole spiral welding Pipe current is stitched in induction level caused by P points, i.e.,

   <mrow> <mi>B</mi> <mo>=</mo> <munder> <mo>&amp;Integral;</mo> <mi>L</mi> </munder> <mi>d</mi> <mi>B</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>

   Wherein L represents the spiral line current of whole spiral weld pipeline, all directions magnetic strength that whole spiral line current is produced in point P The intensity level is answered to be respectively
7. 7. the contactless measurement of the buried pipeline geomagnetic induction current GIC as described in claim 5 or 6, its feature exist In：In described (2) item, determine that the GIC of the specific location of spiral weld pipeline, trend and spiral weld pipeline is most preferably surveyed Point is measured, is specifically comprised the following steps：

   Step 1, the definite spiral weld pipeline in scene is substantially moved towards as z ' directions, by the non-of buried pipeline geomagnetic induction current GIC The yoz faces of contact type measurement are used as mobile face, make horizon rule S2 above, and adjust the axial direction of S2 and bar shaped bar, ensure magnetic The measurement direction of open gate sensor M1 is z ' directions；

   Step 2, a path M is selected, its direction is moved simple component fluxgate sensor M1 along M perpendicular to z ', is protected Card movement track is fully passed over directly over buried spiral weld pipeline, and record Component of magnetic flux door sensor M1 surveys B in real time_z′ Curve；According to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′Curve, judges whether path M passes through symmetric points (x₀, 0,2nb π) or (x₀, 0, (2n+1) b π) and record corresponding symmetric points position；

   Step 3, path M ' parallel and distance l on the right side of it with selected path in above-mentioned steps 2 is selected, by simple component magnetic Open gate sensor M1 is moved along M ', ensures that movement track is fully passed over directly over buried spiral weld pipeline, and remember in real time Record simple component fluxgate sensor M1 surveys B_z′Curve；According to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′ Curve, it is first determined the direction of GIC in spiral weld pipeline, then judge whether path M ' passes through symmetric points (x₀, 0,2nb π) or (x₀, 0, (2n+1) b π) and record corresponding symmetric points position；

   Step 4, path M " parallel and distance l on the right side of it with 3 path of above-mentioned steps is selected, simple component fluxgate is sensed Device M1 is moved along M ", ensures that movement track is fully passed over directly over buried spiral weld pipeline, and record simple component in real time Fluxgate sensor M1 surveys B_z′Curve；According to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′Curve, first Judge whether the direction of GIC in spiral weld pipeline overturns, if upset, rebound step 2 remeasures, if not overturning, sentences Whether broken rail line M " passes through symmetric points (x₀, 0,2nb π) or point (x₀, 0, (2n+1) b π) and record corresponding symmetric points position；

   Step 5,4 are repeated the above steps, according to buried spiral weld pipeline B_z′Curve peak valley positioning mode analyzes B_z′Curve, until Find by symmetric points (x₀, 0,2nb π) and point (x₀, 0, (2n+1) b π) neighbouring two paths, stop repeat step 4, and really Fixed two points substantially position P₁, P₂, 2 lines are new z ' directions；

   Step 6, the distance between two paths l is reduced, 2 is repeated the above steps to step 5, concentrates analysis site P₁, P₂Neighbouring path B_z′Curve, until finding and B_z′Point (the x that the wave crest or trough of curve are completely superposed₀, 0,2nb π) and point (x₀,0,(2n+1)b Untill π), this 2 line, that is, pipelines truly move towards, i.e. z directions；It can use n=0 at this time, then point (x₀, 0,0) and it is optimum measurement Point.
8. 8. the contactless measurement of the buried pipeline geomagnetic induction current GIC as any one of claim 5,6,7, It is characterized in that：Described (3) item filters out influence of the interference field to GIC computational solution precisions, and it is spiral to treat that described (2) item determines Behind the orientation and optimum measuring point of weld seam pipeline, by the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC Heart point O is placed in optimum measuring point (x₀, 0,0) and place, wherein x₀It is existing for the sum of the distance of measurement point to pipeline and pipe radius Field obtains, it is believed that it is known constant, using the non-contact measurement apparatus of buried pipeline geomagnetic induction current GIC, by 2 lists Component of magnetic flux door sensor is placed in above buried spiral weld pipeline, is adjusted the z-axis direction of bar shaped bar, it is truly walked with pipeline To coincidence, horizon rule S2 is adjusted, the measurement direction for making fluxgate sensor is x-axis direction, and fixed compound flux gate probe, then visit Upper two fluxgate sensor points P_mLAnd P_mRCoordinate value be respectivelyWithMeasurement in real time The measurement data B of two sensors_xmL,B_xmR, and measurement data is substituted into formula (56), i.e.,By micro- Processor filters out interference field data, calculates point P_mLSecondary field data B_xL, preserved in real time using memory module, and by remote Journey communication module is uploaded to remote server.
9. 9. the non-contact measurement side of the buried pipeline geomagnetic induction current GIC as any one of claim 5,6,7,8 Method, it is characterised in that：Described (4) item calculates GIC using numerical integrating, specifically includes：

   Step 1, pitch constant b, the space permeability μ of pipe radius a, helix current line are determined₀And point P_mLCoordinate valueDeng field pipes data, and as known parameters；

   Step 2, the non-contact measurement apparatus hair of live buried pipeline geomagnetic induction current GIC remote server received The magnetic induction intensity data B sent_xLAfterwards, the parameter and B step 1 determined_xLFormula (57-3) is substituted into, i.e., Final spiral-selded pipe road geomagnetic induction current GIC is drawn using numerical integration algorithm.