Gasum Now!: Method of calculation

Temporal variations of the magnetic field are associated with an electric field, which in turn drives currents in conductors (earth, pipelines, etc.). To be able to calculate GIC, we need the following input: the geomagnetic field, earth's conductivity model, and a DC model of the pipeline system. Using the magnetic field and the conductivity model, we can calculate the electric field at the earth's surface. Then this allows for the computation of GIC in the pipeline.

There are 9 stations in Finland continuously measuring the geomagnetic field (as a part of the IMAGE magnetometer network), of which the Nurmijärvi Geophysical Observatory is located close to the pipeline system. Magnetic recordings are saved at each second, but in this project we use 10 second averages. Data files are transfered to the server at regular intervals for the calculation of the electric field.

Calculation of the electric field from the magnetic field requires a model of the earth's conductivity. As a starting point, we have used recent models, of which the most recent ones are provided by the BEAR project on the Fennoscandian Shield. The earth's conductivity in Finland varies very much from place to place. However, it is lucky that concerning GIC, small-scale variations are not important. Because the Finnish pipeline covers a relatively small region, we assume a single layered earth model for the whole region, which is also computationally simple. With some experimentation, we have noticed that even a two-layer model yields reasonable results.

When we know the electric field, GIC in the pipeline is determined using an advanced model (DSTL = distributed source transmission line method), which is applicable for any buried pipeline network. Necessary pipeline input parameters are its conductivity and topology. A simplifying fact is that geomagnetic variations are slow so that GIC is nearly DC. Further, it is quite accurate to assume that the electric field is spatially uniform in the area of the pipeline system. Then GIC at each point in the pipeline is a linear function of the two horizontal components of the electric field. The DSTL method allows for using a spatially varying electric field too.

The method has been validated with comparisons to measured values since November 1998 at Mäntsälä. The median error of calculated near-real-time values is about 30%, when GIC of at least 1 A are considered. One error source is that the electric field is not exactly spatially uniform near Mäntsälä. Furthermore, the active cathodic protection system of the pipeline causes some bias to the measurements too. A detailed description of the GIC phenomenon and its modelling is given by Antti Pulkkinen's PhD thesis, whose introduction is publicly available as a PDF file on the e-thesis server of the University of Helsinki.

Example

The starting point is the measured geomagnetic field at Nurmijärvi (BX = north component, BY = east component).

Because the electric field is related to the temporal changes ofnthe magnetic field, the time derivative of the magnetic field is an important quantity. Note especially the rapid variations of the time derivative.

Correspondingly, the electric field behaves in quite a similar manner to the time derivative of the magnetic field.

The same characteristics is finally visible in GIC. At Mäntsälä, GIC along the pipe usually resembles the east component of the electric field (EY).

 

Comparison of modelled and measured values

The nowcasting method is applicable also in post-analysis. We compare here modelled and measured GIC at Mäntsälä. We have used one-minute values in this investigation. The earth consists of two layers: thickness of the upper layer is 150 km and resistivity 38.5 ohmm; the lower layer is infinitely deep and its resistivity is 0.4 ohmm.

Error distribution in November 1998 to May 2003:

Examples of single events

19990113
19990218
19990301
19990417
19990913 - 19990922 - 19990923 - 19990927
19991022
19991107 - 19991113
20000212
20000406 - 20000407
20000517 - 20000524
20000608
20000713 - 20000715 - 20000716
20000811 - 20000812 - 20000813
20000917 - 20000918 - 20000919
20001004 - 20001005 - 20001014
20001106 - 20001107 - 20001110 - 20001127 - 20001129
20010124
20010319 - 20010320 - 20010328 - 20010331
20010401 - 20010408 - 20010411 - 20010412 - 20010413 - 20010418 - 20010428
20010509 - 20010512 - 20010513
20010817
20010923 - 20010925 - 20010926
20011002 - 20011003 - 20011012 - 20011021 - 20011022
20011106 - 20011124
20020417 - 20020418 - 20020419 - 20020420
20020523
20020820
20020907 - 20020908
20021001 - 20021002 - 20021003 - 20021004 - 20021007 - 20021024
20021121
20030202
20030317
20030510 - 20030529 - 20030530