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. ExampleThe 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 valuesThe 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 events1999011319990218 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
|
|