Reconstruction of the surface temperature in the Kura depression (Azerbaijan) by the inversion of borehole data

In order to reconstruct the ground surface temperature history (GSTH) in the Kura depression from the data on borehole temperatures, we selected two thermal logs, which met the requirements of the well temperature inversion into the GSTH. The temperature gradients measured in these wells varied about 20 K/km, which is typical for the region of study. The borehole temperatures were inverted into the ground surface temperatures in the past in accordance with the program developed by Po Yu Shen. It was found that the ground surface temperature had increased by 1–2°C during the last century. The GSTH derived from the well temperature inversion agrees with the climate temperatures measured at the Ganca meteorological station, which has been operating since 1873. These results obtained for the Kura depression are consistent also with the results of similar studies in other regions in the world.     

Climate Change in the Urals,Russia, Inferred from Borehole Temperature Data

Borehole temperatures in the central and south Urals were analysed for the past ground surface temperature (GST) signal. 31 highquality temperature logs were selected for this purpose and inverted with algorithms based on the generalised least squares theory. The signal to noise ratio was improved by averaging the results of individual borehole inversions. No distinct regional trends were found in the studied region except for some indications of more pronounced warming in the south. The mean GST history (GSTH) was characterised by cooling down to –0.6 °C in the 18th century and subsequent warming to 0.5 °C above the longterm mean at the beginning of this century, and to 1 – 1.5 °C by 1980. The stability of the mean GSTH was tested in dependence on the number of holes used for the averaging. It showed that any subset of 15 holes yielded a GSTH similar to that obtained from the whole set. A surface air temperature (SAT) time series comprising the period 1832 – 1989 was combined from 17 meteorological records. Its least squares warming rate of 1.1 °C per 100 years is somewhat higher than that of the GST (0.7 – 0.8°C/100 years) in the same period.     

Pre-Observational Evolution of Surface Temperature in Romania as Inferred from Borehole Temperature Measurements

Temperature data from nine boreholes in the Carpathian orogen in Romania were used to obtain information on the ground surface temperature history (GSTH) in the last 250?years. The temperature measurements were taken with a thermistor probe (sensitivity in the 10 mK range) using the stop-and-go technique, at 10 m intervals, in the depth range of 20–580?m. The least squares inverse modelling approach of Tarantola and Valette (J Geophys 50:159–170, 1982) was used to infer the GSTH. Long-term air temperature records available from the Romanian weather station network were used as a comparison term for the first 100–150?years of the GSTH, and as a forcing function in a POM-SAT model that combines borehole temperature profiles (BTPs) and meteorological time series (surface air temperature, SAT) to produce information on the so-called pre-observational mean (POM). Results from a global circulation model for the Romanian area are incorporated in the discussion as well.     

Quantitative guidance for efficient vertical flow measurements at the sediment–water interface using temperature–depth profiles

Upward discharge to surface water bodies can be quantified using analytical models based on temperature–depth (T-z) profiles. The use of sediment T-z profiles is attractive as discharge estimates can be obtained using point-in-time data that are collected inexpensively and rapidly. Previous studies have identified that T-z methods can only be applied at times of the year when there is significant difference between the streambed–water interface and deeper sediment temperatures (e.g., winter and summer). However, surface water temperatures also vary diurnally, and the influence of these variations on discharge estimates from T-z methods is poorly understood. For this study, synthetic T-z profiles were generated numerically using measured streambed interface temperature data to assess the influence of diurnal temperature variations on discharge estimation and provide insight into the suitable application of T-z methods. Results show that the time of day of data collection can have a substantial influence on vertical flux estimates using T-z methods. For low groundwater discharge fluxes (e.g., 0.1 m d⁻¹), daily transience in streambed temperatures led to relatively large errors in estimated flow magnitude and direction. For higher discharge fluxes (1.5 m d⁻¹), the influence of transient streambed temperatures on discharge estimates was strongly reduced. Discharge estimates from point-in-time T-z profiles were most accurate when the uppermost point in the T-z profile was near the bed interface daily mean (two time periods daily). Where temperature time series data are available, daily averaged T-z profiles can produce accurate discharge estimates across a wide range of discharge rates. Seasonality in shallow groundwater temperature generally had a negligible influence on vertical flow estimates. These findings can be used to plan field campaigns and provide guidance on the optimal application of T-z methods to quantify vertical groundwater discharge to surface water bodies.     

Resistivity inversion of cross-hole and borehole-to-surface EM data using axially symmetric models1

We describe a least-squares inversion approach to estimating the subsurface resistivity structure from cross-hole or borehole-to-surface electromagnetic data. It is assumed that the resistivity distribution is symmetric about the axis of a borehole and that vertical magnetic dipoles are located on the borehole axis. The receivers are placed either in another borehole or on the earth's surface. The inversion scheme uses the finite-element and smoothness-constrained least-squares methods. The computational effort required to obtain partial derivatives is reduced considerably by using the reciprocity principle. Numerical simulations show that the reconstructions are generally in good agreement with the true structures when the assumption of an axisymmetric earth structure holds. An example involving the breakdown of this assumption, which can be obtained by interchanging the source and receiver boreholes, suggests that the inversion result may also be useful for locating a general 3D anomaly although artifacts are present.     

Multiwell study of seismic attenuation at the CO2CRC Otway project geosequestration site: Comparison of amplitude decay,centroid frequency shift and 1D waveform inversion methods

At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO₂-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO₂ plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Q_scat ≈ 50–65.     

Delineation of Ground Water Flow in Fractured Rock in the Southwestern Part of Manhattan, New York, Through Use of Advanced Borehole Geophysical Methods

Advanced borehole-geophysical methods were used to assess the geohydrology of fractured crystalline bedrock at five test boreholes in southwestern Manhattan Island, New York, in preparation for construction of a third water tunnel for New York City. The boreholes penetrated gneiss and other crystalline bedrock that has an overall southwest to northwest dipping foliation with a 60° dip. Most of the fractures encountered are either nearly horizontal or have moderate northwest dip azimuths. Fracture indexes range from 0.25 to 0.44 fracture per foot (0.3 m) of borehole.
Electromagnetic (EM) and heat-pulse flowmeter logs obtained under ambient and pumping conditions, together with other geophysical logs, indicate transmissive fracture zones in each borehole. Pumping tests of each borehole indicated transmissivity ranges from <2 to 360 ft²/day (0.2 to 33 m²/day). Ground water appears to flow within an interconnected fracture network toward the south and west within the study area. No correlation was indicated between the fracture index and the total borehole transmissivity.     

A comparative study of geothermal and meteorological records of climate change in Kamchatka

To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski. Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing 20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature (SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response to the changes occurred on the surface.     

Comparison of formation and fluid-column logs in a heterogeneous basalt aquifer

Deep observation boreholes in the vicinity of active production wells in Honolulu, Hawaii, exhibit the anomalous condition that fluid-column electrical conductivity logs and apparent profiles of pore-water electrical conductivity derived from induction conductivity logs are nearly identical if a formation factor of 12.5 is assumed. This condition is documented in three boreholes where fluid-column logs clearly indicate the presence of strong borehole flow induced by withdrawal from partially penetrating water-supply wells. This result appears to contradict the basic principles of conductivity-log interpretation. Flow conditions in one of these boreholes was investigated in detail by obtaining flow profiles under two water production conditions using the electromagnetic flowmeter. The flow-log interpretation demonstrates that the fluid-column log resembles the induction log because the amount of inflow to the borehole increases systematically upward through the transition zone between deeper salt water and shallower fresh water. This condition allows the properties of the fluid column to approximate the properties of water entering the borehole as soon as the upflow stream encounters that producing zone. Because this condition occurs in all three boreholes investigated, the similarity of induction and fluid-column logs is probably not a coincidence, and may relate to aquifer response under the influence of pumping from production wells.     

Ground surface temperature history in southern Canada: Temperatures at the base of the Laurentide ice sheet and during the Holocene

We use temperature profiles from 7 deep (≈ 2000 m) boreholes located in southern Canada to infer ground surface temperature histories (GSTH) during the Last Glacial Maximum (LGM) and the Holocene. Visual inspection of the heat flow and of the reduced temperature depth profiles reveals significant regional differences with some sites showing conspicuous signs of post glacial warming, and other indicating only very small changes in ground surface temperature. These differences are confirmed by the inversions of the temperature profiles. The most prominent variations in GST are found at the Sudbury, Ontario, sites where the present ground surface temperature is high. With the exception of Sept-Iles, Quebec, the other sites only show moderate or no variation in GST. For all the sites, except possibly Sept-Iles, temperatures at the base of the ice sheet during the LGM were at or slightly below the melting point of ice. Temperatures might have been lower, a few degrees below 0 °C, at Sept-Iles. These results are consistent with field observations and model predictions suggesting high velocity basal flows in the ice sheet above the studied regions. These new data on basal temperatures will provide better quantitative constraints on glacier flow dynamics. The inversions give a chronology for the retreat of the ice sheet comparable to other proxies. Inversion and direct modeling show that, following the ice retreat, there was a warm period between 2 and 5 ka with temperatures 1–2 K higher than present. The inversion yields a time for this episode 1–2 kyr more recent than that inferred by other proxies for the Holocene climate optimum (HCO).     

Predicting horizontal velocities from well data

The Imperial College borehole test site consists of four boreholes with depths lying between 260 and 280 m. The boreholes intersect several cyclical sequences of sandstones, mudstones and limestones. The formations are highly laminated and ultrasonic measurements on preserved core have shown that the mudstones are intrinsically anisotropic. Little or no anisotropy is associated with the sandstones and limestones. A scheme is proposed to predict synthetic vertical and horizontal P- and S-wave logs. Combining (an)isotropic effective medium theories, the Gassmann equation and Backus averaging, the scheme extends previous sand-shale models to transversely isotropic rock formations. The model assumes that the anisotropy is due to layering and due to the preferred horizontal orientation of the clay minerals, pores and cracks within the mudstones. The pores and cracks within the sandstones and limestones are randomly orientated. After fitting the model to the ultrasonic data to obtain the unknown parameters, the model successfully predicts the sonic log and the direct arrival times from a cross-hole survey.     

Impact of Groundwater Flow and Energy Load on Multiple Borehole Heat Exchangers

The effect of array configuration, that is, number, layout, and spacing, on the performance of multiple borehole heat exchangers (BHEs) is generally known under the assumption of fully conductive transport. The effect of groundwater flow on BHE performance is also well established, but most commonly for single BHEs. In multiple‐BHE systems the effect of groundwater advection can be more complicated due to the induced thermal interference between the boreholes. To ascertain the influence of groundwater flow and borehole arrangement, this study investigates single‐ and multi‐BHE systems of various configurations. Moreover, the influence of energy load balance is also examined. The results from corresponding cases with and without groundwater flow as well as balanced and unbalanced energy loads are cross‐compared. The groundwater flux value, 10^?7 m/s, is chosen based on the findings of previous studies on groundwater flow interaction with BHEs and thermal response tests. It is observed that multi‐BHE systems with balanced loads are less sensitive to array configuration attributes and groundwater flow, in the long‐term. Conversely, multi‐BHE systems with unbalanced loads are influenced by borehole array configuration as well as groundwater flow; these effects become more pronounced with time, unlike when the load is balanced. Groundwater flow has more influence on stabilizing loop temperatures, compared to array characteristics. Although borehole thermal energy storage (BTES) systems have a balanced energy load function, preliminary investigation on their efficiency shows a negative impact by groundwater which is due to their dependency on high temperature gradients between the boreholes and surroundings.     

Thermal state of the Campi Flegrei caldera inferred from seismic attenuation tomography

A three-dimensional Qp image of the Campi Flegrei caldera between 0 and 3 km of depth has been inferred by the inversion of P rise time and pulse width data of 87 local earthquakes recorded during the last bradiseismic crisis by a local array deployed in the area by the University of Wisconsin. The availability of both thermal measurements in 5 deep boreholes and of a heat flow surface map of the area allowed us to calibrate the local temperature F vs. Qp relationship. The comparison of Qp, Vp and Vp/Vs images, combined with hydrogeological and geochemical data from deep boreholes, allowed us to distinguish some low-Qp anomalies related to the presence of fluids in the rocks from a deep low-Qp anomaly related to the conductive cooling of a magma reservoir. The deep anomaly is located in the same zone where several authors believe that the volcanic and magmatic activity migrated after the Neapolitan Yellow Tuff eruption. Moreover this anomaly includes the area where the existence of a magma chamber at depth between 4 and 5 km was inferred by an active seismic experiment.     

Ground Surface Warming History in Northern Canada Inferred from Inversions of Temperature Logs and Comparison with Other Proxy Climate Reconstructions

— Well temperature logs from 61 sites located in discontinuous and continuous permafrost regions of northern Canada are analyzed. The method of functional space inversion (FSI) is applied to the set of precise temperature logs from wells located between 60 ° and 82 °N. There is strong evidence of extensive ground surface temperature (GST) warming beginning in the late 18^th century and lasting until the 20^th century. This was preceded by a lengthy period of cooling. The approximate average increase of the surface temperature of Canadian Arctic, based on all individual GST histories, is > 1.3 °C for the last 200 years. Simultaneous inversion of all well temperature data suggests an even higher warming (approximately 2 °C). There has been no strong south-to-north gradient in the ground warming magnitude when northern Canadian data are compared with eastern and central Canadian data south of 60 °N which also shows warming close to 2 °C. In Alaska, warming of some 2 °C has been restricted mainly to the 20^th century. In general, however, a high warming magnitude is common for most of Canada and Alaska for the previous century. The averaged GST history (GSTH) for the Canadian Arctic is calibrated with and compares visually with a variety of recently published regional and hemispheric proxy climate reconstructions. These show that GST warming derived from well temperature logs is generally higher than one shown by other proxy (mainly tree-ring reconstructions).Received April, 2003     

Broadband synthetic aperture borehole radar interferometry

Trials in mines have established that wideband VHF borehole radars (BHR), working in the 10–100-MHz band, can be used to probe the rockmass between boreholes over ranges from <5 m to as much as 150 m with submeter resolution. There is evidence that ore bodies reflect these radar signals both specularly and diffusely, much as the ground/air interface does when overflown by synthetic aperture radar (SAR). In both SAR and BHR, multiple flight lines, together with diffuse reflections admit the possibility of developing interferometric 3D images of the object. This paper examines the possibility of imaging buried objects in three dimensions by interferometrically combining broadband VHF borehole radar profiles shot in adjacent pairs of boreholes. Broadbanding in BHR has the advantage of releasing the image from 2nπ phase ambiguities, but practically, interferometric borehole radar (InBHR) needs high signal-to-noise ratios (SNR) to avoid noise capture. This means that 3D InBHR is limited to ranges in wavelengths which are less than the rock's attenuation factor Q. Interferometric methods are developed which are capable of mapping ore bodies and other structures in three dimensions. Tangent plane migration methods are developed here in order to reconstruct surfaces that lie in the near-field of sparse interferometric arrays.     

Physical property anisotropy of foliated fault rocks: Study from the Nobeoka Thrust,Shimanto Belt,southwest Japan

To investigate the physical property anisotropies of foliated fault rocks in subduction zones, the hanging wall phyllites and footwall cataclasites exhumed along the Nobeoka Thrust, a fossilized out‐of‐sequence‐thrust in the Shimanto Belt, Japan, was focused. Discrete physical property (electric resistivity, P‐ and S‐wave velocities, and porosity) measurements were conducted employing geologic coordinates (depth‐parallel direction, strike direction, and maximum dip direction of foliation), using the core samples obtained from the Nobeoka Thrust Drilling Project and compared the data to borehole geophysical logs. A higher sample P‐wave velocity (Vp), lower S‐wave velocity (Vs), higher Vp/Vs, and lower sample porosity and resistivity compared to the logs, are inferred to have been caused by the larger sampling scale of the logs and lower fluid saturation of the borehole. The phyllites and cataclasites exhibited substantial vertical and horizontal anisotropy of Vp (0.4–17.3 % and 2.7–13.8 %, respectively), Vs (0.5–56 % and 7.7–43 %, respectively), and resistivity (0.9–119 % and 2.0–65.9 %, respectively). The physical property anisotropies are primarily affected by the dip angles of foliation. The fault rocks that have gentler dip angles exhibit a higher Vp in the strike and maximum dip direction and a lower Vp in the depth‐parallel direction. In contrast, the fault rocks that have steeply dipping structures show a higher Vp in the strike and depth‐parallel directions with a lower velocity in the maximum dip direction. Resistivity anisotropy show a trend opposite to that of the Vp in relation to the dip angles. Our results show lower Vp anisotropy than those obtained in previous studies, which measured wave speeds perpendicular or parallel to foliation under confining pressure. This study highlights the significance of dip angles on vertical properties in geophysical surveys across foliated fault rocks.     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

New combined log and tracer test interpretation method for identifying transfers in fissured aquifers

Tracer tests represent the most appropriate approach for assessing hydrodispersive parameters such as transversal and longitudinal dispersivities or kinematic porosity on an aquifer scale. They are generally carried out by injecting a tracer in a borehole and measuring its concentration over time in neighboring boreholes by extracted volume sampling or downhole measurements. Logging is one of the most suitable methods for evaluating fissured reservoirs. But short circuits between fractures with different hydraulic potential through boreholes induce mixing phenomena that cannot be avoided without packers. This mixing can shift the breakthrough curves deduced from the logs for each producing fracture and distort determination of their location.
The method proposed in this paper aims at measuring the flow rate and the solute breakthrough for hydraulically active fractures, in open boreholes. It involves estimating a velocity profile along the borehole column by the analysis of two successive logs: a shift function according to depth is thus determined by comparison between log portions on each successive one. The velocity gradients reflect the inward or outward flow rates produced by each fracture. On the basis of these flow rates, it is possible to determine the mixing effects inside the borehole and then to plot unbiased breakthrough curves for each producing fracture.
This method was applied at a granitic site in the eastern Pyrenees. In spite of some questionable limitations, the results showed that the method seems adapted to situations with many fractures. The precise hydraulic pattern which is obtained at the borehole scale is discussed in terms of a dual porosity model. Furthermore, interpretation of the breakthrough curves for fractures corrected for mixing effects revealed that Peclet numbers are strongly underestimated if this phenomenon is not considered.     

Fracture-frequency prediction from borehole wireline logs using artificial neural networks

Borehole-wall imaging is currently the most reliable means of mapping discontinuities within boreholes. As these imaging techniques are expensive and thus not always included in a logging run, a method of predicting fracture frequency directly from traditional logging tool responses would be very useful and cost effective. Artificial neural networks (ANNs) show great potential in this area. ANNs are computational systems that attempt to mimic natural biological neural networks. They have the ability to recognize patterns and develop their own generalizations about a given data set. Neural networks are trained on data sets for which the solution is known and tested on data not previously seen in order to validate the network result. We show that artificial neural networks, due to their pattern recognition capabilities, are able to assess the signal strength of fracture-related heterogeneity in a borehole log and thus fracture frequency within a borehole. A combination of wireline logs (neutron porosity, bulk density, P-sonic, S-sonic, deep resistivity and shallow resistivity) were used as input parameters to the ANN. Fracture frequency calculated from borehole televiewer data was used as the single output parameter. The ANN was trained using a back-propagation algorithm with a momentum learning function. In addition to fracture frequency within a single borehole, an ANN trained on a subset of boreholes in an area could be used for prediction over the entire set of boreholes, thus allowing the lateral correlation of fracture zones.     

3D waveform inversion of downhole microseismic data for transversely isotropic media

3D anisotropic waveform inversion could provide high-resolution velocity models and improved event locations for microseismic surveys. Here we extend our previously developed 2D inversion methodology for microseismic borehole data to 3D transversely isotropic media with a vertical symmetry axis. This extension allows us to invert multicomponent data recorded in multiple boreholes and properly account for vertical and lateral heterogeneity. Synthetic examples illustrate the performance of the algorithm for layer-cake and ‘hydraulically fractured’ (i.e. containing anomalies that simulate hydraulic fractures) models. In both cases, waveform inversion is able to reconstruct the areas which are sufficiently illuminated for the employed source-receiver geometry. In addition, we evaluate the sensitivity of the algorithm to errors in the source locations and to band-limited noise in the input displacements. We also present initial inversion results for a microseismic data set acquired during hydraulic fracturing in a shale reservoir.