Analysis of distrubances in large interconnected power systems

World economies increasingly demand reliable and economical power supply and distribution. To achieve this aim the majority of power systems are becoming interconnected, with several power utilities supplying the one large network. One problem that occurs in a large interconnected power system is the regular occurrence of system disturbances which can result in the creation of intra-area oscillating modes. These modes can be regarded as the transient responses of the power system to excitation, which are generally characterised as decaying sinusoids. For a power system operating ideally these transient responses would ideally would have a “ring-down” time of 10-15 seconds. Sometimes equipment failures disturb the ideal operation of power systems and oscillating modes with ring-down times greater than 15 seconds arise. The larger settling times associated with such “poorly damped” modes cause substantial power flows between generation nodes, resulting in significant physical stresses on the power distribution system. If these modes are not just poorly damped but “negatively damped”, catastrophic failures of the system can occur. To ensure system stability and security of large power systems, the potentially dangerous oscillating modes generated from disturbances (such as equipment failure) must be quickly identified. The power utility must then apply appropriate damping control strategies. In power system monitoring there exist two facets of critical interest. The first is the estimation of modal parameters for a power system in normal, stable, operation. The second is the rapid detection of any substantial changes to this normal, stable operation (because of equipment breakdown for example). Most work to date has concentrated on the first of these two facets, i.e. on modal parameter estimation. Numerous modal parameter estimation techniques have been proposed and implemented, but all have limitations [1-13]. One of the key limitations of all existing parameter estimation methods is the fact that they require very long data records to provide accurate parameter estimates. This is a particularly significant problem after a sudden detrimental change in damping. One simply cannot afford to wait long enough to collect the large amounts of data required for existing parameter estimators. Motivated by this gap in the current body of knowledge and practice, the research reported in this thesis focuses heavily on rapid detection of changes (i.e. on the second facet mentioned above). This thesis reports on a number of new algorithms which can rapidly flag whether or not there has been a detrimental change to a stable operating system. It will be seen that the new algorithms enable sudden modal changes to be detected within quite short time frames (typically about 1 minute), using data from power systems in normal operation. The new methods reported in this thesis are summarised below. The Energy Based Detector (EBD): The rationale for this method is that the modal disturbance energy is greater for lightly damped modes than it is for heavily damped modes (because the latter decay more rapidly). Sudden changes in modal energy, then, imply sudden changes in modal damping. Because the method relies on data from power systems in normal operation, the modal disturbances are random. Accordingly, the disturbance energy is modelled as a random process (with the parameters of the model being determined from the power system under consideration). A threshold is then set based on the statistical model. The energy method is very simple to implement and is computationally efficient. It is, however, only able to determine whether or not a sudden modal deterioration has occurred; it cannot identify which mode has deteriorated. For this reason the method is particularly well suited to smaller interconnected power systems that involve only a single mode. Optimal Individual Mode Detector (OIMD): As discussed in the previous paragraph, the energy detector can only determine whether or not a change has occurred; it cannot flag which mode is responsible for the deterioration. The OIMD seeks to address this shortcoming. It uses optimal detection theory to test for sudden changes in individual modes. In practice, one can have an OIMD operating for all modes within a system, so that changes in any of the modes can be detected. Like the energy detector, the OIMD is based on a statistical model and a subsequently derived threshold test. The Kalman Innovation Detector (KID): This detector is an alternative to the OIMD. Unlike the OIMD, however, it does not explicitly monitor individual modes. Rather it relies on a key property of a Kalman filter, namely that the Kalman innovation (the difference between the estimated and observed outputs) is white as long as the Kalman filter model is valid. A Kalman filter model is set to represent a particular power system. If some event in the power system (such as equipment failure) causes a sudden change to the power system, the Kalman model will no longer be valid and the innovation will no longer be white. Furthermore, if there is a detrimental system change, the innovation spectrum will display strong peaks in the spectrum at frequency locations associated with changes. Hence the innovation spectrum can be monitored to both set-off an “alarm” when a change occurs and to identify which modal frequency has given rise to the change. The threshold for alarming is based on the simple Chi-Squared PDF for a normalised white noise spectrum [14, 15]. While the method can identify the mode which has deteriorated, it does not necessarily indicate whether there has been a frequency or damping change. The PPM discussed next can monitor frequency changes and so can provide some discrimination in this regard. The Polynomial Phase Method (PPM): In [16] the cubic phase (CP) function was introduced as a tool for revealing frequency related spectral changes. This thesis extends the cubic phase function to a generalised class of polynomial phase functions which can reveal frequency related spectral changes in power systems. A statistical analysis of the technique is performed. When applied to power system analysis, the PPM can provide knowledge of sudden shifts in frequency through both the new frequency estimate and the polynomial phase coefficient information. This knowledge can be then cross-referenced with other detection methods to provide improved detection benchmarks.

Planar difference functions

In 1980 Alltop produced a family of cubic phase sequences that nearly meet the Welch bound for maximum non-peak correlation magnitude. This family of sequences were shown by Wooters and Fields to be useful for quantum state tomography. Alltop’s construction used a function that is not planar, but whose difference function is planar. In this paper we show that Alltop type functions cannot exist in fields of characteristic 3 and that for a known class of planar functions, x^3 is the only Alltop type function.

Phase composition of the rapidly quenched melt of YBa2Cu3O7-y+20 mol% Y2BaCuO5

Samples of YBa2Cu3O7-y+20 mol% Y2BaCuO5, with thicknesses ranging between 50-250 μm, have been melt processed and rapidly quenched from temperatures between 985 and 1100°C by immersing them in liquid nitrogen. The phase composition and microstructures of these samples have been characterised using a combination of X-ray diffractometry, optical microscopy and scanning electron microscopy with energy dispersive X-ray spectroscopy. The quenched melt of samples quenched from temperatures greater than 985°C appears relatively homogeneous but consists of Ba2Cu3Ox (BC1.5) and BaCu2O2 (BC2) regions. At about 985°C, BaCuO2 (BC1) crystallises from the melt and most of the BC1.5 decomposes into BC1 and CuO or into BC1 and BC2. The crystallisation of BC1 induces segregation of elements in the melt and this is very significant for the melt texturing of YBCO.

Sol-gel synthesis and characterization of cubic bismuth zinc niobium oxide nanopowders

Bismuth zinc niobium oxide (BZN) was successfully synthesized by a diol-based sol-gel reaction utilizing metal acetate and alkoxide precursors. Thermal analysis of a liquid suspension of precursors suggests that the majority of organic precursors decompose at temperatures up to 150°C, and organic free powders form above 350°C. The experimental results indicate that a homogeneous gel is obtained at about 200°C and then converts to a mixture of intermediate oxides at 350–400°C. Finally, single-phased BZN powders are obtained between 500 and 900°C. The degree of chemical homogeneity as determined by X-ray diffraction and EDS mapping is consistent throughout the samples. Elemental analysis indicates that the atomic ratio of metals closely matches a Bi1.5ZnNb1.5O7 composition. Crystallite sizes of the BZN powders calculated from the Scherrer equation are about 33–98 nm for the samples prepared at 500–700°C, respectively. The particle and crystallite sizes increase with increased sintering temperature. The estimated band gap of the BZN nanopowders from optical analysis is about 2.60–2.75 eV at 500-600°C. The observed phase formations and measured results in this study were compared with those of previous reports.

Influences of molecular structure on co-crystallization of polypyridyl metal complexes.

Heteroleptic complexes of the type \[RuL2L′](PF6)2 (L, L′ = combinations of 1,10-phenanthroline (phen) and 2,2′-bipyridine (bipy)) were found to cocrystallize with \[Ni(phen)3](PF6)2 to produce cocrystals of \[Ni(phen)3]x\[RuL2L′]1–x(PF6)2. In this report we show that the ability of the complexes to cocrystallize is influenced by the number of common ligands between complexes in solution. Supramolecular selection is a phenomenon caused by molecular recognition through which cocrystals can grow from the same solution but contain different ratios of the molecular components. It was found that systems where L = phen displayed less supramolecular selection than systems where L = bipy. With increasing supramolecular selection, the composition of cocrystals was found to vary significantly from the initial relative concentration in the cocrystallizing solution, and therefore it was increasingly difficult to control the final composition of the resultant cocrystals. Consequently, modulation of concentration-dependent properties such as phase was also found to be less predictable with increasing supramolecular selection. Notwithstanding the complication afforded by the presence of supramolecular selection, our results reaffirm the robustness of the \[M(phen)3](PF6)2 structure because it was maintained even when ca. 90% of the complexes in the cocrystals were \[Ru(phen)(bipy)2](PF6)2, which in its pure form is not isomorphous with \[M(phen)3](PF6)2. Experiments between complexes without common ligands, i.e., \[Ru(bipy)3](PF6)2 cocrystallized with \[Ni(phen)3](PF6)2, were found to approach the limit to which molecular recognition processes can be confused into cocrystallizing different molecules to form single cocrystals. For these systems the result was the formation of block-shaped crystals skewered by a needle-shaped crystals.

Inhibition by poly(acrylic acid) and morphological changes in calcium carbonate and calcium carbonate/calcium sulfate crystallization on silica fibers

An intrinsic exposed core optical fiber sensor (IECOFS) made from fused silica was used to monitor the crystallization of calcium carbonate (CaCO3) and CaCO3/calcium sulfate (CaSO4) composite at 100 and 120 °C in the absence and presence of low-molar-mass (Mn ≤ 2000) poly(acrylic acid) (PAA) with different end groups. The IECOFS responded only to deposition and growth processes on the fiber surface rather than changes occurring in the bulk of the solution. Hexyl isobutyrate-terminated PAA (Mn = 1400) and hexadecyl isobutyrate-terminated PAA (Mn = 1700) were the most effective species in preventing CaCO3 deposition. Phase transformation from vaterite to aragonite/calcite decreased with increasing hydrophobicity of the PAA end group. Low-molar-mass PAA at 10 ppm showed very significant inhibition of CaCO3/CaSO4 composite formation for all end groups investigated.

In situ investigation of explosive crystallization in a-Ge: Experimental determination of the interface response function using dynamic transmission electron microscopy

The crystallization of amorphous semiconductors is a strongly exothermic process. Once initiated the release of latent heat can be sufficient to drive a self-sustaining crystallization front through the material in a manner that has been described as explosive. Here, we perform a quantitative in situ study of explosive crystallization in amorphous germanium using dynamic transmission electron microscopy. Direct observations of the speed of the explosive crystallization front as it evolves along a laser-imprinted temperature gradient are used to experimentally determine the complete interface response function (i.e., the temperature-dependent front propagation speed) for this process, which reaches a peak of 16 m/s. Fitting to the Frenkel-Wilson kinetic law demonstrates that the diffusivity of the material locally/immediately in advance of the explosive crystallization front is inconsistent with those of a liquid phase. This result suggests a modification to the liquid-mediated mechanism commonly used to describe this process that replaces the phase change at the leading amorphous-liquid interface with a change in bonding character (from covalent to metallic) occurring in the hot amorphous material.

Synthesis, characterization, and electronic structure studies of cubic Bi1.5ZnTa1.5O7 for photocatalytic applications

Bi1.5ZnTa1.5O7 (BZT) has been synthesized using an alkoxide based sol-gel reaction route. The evolution of the phases produced from the alkoxide precursors and their properties have been characterized as function of temperature using a combination of thermogravimetric analysis (TGA) coupled with mass spectrometry (MS), infrared emission spectrometry (IES), X-ray diffraction (XRD), ultraviolet and visible (UV-Vis) spectroscopy, Raman spectroscopy, and N2 adsorption/desorption isotherms. The lowest sintering temperature (600∘C) to obtain phase pure BZT powders with high surface area (14.5m2/g) has been determined from the thermal decomposition and phase analyses.The photocatalytic activity of the BZT powders has been tested for the decolorization of organic azo-dye and found to be photoactive under UV irradiation.The electronic band structure of the BZT has been investigated using density functional theory (DFT) calculations to determine the band gap energy (3.12 eV) and to compare it with experimental band gap (3.02 eV at 800∘C) from optical absorptionmeasurements. An excellent match is obtained for an assumption of Zn cation substitutions at specifically ordered sites in the BZT structure.

Condensed phase combustion travelling waves with sequential exothermic or endothermic reactions

The one-dimensional propagation of a combustion wave through a premixed solid fuel for two-stage kinetics is studied. We re-examine the analysis of a single reaction travelling-wave and extend it to the case of two-stage reactions. We derive an expression for the travelling wave speed in the limit of large activation energy for both reactions. The analysis shows that when both reactions are exothermic, the wave structure is similar to the single reaction case. However, when the second reaction is endothermic, the wave structure can be significantly different from single reaction case. In particular, as might be expected, a travelling wave does not necessarily exist in this case. We establish conditions in the limiting large activation energy limit for the non-existence, and for monotonicity of the temperature profile in the travelling wave.