Estimação de Fasores para Proteção de Sistemas Elétricos Baseada em Mínimos Quadrados e Morfologia Matemática

This work proposes a new technique for phasor estimation applied in microprocessor numerical relays for distance protection of transmission lines, based on the recursive least squares method and called least squares modified random walking. The phasor estimation methods have compromised their performance, mainly due to the DC exponential decaying component present in fault currents. In order to reduce the influence of the DC component, a Morphological Filter (FM) was added to the method of least squares and previously applied to the process of phasor estimation. The presented method is implemented in MATLABr and its performance is compared to one-cycle Fourier technique and conventional phasor estimation, which was also based on least squares algorithm. The methods based on least squares technique used for comparison with the proposed method were: forgetting factor recursive, covariance resetting and random walking. The techniques performance analysis were carried out by means of signals synthetic and signals provided of simulations on the Alternative Transient Program (ATP). When compared to other phasor estimation methods, the proposed method showed satisfactory results, when it comes to the estimation speed, the steady state oscillation and the overshoot. Then, the presented method performance was analyzed by means of variations in the fault parameters (resistance, distance, angle of incidence and type of fault). Through this study, the results did not showed significant variations in method performance. Besides, the apparent impedance trajectory and estimated distance of the fault were analysed, and the presented method showed better results in comparison to one-cycle Fourier algorithm

Avaliação e modelagem da absorção de H2S do gás natural em coluna de leito estagnado

Natural gas, although basically composed by light hydrocarbons, also presents contaminant gases in its composition, such as CO2 (carbon dioxide) and H2S (hydrogen sulfide). The H2S, which commonly occurs in oil and gas exploration and production activities, causes damages in oil and natural gas pipelines. Consequently, the removal of hydrogen sulfide gas will result in an important reduction in operating costs. Also, it is essential to consider the better quality of the oil to be processed in the refinery, thus resulting in benefits in economic, environmental and social areas. All this facts demonstrate the need for the development and improvement in hydrogen sulfide scavengers. Currently, the oil industry uses several processes for hydrogen sulfide removal from natural gas. However, these processes produce amine derivatives which can cause damage in distillation towers, can cause clogging of pipelines by formation of insoluble precipitates, and also produce residues with great environmental impact. Therefore, it is of great importance the obtaining of a stable system, in inorganic or organic reaction media, able to remove hydrogen sulfide without formation of by-products that can affect the quality and cost of natural gas processing, transport, and distribution steps. Seeking the study, evaluation and modeling of mass transfer and kinetics of hydrogen removal, in this study it was used an absorption column packed with Raschig rings, where the natural gas, with H2S as contaminant, passed through an aqueous solution of inorganic compounds as stagnant liquid, being this contaminant gas absorbed by the liquid phase. This absorption column was coupled with a H2S detection system, with interface with a computer. The data and the model equations were solved by the least squares method, modified by Levemberg-Marquardt. In this study, in addition to the water, it were used the following solutions: sodium hydroxide, potassium permanganate, ferric chloride, copper sulfate, zinc chloride, potassium chromate, and manganese sulfate, all at low concentrations (»10 ppm). These solutions were used looking for the evaluation of the interference between absorption physical and chemical parameters, or even to get a better mass transfer coefficient, as in mixing reactors and absorption columns operating in counterflow. In this context, the evaluation of H2S removal arises as a valuable procedure for the treatment of natural gas and destination of process by-products. The study of the obtained absorption curves makes possible to determine the mass transfer predominant stage in the involved processes, the mass transfer volumetric coefficients, and the equilibrium concentrations. It was also performed a kinetic study. The obtained results showed that the H2S removal kinetics is greater for NaOH. Considering that the study was performed at low concentrations of chemical reagents, it was possible to check the effect of secondary reactions in the other chemicals, especially in the case of KMnO4, which shows that your by-product, MnO2, acts in H2S absorption process. In addition, CuSO4 and FeCl3 also demonstrated to have good efficiency in H2S removal

Evolução da biodegradabilidade da matéria orgânica em um sistema de lagoas de estabilização

Waste stabilization ponds (WSP) have been widely used for sewage treatment in hot climate regions because they are economic and environmentally sustainable. In the present study a WSP complex comprising a primary facultative pond (PFP) followed by two maturation ponds (MP-1 and MP-2) was studied, in the city of Natal-RN. The main objective was to study the bio-degradability of organic matter through the determination of the kinetic constant k throughout the system. The work was carried out in two phases. In the first, the variability in BOD, COD and TOC concentrations and an analysis of the relations between these parameters, in the influent raw sewage, pond effluents and in specific areas inside the ponds was studied. In the second stage, the decay rate for organic matter (k) was determined throughout the system based on BOD tests on the influent sewage, pond effluents and water column samples taken from fixed locations within the ponds, using the mathematical methods of Least Squares and the Thomas equation. Subsequently k was estimated as a function of a hydrodynamic model determined from the dispersion number (d), using empirical methods and a Partial Hydrodynamic Evaluation (PHE), obtained from tracer studies in a section of the primary facultative pond corresponding to 10% of its total length. The concentrations of biodegradable organic matter, measured as BOD and COD, gradually reduced through the series of ponds, giving overall removal efficiencies of 71.95% for BOD and of 52.45% for COD. Determining the values for k, in the influent and effluent samples of the ponds using the mathematical method of Least Squares, gave the following values respectively: primary facultative pond (0,23 day-1 and 0,09 day-1), maturation 1 (0,04 day-1 and 0,03 day-1) and maturation 2 (0,03 day-1 and 0,08 day-1). When using the Thomas method, the values of k in the influents and effluents of the ponds were: primary facultative pond (0,17 day-1 and 0,07 day-1), maturation 1 (0,02 day-1 and 0,01 day-1) and maturation 2 (0,01 day-1 and 0,02 day-1). From the Partial Hydrodynamic Evaluation, in the first section of the facultative pond corresponding to 10% of its total length, it can be concluded from the dispersion number obtained of d = 0.04, that the hydraulic regime is one of dispersed flow with a kinetic constant value of 0.20 day-1

Integração de redes GNSS: uma proposta metodológica de densificação da rede SIRGAS na América do Sul

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Análise de séries temporais de coordenadas estimadas com GPS: uma proposta metodológica para detecção, remoção e recuperação de efeitos sazonais

GPS active networks are more and more used in geodetic surveying and scientific experiments, as water vapor monitoring in the atmosphere and lithosphere plate movement. Among the methods of GPS positioning, Precise Point Positioning (PPP) has provided very good results. A characteristic of PPP is related to the modeling and/or estimation of the errors involved in this method. The accuracy obtained for the coordinates can reach few millimeters. Seasonal effects can affect such accuracy if they are not consistent treated during the data processing. Coordinates time series analyses have been realized using Fourier or Harmonics spectral analyses, wavelets, least squares estimation among others. An approach is presented in this paper aiming to investigate the seasonal effects included in the stations coordinates time series. Experiments were carried out using data from stations Manaus (NAUS) and Fortaleza (BRFT) which belong to the Brazilian Continuous GPS Network (RBMC). The coordinates of these stations were estimated daily using PPP and were analyzed through wavelets for identification of the periods of the seasonal effects (annual and semi-annual) in each time series. These effects were removed by means of a filtering process applied in the series via the least squares adjustment (LSQ) of a periodic function. The results showed that the combination of these two mathematical tools, wavelets and LSQ, is an interesting and efficient technique for removal of seasonal effects in time series.

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Improving the GNSS positioning stochastic model in the presence of ionospheric scintillation

Ionospheric scintillations are caused by time-varying electron density irregularities in the ionosphere, occurring more often at equatorial and high latitudes. This paper focuses exclusively on experiments undertaken in Europe, at geographic latitudes between similar to 50 degrees N and similar to 80 degrees N, where a network of GPS receivers capable of monitoring Total Electron Content and ionospheric scintillation parameters was deployed. The widely used ionospheric scintillation indices S4 and sigma(phi) represent a practical measure of the intensity of amplitude and phase scintillation affecting GNSS receivers. However, they do not provide sufficient information regarding the actual tracking errors that degrade GNSS receiver performance. Suitable receiver tracking models, sensitive to ionospheric scintillation, allow the computation of the variance of the output error of the receiver PLL (Phase Locked Loop) and DLL (Delay Locked Loop), which expresses the quality of the range measurements used by the receiver to calculate user position. The ability of such models of incorporating phase and amplitude scintillation effects into the variance of these tracking errors underpins our proposed method of applying relative weights to measurements from different satellites. That gives the least squares stochastic model used for position computation a more realistic representation, vis-a-vis the otherwise 'equal weights' model. For pseudorange processing, relative weights were computed, so that a 'scintillation-mitigated' solution could be performed and compared to the (non-mitigated) 'equal weights' solution. An improvement between 17 and 38% in height accuracy was achieved when an epoch by epoch differential solution was computed over baselines ranging from 1 to 750 km. The method was then compared with alternative approaches that can be used to improve the least squares stochastic model such as weighting according to satellite elevation angle and by the inverse of the square of the standard deviation of the code/carrier divergence (sigma CCDiv). The influence of multipath effects on the proposed mitigation approach is also discussed. With the use of high rate scintillation data in addition to the scintillation indices a carrier phase based mitigated solution was also implemented and compared with the conventional solution. During a period of occurrence of high phase scintillation it was observed that problems related to ambiguity resolution can be reduced by the use of the proposed mitigated solution.

Orientação indireta de imagens digitais usando telhados de duas águas como apoio de campo

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Orbit determination using GPS including perturbations due to geopotential coefficients of high degree and order and solar radiation pressure

The main target here is to determine the orbit of an artificial satellite, using signals of the GPS constellation and least squares algorithms implemented through sequential Givens rotations as a method of estimation, with the aim of improving the performance of the orbit estimation process and, at the same time, minimizing the computational procedure cost. Geopotential perturbations up to high order and direct solar radiation pressure were taken into account. It was also considered the position of the GPS antenna on the satellite body that, lately, consists of the influence of the satellite attitude motion in the orbit determination process. An application has been done, using real data from the Topex/Poseidon satellite, whose ephemeris is available at Internet. The best accuracy obtained in position was smaller than 5 meters for short period (2 hours) and smaller than 28 meters for long period (24 hours) orbit determination. In both cases, the perturbations mentioned before were taken into consideration and the analysis occurred without selective availability on the signals measurements.

Ferramentas quimiométricas aplicadas a classificação de amostras de leite em pó e quantificação de proteínas

In this work we used chemometric tools to classify and quantify the protein content in samples of milk powder. We applied the NIR diffuse reflectance spectroscopy combined with multivariate techniques. First, we carried out an exploratory method of samples by principal component analysis (PCA), then the classification of independent modeling of class analogy (SIMCA). Thus it became possible to classify the samples that were grouped by similarities in their composition. Finally, the techniques of partial least squares regression (PLS) and principal components regression (PCR) allowed the quantification of protein content in samples of milk powder, compared with the Kjeldahl reference method. A total of 53 samples of milk powder sold in the metropolitan areas of Natal, Salvador and Rio de Janeiro were acquired for analysis, in which after pre-treatment data, there were four models, which were employed for classification and quantification of samples. The methods employed after being assessed and validated showed good performance, good accuracy and reliability of the results, showing that the NIR technique can be a non invasive technique, since it produces no waste and saves time in analyzing the samples

Métodos estatísticos recursivos aplicados ao problema de estimação de vazão de gás em plantas de Plunger-lift

This work has as main objective to find mathematical models based on linear parametric estimation techniques applied to the problem of calculating the grow of gas in oil wells. In particular we focus on achieving grow models applied to the case of wells that produce by plunger-lift technique on oil rigs, in which case, there are high peaks in the grow values that hinder their direct measurement by instruments. For this, we have developed estimators based on recursive least squares and make an analysis of statistical measures such as autocorrelation, cross-correlation, variogram and the cumulative periodogram, which are calculated recursively as data are obtained in real time from the plant in operation; the values obtained for these measures tell us how accurate the used model is and how it can be changed to better fit the measured values. The models have been tested in a pilot plant which emulates the process gas production in oil wells

Estudo sísmico no açude Tucunduba, Senador Sá, CE

The Tucunduba Dam, is west of Fortaleza, Ceará State. The seismic monitoring of the area, with an analogical station and seven digital stations, had beginning on June 11, 1997. The digital stations, operated from June to November 1997. The data collected in the period of digital monitoring was analyzed for determination of hypocenters, focal mechanisms, and shear-wave anisotropy analysis. For determination of hypocenters, it was possible to find an active zone of nearly 1 km in length, with depth between 4.5 and 5.2 km. A 60AZ/88SE fault plane was determined using the least-squares method and hypocenters of a selected set of 16 earthquakes recorded. Focal mechanisms were determined, in the composite fault plane solution, a strike-slip fault, trending nearly E-W, was found. Single fault plane solutions were obteined to some earthquakes presented mean values of 65 (azimuth), and 80 (dip). Shear-wave anisotropy was found in the data. Polarization directions and travel time delays, between S spliting waves, were determined. It was not possible to obtain any conclusion on the cause of the observed anisotropy. It is not clear if there is correlation between seismicity and mapped faults in the area, although the directions obtained starting from the hipocentros and focal mechanism are they are consistent with directions, observed in the area, photo, topographic and fractures directions observed in the area

Effects of GHR gene polymorphisms on growth and carcass traits in Zebu and crossbred beef cattle

The growth hormone receptor (GHR) is the cell surface receptor for growth hormone (GH) and is required for GH to carry out its effects on target tissues. The objectives of the present study were to estimate the allele and genotype frequencies of the GHR/Alu I gene polymorphism located in the regulatory region in beef cattle belonging to different genetic groups and to determine associations between this polymorphism and growth and carcass traits. Genotyping was performed on 384 animals, including 79 Nellore (Zebu), 30 Canchim (5/8 Charolais+3/8 Zebu), 30 Simmental X Nellore crossbred and 245 Angus x Nellore crossbred cattle. Alleles Alu I(+), Alu I(-) and Alu I(N)-null allele-were evidenced for the GHR/Alu I polymorphism and the frequency of the Alu I(N) allele was significantly higher than the frequency of the Alu I(+) and Alu I(-) alleles in all genetic groups. Genotype Alu I(N/N) of the GHRIAlu I predominated in Nellore animals, while the Alu I(N/+) and Alu I(N/-) predominated in the other genetic groups. In the association studies, traits of interest were analyzed using the General Linear Model (GLM) procedure of the SAS program and least squares means of the genotypes were compared by the Tukey test. Significant associations (P < 0.05) were observed between the Alu I(N/N) genotype of the GHRIAlu I polymorphism and lower weight gain and body weight at slaughter, although a confounding between genotypes and genetic groups may have occurred. (c) 2005 Elsevier B.V. All rights reserved.

Parâmetros e tendência genética da produção de leite de cabra no Brasil

A caprinocultura leiteira no Brasil, apesar de ser uma atividade rural consolidada há algumas décadas, tem se mostrado totalmente dependente de outros países no que se refere ao melhoramento genético. A maioria dos plantéis existentes atualmente tem como base animais importados, e a renovação do material genético é feita por meio da importação de sêmen. Inexistem informações sobre o valor genético dos animais e sua evolução no decorrer dos anos. No presente trabalho, foram estimadas a herdabilidade e a repetibilidade da produção de leite utilizando o REML. Os valores obtidos foram 0,21557 e 0,21564, respectivamente. Para a predição do valor gênico dos animais, foi usado o procedimento BLUP com modelo animal. A mudança na tendência genética anual estimada por um modelo quadrático foi -0,8109 kg/ano², indicando desaceleração no ganho genético. A correlação de Pearson entre os valores gênicos dos bodes estimados com base na média da capacidade provável de produção das filhas obtida pelo método de mínimos quadrados com as estimadas pelas equações do modelo misto foi de 0,5751. A correlação de SPEARMAN entre as classificações dos bodes obtidos pelos dois métodos foi de 0,5813.

Fatores genéticos e de meio na produção de leite de caprinos leiteiros

Objetivou-se, no presente estudo, avaliar a produção de leite de caprinos leiteiros da região sudeste do Brasil, com intuito de verificar os fatores de meio e estimar os parâmetros genéticos pelo método dos mínimos quadrados (MMQ). Os controles de 1336 lactações foram inicialmente ajustados pela função multifásica (difásica) e calculou-se a produção de leite total (PLT). Os dados foram provenientes de sete propriedades e três raças (Parda Alpina, Saanen e Toggenburg). A média e o erro-padrão da PLT estimados pelo MMQ foram de 635,31 ±39,75 kg. A interação ano x estação do parto influenciou a PLT. em um dos anos estudados, a PLT foi menor para as cabras paridas no final da estação. Nas três estações de parto, observou-se comportamento quadrático da PLT, em função dos anos de parto. Para as três estações, a PLT aumentou de 1986 até meados de 1990, decrescendo em seguida. A idade de máxima PLT foi observada aos 46,65 meses. Das três raças estudadas, observou-se que as raças Parda Alpina e Saanen apresentaram alternância de superioridade na PLT em algumas fazendas, porém maiores que a Toggenburg. Os coeficientes de herdabilidade e repetibilidade da PLT estimados pelo MMQ foram de 0,296 ± 0,079 e 0,277 ± 0,033, respectivamente. Estes resultados revelam baixa confiabilidade em poucas observações dessas características ou na inconsistência das estimativas da função multifásica. Como a PLT é uma característica limitada ao sexo, sugere-se o teste de progênie como método de seleção mais eficiente para os reprodutores e uso de inseminação artificial como processo de disseminação do material genético selecionado.