Parameter estimation for transformer modeling

Large Power transformers, an aging and vulnerable part of our energy infrastructure, are at choke points in the grid and are key to reliability and security. Damage or destruction due to vandalism, misoperation, or other unexpected events is of great concern, given replacement costs upward of $2M and lead time of 12 months. Transient overvoltages can cause great damage and there is much interest in improving computer simulation models to correctly predict and avoid the consequences. EMTP (the Electromagnetic Transients Program) has been developed for computer simulation of power system transients. Component models for most equipment have been developed and benchmarked. Power transformers would appear to be simple. However, due to their nonlinear and frequency-dependent behaviors, they can be one of the most complex system components to model. It is imperative that the applied models be appropriate for the range of frequencies and excitation levels that the system experiences. Thus, transformer modeling is not a mature field and newer improved models must be made available. In this work, improved topologically-correct duality-based models are developed for three-phase autotransformers having five-legged, three-legged, and shell-form cores. The main problem in the implementation of detailed models is the lack of complete and reliable data, as no international standard suggests how to measure and calculate parameters. Therefore, parameter estimation methods are developed here to determine the parameters of a given model in cases where available information is incomplete. The transformer nameplate data is required and relative physical dimensions of the core are estimated. The models include a separate representation of each segment of the core, including hysteresis of the core, λ-i saturation characteristic, capacitive effects, and frequency dependency of winding resistance and core loss. Steady-state excitation, and de-energization and re-energization transients are simulated and compared with an earlier-developed BCTRAN-based model. Black start energization cases are also simulated as a means of model evaluation and compared with actual event records. The simulated results using the model developed here are reasonable and more correct than those of the BCTRAN-based model. Simulation accuracy is dependent on the accuracy of the equipment model and its parameters. This work is significant in that it advances existing parameter estimation methods in cases where the available data and measurements are incomplete. The accuracy of EMTP simulation for power systems including three-phase autotransformers is thus enhanced. Theoretical results obtained from this work provide a sound foundation for development of transformer parameter estimation methods using engineering optimization. In addition, it should be possible to refine which information and measurement data are necessary for complete duality-based transformer models. To further refine and develop the models and transformer parameter estimation methods developed here, iterative full-scale laboratory tests using high-voltage and high-power three-phase transformer would be helpful.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

Gotland Deep (Baltic Sea) independent age model based on Palaeomagnetic Secular Variation (PSV) and lead (Pb) deposition profiles

Dating of sediment cores from the Baltic Sea has proven to be difficult due to uncertainties surrounding the 14C reservoir age and a scarcity of macrofossils suitable for dating. Here we present the results of multiple dating methods carried out on cores in the Gotland Deep area of the Baltic Sea. Particular emphasis is placed on the Littorina stage (8 ka ago to the present) of the Baltic Sea and possible changes in the 14C reservoir age of our dated samples. Three geochronological methods are used. Firstly, palaeomagnetic secular variations (PSV) are reconstructed, whereby ages are transferred to PSV features through comparison with varved lake sediment based PSV records. Secondly, lead (Pb) content and stable isotope analysis are used to identify past peaks in anthropogenic atmospheric Pb pollution. Lastly, 14C determinations were carried out on benthic foraminifera (Elphidium spec.) samples from the brackish Littorina stage of the Baltic Sea. Determinations carried out on smaller samples (as low as 4 µg C) employed an experimental, state-of-the-art method involving the direct measurement of CO2 from samples by a gas ion source without the need for a graphitisation step - the first time this method has been performed on foraminifera in an applied study. The PSV chronology, based on the uppermost Littorina stage sediments, produced ten age constraints between 6.29 and 1.29 cal ka BP, and the Pb depositional analysis produced two age constraints associated with the Medieval pollution peak. Analysis of PSV data shows that adequate directional data can be derived from both the present Littorina saline phase muds and Baltic Ice Lake stage varved glacial sediments. Ferrimagnetic iron sulphides, most likely authigenic greigite (Fe3S4), present in the intermediate Ancylus Lake freshwater stage sediments acquire a gyroremanent magnetisation during static alternating field (AF) demagnetisation, preventing the identification of a primary natural remanent magnetisation for these sediments. An inferred marine reservoir age offset (deltaR) is calculated by comparing the foraminifera 14C determinations to a PSV & Pb age model. This deltaR is found to trend towards younger values upwards in the core, possibly due to a gradual change in hydrographic conditions brought about by a reduction in marine water exchange from the open sea due to continued isostatic rebound.

Acute survivorship of the deep-sea coral Lophelia pertusa from the Gulf of Mexico under acidification,warming,and deoxygenation

Changing global climate due to anthropogenic emissions of CO2 are driving rapid changes in the physical and chemical environment of the oceans via warming, deoxygenation, and acidification. These changes may threaten the persistence of species and populations across a range of latitudes and depths, including species that support diverse biological communities that in turn provide ecological stability and support commercial interests. Worldwide, but particularly in the North Atlantic and deep Gulf of Mexico, Lophelia pertusa forms expansive reefs that support biological communities whose diversity rivals that of tropical coral reefs. In this study, L. pertusa colonies were collected from the Viosca Knoll region in the Gulf of Mexico (390 to 450 m depth), genotyped using microsatellite markers, and exposed to a series of treatments testing survivorship responses to acidification, warming, and deoxygenation. All coral nubbins survived the acidification scenarios tested, between pH of 7.67 and 7.90 and aragonite saturation states of 0.92 and 1.47. However, calcification generally declined with respect to pH, though a disparate response was evident where select individuals net calcified and others exhibited net dissolution near a saturation state of 1. Warming and deoxygenation both had negative effects on survivorship, with up to 100% mortality observed at temperatures above 14ºC and oxygen concentrations of approximately 1.5 ml·l-1. These results suggest that, over the short-term, climate change and OA may negatively impact L. pertusa in the Gulf of Mexico, though the potential for acclimation and the effects of genetic background should be considered in future research.

Bulk rock magnetic measurements and thermomagnetic analyses of sediment core GeoB4901-8 from the southeastern flank of the Niger deep-sea fan

Diagenesis has extensively affected the magnetic mineral inventory of organic-rich late Quaternary sediments in the Niger deep-sea fan. Changes in concentration, grain size, and coercivity document modifications of the primary magnetic mineral assemblages at two horizons. The first front, the modern iron redox boundary, is characterized by a drastic decline in magnetic mineral content, coarsening of the grain size spectrum, and reduction in coercivity. Beneath a second front, the transition from the suboxic to the sulfidic anoxic domain, a further but less pronounced decrease in concentration and bulk grain size occurs. Finer grains and higher coercive magnetic constituents substantially increase in the anoxic environment. Low- and high-temperature experiments were performed on bulk sediments and on extracts which have also been examined by X-ray diffraction. Thermomagnetic analyses proved ferrimagnetic titanomagnetites of terrigenous provenance as the principal primary magnetic mineral components. Their broad range of titanium contents reflects the volcanogenic traits of the Niger River drainage areas. Diagenetic alteration is not only a grain size selective process but also critically depends on titanomagnetite composition. Low-titanium compounds are less resistant to diagenetic dissolution. Intermediate titanium content titanomagnetite thus persists as the predominant magnetic mineral fraction in the sulfidic anoxic sediments. At the Fe redox boundary, precipitation of authigenic, possibly bacterial, magnetite is documented. The presence of hydrogen sulfide in the pore water suggests a formation of secondary magnetic iron sulfides in the anoxic domain. Grain size-specific data argue for a gradual development of a superparamagnetic and single-domain iron sulfide phase in this milieu, most likely greigite.

Low-temperature magnetic remanence and hysteresis measurements on sediment core GeoB4901-8 from the Niger deep-sea fan

Low-temperature (LT) magnetic remanence and hysteresis measurements, in the range 300-5 K, were combined with energy dispersive spectroscopy (EDS) in order to characterize the magnetic inventory of strongly diagenetically altered sediments originating from the Niger deep-sea fan. We demonstrate the possibility of distinguishing between different compositions of members of the magnetite-ulvöspinel and ilmenite-hematite solid solution series on a set of five representative samples, two from the upper suboxic and three from the lower sulfidic anoxic zone of gravity core GeoB 4901. Highly sensitive LT magnetic measurements were performed on magnetic extracts resulting in large differences in the magnetic behavior between samples from the different layers. This emphasizes that both Fe-Ti oxide phases occur in different proportions in the two geochemical environments. Most prominent are variations in the coercivity sensitive parameter coercive field (BC). At room-temperature (RT) hysteresis loops for all extracts are narrow and yield low coercivities (6-13 mT). With decreasing temperature the loops become more pronounced and wider. At 5 K an approximately 5-fold increase in BC for the suboxic samples contrasts a 20-25-fold increase for the samples from the anoxic zone. We demonstrate that this distinct increase in BC at LT corresponds to the increasing proportion of the Ti-rich hemoilmenite phase, while Fe-rich (titano-)magnetite dominates the magnetic signal at RT. This trend is also seen in the room-temperature saturation isothermal remanent magnetization (RT-SIRM) cycles: suboxic samples show remanence curves dominated by Fe-rich mineral phases while anoxic samples display curves clearly dominated by Ti-rich particles. We show that the EDS intensity ratios of the characteristic Fe Kalpha and Ti Kalpha lines of the Fe-Ti oxides may be used to differentiate between members of the magnetite-ulvöspinel and ilmenite-hematite solid solution series. Furthermore it is possible to calculate an approximate composition for each grain if the intensity ratios of natural particles are linked to well-known standards. Thus, element spectra with high Fe/Ti intensity ratios were found to be rather typical of titanomagnetite while low Fe/Ti ratios are indicative of hemoilmenite. The EDS analyses confirm the LT magnetic results, Fe-rich magnetic phases dominate in the upper suboxic environment whereas Ti-rich magnetic phases comprise the majority of particles in the lower anoxic domain: The mineral assemblage of the upper suboxic environments is composed of magnetite (~19%), titanomagnetite (~62%), hemoilmenite (~17%) and ~2% other particles. In the lower anoxic sediments, reductive diagenetic alteration has resulted in more extensive depletion of the (titano-)magnetite phase, resulting in a relative enrichment of the hemoilmenite phase (~66%). In these strongly anoxic sediments stoichiometric magnetite is barely preserved and only ~5% titanomagnetite was detected. The remaining ~28% comprises Ti-rich particles such as pseudobrookite or rutile.

CO2-driven decrease in pH disrupts olfactory behaviour and increases individual variation in deep-sea hermit crabs

Deep-sea species are generally thought to be less tolerant of environmental variation than shallow-living species due to the relatively stable conditions in deep waters for most parameters (e.g. temperature, salinity, oxygen, and pH). To explore the potential for deep-sea hermit crabs (Pagurus tanneri) to acclimate to future ocean acidification, we compared their olfactory and metabolic performance under ambient (pH 7.6) and expected future (pH 7.1) conditions. After exposure to reduced pH waters, metabolic rates of hermit crabs increased transiently and olfactory behaviour was impaired, including antennular flicking and prey detection. Crabs exposed to low pH treatments exhibited higher individual variation for both the speed of antennular flicking and speed of prey detection, than observed in the control pH treatment, suggesting that phenotypic diversity could promote adaptation to future ocean acidification.

Sr/Ca and Cd/Ca ratios of benthic foraminifera in surface deep-sea sediments

Measurements of Sr/Ca of benthic foraminifera show a linear decrease with water depth which is superimposed upon significant variability identified by analyses of individual foraminifera. New data for Cd/Ca support previous work in defining a contrast between waters shallower and deeper than ~2500 m. Measured element partition coefficients in foraminiferal calcium carbonate relative to sea water (D) have been described by means of a one-box model in which elements are extracted by Rayleigh distillation from a biomineralization reservoir that serves for calcification with a constant fractionation factor (alpha), such that D = (1 - f**alpha)/(l - f), where f is the proportion of Ca remaining after precipitation. A modification to the model recognises differences in element speciation. The model is consistent with differences between D[Sr], D[Ba], and D[Cd] in benthic but not planktonic foraminifera. Depth variations in D for Sr and Ba are consistent with the model, as are differences in depth variation of D[Cd] in calcitic and aragonitic benthic foraminifera. The shallower depth variations may reflect increasing calcification rates with increasing water depth to an optimum of about 2500 m. Observations of unusually lower DCd for some deep waters, not accompanied by similar [Sr], or D[Ba] may be because of dissolution or a calcification response to a lower carbonate saturation state.

In search of a hidden long-term isolated sub-chondritic 142Nd/144Nd reservoir in the deep mantle : implications for the Nd isotope systematics of the Earth

Here we search for evidence of the existence of a sub-chondritic 142Nd/144Nd reservoir that balances the Nd isotope chemistry of the Earth relative to chondrites. If present, it may reside in the source region of deeply sourced mantle plume material. We suggest that lavas from Hawai’i with coupled elevations in 186Os/188Os and 187Os/188Os, from Iceland that represent mixing of upper mantle and lower mantle components, and from Gough with sub-chondritic 143Nd/144Nd and high 207Pb/206Pb, are favorable samples that could reflect mantle sources that have interacted with an Early-Enriched Reservoir (EER) with sub-chondritic 142Nd/144Nd. High-precision Nd isotope analyses of basalts from Hawai’i, Iceland and Gough demonstrate no discernable 142Nd/144Nd deviation from terrestrial standards. These data are consistent with previous high-precision Nd isotope analysis of recent mantle-derived samples and demonstrate that no mantle-derived material to date provides evidence for the existence of an EER in the mantle. We then evaluate mass balance in the Earth with respect to both 142Nd/144Nd and 143Nd/144Nd. The Nd isotope systematics of EERs are modeled for different sizes and timing of formation relative to ε143Nd estimates of the reservoirs in the μ142Nd = 0 Earth, where μ142Nd is ((measured 142Nd/144Nd/terrestrial standard 142Nd/144Nd)−1 * 10−6) and the μ142Nd = 0 Earth is the proportion of the silicate Earth with 142Nd/144Nd indistinguishable from the terrestrial standard. The models indicate that it is not possible to balance the Earth with respect to both 142Nd/144Nd and 143Nd/144Nd unless the μ142Nd = 0 Earth has a ε143Nd within error of the present-day Depleted Mid-ocean ridge basalt Mantle source (DMM). The 4567 Myr age 142Nd–143Nd isochron for the Earth intersects μ142Nd = 0 at ε143Nd of +8 ± 2 providing a minimum ε143Nd for the μ142Nd = 0 Earth. The high ε143Nd of the μ142Nd = 0 Earth is confirmed by the Nd isotope systematics of Archean mantle-derived rocks that consistently have positive ε143Nd. If the EER formed early after solar system formation (0–70 Ma) continental crust and DMM can be complementary reservoirs with respect to Nd isotopes, with no requirement for significant additional reservoirs. If the EER formed after 70 Ma then the μ142Nd = 0 Earth must have a bulk ε143Nd more radiogenic than DMM and additional high ε143Nd material is required to balance the Nd isotope systematics of the Earth.