Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

Front Speed of Language Replacement

We use two coupled equations to analyze the space-time dynamics of two interacting languages. Firstly, we introduce a cohabitation model, which is more appropriate for human populations than classical (non-cohabitation) models. Secondly, using numerical simulations we nd the front speed of a new language spreading into a region where another language was previously used. Thirdly, for a special case we derive an analytical formula that makes it possible to check the validity of our numerical simulations. Finally, as an example, we nd that the observed front speed for the spread of the English language into Wales in the period 1961-1981 is consistent with the model predictions. We also nd that the e¤ects of linguistic parameters are much more important than those of parameters related to population dispersal and reproduction. If the initial population densities of both languages are similar, they have no e¤ect on the front speed. We outline the potential of the new model to analyze relationships between language replacement and genetic replacement

Interdisciplinary Medication Adherence Program: The Example of a University Community Pharmacy in Switzerland.

The Community Pharmacy of the Department of Ambulatory Care and Community Medicine (Policlinique Médicale Universitaire, PMU), University of Lausanne, developed and implemented an interdisciplinary medication adherence program. The program aims to support and reinforce medication adherence through a multifactorial and interdisciplinary intervention. Motivational interviewing is combined with medication adherence electronic monitors (MEMS, Aardex MWV) and a report to patient, physician, nurse, and other pharmacists. This program has become a routine activity and was extended for use with all chronic diseases. From 2004 to 2014, there were 819 patient inclusions, and 268 patients were in follow-up in 2014. This paper aims to present the organization and program's context, statistical data, published research, and future perspectives.

Learning Preferences with Kernel-Based Methods

Learning of preference relations has recently received significant attention in machine learning community. It is closely related to the classification and regression analysis and can be reduced to these tasks. However, preference learning involves prediction of ordering of the data points rather than prediction of a single numerical value as in case of regression or a class label as in case of classification. Therefore, studying preference relations within a separate framework facilitates not only better theoretical understanding of the problem, but also motivates development of the efficient algorithms for the task. Preference learning has many applications in domains such as information retrieval, bioinformatics, natural language processing, etc. For example, algorithms that learn to rank are frequently used in search engines for ordering documents retrieved by the query. Preference learning methods have been also applied to collaborative filtering problems for predicting individual customer choices from the vast amount of user generated feedback. In this thesis we propose several algorithms for learning preference relations. These algorithms stem from well founded and robust class of regularized least-squares methods and have many attractive computational properties. In order to improve the performance of our methods, we introduce several non-linear kernel functions. Thus, contribution of this thesis is twofold: kernel functions for structured data that are used to take advantage of various non-vectorial data representations and the preference learning algorithms that are suitable for different tasks, namely efficient learning of preference relations, learning with large amount of training data, and semi-supervised preference learning. Proposed kernel-based algorithms and kernels are applied to the parse ranking task in natural language processing, document ranking in information retrieval, and remote homology detection in bioinformatics domain. Training of kernel-based ranking algorithms can be infeasible when the size of the training set is large. This problem is addressed by proposing a preference learning algorithm whose computation complexity scales linearly with the number of training data points. We also introduce sparse approximation of the algorithm that can be efficiently trained with large amount of data. For situations when small amount of labeled data but a large amount of unlabeled data is available, we propose a co-regularized preference learning algorithm. To conclude, the methods presented in this thesis address not only the problem of the efficient training of the algorithms but also fast regularization parameter selection, multiple output prediction, and cross-validation. Furthermore, proposed algorithms lead to notably better performance in many preference learning tasks considered.

Catalase in fluvial biofilms: a comparison between different extraction methods and example of application in a metal-polluted river

Antioxidant enzymes are involved in important processes of cell detoxification during oxidative stress and have, therefore, been used as biomarkers in algae. Nevertheless, their limited use in fluvial biofilms may be due to the complexity of such communities. Here, a comparison between different extraction methods was performed to obtain a reliable method for catalase extraction from fluvial biofilms. Homogenization followed by glass bead disruption appeared to be the best compromise for catalase extraction. This method was then applied to a field study in a metal-polluted stream (Riou Mort, France). The most polluted sites were characterized by a catalase activity 4–6 times lower than in the low-polluted site. Results of the comparison process and its application are promising for the use of catalase activity as an early warning biomarker of toxicity using biofilms in the laboratory and in the field

Numerical Dominance of the Argentine Ant vs Native Ants and Consequences on Soil Resource Searching in Mediterranean Cork-Oak Forests (Hymenoptera: Formicidae)

This study is focused on the dominance exerted by the invasive Argentine ant over native ants in a coastal Mediterranean area. Theimpact of this invasive ant on native ant assemblages and its consequences on total ant biomass and on the intensity of habitat explorationwere evaluated. Foraging ants were observed and their trajectories recorded during 5-minute periods in two study zones, one invaded andthe other non-invaded. Ant species detected, ant worker abundance, ant biomass and the intensity of soil surface searching done by antswere compared between the two zones. The Argentine ant invasion provoked a drastic reduction of the ant species richness. Apparentlyonly one native ant species is able to coexist with the Argentine ant, the cryptic Plagiolepis pygmaea. Ant worker abundance was also modified after the invasion: the number of Argentine ant workers detected, which represented 92% of the invaded zone, was two times higher than the number of native ant workers detected in the non-invaded zone. The total ant biomass was inversely affected, becoming four times lower in the invaded zone highly dominated by Linepithema humile. The higher number of Argentine ant workers and their fast tempo of activity implied an alteration of the intensity of soil surface searching: scanning by the Argentine ants in the invaded zone was higher than that done by the native ants in the non-invaded zone, and the estimated time for a complete soil surface scan was 64 minutes in the invaded zone and 108 minutes in the non-invaded zone. Consequently, resources will be discovered faster by ants in the invaded zone than in the non-invaded zone. The increase of the mean temperature and the decrease of the relative humidity from May to August reduced the ant activity in the two study zones but this reduction was greater in the invaded zone

Performance Improvement of Centrifugal Compressor Stage with Pinched Geometry or Vaned Diffuser

Centrifugal compressors are widely used for example in refrigeration processes, the oil and gas industry, superchargers, and waste water treatment. In this work, five different vaneless diffusers and six different vaned diffusers are investigated numerically. The vaneless diffusers vary only by their diffuser width, so that four of the geometries have pinch implemented to them. Pinch means a decrease in the diffuser width. Four of the vaned diffusers have the same vane turning angle and a different number of vanes, and two have different vane turning angles. The flow solver used to solve the flow fields is Finflo, which is a Navier-Stokes solver. All the cases are modeled with the Chien's k – έ- turbulence model, and selected cases are modeled also with the k – ώ-SST turbulence model. All five vaneless diffusers and three vaned diffusers are investigated also experimentally. For each construction, the compressor operating map is measured according to relevant standards. In addition to this, the flow fields before and after the diffuser are measured with static and total pressure, flow angle and total temperature measurements. When comparing the computational results to the measured results, it is evident that the k – ώ-SST turbulence model predicts the flow fields better. The simulation results indicate that it is possible to improve the efficiency with the pinch, and according to the numerical results, the two best geometries are the ones with most pinch at the shroud. These geometries have approximately 4 percentage points higher efficiency than the unpinched vaneless diffusers. The hub pinch does not seem to have any major benefits. In general, the pinches make the flow fields before and after the diffuser more uniform. The pinch also seems to improve the impeller efficiency. This is down to two reasons. The major reason is that the pinch decreases the size of slow flow and possible backflow region located near the shroud after the impeller. Secondly, the pinches decrease the flow velocity in the tip clearance, leading to a smaller tip leakage flow and therefore slightly better impeller efficiency. Also some of the vaned diffusers improve the efficiency, the increment being 1...3 percentage points, when compared to the vaneless unpinched geometry. The measurement results confirm that the pinch is beneficial to the performance of the compressor. The flow fields are more uniform with the pinched cases, and the slow flow regions are smaller. The peak efficiency is approximately 2 percentage points and the design point efficiency approximately 4 percentage points higher with the pinched geometries than with the un- pinched geometry. According to the measurements, the two best geometries are the ones with the most pinch at the shroud, the case with the pinch only at the shroud being slightly better of the two. The vaned diffusers also have better efficiency than the vaneless unpinched geometries. However, the pinched cases have even better efficiencies. The vaned diffusers narrow the operating range considerably, whilst the pinch has no significant effect on the operating range.

Vivianite and siderite in lateritic iron crust: an example of bioreduction

Lateritic iron crust (LIC) samples from Padauari (AM) were analysed by XRD, optical microscopy and SEM-EDS. The equilibrium of iron minerals (IM) was studied using Eh-pH diagram. It was shown that the minerals of the LIC are goethite (alpha-FeOOH), vivianite [Fe3(PO4)2.8H2O] and siderite (FeCO3). Carbonate grains are a solid solution of FeCO3-MnCO3. The LIC presents textures and structures of dissolution of IM. The siderite and vivianite are stable from Eh =-0.3 to 0.0 V and pH=5.0-7.5. These results indicate that vivianite and siderite are products of bioreduction through biogenic dissolution of IM, the new conditions of ecosystems of the Amazon region.

Development of beam and plate finite elements based on the absolute nodal coordinate formulation

The focus of this dissertation is to develop finite elements based on the absolute nodal coordinate formulation. The absolute nodal coordinate formulation is a nonlinear finite element formulation, which is introduced for special requirements in the field of flexible multibody dynamics. In this formulation, a special definition for the rotation of elements is employed to ensure the formulation will not suffer from singularities due to large rotations. The absolute nodal coordinate formulation can be used for analyzing the dynamics of beam, plate and shell type structures. The improvements of the formulation are mainly concentrated towards the description of transverse shear deformation. Additionally, the formulation is verified by using conventional iso-parametric solid finite element and geometrically exact beam theory. Previous claims about especially high eigenfrequencies are studied by introducing beam elements based on the absolute nodal coordinate formulation in the framework of the large rotation vector approach. Additionally, the same high eigenfrequency problem is studied by using constraints for transverse deformation. It was determined that the improvements for shear deformation in the transverse direction lead to clear improvements in computational efficiency. This was especially true when comparative stress must be defined, for example when using elasto-plastic material. Furthermore, the developed plate element can be used to avoid certain numerical problems, such as shear and curvature lockings. In addition, it was shown that when compared to conventional solid elements, or elements based on nonlinear beam theory, elements based on the absolute nodal coordinate formulation do not lead to an especially stiff system for the equations of motion.

A systematic review of Twitter’s hashtags in public health: an example of a globally adopted standard

The objective is to find evidence and make recommendations on the use of Twitter in public health,particularly through the study of hashtags(#). This systematic review shows the use of Twitter in different areas of public health: epidemiological surveillance, health promotion, health protection and disease prevention. Articles on this subject published in indexed journals with impact factor show the importance of conversations to engage the attention of Twitter users by using citations (@ user) and retweets(RT); however, not much importance seems to be given to the use of hashtags(#), which are often assimilated to the concept of keywords. Although tracking recurring hashtagsshould be less expensive than computing Twitter content, the potential of hashtaggeddata has not been properly exploited or recognized over the past years, probably due to the lack of efficient tools.

De los contenidos a las competencias numéricas en la enseñanza obligatoria = From contents to numerical competence in Compulsory Secondary Education

Mathematics has formed part of all our daily lives since the most remote origins of our civilization, although on too many occasions schools have done little for this functional view of mathematics reclaiming the formal role of this discipline in detriment of its more practical and applied element. From this perspective the curriculum organized by competence arises due to the need to fill this vacuum and allow our students to function better in the constant elements of real situations which they are going to have to resolve throughout their lives. The specifying of this general view of mathematic competence in the numeric competence should allow all our students to progressively acquire numeric sense, that is, that they have the capacity to apply good quantitative ideas in real situations

Cálculo do volume na equação de van der Waals pelo método de cardano

Analytical solutions of a cubic equation with real coefficients are established using the Cardano method. The method is first applied to simple third order equation. Calculation of volume in the van der Waals equation of state is afterwards established. These results are exemplified to calculate the volumes below and above critical temperatures. Analytical and numerical values for the compressibility factor are presented as a function of the pressure. As a final example, coexistence volumes in the liquid-vapor equilibrium are calculated. The Cardano approach is very simple to apply, requiring only elementary operations, indicating an attractive method to be used in teaching elementary thermodynamics.

Nanoporous materials: pillared clays and regular silicas as an example of synthesis and their porosity characterization by X-ray diffraction

Because of their practical applications, porous materials attract the attention of undergraduate students in a way that can be used to teach techniques and concepts in various chemistry disciplines. Porous materials are studied in various chemistry disciplines, including inorganic, organic, and physical chemistry. In this work, the syntheses of a microporous material and a mesoporous material are presented. The porosity of the synthesized materials is characterized by X-ray diffraction analysis. We show that this technique can be used to determine the pore dimensions of the synthesized materials.

Polymeeripinnoitteen aiheuttamat itseherätteiset värähtelyt nipillisessä telasysteemissä

Lukuisissa teollisuussovelluksissa materiaalien, kuten paperin ja teräslevyjen, muokkaamiseen käytettävät pyörivät nippitelat kärsivät aina erilaisten herätteiden synnyttämistä mekaanisista värähtelyistä, jotka voivat aiheuttaa virheitä valmistettaviin tuotteisiin. Tässä työssä tutkittiin viskoelastisia polymeerejä ja polymeeripinnoitteen nipilliseen telasysteemiin synnyttämiä haitallisia itseherätteisiä värähtelyjä. Työn polymeerejä käsittelevässä kirjallisuusosassa luotiin katsaus amorfisten polymeerien fysikaalisiin ominaisuuksiin. Kokeellisessa osuudessa tutkittiin tarkemmin kahden amorfisen telapinnoitepolymeerin termoreologisia ja mekaanisia ominaisuuksia suoritettujen DMTA-mittausten perusteella. Sovittamalla toisen polymeerin master-käyrään yleistetty lineaarisen standardiaineen malli saatiin selville polymeerin mekaaniset parametrit ja approksimaatio sen relaksaatiospektrille. Telapinnoitteen nipilliseen systeemiin synnyttämiä itseherätteisiä värähtelyjä ja niiden seurauksia tarkasteltiin kahdelle telalle ja polymeeripinnoitteelle kehitetyn analyyttisen mallin ja numeeristen laskujen avulla. Pinnoite mallinnettiin lineaarisen standardiaineen mukaisesti. Telasysteemin parametrit määritettiin DMTA-mittaustuloksista ja systeemiä vastaavasta koelaitteesta kokeellisella moodianalyysillä ja elementtimenetelmällä. Numeerisesta stabiilisuusanalyysistä ja liikeyhtälöiden integroinneista saadut tulokset kertovat telapinnoitteen aaltomaisista deformaatiomuodoista ja niiden synnyttämistä taajuusalueittain esiintyvistä epästabiileista värähtelyistä. Telasysteemi on epästabiili pinnoitedeformaatiokuvion systeemiin aiheuttaman herätevoiman taajuuden ollessa lähellä systeemin korkeampaa ominaistaajuutta. Numeerisista tuloksista voitiin ennustaa nopean ja hitaan barringin olemassaolo.