This thesis investigates the validity of the standard thermal response test (TRT) results when performed on a series of boreholes (string). The typical TRT consists of subjecting a single borehole to a constant heat injection rate to obtain the temperature response in the ground which can then be used to determine the ground thermal conductivity. When completed on a single borehole, the results may be analyzed with the line source theory, since the assumption of a single line heat source is valid. For multiple boreholes, the assumption of a single line source becomes invalid if the spacing between the boreholes is small enough for borehole thermal interaction to occur. Moreover, for boreholes that are charged in series, heat transfer from the horizontal pipes that connect the vertical boreholes may also influence the ground thermal response. This thesis takes an in-depth look at the different factors that affect the results of TRTs performed on borehole strings. Different analysis methods are implemented to determine areas of improvement for determining the thermal conductivity of the soil surrounding the borehole string. For the analysis, the infinite line source (ILS) model and a model developed using TRNSYS 18 were used to determine the effective thermal conductivity. The results show that TRNSYS is unable to accurately model a TRT performed on a borehole string. The horizontal pipe model within TRNSYS proved to have significant fundamental issues, as the effective thermal conductivity is greatly underestimated with values of 1.2±0.1W/mK and the results of increasing the horizontal length both increased and decreased the effective thermal conductivities. The results from the ILS demonstrate that an effective thermal conductivity of 1.7±0.2W/mK is an appropriate estimate of the soil at the BTES field tested, as the borehole string with the furthest spacing between boreholes gave an effective thermal conductivity of 1.7W/mK. Performing multiple thermal response tests within the same BTES field also provided evidence of the need to implement multiple TRTs as common practise. The testing presented shows that the effective thermal conductivity can vary within ±0.2W/mK within the same relative location. With better knowledge of the thermal properties within the BTES field location comes the opportunity for improved planning of operation and control of thermal distribution within the field. This would be especially beneficial when dealing with seasonal BTES fields / Thesis / Master of Applied Science (MASc)

No description available.

Energy demand for buildings has been rising during recent years. Increasing building energy consumption has caused many energy-related problems and environmental issues. The on-site community energy system application is a promising way of providing energy for buildings. Community energy system usage reduces the primary energy consumption and environmental effects of greenhouse gas (GHG) emissions compared to the implementation of the stand-alone energy systems. Furthermore, due to the increase in electricity price and shortage of fossil fuel resources, renewable energies and energy storage technologies could be great alternative solutions to solve energy-related problems. Generally, the energy system might include various technologies such as internal combustion engine, heat recovery system, boiler, thermal storage tank, battery, absorption chiller, ground source heat pump, heating coil, electric chiller, solar photovoltaics (PV) and solar thermal collectors, and seasonal thermal energy storage. The economic, technical and environmental impacts of energy systems depend on the system design and operational strategy. The focus of this thesis is to propose unified frameworks, including the mathematical formulation of all of the components to determine the optimal energy system configuration, the optimal size of each component, and optimal operating strategy. The proposed methodologies address the problems related to the optimal design of the energy system for both deterministic and stochastic cases. By the use of the proposed frameworks, the design of the energy system is investigated for different specified levels of GHG emissions ratio, and the purpose is to minimize the annual total cost. To account for uncertainties and to reduce the computational times and maintain accuracy, a novel strategy is developed to produce scenarios for the stochastic problem. System design is carried out to minimize the annual total cost and conditional value at risk (CVaR) of emissions for the confidence level of 95%. The results demonstrate how the system size changes due to uncertainty and as a function of the operational GHG emissions ratio. It is shown that with the present-day technology (without solar technologies and seasonal storage), the lowest amount of GHG emissions ratio is 37%. This indicates the need for significant technological development to overcome that ratio to be 10% of stand-alone systems. This thesis introduces novel performance curves (NPC) for determining the optimal operation of the energy system. By the use of this approach, it is possible to identify the optimal operation of the energy system without solving complex optimization procedures. The application of the proposed NPC strategy is investigated for various case studies in different locations. The usage of the proposed strategy leads to the best-operating cost-saving and operational GHG savings when compared to other published approaches. It has shown that other strategies are special (not always optimal) cases of the NPC strategy. Based on the extensive literature review, it is found that it is exceptionally complicated to apply the previously proposed models of seasonal thermal energy storage in optimization software. Besides, the high computational time is required to obtain an optimum size and operation of storage from an optimization software. This thesis also proposes a new flexible semi-analytical, semi-numerical methodology to model the heat transfer process of the borehole thermal energy storage to solve the above challenges. The model increases the flexibility of the storage operation since the model can control the process of the storage by also deciding the appropriate storage zone for charging and discharging. / Thesis / Doctor of Engineering (DEng)

Fracture systems or networks always control the stability, deformability, fluid and gas storage capacity and permeability, and other mechanical and hydraulic behavior of rock masses. The characterization of fracture systems is of great significance for understanding and analyzing the impact of fractures to rock mass behavior. Fracture trace data have long been used by engineers and geologists to character fracture system. For subsurface fractures, however, boreholes, wells, tunnels and other cylindrical samplings of fractures often provide high quality fracture trace data and have not been sufficiently utilized. The research work presented herein is intended to interpret fracture traces on borehole walls and other cylindrical surfaces by using stereology. The relationships between the three-dimension fracture intensity measure, P32, and the lower dimension fracture intensity measures are studied. The analytical results show that the conversion factor between the three-dimension fracture intensity measure and the two-dimension intensity measure on borehole surface is not dependent on fracture size, shape or circular cylinder radius, but is related to the orientation of the cylinder and the orientation distribution of fractures weight by area. The conversion factor between the two intensity measures is determined to be in the range of [1.0, π/2]. The conversion factors are also discussed when sampling in constant sized or unbounded fractures with orientation of Fisher distribution. At last, the author proposed estimators for mean fracture size (length and width) with borehole/shaft samplings in sedimentary rocks based on a probabilistic model. The estimators and the intensity conversion factors are tested and have got satisfactory results by Monte Carlo simulations. / Ph. D.

Multiple saline tracers were used to explore the role of geologic structure on groundwater recharge at the Fractured Rock Research Site in Floyd County, Virginia. Tracer migration was monitored through soil, saprolite, and fractured crystalline bedrock for a period of 3 months with chemical, physical, and geophysical techniques. Potassium chloride (KCl) and potassium bromide (KBr) tracers were applied at specific locations on the ground surface to directly test flow pathways in a shallow saprolite and deep fractured rock aquifer. Previous work at the Fractured Rock Research Site have identified an ancient thrust fault complex that is present in the otherwise competent metamorphic bedrock; fracturing along this fault plane has resulted in a highly transmissive aquifer that receives recharge along the vertically oriented portion of the fault zone. A shallow aquifer has been located above the thrust fault aquifer in a heterogeneous saprolite layer that rapidly transmits precipitation to a downgradient spring. Tracer monitoring was accomplished with differential electrical resistivity, chemical sampling, and physical monitoring of water levels and spring discharge. Tracer concentrations were monitored quantitatively with ion chromatography and qualitatively with differential resistivity surveys. KCl, applied at a concentration of 10,000 mg/L, traveled 160 meters downgradient through the thrust fault aquifer to a spring outlet in 24 days. KBr, applied at a concentration of 5,000 mg/L, traveled 90m downgradient through the saprolite aquifer in 19 days. KCl and KBr were present at the sampled springheads for 30 days and 33 days, respectively. Tracer breakthrough curves indicate diffuse flow through the saprolite aquifer and fracture flow through the crystalline thrust fault aquifer. Heterogeneities in the saprolite aquifer had a large effect on tracer transport, with breakthrough peaks varying several days over vertical distances of several meters. Monitoring saline tracer migration through soil, saprolite, and fractured rock provided data on groundwater recharge that would not have been available using other traditional hydrologic methods. Travel times and flowpaths observed during this study support preferential groundwater recharge controlled by geologic structure. Geologic structure, which is not currently considered an important factor in current models of Blue Ridge hydrogeology, should be evaluated on a local or regional scale for any water resources investigation, wellhead protection plan, or groundwater remediation project. / Master of Science

Thesis (PhD (Electrical and Electronic Engineering)--University of Stellenbosch, 2007. / ENGLISH ABSTRACT: A monostatic baseband radar is required with the capability of detecting close range targets that appear at distances comparable to the system’s resolution, without compromising the radar’s maximum range. The application in borehole radar imposes further constraints associated with the physical limitations and variable electromagnetic environment of different borehole diameters and conditions. This dissertation discusses the complete design process of the analog section of a monostatic radar that successfully addresses these issues. The proposed transceiver employs a series duplexing arrangement consisting of an antenna, transmitter, receiver and an isolation switch. An exponentially decaying tail is observed in the current flowing on a borehole radar antenna when excited by pulse waveforms. The characteristics of this tail depend strongly on the borehole environment. A measurement technique is developed that accurately quantifies this exponential decay by digitizing a logarithmic representation of the antenna current while it is operating in various boreholes. Transmitters are then designed to drive these antennas with waveforms that prevent the formation of current tails. This is achieved through the use of pole-zero networks or alternatively by generating certain asymmetric, bipolar waveforms. The transmitters are simultaneously designed to have an output impedance approximating a short circuit after the transient is generated. In the series configuration proposed here, the duplexing of the antenna between transmitter and receiver is then reduced to simply isolating the receiver during transmit-mode. The switch responsible for this isolation disconnects the receiver and presents a short circuit between antenna and transmitter during transmit-mode, while connecting the receiver terminals between the antenna and the short circuited transmitter terminals in receive-mode. The required close-in performance of the transceiver dictates that the transition between these two states of the isolation switch occur in a time similar to the duration of the transmitter waveform. The switching artefacts generated by the switch are consequently similar to the radar data signal. The isolation switch employs an innovative configuration (using both transistors and diodes) which accepts a single control signal and causes the switching artefacts to be generated as a common mode signal, while a differential path is created for the radar data signal which is being switched. This leads to effective suppression of the switching signal in the signal passed to the receiver. Dissipative filtering is advocated as a fundamental design principle for high fidelity receivers and it is shown how it can be applied by using constant impedance equalizers and diplexers as basic building blocks. This principle is used as the basis for the design of this transceiver's receivers, which incorporate both standard gain blocks and operational amplifiers. A complete borehole radar system, based on the transceiver developed here, was built and tested; resulting in the first known practical monostatic borehole radar system. Data obtained in field trials are presented and suggest that the monostatic system compares well with current state of the art bi-static systems. / AFRIKAANSE OPSOMMING: Die behoefte is geïdentifiseer vir 'n monostatiese basisbandradar wat oor die vermoë beskik om nabygeleë teikens op 'n afstand soortgelyk aan die resolusie van die stelsel waar te neem, sonder om die maksimum bereik van die stelsel in te kort. Die toepassing daarvan in 'n boorgatradarstelsel lei tot verdere vereistes vanweë die fisiese beperkings en veranderende elektromagnetiese omgewing van boorgate met verskillende deursnitte en toestande. Hierdie proefskrif is gemoeid met die volledige ontwerpsprosedure van die analoog gedeelte van 'n monostatiese radar wat al hierdie kwessies aanspreek. 'n Serie verbinding van antenne, sender, ontvanger en isolasieskakelaar word ingespan vir hierdie ontwerp. Eksponensieel wegsterwende stertjies word waargeneem in die antennestroom van 'n boorgatradarantenne wanneer dit aangedryf word deur puls golfvorms. 'n Meettegniek word ontwikkel wat hierdie eksponensiële verslapping noukeurig kan monitor deur 'n logaritmiese voorstelling van die antennastroom te versyfer terwyl dit ontplooi word in verskillende boorgate. Senders word dan ontwikkel om hierdie antennes aan te dryf met golfvorms wat juis die vorming van hierdie stertjies voorkom. Dit word bewerkstellig deur die gebruik van pool-zero netwerke of andersins deur die opwek van sekere asimmetriese, bipolêre golfvorms. Die senders se uittree-impedansies moet egter terselfdertyd ontwerp word om 'n kortsluiting te benader sodra die oorgang klaar opgewek is. Met die serie verbinding wat hier gebruik word, raak die vereiste tyddeling van die antenna tussen die sender en ontvanger dan bloot 'n geval van ontvanger-isolasie gedurende uitsaai-modus. Die skakelaar wat verantwoordelik is vir hierdie isolasie ontkoppel die ontvanger en vertoon soos 'n kortsluiting tussen sender en antenne tydens uitsaai-modus, maar verbind weer die terminale van die ontvanger tussen die antenne en kortgeslote senderterminale tydens ontvang-modus. Die vereiste kortafstand vermoë van die stelsel veroorsaak dat die tysduur van die oorgang tussen hierdie twee modusse soortgelyk is aan dié van die sender golfvorm en enige skakelverskynsels wat opgewek word deur die skakelaar is gevolglik soortgelyk aan die radardatasein self. Die isolasieskakelaar gebruik egter 'n innoverende konfigurasie (met transistors sowel as diodes) wat funksioneer met 'n enkele beheersein en die skakelverskynsels as gemene modus seine opwek, terwyl 'n differensiële seinpad geskep word vir die radardatasein wat geskakel word. Die skakelseine word gevolglik effektief onderdruk in die sein wat oorgedra word aan die ontvanger. Die gebruik van verkwistende filters word voorgestel as 'n fundamentele ontwerpsbeginsel vir hoëtrou ontvangers en daar word getoon hoe dit toegepas kan word met konstante impedansie vereffeningsbane en dipleksers. Hierdie beginsel is dan ook gebruik as basis vir die ontwerp van hierdie stelsel se ontvangers, wat gebruik maak van beide standard aanwinsblokke sowel as operasionel versterkers. 'n Volledige boorgatradarstelsel, gebaseer op die stelsel wat hier ontwikkel is, is gebou en getoets. Die gevolg is die eerste bekende, praktiese monostatiese boorgatradarstelsel. Data wat hiermee verwerf is word aangebied en dui daarop dat die monostatiese stelsel baie goed opweeg teen huidige bi-statiese stelsels.

This study includes a systematic investigation of the bond strength of cementitious borehole plugs in welded tuff. Analytical and numerical analysis of borehole plug-rock stress transfer mechanics is performed. The interface strength and deformation are studied as a function of Young's modulus ratio of plug and rock, plug length and rock cylinder outside-to-inside radius ratio. The tensile stresses in and near an axially loaded plug are analyzed. The frictional interface strength of an axially loaded borehole plug, the effect of axial stress and lateral external stress, and thermal effects are also analyzed. Implications for plug design are discussed. Push-out tests are used to determine the bond strength by applying an axial load to the cement plugs. A total of 130 push-out tests are performed as a function of borehole size, plug length, temperature, and degree of saturation of the tuff cylinder. The use of four different borehole radii enables evaluation of size effects. A well-defined exponential strength decrease with increasing plug diameter results.

Les capteurs distribués à fibre optique (aussi nommés DAS) sont une nouvelle technologie d'acquisition sismique qui utilise des câbles traditionnels à fibre optique pour fournir une mesure de la déformation le long du câble. Ce système d'acquisition est largement utilisé dans les profils sismiques verticaux (PSV). Le couplage est un facteur clé qui a une grande influence sur la qualité des données. Alors que, pour les acquisitions PSV, les géophones sont attachés à la paroi du puits, le câble de fibre optique est soit cimenté derrière le tubage, soit attaché avec des pinces rigides au tubage ou simplement descendu dans le puits. Cette dernière stratégie de déploiement donne généralement le plus petit rapport signal sur bruit, mais est considérée comme la plus rentable en particulier pour les installations dans des puits existants. Cette thèse porte sur la problématique du couplage du DAS quand le câble est simplement descendu dans le puits. Nous développons des modèles numériques pour analyser les données réelles. L'interprétation de ces résultats nous permet de conclure qu'un contact immédiat du câble avec la paroi du puits avec une force de contact calculée est nécessaire pour fournir des bonnes conditions de couplage. Sur la base de ces résultats, nous proposons des solutions pour optimiser davantage les acquisitions avec le système DAS. Nous modifions numériquement la force de contact et les propriétés élastiques du câble DAS et démontrons comment ces modifications peuvent améliorer mais aussi détériorer la qualité des données. Enfin, nous proposons un algorithme de détection du couplage qui permet d'assurer l'acquisition de données réelles avec un rapport signal / bruit élevé. / Distributed Acoustic Sensing (DAS) is a new technology of seismic acquisition that relies on traditional fibre-optic cables to provide inline strain measurement. This acquisition system is largely used in vertical seismic profiling (VSP) surveys. Coupling is a key factor influencing data quality. While geophones and accelerometers are clamped to the borehole wall during VSP surveys, the fibre cable is either clamped and then cemented behind the casing, or attached with rigid clamps to the tubing, or loosely lowered into the borehole. The latter deployment strategy, also called wireline deployment, usually acquires the lowest level of signal but is regarded as the most cost-effective in particular for existing well installations. This PhD thesis addresses the problematic of coupling of DAS using wireline deployment. We develop numerical models that are used to analyse real data. The interpretation of these results allows us concluding that an immediate contact of the cable with the borehole wall with a computed contact force is required to provide good coupling conditions. Based on those findings, we propose solutions to further optimise DAS acquisitions. We numerically modify the contact force and the elastic properties of the DAS cable and show how these modifications can improve but also deteriorate data quality. Finally, we propose a coupling detection algorithm that is applied to real datasets and allows ensuring the acquisition of data with a high signal-to-noise ratio.

A model has been developed which simulates a system consisting of a horizontal and vertical ground-coupled heat exchanger connected in parallel to the same heat pump. The model was used in computer simulations to investigate how the annual minimum and mean fluid temperatures at the heat pump varied as several parameters of the combined system were changed. A comparison was also made between different control settings for fluid flow rate distribution between the two exchangers. For the case when the flow rate distribution was not controlled, the effect of viscosity differences between a colder and warmer exchanger was investigated. The short term effects of letting the vertical heat source rest during the warm summer months was then tested. Lastly, the results of the model was compared to a simple 'rule of thumb' that have been used in the industry for this kind of combined system. The results show that using a combined system might not always result in increased performance, if the previously existing exchanger is a vertical ground-coupled heat exchanger. The effects of viscosity differences on the flow distribution seems to be negligible, especially for high net flows. Controlling the fluid flow rates seems to only be worth the effort if the the pipe lengths of the two combined exchangers differ heavily. Letting the vertical ground-coupled heat exchanger rest during summer was shown to in some cases yield an increased short-term performance in addition to the already known positive long term effects. The rule of thumb was shown to recommend smaller dimensions for combination systems than the more realistic analytical model.

Local soil conditions influence the characteristics of earthquake ground shaking and these effects must be taken into account when specifying ground shaking levels for seismic design. These effects are quantified via site response analysis, which involves the propagation of earthquake motions from the base rock through the overlying soil layers to the ground surface. Site response analysis provides surface acceleration-time series, surface acceleration response spectra, and/or spectral amplification factors based on the dynamic response of the local soil conditions. This dissertation investigates and compares the results from different site response methods. Specifically, equivalent-linear time series analysis, equivalent-linear random vibration theory analysis, and nonlinear time series analysis are considered. In the first portion of this study, hypothetical sites and events are used to compare the various site response methods. The use of hypothetical events at hypothetical sites allowed for the seismic evaluation process used in engineering practice to be mimicked. The hypothetical sites were modeled after sites with characteristics that are representative of sites in the Eastern and Western United States. The input motions selected to represent the hypothetical events were developed using the following methods: stochastically-simulated time series, linearly-scaled recorded time series, and spectrally-matched time series. The random vibration theory input motions were defined using: seismological source theory, averaging of the Fourier amplitude spectra computed from scaled time series, and a response spectrum compatible motion. All of the different input motions were then scaled to varying intensity levels and propagated through the sites to evaluate the relative differences between the methods and explain the differences. Data recorded from borehole arrays, which consist of instrumentation at surface and at depth within the soil deposit, are used to evaluate the absolute bias of the site response methods in the second portion of this study. Borehole array data is extremely useful as it captures both the input motion and the surface motion, and can be used to study solely the wave propagation process within the soil deposit. However, comparisons using the borehole data are complicated by the assumed wavefield at the base of the array. In this study, sites are selected based on site conditions and the availability of high intensity input motions. The site characteristics are then developed based on site specific information and data from laboratory soil testing. Comparisons between the observed and computed response are used to first assess the wavefield at the base of the array, and then to evaluate the accuracy of the site response methods. / text