Integration of 3D geological and numerical models based on tetrahedral meshes for hydrogeological simulations in fractured porous media

Une nouvelle approche de modélisation des milieux géologiques fracturés représentés par un modèle conceptuel de fractures discrètes et déterministes est présentée dans cette thèse. L'objectif principal de l'étude est de reproduire l'hétérogénéité et la complexité des milieux poreux fracturés dans un modèle géométrique tridimensionnel afin d'effectuer des simulations numériques dans le but d'améliorer les capacités de modélisation en hydrogéologie. Ceci est réalisé à travers le couplage entre une plateforme de modélisation géologique (GOCAD) et un code numérique (HydroGeoSphere). Les principaux défis à relever sont: la représentation géométrique du réseau de fractures, la sélection d'un maillage approprié pour la discrétisation spatiale du domaine de simulation et l'adaptation du code numérique à ce maillage. La nouvelle approche est basée sur une première phase de modélisation géologique 3D, suivie par la génération d'un maillage tétraédrique 3D et par la simulation numérique de l'écoulement souterrain en conditions saturées et du transport de solutés. En général, le maillage tétraédrique s'avère plus adéquat que les maillages de blocs ou de prismes pour discrétiser les geometries complexes telles que les milieux fracturés. De plus, une définition alternative du maillage "dual", qui est essentiel pour appliquer la méthode numérique élément finis - volume de contrôle utilisée par HydroGeoSphere, est analysée et intégrée dans le code numérique. Le code numérique proposé est d'abord vérifié par l'intermédiaire de simples scénarios de simulation dont les solutions, analytiques et numériques, sont déjà connues. La complexité des simulations est augmentée de façon graduelle. L'approche de modélisation est finalement appliquée au site Olkiluoto (Finlande) où un laboratoire de recherche souterrain est en construction afin d'évaluer la faisabilité du stockage géologique profond de déchets nucléaires à haute activité. Les techniques de modélisation géologique mises au point permettent de modéliser facilement la géométrie des fractures identifiées à travers la caractérisation géologique in situ. De plus, le modèle numérique s'avère adéquat pour la simulation de l'écoulement et du transport de solutés dans ce site complexe. Ce travail de recherche présente une contribution au développement des techniques de modélisation hydrogéologique des milieux fracturés

Capturing the two-way hydromechanical coupling effect on fluid-driven fracture in a dual-graph lattice beam model

Fluid-driven fractures of brittle rock is simulated via a dual-graph lattice model. The new discrete hydromechanical model incorporates a two-way coupling mechanism between the discrete element model and the flow network.By adopt- ing an operator-split algorithm, the coupling model is able to replicate the transient poroelasticity coupling mechanism and the resultant Mandel-Cryer hydromechanical coupling effect in a discrete mechanics framework. As crack propagation, coalescence and branching are all path-dependent and irreversible processes, capturing this transient coupling effect is important for capturing the essence of the fluid-driven fracture in simulations. Injection simulations indicate that the onset and propagation of fractures are highly sensitive to the ratio between the injection rate and the effective permeability. Furthermore, we show that in a permeable rock, the borehole breakdown pressure, the pressure at which fractures start to grow from the borehole, depends on both the given ratio between injection rate and permeability and the Biot coefficient

Notranja oksidacija Cu-C in Ag-C kompozitov

The internal oxidation in copper-carbon and silver-carbon composites occurs when they are exposed to air or oxygen at high temperature. Solubility of carbon in copper or in silver is very low. The kinetics of oxidation at high temperature and activation energy were determined and the mechanism of internal oxidation was analysed. The kinetics of internal oxidation was determined for both cases and it is depended from the diffusion of oxygen following parabolic time dependence according to Wagner\u27s theory. The activation energy for Cu-C composite is 70.5 kJ/mol, and for Ag-C composite is 50.1 kJ/mol, what is in both cases close to the activation energy for the volume diffusion of oxygen in copper or in silver. In both cases gas products are formed during the internal oxidation of composites. In the internal oxidation zone pores, bubbles occur. The carbon oxidates directly with the oxygen from solid solution as long there is a contact, which breaks down with the presence of gas products. Then the oxidation occurs over the gas mixture of CO and CO2.Pri visokih temperaturah kompoziti bakra in srebra z ogljikom na zraku ali v kisiku reagirajo po mehanizmu notranje oksidacije. Topnost ogljika v trdnem bakru in trdnem srebru je zelo majhna. Analizirali smo kinetiko oksidacije kompozitov, določili aktivacijsko energijo in mehanizem notranje oksidacije. Kinetika oksidacije je pri obeh skupinah materialov odvisna od difuzije kisika in sledi parabolični odvisnosti od časa v skladu z Wagnerjevo teorijo. Aktivacijska energija procesa je za kompozit Cu-C enaka 70,5 kJ/mol, za kompozit Ag-C pa 50,1 kJ/mol, kar je blizu aktivacijski energiji za volumsko difuzijo kisika v trdnem bakru oziroma srebru. Pri oksidaciji kompozita nastajajo plinski produkti. Oksidacija ogljika poteka neposredno s kisikom iz trdne raztopine, ko pa se zaradi nastanka plinske faze stik prekine, pa preko plinske zmesi CO in CO2

Enabling the Development and Implementation of Digital Twins : Proceedings of the 20th International Conference on Construction Applications of Virtual Reality

Welcome to the 20th International Conference on Construction Applications of Virtual Reality (CONVR 2020). This year we are meeting on-line due to the current Coronavirus pandemic. The overarching theme for CONVR2020 is "Enabling the development and implementation of Digital Twins". CONVR is one of the world-leading conferences in the areas of virtual reality, augmented reality and building information modelling. Each year, more than 100 participants from all around the globe meet to discuss and exchange the latest developments and applications of virtual technologies in the architectural, engineering, construction and operation industry (AECO). The conference is also known for having a unique blend of participants from both academia and industry. This year, with all the difficulties of replicating a real face to face meetings, we are carefully planning the conference to ensure that all participants have a perfect experience. We have a group of leading keynote speakers from industry and academia who are covering up to date hot topics and are enthusiastic and keen to share their knowledge with you. CONVR participants are very loyal to the conference and have attended most of the editions over the last eighteen editions. This year we are welcoming numerous first timers and we aim to help them make the most of the conference by introducing them to other participants

Time-lapse seismic and electrical resistivity tomography combined for monitoring of the CO2 storage site, Ketzin, Germany

Imaging of subsurface fluid migration is a longstanding geophysical research task, in which CO2 storage monitoring is a relatively recent application. It is well known that such monitoring can benefit from a combination of different geophysical methods, which complement each other with regard to imaging characteristics and sensitivity to variable saturation of CO2. At the Ketzin site, Germany, pilot-scale CO2 injection is performed with then aim to improve the understanding of in-situ physical, chemical and biological processes, and to provide practical and operational experience for future geological storage of CO2. This thesis addresses the geophysical monitoring of the CO2 that was injected at the Ketzin site. Therefore, time-lapse analyses are carried out on the basis of repeated seismic experiments and electrical resistivity tomographies (ERT), both of which are presented separately in the first place. Finally, a methodical combination of these methods is proposed and applied to image the injected CO2. The thesis gathers these investigations as follows: First, a synthetic modelling of crosshole ERT and surface seismic experiments is performed. This study is carried out in order to identify key sensitivities of the time-lapse images from the Ketzin site and to provide assistance for the subsequent processing of real data. Therefore, a static reservoir model is derived from well logs and the structural interpretation of the baseline 3D seismic data. Further, multiphase flow simulations are used to establish dynamic CO2 distributions in this model. The study indicates that a CO2 quantification may not be possible on the basis of the seismic amplitude information only, since the thicknesses of the CO2 distributions do not exceed the seismic resolution limit with regard to typical wavelengths that are provided by surface seismic data. Secondly, the time-lapse processing of the baseline and first repeat 2D seismic surveys from the Ketzin site is presented. The first survey was performed before the start of the injection in 2005 and the second survey in 2009, after approximately 22 kilotons of CO2 had been injected. Although the datasets show a good repeatability, near-surface velocity changes are found to have a degrading impact on the quality of the time-lapse images. As these changes could only be imperfectly resolved by refraction static corrections, a pre-stack static correction is proposed, which decomposes the timing delays of baseline and repeat traces in a surface-consistent manner. Together with a test of a post-stack static correction, this pre-stack static correction is shown to improve the quality of the time-lapse images considerably. In the subsequent interpretation of the time-lapse images, this provided evidence that in 2009 no CO2 related amplitude change is observable where the 2D lines allow for monitoring of the reservoir. This finding is further confirmed by an AVO analysis and a comparison with the corresponding 3D surveys. Thirdly, large-scale ERT surveys are presented, which have been carried out repeatedly during the site startup and the first year of CO2 injection. The experimental setup of these surveys combines surface measurements with downhole measurements by utilizing a permanent electrode array which has been deployed in the three Ketzin wells. One baseline and three repeat experiments are presented, which show a resistivity increase over time at the CO2 injector and indicate a preferential CO2 migration towards the northwest. Using an experimental resistivity-saturation relationship, CO2 saturations of up to 70% are mapped near the injection well, which is consistent with CO2 saturations obtained from pulsed neutron-gamma logging. Finally, a combination of seismic and geoelectric investigations is presented using a structurally constrained inversion approach. For this purpose, lithostructural constraints are interpreted from the seismic reflection data and implemented in the geoelectric inversion by means of a local regularization technique. Consequently, seismics and constrained resistivity inversion are arranged in a sequential workflow which is based on a structural similarity of seismic parameters and electric resistivity. In other words, a change in elastic impedance is expected to occur together with a change in resistivity and vice versa. Prior to an application to the Ketzin datasets, this approach was tested through a synthetic data example. In consistency, both the synthetic and the real data example demonstrated that the constrained inversion allows for an enhanced resistivity imaging along the caprock-reservoir boundary than a conventional ERT inversion. The practical demonstration for the Ketzin datasets shows, that the approach has a potential for an integrated geophysical monitoring of CO2 storages and is also significant for imaging of other subsurface processes, which trigger changes in elastic parameters and electric resistivity

Sestava in morfologija diboridov v zlitinah Al-Ti-B po žarjenju na 1873 K

In the investigation, Al-Ti-B alloys with different Ti/B ratios were annealed at 1873 K for 10 hours. Under this condition diboride particles were in contact with molten aluminium. During annealing considerable growth of remaining unmelted diboride particles took place. The composition of diboride particles shifted extremely close to the composition of pure TiB2 regardless of the Ti/B ratios in the alloys, indicating that in equilibrium almost pure TiB2 coexisted with a melt. Diboride particles were bound with well-defined crystal facets. Types and relative size of the facets depended on the Ti/B ratio. It is strongly indicated that an excess of Ti (B) in the melt over that needed to form TiB2 stabilizes crystal planes having planar arrangement of Ti (B) atoms.V tej raziskavi smo več zlitin Al-Ti-B, ki so imele različna razmerja Ti/B, žarili 10 ur pri 1873 K. V teh razmerah so bili delci diborida v kontaktu z aluminijevo talino. Med žarjenjem so diboridni delci močno zrasli. Njihova sestava se je v vseh raziskanih zlitinah močno približala sestavi čistega TiB2, kar kaže, da je v termodinamskem ravnotežju s talino praktično čisti TiB2. Diboridni delci so bili obdani s kristalnimi fasetami. Vrsta in relativna velikost faset je bila močno odvisna od razmerja Ti/B. Obstajajo trdni dokazi, da presežek Ti (B) nad vrednostjo, ki je potrebna za tvorbo TiB2, stabilizira kristalne ravnine, ki imajo ravninsko razporeditev Ti (B) atomov

3-D Electrical Resistivity Tomography for Cliff Stability Assessment at Pointe du Hoc in Normandy, France

Pointe du Hoc overlooking the English Channel in Normandy, France was host to one of the most important military engagements of World War II. While the site is a valuable historic cultural resource, it is vulnerable to cliff collapses that already have endangered the observation post and Rudder?s command post. The observation post has been closed to visitors for some time due to safety concerns. Geophysical techniques have been used increasingly in recent years for slope stability investigation purposes. The objective of this study is to apply advanced 3?D resistivity tomography toward a detailed site stability assessment with special attention to the two at-risk buildings. 3?D resistivity tomography datasets at Pointe du Hoc in the presence of extreme topography and dense cultural clutter have been successfully acquired, inverted, and interpreted. The cliff stability in the areas around the two at?risk buildings has been analyzed. A hazard assessment scheme has been designed in which regions of high resistivity are interpreted as zones of open, dry fractures with a moderate mass movement potential. Regions of low resistivity are zones of wet, clay?filled fractures with a high mass movement potential. The observation post tomography results indicate that the highest mass movement hazard appears to be associated with the marine caverns at the base of the cliff that are positioned at the point of strongest wave attack. These caverns likely occupy the future site of development of a sea arch which will definitely threaten the observation post building. A high probability of a soil wedge failure is on the east?facing cliff edge close to the observation post that could damage or destroy the building. The mass movement potential at the Rudder?s command post area is low to moderate. The greatest risk is associated with soil wedge failures at the top of the cliffs. The resistivity geophysical data add great value to the natural geohazard assessment at Pointe du Hoc and constitute an integral component of an interdisciplinary approach to the problem of cultural resource preservation at the site. Geophysics is a non?invasive and relatively inexpensive technology that provides unique constraints which are unobtainable using traditional engineering geology methods for site characterization. However, the technology is difficult to master and the inherent limitations must be carefully understood to ensure a reliable geotechnical interpretation