Anisotropic seismic inversion using a multigrid Monte Carlo approach

We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (epsilon, delta). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (epsilon, delta) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen epsilon parameter is particularly robust compared to that of slowness and the Thomsen delta parameter, even in the face of complex subsurface anomalies. The Thomsen epsilon and delta parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that epsilon and delta in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.

How Porosity and Permeability Vary Spatially With Grain Size, Sorting, Cement Volume, and Mineral Dissolution In Fluvial Triassic Sandstones: The Value of Geostatistics and Local Regression

Although it is well known that sandstone porosity and permeability are controlled by a range of parameters such as grain size and sorting, amount, type, and location of diagenetic cements, extent and type of compaction, and the generation of intergranular and intragranular secondary porosity, it is less constrained how these controlling parameters link up in rock volumes (within and between beds) and how they spatially interact to determine porosity and permeability. To address these unknowns, this study examined Triassic fluvial sandstone outcrops from the UK using field logging, probe permeametry of 200 points, and sampling at 100 points on a gridded rock surface. These field observations were supplemented by laser particle-size analysis, thin-section point-count analysis of primary and diagenetic mineralogy, quantitiative XRD mineral analysis, and SEM/EDAX analysis of all 100 samples. These data were analyzed using global regression, variography, kriging, conditional simulation, and geographically weighted regression to examine the spatial relationships between porosity and permeability and their potential controls. The results of bivariate analysis (global regression) of the entire outcrop dataset indicate only a weak correlation between both permeability porosity and their diagenetic and depositional controls and provide very limited information on the role of primary textural structures such as grain size and sorting. Subdividing the dataset further by bedding unit revealed details of more local controls on porosity and permeability. An alternative geostatistical approach combined with a local modelling technique (geographically weighted regression; GWR) subsequently was used to examine the spatial variability of porosity and permeability and their controls. The use of GWR does not require prior knowledge of divisions between bedding units, but the results from GWR broadly concur with results of regression analysis by bedding unit and provide much greater clarity of how porosity and permeability and their controls vary laterally and vertically. The close relationship between depositional lithofacies in each bed, diagenesis, and permeability, porosity demonstrates that each influences the other, and in turn how understanding of reservoir properties is enhanced by integration of paleoenvironmental reconstruction, stratigraphy, mineralogy, and geostatistics.

Weak interactions, omnivory and emergent food-web properties

Empirical studies have shown that, in real ecosystems, species-interaction strengths are generally skewed in their distribution towards weak interactions. Some theoretical work also suggests that weak interactions, especially in omnivorous links, are important for the local stability of a community at equilibrium. However, the majority of theoretical studies use uniform distributions of interaction strengths to generate artificial communities for study. We investigate the effects of the underlying interaction-strength distribution upon the return time, permanence and feasibility of simple Lotka-Volterra equilibrium communities. We show that a skew towards weak interactions promotes local and global stability only when omnivory is present. It is found that skewed interaction strengths are an emergent property of stable omnivorous communities, and that this skew towards weak interactions creates a dynamic constraint maintaining omnivory. Omnivory is more likely to occur when omnivorous interactions are skewed towards weak interactions. However, a skew towards weak interactions increases the return time to equilibrium, delays the recovery of ecosystems and hence decreases the stability of a community. When no skew is imposed, the set of stable omnivorous communities shows an emergent distribution of skewed interaction strengths. Our results apply to both local and global concepts of stability and are robust to the definition of a feasible community. These results are discussed in the light of empirical data and other theoretical studies, in conjunction with their broader implications for community assembly.

Fluid Shear Induces Conformation Change in Human Blood Protein von Willebrand Factor in Solution

Many of the physiological functions of von Willebrand Factor (VWF), including its binding interaction with blood platelets, are regulated by the magnitude of applied fluid/hydrodynamic stress. We applied two complementary strategies to study the effect of fluid forces on the solution structure of VWF. First, small-angle neutron scattering was used to measure protein conformation changes in response to laminar shear rates (G) up to 3000/s. Here, purified VWF was sheared in a quartz Couette cell and protein conformation was measured in real time over length scales from 2-140 nm. Second, changes in VWF structure up to 9600/s were quantified by measuring the binding of a fluorescent probe 1,1'-bis(anilino)-4-,4'-bis(naphtalene)-8,8'-disulfonate (bis-ANS) to hydrophobic pockets exposed in the sheared protein. Small angle neutron scattering studies, coupled with quantitative modeling, showed that VWF undergoes structural changes at G < 3000/s. These changes were most prominent at length scales <10 nm (scattering vector (q) range >0.6/nm). A mathematical model attributes these changes to the rearrangement of domain level features within the globular section of the protein. Studies with bis-ANS demonstrated marked increase in bis-ANS binding at G > 2300/s. Together, the data suggest that local rearrangements at the domain level may precede changes at larger-length scales that accompany exposure of protein hydrophobic pockets. Changes in VWF conformation reported here likely regulate protein function in response to fluid shear.

Modified Lipoproteins in Diabetic Retinopathy: A Local Action in the Retina

Clinical epidemiological studies have revealed relatively weak, yet statistically significant, associations between dyslipidemia/dyslipoproteinemia and diabetic retinopathy (DR). Recent large interventional studies, however, demonstrated an unexpectedly robust efficacy of fenofibrate on the development of DR, possibly independent of plasma lipids. To unify the apparent discrepancies, we hypothesize that plasma lipoproteins play an indirect but important role in DR, contingent on the integrity of the blood-retina-barrier (BRB). In retinas with an intact BRB, plasma lipoproteins may be largely irrelevant; however, important effects become operative after the BRB is impaired in diabetes, leading to lipoprotein extravasation and subsequent modification, hence toxicity to the neighbouring retinal cells. In this hypothesis, BRB leakage is the key, plasma lipoprotein concentrations mainly modulate its consequences, and fenofibrate has intra-retinal actions. This review summarizes our current knowledge of the direct effects and mechanisms of modified lipoproteins on retinal cells and their potential contribution to the pathogenesis of DR.

Effect of Post-weld Heat-treatment on the Stress Corrosion Cracking Behaviour of Laser-welded Shape Memory NiTi Wires in Hanks’ Solution

In this study, the stress-corrosion cracking (SCC) behaviour of laser-welded NiTi wires before and after post-weld heat-treatment (PWHT) was investigated. The samples were subjected to slow strain rate testing (SSRT) under tensile loading in Hanks’ solution at 37.5 °C (or 310.5 K) at a constant anodic potential (200 mVSCE). The current density of the samples during the SSRT was captured by a potentiostat, and used as an indicator to determine the susceptibility to SCC. Fractography was analyzed using scanning-electron microscopy (SEM). The experimental results showed that the laser-welded sample after PWHT was immune to the SCC as evidenced by the stable current density throughout the SSRT. This is attributed to the precipitation of fine and coherent nano-sized Ni4Ti3 precipitates in the welded regions (weld zone, WZ and heat-affected zone, HAZ) after PWHT, resulting in (i) enrichment of TiO2 content in the passive film and (ii) higher resistance against the local plastic deformation in the welded regions.

Minimal intervention dentistry and older patients Part 1: Risk assessment and caries prevention

Ten million people in the UK today are aged over 65. The latest projections estimate that there will be 5 1/2 million more people aged 65 and older in the next 20 years. This projected pattern of population ageing will have profound consequences for dentistry. Minimal intervention dentistry (MID) is a modern evidence-based approach to caries management in dentate patients that uses the 'medical model' whereby disease is controlled by the 'oral physician'. This approach offers considerable benefits over conventional dentistry for older patients. It encourages patients to be responsible for their oral health through the provision of both knowledge and motivation. MID encompasses risk assessment for dental disease, early detection and control of disease processes, and minimally invasive treatment. 

Clinical Relevance: Risk assessment tools can aid the general dental practitioner and the patient to develop a suitable caries prevention programme for that individual and reduce the need for future operative intervention.

Partial Palmprint Matching Using Invariant Local Minutiae Descriptors

In forensic investigations, it is common for forensic investigators to obtain a photograph of evidence left at the scene of crimes to aid them catch the culprit(s). Although, fingerprints are the most popular evidence that can be used, scene of crime officers claim that more than 30% of the evidence recovered from crime scenes originate from palms. Usually, palmprints evidence left at crime scenes are partial since very rarely full palmprints are obtained. In particular, partial palmprints do not exhibit a structured shape and often do not contain a reference point that can be used for their alignment to achieve efficient matching. This makes conventional matching methods based on alignment and minutiae pairing, as used in fingerprint recognition, to fail in partial palmprint recognition problems. In this paper a new partial-to-full palmprint recognition based on invariant minutiae descriptors is proposed where the partial palmprint’s minutiae are extracted and considered as the distinctive and discriminating features for each palmprint image. This is achieved by assigning to each minutiae a feature descriptor formed using the values of all the orientation histograms of the minutiae at hand. This allows for the descriptors to be rotation invariant and as such do not require any image alignment at the matching stage. The results obtained show that the proposed technique yields a recognition rate of 99.2%. The solution does give a high confidence to the judicial jury in their deliberations and decision.

Problems of minimal aerodynamic resistance

Nesta tese estamos preocupados com o problema da resistência mínima primeiro dirigida por I. Newton em seu Principia (1687): encontrar o corpo de resistência mínima que se desloca através de um médio. As partículas do médio não interagem entre si, bem como a interação das partículas com o corpo é perfeitamente elástica. Diferentes abordagens desse modelo foram feitas por vários matemáticos nos últimos 20 anos. Aqui damos uma visão geral sobre estes resultados que representa interesse independente, uma vez que os autores diferentes usam notações diferentes. Apresentamos uma solução do problema de minimização na classe de corpos de revolução geralmente não convexos e simplesmente conexos. Acontece que nessa classe existem corpos com resistência menor do que o mínimo da resistência na classe de corpos convexos de revolução. Encontramos o infimum da resistência nesta classe e construimos uma sequência regular de corpos que aproxima este infimum. Também apresentamos um corpo de resistência nula. Até agora ninguém sabia se tais corpos existem ou não, evidentemente o nosso corpo não pertence a nenhuma classe anteriormente analisado. Este corpo é não convexo e não simplesmente conexo; a forma topológica dele é um toro, parece um UFO extraterrestre. Apresentamos aqui várias famílias de tais corpos e estudamos as suas propriedades. Também apresentamos um corpo que é natural de chamar um corpo "invisíveis em uma direção", uma vez que a trajectória de cada partícula com a certa direcção coincide com a linha recta fora do invólucro convexo do corpo. ABSTRACT: In this thesis we are concerned with the problem of minimal resistance first addressed by I. Newton in his Principia (1687): find the body of minimal resistance moving through a medium. The medium particles do not mutually interact, and the interaction of particles with the body is perfectly elastic. Different approaches to that model have been tried by several mathematicians during the last 20 years. Here we give an overview of these results that represents interest in itself since all authors use different notations. We present a solution of the minimization problem in the class of generally non convex, simply connected bodies of revolution. It happens that in this class there are bodies with smaller resistance than the minimum in the class of convex bodies of revolution. We find the infimum of the resistance in this class, and construct a sequence of bodies which approximates this infimum. Also we present a body of zero resistance. Since earlier it was unknown if such bodies exists or not, evidently our body does not belong to any class previously examined. The zero resistance body found by us is non-convex and non-simply connected; topologically it is a torus, and it looks like an extraterrestrial UFO. We present here several families of such bodies and study their properties. We also present a body which is natural to call a body "invisible in one direction", since the trajectory of each particle with the given direction, outside the convex hull of the body, coincides with a straight line.

Heurísticas de pesquisa local para problemas de máquina única

O escalonamento é uma das decisões mais importantes no funcionamento de uma linha de produção. No âmbito desta dissertação foi realizada uma descrição do problema do escalonamento, identificando alguns métodos para a optimização dos problemas de escalonamento. Foi realizado um estudo ao caso do problema de máquina única através do teste de várias instâncias com o objectivo de minimizar o atraso pesado, aplicando uma Meta-Heurística baseada na Pesquisa Local e dois algoritmos baseados no SB. Os resultados obtidos reflectem que os algoritmos baseados no SB apresentaram resultados mais próximos do óptimo, em relação ao algoritmo baseado na PL. Os resultados obtidos permitem sustentar a hipótese de não existirem algoritmos específicos para os problemas de escalonamento. A melhor forma de encontrar uma solução de boa qualidade em tempo útil é experimentar diferentes algoritmos e comparar o desempenho das soluções obtidas.

Eco-wall modular solutions for buildings

This paper discusses the development of modular solutions for eco low-cost houses based on a pre-fabricated modular wall system environmentally sustainable, socioeconomically competitive and geared towards developing African nations with a housing deficit. The key point to the research of a modular wall solution is a structural layer complemented with local and materials made by non-specialized workforce. This wall also meets also hydrothermal acoustic and mechanical properties. Thus,the solution also offers good safety and interior comfort conditions to its users while maintaining the flexibility to expand the size of the house. Parameters as dimensions, materials and constructive processes of the existing housing stock were studied. Features such as the family size, typology, different uses, common materials, existing regulations, minimal living conditions, safety and comfort have also been considered to achieve the most efficient solution.

A videovigilância no local de trabalho

The emergence of new technologies has introduced significant changes in the citizens life’s. There is a constant evolution of technological means and profound impact of their use in the habits of life of the human being. These new technological media are important tools in labor relations. The working and businesses worlds are increasingly turning to these new technologies, so that the use of video surveillance in the workplace is nowadays common. New technologies in general and the use of video surveillance in workplace in particular are providing ways to allow control of the work performance that are desired by most employers. However, the collection of images in the workplace often collides with the fundamental rights and freedoms of workers, in particular, with the right to privacy. The subject concerns the question of investigating is whether the images collected in workplace can be used as evidence in disciplinary proceedings. In fact, this issue is controversial. Doctrine and jurisprudence defend, at least, two responses for the same question. Those who understand that the evidence may be admitted for not violate any right of the worker, and others who argue that the evidence should not be admitted in disciplinary office. In the Portuguese legal system, there is, even, a new intermediate theory that begins to be defended, that only on certain occasions the evidence may be admitted. The solution to this problem involves the study of employment law and data video surveillance processing. Analysis of workers fundamental rights is fundamental to come to a grounded conclusion.

Semi-Parametric Weak Instrument Regressions with an Application to the Risk-Return Trade-off

Recent work shows that a low correlation between the instruments and the included variables leads to serious inference problems. We extend the local-to-zero analysis of models with weak instruments to models with estimated instruments and regressors and with higher-order dependence between instruments and disturbances. This makes this framework applicable to linear models with expectation variables that are estimated non-parametrically. Two examples of such models are the risk-return trade-off in finance and the impact of inflation uncertainty on real economic activity. Results show that inference based on Lagrange Multiplier (LM) tests is more robust to weak instruments than Wald-based inference. Using LM confidence intervals leads us to conclude that no statistically significant risk premium is present in returns on the S&P 500 index, excess holding yields between 6-month and 3-month Treasury bills, or in yen-dollar spot returns.

Projection-Based Statistical Inference in Linear Structural Models with Possibly Weak Instruments

It is well known that standard asymptotic theory is not valid or is extremely unreliable in models with identification problems or weak instruments [Dufour (1997, Econometrica), Staiger and Stock (1997, Econometrica), Wang and Zivot (1998, Econometrica), Stock and Wright (2000, Econometrica), Dufour and Jasiak (2001, International Economic Review)]. One possible way out consists here in using a variant of the Anderson-Rubin (1949, Ann. Math. Stat.) procedure. The latter, however, allows one to build exact tests and confidence sets only for the full vector of the coefficients of the endogenous explanatory variables in a structural equation, which in general does not allow for individual coefficients. This problem may in principle be overcome by using projection techniques [Dufour (1997, Econometrica), Dufour and Jasiak (2001, International Economic Review)]. AR-types are emphasized because they are robust to both weak instruments and instrument exclusion. However, these techniques can be implemented only by using costly numerical techniques. In this paper, we provide a complete analytic solution to the problem of building projection-based confidence sets from Anderson-Rubin-type confidence sets. The latter involves the geometric properties of “quadrics” and can be viewed as an extension of usual confidence intervals and ellipsoids. Only least squares techniques are required for building the confidence intervals. We also study by simulation how “conservative” projection-based confidence sets are. Finally, we illustrate the methods proposed by applying them to three different examples: the relationship between trade and growth in a cross-section of countries, returns to education, and a study of production functions in the U.S. economy.

Characterization of the unfolding of a weak focus and modulus of analytic classification

La thèse présente une description géométrique d’un germe de famille générique déployant un champ de vecteurs réel analytique avec un foyer faible à l’origine et son complexifié : le feuilletage holomorphe singulier associé. On montre que deux germes de telles familles sont orbitalement analytiquement équivalents si et seulement si les germes de familles de difféomorphismes déployant la complexification de leurs fonctions de retour de Poincaré sont conjuguées par une conjugaison analytique réelle. Le “caractère réel” de la famille correspond à sa Z2-équivariance dans R^4, et cela s’exprime comme l’invariance du plan réel sous le flot du système laquelle, à son tour, entraîne que l’expansion asymptotique de la fonction de Poincaré est réelle quand le paramètre est réel. Le pullback du plan réel après éclatement par la projection monoidal standard intersecte le feuilletage en une bande de Möbius réelle. La technique d’éclatement des singularités permet aussi de donner une réponse à la question de la “réalisation” d’un germe de famille déployant un germe de difféomorphisme avec un point fixe de multiplicateur égal à −1 et de codimension un comme application de semi-monodromie d’une famille générique déployant un foyer faible d’ordre un. Afin d’étudier l’espace des orbites de l’application de Poincaré, nous utilisons le point de vue de Glutsyuk, puisque la dynamique est linéarisable auprès des points singuliers : pour les valeurs réels du paramètre, notre démarche, classique, utilise une méthode géométrique, soit un changement de coordonée (coordonée “déroulante”) dans lequel la dynamique devient beaucoup plus simple. Mais le prix à payer est que la géométrie locale du plan complexe ambiante devient une surface de Riemann, sur laquelle deux notions de translation sont définies. Après avoir pris le quotient par le relèvement de la dynamique nous obtenons l’espace des orbites, ce qui s’avère être l’union de trois tores complexes plus les points singuliers (l’espace résultant est non-Hausdorff). Les translations, le caractère réel de l’application de Poincaré et le fait que cette application est un carré relient les différentes composantes du “module de Glutsyuk”. Cette propriété implique donc le fait qu’une seule composante de l’invariant Glutsyuk est indépendante.