Wire-line log-based stratigraphy of flood basalts from the Lopra-1/1A well, Faroe Islands 7 Wire-line log-based stratigraphy of flood basalts from the Lopra-1/1A well, Faroe Islands Lars O. Boldreel The present study shows that it is possible to use conventional borehole logs to perform a detailed lithological/stratigraphical division of a column of subaerially extruded basalt. A stratigraphical divi- sion of the subaerial flood basalts penetrated by the Lopra-1/1A well has been carried out using new wire-line logging data measured in 1996 in the interval 200–2489 m depth. Resistivity data acquired in the interval 200–2178 m depth during 1981 after the initial drilling of the Lopra-1 well have also been incorporated. Eighty-six individual flow units, 18 compound flows and two dolerite dykes have been identified by combining the NPHI porosity, RHOB density, P-, S- and Stonely-sonic transit time, calliper and resistivity logs. Fifty-two sedimentary/tuffaceous layers have also been identified using the CGR and SGR gamma ray and potassium logs in combination with the aforementioned logs. Within the flow units, sonic velocity, density and resistivity are highest in the core where porosity is lowest. This relation is reversed in the uppermost and basal zones of the flow units. The sonic velocity in the core seems to be independent of the thickness of the flow unit. Porous zones seem abundant in some cores and the total section of cores containing porous zones constitutes more than 70% of the thickness of its flow unit, but where porous zones are absent the core makes up only roughly 50% of the thickness of the flow. It is suggested that the flow units with porous cores repre- sent aa flows (88% of the flow units) and the others pahoehoe flows (12% of the flow units). The log pattern of the flow units (crust, core and basal zone) is similar to log patterns reported from other basalt plateaux. However the patterns in Lopra-1/1A show a larger variation than else- where, suggesting that the flow units are more complex vertically than previously thought. Statistical analysis of P-, S- and Stonely-waves, RHOB, NPHI, resistivity, gamma and calliper logs has been carried out. Cross-plots based on the lithological divisions have been produced that show a pro- nounced reduction in scattering versus P-sonic transit time and P- and S-sonic, RHOB and NPHI logs correlate with depth. The geochemical logs do not reflect the cyclic structure of the flow units and probably represent the primary composition of the basalt. The thorium log especially indicates flow units with high and low radioactivity and it is suggested that a minimum of 36 flow fields form the logged part of the lower basalt series. Dolerite units described in previous works have been con- firmed based on the combined interpretation of wire-line logs. The log data suggest that the subaeri- ally extruded basalt has its base at a depth of approximately 2490 m and that a hyaloclastite succession is found below that depth. The transition from subaerially extruded basalt to hyaloclastite produces a negative acoustic impedance and it is found that the transition corresponds to a negative reflection interpreted on VSP surveys from 1988 and 1994. Keywords: Flood basalts, flood basalt petrophysics, North Atlantic, Faroe Islands, Lopra-1/1A well, Faroese basalt lithology, wire-line logs _________________________________________________________________________________________________________________ Geological Institute, University of Copenhagen, Øster Voldgade 10, DK-1350 Copenhagen K, Denmark. E-mail: lob@geol.ku.dk © GEUS, 2006. Geological Survey of Denmark and Greenland Bulletin 9, 7–22. Available at: www.geus.dk/publications/bull GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:197 8 The Lopra-1 well was deepened in 1996 to drill through a pronounced negative seismic reflection at an estimated depth of 2366 m and an additional reflection at a depth of 3486 m found from the analysis of the data acquired in Vertical Seismic Profiling (VSP) surveys carried out in 1988 (Kiørboe & Petersen 1995) and 1994. Several wire- line logs were run in the upper part of the well (0–2184 m) that had been drilled in 1981 to supplement the previous logging (Nielsen et al. 1984). A more extensive logging programme was carried out in the deepened part of Lopra- 1/1A and the new wire-line data for the entire well have been analysed. The wire-line logs were used to investigate the lower basalt series of the Faroe Islands, to construct a lithologi- cal and stratigraphical division of the lava sequence and to test whether a subdivision of the basalt column would result in better statistics of its physical properties. The re- sults of the new logging have been both new understand- ing and results in greater detail than was possible from the data available from the 1981 drilling and logging, and they supplement studies carried out in other basalt cov- ered regions in the world. It was intended that analysis of the wire-line logs ac- quired in 1996 should include the logs acquired in 1981, but this turned out to be complicated. It was possible to use the older resistivity data only after they had been cor- rected manually within individual flow units. Consequent- ly this paper is concerned mainly with interpretation of the wire-line log data measured in 1996 in the subaerially extruded lower basalt series. Geology of the Faroe Islands The Faroe Islands consist of subaerially extruded tholeiit- ic basalt flows erupted in the late Paleocene in connection with the opening of the northern North Atlantic Ocean. The basalts consist of three lava formations called the upper, middle and lower basalt formations (or series). A tuff-agglomerate zone and a coal-bearing horizon sepa- rate the lower and the middle formations (e.g. Rasmussen & Noe-Nygaard 1969, 1970; 1990). The three lava for- mations are all exposed on the islands and the lowest ex- posed stratigraphical level is found in the southern part of the southernmost island, Suðuroy, where the original Lo- pra-1 well was drilled. Geological and geophysical inves- tigations based on data from the Lopra-1 well resulted in six papers published in Berthelsen et al. (1984). The geol- ogy of the well was described by Hald & Waagstein (1984) who found the lava to consist of tholeiitic basalt. A strati- graphical division of the well was established on the basis of drill cuttings collected every 2 m, five cores, the wire- line logs (acquired in 1981) and geochemical rock analy- sis. This resulted in a division of the sequence into 87 flow units, 27 sedimentary units and two dolerite dykes within the interval 200–2180 m depth plus an additional 10 flow units and five sedimentary beds in the cased in- terval 0–200 m. The bedding in the well was established using chilled vesicular (amygdaloidal) basalt as a marker of flow tops. The observed transition between vesicular and non-vesicular basalt was afterwards correlated to the neutron-neutron and resistivity logs (Waagstein et al. 1982; Nielsen et al. 1984). The flow units range in thickness from a few metres to about 50 m with the majority of flows having thicknesses of 5–10 m or 10–15 m. Most of the flows were interpreted as aa type based on the charac- teristics of the topmost cuttings. Sediments were observed on the top of about one third of the drilled lava flows. Two dolerite dykes were found at 834–738 m and 616– 508 m, and Hald & Waagstein (1984) suggested that the real thickness of both dykes is about 9 m assuming an inclination of 85°. The volcaniclastic sediments were dis- tributed more or less randomly throughout the drilled sequence, and no correlation was found either with flow thickness or with chemical composition of the flows. It was concluded that there were no major hiatuses during eruption of the drilled sequence or the overlying 700 m of exposed lower basalt formation. Data from the Lopra-1 well The original Lopra-1 stratigraphic well reached a total depth (TD) of 2184 m in 1981. It was deepened in 1996 Table 1. Wire-line log suites 1, 2 and 5 used in this study Suite 1: 196–2184 m (in the existing hole drilled in 1981). Two log runs were carried out: Run No. 1: LDL-CNL-NGS (recording a maximum hole deviation of 1º from the vertical). GR was run to the surface Run No. 2 pass 1: DSI-Cross and U&L Dipole Run No. 2 pass 2: DSI-P&S and Stoneley. Suite 2: 2185–3158 m Three runs were carried out: Run No. 1: LDL-CNL-NGS-DSI. Run No. 2: SHDT-GR-AMS recording a maximum hole deviation from the vertical of 14° Run No. 3: GHMT-GR. Suite 5: 2204–3531 m Run No 1: MSCT-GR GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:198 9 to a TD of 3158 m and was sidetracked from 3091 m after a partly unsuccessful fishing operation. The sidetrack (Lo- pra-1A) well reached TD at 3565 m and the Lopra-1/1A well sections were plugged and abandoned after final wire-line logging. The Lopra-1/1A well has been logged on various occasions. Caliper, gamma ray, resistivity, neutron-neu- tron porosity and sonic logs were acquired in 1981 (Nielsen et al. 1984), but the sonic log and part of the resistivity log turned out to be unsuccessful (Nielsen et al. 1984). Suc- cessful VSP surveys were carried out in 1988 (Kiørboe & Pedersen 1995) and 1994, and ahead-of-tool predictions were carried out in order to evaluate the remaining thick- ness of the basalt. The deepened well was logged in 1996 and a large number of logs, including a VSP, were ob- tained. A total of five suites were run between 196 and 3565 m depth. Suites 1 to 5 were acquired in the intervals 196–2184 m, 2185–3158 m, 3091–3516 m, 3516–3561 m and 2204–3531 m, respectively. Three of the suites have been used in this study, suites 1, 2 and 5 (Table 1). From the suite of wire-line logs, the calliper, SGR (total gamma ray), CGR gamma ray (computed gamma ray, including measurements of thorium and potassium), P-, S- and Stonely-wave transit time, RHOB (bulk density), NPHI (neutron porosity hydrogen index), resistivity (short and long from 1981), potassium, thorium and uranium logs were selected for this study. After processing (apart from the resistivity logs) by Schlumberger, the data were interpreted on an SGI UNIX workstation using Land- mark software (Stratworks® and Petroworks®). There were problems with the logs acquired by the spectral gamma log tool in the deepened part of the Lopra-1 well, probably because calibration had not been carried out properly. The data were corrected by shifting the thorium, potassium and CGR logs to the responses measured in the previously drilled hole. An attempt was made to correlate the wire- line logs from 1981 and 1996, as both suites include a calliper. That task proved not to be simple, probably be- cause of different cable tensions and depth annotations. It is felt that the 1981 resistivity data have been correctly adjusted by manual editing and adjusting the 1981 data to the interpretation of the wire-line suite obtained in 1996. Characteristics of subaerially extruded basalt It is known from various workers (e.g. Cashman & Kaua- hikaua 1997; Self et al. 1997; Thordarson & Self 1998; Waagstein 1999) that subaerially extruded lava flows have some characteristics in common and the following dis- cussion focuses on those characteristics which may be ex- pected to be confirmed from the logging suite run in the Lopra-1 well. A lava formed in a single volcanic eruption is called a flow field and consists of one or several lava flows erupted more or less continuously from the same vent area or fissure. The commonest diagnostic criterion likely to be used is the chemical composition of the basalt magma that is assumed to be fairly constant during an eruption. A lava flow may consist of several flow units, partly or completely surrounded by chilled crust. A flow unit con- sists of a lava crust, a lava core and a basal zone. The appea- rance of the crust has been used to define two morpholo- gical types, aa and pahoehoe, in both of which the crust is usually much thicker than the basal zone. The core typi- cally consists of massive basalt with few or no vesicles, and sometimes displays flow banding. Thick cores often show a more or less columnar jointing in the lower part which propagated upwards perpendicular to the base of the flow during cooling. Thinner and less regular columns in the upper part may have propagated downwards from the top of the flow. Vesicular zones in the core have been reported by McMillan et al. (1987) who showed a two-stage vesiculation of basalt from the Colombia River Basalt Pla- teau. Walker (1989) studied vesicle distribution in pro- files across flow units, on Hawaii primarily of pahoehoe type but also of aa type, and concluded that the distribu- tion of vesicles does not need to be symmetrical about a horizontal median plane. The rubble near the surface of a flow unit in aa lava may be unconsolidated with fractures and voids filled by tuff, and sometimes the rubble fragments are welded. The fragmentation decreases downwards, and a complete grada- tion into unbrecciated basalt in the core of the flow unit is commonly seen. Thick flow units in the lower forma- tion on the Faroes are typically from several metres to sev- eral tens of metres thick, as seen on Suðuroy (Waagstein 1999). Thick flow units usually have an upper rubbly crust several metres thick whereas the lower crust or basal zone is generally less than one metre. The vesicles and fractures are usually filled with secondary minerals. Aa lavas on Hawaii normally consist of a large number of thin flow units, 0.2–2.0 m thick. Pahoehoe type lava has a continuous vesicular crust that is distinguished from its core mainly by its vesicularity, and roughly half of the total thickness of a pahoehoe flow consists of crust. The crust is highly vesicular near its top surface with abundant small vesicles that often exceed 50% of the bulk rock volume. The vesicles show an overall de- crease in abundance and increase in size downwards re- flecting the increase in gas pressure exerted by the thick- ening lava crust. The vesicles are often arranged in hori- GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:199 10 zontal layers with a vertical spacing of about 10 cm, espe- cially in the upper part of the crust. Horizontal vesicular veins or sheets of basalt up to about 10 cm thick may occur near the top of the core, but otherwise the core of the flow is usually massive with no or only a few relatively large vesicles. The basal vesicular zone is usually a few tens of centimetres thick at most and sometimes sub-vertical pipes are found within it. Sediments Thin beds of tuff with a vertical extent of between a few centimetres and few tens of metres are found within the Faroes lava succession. The beds are typically very fine- grained, clayey, fissile with vivid colours. The tuff beds are buried soils (palaeosols) formed by disintegration and che- mical break down of volcanic materials. The beds may originate from either airfall volcanic ash or from in situ weathering of the lava flows (Hald & Waagstein 1984; Waagstein 1999). Microfossils suggest that the latter formed slowly enough to allow immigration of plants and other living organisms before being covered by the next lava flow, thus representing quiescence in the erupting environment (Lund 1983). Previous work on the relationship between lithology and logging response The geophysical response of interbedded sediment/tuff layers and subaerially extruded basalt, characterised both by the rhythmic succession of flow units and the subdivi- sion of the flow unit into three parts, may be examined from wire-line logging data and drilling penetration rate. This has been done by a number of authors from various places around the world: the Faroe Islands (Kern & Rich- ter 1979; Nielsen et al. 1984), ODP Hole 642 (Planke & Flovenz 1996), ODP Hole 642E (Planke 1994), Newark Basin (Goldberg et al. 1994), North Sea basalt flows from the Middle Jurassic (Rider 1996), Deccan Trap area (Buck- ley & Oliver 1990; Singh 1996), Japan and Chile (Berlitz et al. 1988) and a geothermal well with no cited geograph- ical location (Sanyal et al. 1980). Physical parameters ob- tained from the analyses of wire-line data by these workers are presented in Table 2, and the results from the present study are found in Table 3. In two studies where the lay- ered basalt flows represent the Seaward Dipping Reflector Sequence, from ODP Hole 917 (which also had forma- tion microscanner images) (Planke & Flovenz 1996) and from ODP Hole 642E (Planke 1994), the interpretation of wire-line logs was integrated with continuous coring. The above studies suggest that the following features are of special interest in relation to the wire-line logs ob- tained from Lopra-1/1A. The often vesicular, fractured and altered crust and flow tops are characterised by gradients in the values represented * Chile; † Japan; # uncorrected values; § average cps: counts per second; p.u.: neutron porosity units; GAPI: Gamma log in API units; mas: The massive core of a flow unit; por: The porous crust of a flow unit; sed: sediments between flow units. Kern & Richter (1979) Faroese basalt Berlitz et al. (1988) Chile Japan Goldberg et al. (1994) USA Planke (1994) ODP Hole 642E Planke & Flovenz (1996) ODP Hole 917 Singh (1996) Deccan traps 10 p.u. 3.06*; 2.8† 15 3–5.5 4§ 2.5§ 6–8 m 4.79–6.21 2.73–3.52 2.73–3.52 Christensen (1996) 5.91§ 3.03§ 2.882§ 6.2§ 10%# 3.05§ 75 ± 25 5.0–6.0 3.0–5.3 < 5% corrected 15–30% corrected 10–20% corrected 2.3–2.8 5 GAPI 20 GAPI 15–40 cps 40–100 cps 100–300 cps Table 2. Physical properties of basalt based on wire-line logs Vp (km/sec) Vp massive (km/sec) Vp sediments (km/sec) Vs (km/sec) Porosity Porosity massive Porosity porous Porosity sediment Density (g/cm3) Gamma-ray (GAPI) Gamma-ray (mas) Gamma-ray (por) Gamma-ray (sed) Unit thickness 6–8 m GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1910 11 Vp flow unit Vp core Vp crust Vp compound flow Vp dolerite dyke Vp sediments Vs flow unit Vs core Vs crust Vs compound flow Vs dolerite dyke Vs sediments Vst flow unit Vst core Vst crust Vst compound flow Vst dolerite dyke Vst sediments Porosity flow unit Porosity core Porosity crust Porosity compound flow Porosity dolerite dyke Porosity sediments Density flow unit Density core Density crust Density compound flow Density dolerite dyke Density sediments Gamma ray flow unit Gamma ray core Gamma ray crust Gamma ray compound flow Gamma ray dolerite dyke Gamma ray sediments Calliper flow unit Calliper core Calliper crust Calliper compound flow Calliper dolerite dyke Calliper sediments Resistivity (lld) flow unit Resistivity (lld) core Resistivity (lld) crust Resistivity (lld) compound flow Resistivity (lld) dolerite dyke Resistivity (lld) sediments Resistivity (Ild) flow unit Resistivity (lls) core Resistivity (lls) crust Resistivity (lls) compound flow Resistivity (lls) dolerite dyke Resistivity (lls) sediments 5.53 5.83 4.69 4.92 6.08 4.26 2.97 3.14 2.52 2.59 3.28 2.34 1.46 1.47 1.42 1.44 1.45 1.40 0.19 0.15 0.31 0.30 0.12 0.31 2.85 2.90 2.66 2.74 2.90 2.58 10.40 10.23 9.75 10.32 4.74 27.09 10.47 10.52 10.15 10.26 11.81 11.84 11.96 11.96 11.73 11.80 12.12 12.08 Average Table 3. Physical properties of basalt based on wire-line logs from Lopra-1 317.16 369.51 165.68 167.66 515.80 91.91 558.39 661.47 283.45 (37.35–1777.40) 1038.07 288.27 356.15 122.45 121.70 * limestone porosity units; om: ohms m2/m 1055.64 601.04 187.91 144.09 9.50 (0.7–61.1) API 9.39 (0.0–33.77) 9.43 (0.0–33.76) 9.18 (1.34–36.53) 4.57 (1.36–26.95) (1.66–78.28)25.21 2.96 2.93 2.66 2.76 2.95 2.64 0.32 0.10 0.30 0.32 0.14 0.16 1.41 1.47 1.44 1.43 1.48 1.47 2.34 3.39 2.60 (1.90–3.34) (1.90–4.03) (1.82–3.44) 2.54 (1.74–3.82) 3.19 3.11 4.30 6.29 5.02 4.73 5.90 5.78 473.36 527.20 76.63 (24.07–1777.40) om (6.64–15.83) inches (6.64–15.83) (8.24–15.62) (10.10–14.61) (9.49–18.58) (10.09–18.41) (2.03–3.01) (1.51–3.04) (2.53–2.95) (1.55–3.15) (1.57–3.15) (1.55–3.15) g/cm3 (0.11–0.50) (0.04–0.45) (0.15–0.49) (0.15–0.44) (0.04–0.47) (0.04–0.50) % lspu* (1.19–1.49) (1.28–1.48) (1.19–1.49) (1.11–1.50) (1.22–1.52) (1.11–1.52) km/sec (2.00–4.01) (1.74–4.01) km/sec (2.86–6.31) (3.37–7.58) (2.43–6.57) (2.94–6.63) (2.94–7.47) (2.94–7.47) km/sec Median (Min–Max) 164.21255.35 (22.37–309.21) (22.37–760.66) (26.32–596.06) (18.52–731.21) (32.43–766.68) (18.52–766.66) om (26.74–890.36) (145.01–1773.44) (52.07–1135.72) (24.67–1197.98) 178.05 GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1911 12 by the porosity dependent wire-line logs. The P-wave velo- city, density and resistivity values are rather low but in- crease into the flow, in contrast to the porosity log that shows rather high values that decrease into the flow (Planke & Flovenz 1996). Natural gamma counts are intermedi- ate (Goldberg et al. 1994; Planke 1994; Rider 1996; Singh 1996). The crust is often soft and weathered and can be recognised from an increase in the rate of drilling pentra- tion (Buckley & Oliver 1990). The massive core is characterised by rather high values of P-wave velocity, density and resistivity in contrast to rather low values of porosity and low natural gamma counts (Goldberg et al. 1994; Planke 1994; Planke & Flovenz 1996; Rider 1996; Singh 1996). A number of Faroes basalt samples, selected to be free of fractures and secondary alteration, showed in laboratory experiments that the P- and S-wave velocities increase as a function of pressure (Kern & Richter 1979). Two distinct gamma ray count rates of primary potassium and thorium show that high-gamma basalts exist (Buckley & Oliver 1990). The basal zone is characterised by gradients in values from the wire-line logs representing porosity. P-wave velo- city, density and resistivity decrease rapidly downwards whereas porosity increases downwards near the base of the unit (Goldberg et al. 1994; Planke 1994; Planke & Flovenz 1996). The sediments interbedding the flow units show very low resistivity and very high gamma counts (Nielsen et al. 1984; Singh 1996) and intermediate porosity and density values (Planke 1994). Nielsen et al. (1984) showed that different series with- in a large vertical section of basalt flows could be identi- fied primarily from the results of the natural gamma ray log. Sanyal et al. (1980) showed that lithology could be discriminated by analysing wire-line logs and making cross- plots of especially the gamma-ray and neutron responses. Some of the results are of special interest in the context of the present work on the Lopra-1 logs. A positive corre- lation between flow thickness and mean velocity simply reflects increased proportions of massive basalt in thicker flows (Planke 1994); velocity and density gradients at the top and bottom of a flow unit are largely independent of flow thickness (Planke & Flovenz 1996); Goldberg et al. (1994) and Kern & Richter (1979) found that travel time correlates directly with porosity and showed that it also decreases with depth in sediments. Planke (1994) found flows between 0.6 and 18.5 m thick while most sediment was less than 1 m thick. Radioactive gamma ray and po- tassium logs are asymmetrical, having high values at flow tops, decrease gradually towards a central wide minimum and increase rapidly again near the flow base. Principles of interpretation of the Lopra-1 well wire-line data It is assumed that the principles and observations about the wire-line interpretation of subaerial basalt successions discussed in the previous section can be applied to the interpretation of the wire-line logging suite carried out in the Lopra-1 well. The physical properties of a basalt flow are related primarily to active mechanical processes (i.e. growth through inflation, cooling history, the environ- ment into which the flow was extruded, stress, faulting and whether it is a pahoehoe or aa flow type) and to secon- dary processes such as leaching and mineral growth. Thus it is expected that, within subaerially extruded flows, the responses of the sonic, porosity, density, resistivity and calliper logs, all of which record physical properties, will follow a generalised cyclic pattern. On the other hand the geochemical properties reflect a combination of magma composition, primary mineralogy, weathering, leaching of various elements and growth of secondary minerals so the response from the gamma (SGR and CGR), thorium, potassium and uranium wire-line logs, all of which repre- sent geochemical properties, will not necessarily show simi- lar cyclic patterns. As a first approach, an idealised cross-section of a flow is envisaged as follows: a porous crust at the top of the flow grades into a massive core. A basal zone below the massive core appears abruptly and is located on top of the preceding tuff/sedimentary layer or, in case such a layer is absent, on top of the crust of the preceding flow. The porous crust represents a transition zone between the tuff/sediment layer above it and the massive core ben- eath it. Its physical properties thus change gradually, rep- resented by the readings of the physical logs changing grad- ually downwards from those observed at the very top of the flow or in the overlying tuff/sediment layer to the val- ues of the underlying massive core. It is thus expected that values for the transit time of both the P- and S- waves, porosity and the diameter of the well will decrease down- wards, whereas the density and resistivity will increase downwards. The interpretation of the geochemical logs is ambiguous because the crust could either be leached or enriched in elements compared to the original composi- tion of the flow. The massive core is often found to be characterised by uniform values of the sonic, porosity, density, resistivity and calliper logs. More specifically, the core is recognised by the low transit time (high velocity) of the P-, S- and Stonely-waves, low values of porosity, the well diameter and drilling rate, but high values of density and resistivity. The chemical logs are expected to show a concentration GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1912 13 for the core and crust that more or less reflect the original chemical composition of the flow. The thin basal zone is recognised physically as a sharp transition above the underlying lava flow/tuff horizon. An abrupt change downwards of the log patterns is expected. Values of transit time of P-, S- and Stonely-waves are ex- pected to increase downwards (i.e. a downward fall in ve- locities) as are those of porosity and well diameter, but resistivity and density should decrease downwards. Sediment deposited on top of a lava flow should appear on the petrophysical logs as follows: high values of transit times of P-, S- and Stonely-wave sonic logs (low velocity), high values of porosity, low values of density, resistivity, distinctly high values of the gamma ray, especially the spec- tral gamma ray responses (particularly for the compen- sated gamma ray and potassium curves) and enlarged dia- meter of the well. Lopra-1 results The wire-line-logs are shown on Enclosure 1. Individual curves are arranged such that their visual impression can be correlated easily with the inferred lithology. Inspection shows that while the basal zone of many flow units can be recognised, it is too thin to be treated statistically as has been done for the core, crust, compound flows, dolerite dykes and sediments. Analysis of the basal zone has there- fore been included in that of the core. Calliper The calliper log shows that the diameter of the borehole varies cyclically in the range of 9–18″ in the upper part of the hole that was drilled in 1981 with an 8.5″ bit (Enclo- sure 1). The diameter is within the interval 6.5″ to 13″ in the deepened hole drilled with a 6.5″ bit. The diameter tends to decrease with depth in the upper hole but it is close to being gauge only in a few places. The measured calliper was expected to have been close to gauge in size; that the final diameter of the hole is irregular is probably due to the method of drilling (Nielsen et al. 1984). In the lower section drilled in 1996 with a different technique, the calliper is much closer to the bit size within cores (e.g. flow units 83, 86). The hole is narrowest in massive basalt and is usually wider in flow crusts and basal zones and is usually widest of all in sediment. In a few places, the larg- est diameter is found in porous crust rather than in over- lying sediment, possibly indicating that the sediment is more compact or less friable than the lava crust. Where the sediment beds are closely spaced, they do not seem to affect the condition of the hole, suggesting that the pres- ence of nearby basalt stabilises the hole and that the sedi- ment beds are rather competent (Enclosure 1). In places, the massive cores of the flows show minor deviations from stable hole condition that might indicate that the core is fractured or not as massive as presumed (e.g. flow units 42, 46 and 58). P-, S- and Stonely-wave transit time The log curves in Enclosure 1 have been plotted at the same scale but using different intersects, so that S-transit times are within the range 50–170 µsec/ft and P-transit times within 10–130 µsec/ft. Plotting the curves in this way ensures that they plot on top of one another in mas- sive cores where the hole is at or close to gauge as inferred from the combined interpretation of P-, S- and Stonely transit time, RHOB and NPHI in combination. The steady relationship between the two transit times probably indi- cates that the physical properties of the cores are fairly stable with a P- to S-wave velocity ratio of 1.31. Both P- and S-wave transit time logs show a cyclic pat- tern (Enclosure 1). The sediment beds are characterised by the highest transit time (70–100 µsec/ft sonic-P) fol- lowed by porous crusts and massive cores. Maximum sep- aration of the curves is a strong indicator of sediments, and intermediate separation is considered to represent crust. The velocity gradient in the upper part of flow units 24, 26, 28, 40, 41 and 42 seems to be independent of the thickness of the flow. Transit times within the core units seem to be independent of the thickness of the flow (e.g. flow units 21, 29, 52, 59 and 79) (Fig. 1). The Stonely-wave transit time is a part of the full-wave sonic survey. Attenuation of the Stonely-waves occurs where fractures are open and may also be caused by changes in permeability. Such attenuation will be more marked in hard formations where the acoustic contrast is greatest between the formation and mud-filled voids (Rider 1996). The cyclic pattern outlined by the P- and S-wave transit time logs is also shown on the Stonely-wave curve, and the smallest amount of oscillation of the curve occurs in the dolerite dykes and in the core of the flow units (En- closure 1). Density-porosity The bulk density and NPHI porosity curves (measured in limestone porosity units) also show a cyclic nature. The GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1913 14 highest ‘porosity’ and lowest density values are recorded within the sediments, intermediate values in lava crusts and lowest ‘porosity’ and highest density values in the lava cores (Enclosure 1). The curves have been plotted at com- parable scales, which have been chosen so that they over- lap within a typical massive basalt flow, which has been ascribed a density of 2.85–3.00 g/cm3 and a porosity (NPHI) of 0.1–0.15. Plotted in this way, a positive sepa- ration (where the porosity curve is located to the left of the density curve) is diagnostic of sediment or lava crust whereas no separation of the curves is diagnostic of lava core. Negative separations occur within lava cores that have a lower than average porosity. Resistivity logs The shallow (lls) and deep (lld) resistivity curves were ob- tained during the first drilling operations in the Lopra-1 well. The log pattern shows a cyclic behaviour compara- ble to the other logs representing physical properties, re- flecting the separation of flow units into crust, core and basal zones. The resistivity values are highest in the core of the flow units where the lld-log reads higher than the lls-log, showing an increase in resistivity away from the hole and into the formation. That could be due to either the formation water being fresher than the mud water and/or that the basalt becomes more massive away from the well. The values in the flow unit crusts are lower than those in the cores, and the separation between the curves is small. The composite flows have a behaviour similar to that of the crust. Sediments have lowest values but the lld-values are still larger than those of the lls-log. Dolerite dykes have the largest values of resistivity, and a large sep- aration of the logs is found. Geochemical logs Five geochemical logs are presented: SGR and CGR gam- ma ray, potassium, thorium and uranium (Enclosure 1). Basalt that contains both consistently high and consist- ently low amounts of radioactive material is seen. The shifts from high to low radioactive basalt are remarkably abrupt, and there seems to be no gradual changes, especially for thorium, indicating that thorium represents the original magma (e.g. flow units 7–8, 12–13, 22–23, 23–24, 30– 31, 31–32, 46–47, 55–56, Enclosure 1). In general terms, high thorium values are associated with high potassium and thus high CGR gamma ray, as well as high uranium values. There is a close relationship between the CGR gam- ma ray and potassium spikes and a less pronounced one between CGR gamma ray and uranium spikes. Sediment is characterised by high localised potassium spikes, often recognised in the CGR gamma ray, that correlate to anom- alies in the physical properties logs. In some flow unit Flow unit 21 Flow unit 29 Flow unit 52 Flow unit 59 Flow unit 75 P-sonic Cali_merged 140 40 8 23 NPHI 0 0.5 2.2 RHOB 3.2 CGR 0 30 THOR 30 10 15 20 25 30 35 5 0 5 10 M fr om c ru st /c or e bo un da ry C ru st C or e C al ip er P- so ni c Fig. 1. Variation with depth in five selected flow units. The bimodal behav- iour within the flow units is found in the logs representing physical properties but not in the geochemical logs. Variations with depth in the core are most pro- nounced for the NPHI porosity and RHOB density as compared to the P- sonic transit time and calliper. The P- sonic transit time appears rather constant. The CGR gamma-ray and thorium curves show the large variations present between different flow units. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1914 15 crusts, a decrease in potassium values with depth was found. Lithofacies analysis The cyclicity of the wire-line log curves has been used to distinguish a number of lithological units including (1) sedimentary layers, (2) flow units of simple type composed of crust, core and basal zone, (3) groups of thin flow units presumably defining composite flows in the sense of Walker (1991), (4) dolerite dykes. Fifty-two units have been interpreted to be sedimenta- ry beds, identified from a combined interpretation of the potassium, SGR and CGR gamma ray, sonic P-, S- and Stonely-waves, NPHI porosity, RHOB density, calliper and resistivity logs (Enclosure 1). The sediment layers are rec- ognised by gamma ray and potassium spikes relative to the overlying core (and basal zone) and the underlying crust. The gamma ray spike has to be larger than 25 GAPI and be localised within a short thickness interval less than 3.0 m (98% of the sediment layers; Fig. 2) and the pro- nounced potassium spike must contain more than 1% K. The uranium values are often high too, but are less diag- nostic. In addition, the calliper shows a localised maxi- mum deviation from the bit size, the sonic and NPHI logs record high values whereas the density is low. The thick- nesses of the sedimentary layers are generally small, and most of the bed thicknesses fall within the intervals 1.0– 1.5 m (38%) and 0.5–1.0 m (25%) (Fig. 2). The total thickness of sediments is estimated to be 75.8 m as com- pared to 40 ± 4 m reported by Nielsen et al. (1984). Within the depth intervals displayed in Enclosure 1, the NPHI- porosity and density curves often show a large negative separation. The NPHI-porosity index is large, but the NPHI-porosity often shows a higher value within the cru- stal part of the underlying flow unit. The density value is low, and generally the lowest value is found in the sedi- ments rather than in the uppermost part of the underly- ing flow unit (e.g. flow units 12, 20, 29 and 61, Enclo- sure 1). There are, however, examples where the density values are lower within flow units than within sediments. Thus porosity and density alone are not diagnostic of sedi- ment beds. The P-, S- and Stonely-wave transit times are high within small depth intervals. The shifted transit time curves of the P- and S-waves show large negative separa- tion, but the values are not diagnostic enough to separate sediment layers from the porous crust of the lava flow. The resistivity values are extremely low. Hald & Waag- stein (1984) and Waagstein et al. (1982) recognised a to- tal of 36 sedimentary beds from analysis of cuttings sam- pled every 2 m and, among these, 32 are located below the casing of the well and thus within the 1996 log runs. It has been possible to establish a correlation between the sediment layers interpreted from cuttings and the 1996 log curves (Enclosure 1), and 27 out of the 32 sedimenta- ry beds can be confirmed by the criteria put forward here. Thus 22 additional sedimentary layers have been identi- fied within the shallower part of the Lopra-1 well. A total of 87 flow units made up of crust, core and basal zone have been identified using a combination of P-, S- and Stonely-wave transit times, porosity and density values, resistivity values and the calliper log. From the top of each unit downwards to the core, the lava crust is indi- cated by a general gentle decrease in the sonic transit times combined with a decrease in porosity and increase in den- sity and resistivity (Enclosure 1). Below this crustal zone, the log patterns are replaced by fairly constant low values of transit times and porosity combined with a high value of the density and resistivity, defining the lava core. The core is characterised by the calliper being closely in gauge and low values of the shifted sonic transit time curves, low and small oscillations of the Stonely-transit time curve, no separation of the density and porosity curves, and high values of the resistivity values (Enclosure 1). This fairly uniform behaviour of the log curves is abruptly different near the base of the unit where a pronounced increase in sonic transit time, porosity and calliper values combined with a pronounced decrease in density and resistivity val- ues marks the basal zone. The thickness of the basal zone is generally small, and it has not been possible to differen- tiate it from the overlying core and the underlying sedi- mentary layer/crust in a consistent manner, so here the 0 5 10 15 20 Fr eq ue nc y (% ) 1.0 2.0 3.00 4.0 Sediment thickness in metres Fig. 2. Frequency diagram illustrating the distribution of sedi- ment thicknesses interpreted in the present study. The sediment layers have been grouped into 0.5 m thickness intervals. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1915 16 basal zones have been included in the cores. The behav- iour of the log curves representing the basal zone seems independent of whether or not a sediment layer is present beneath it. The individual flow units are generally reflected in the geochemical logs, primarily in the thorium curve and sec- ondarily in the potassium curve, whereas the SGR, CGR gamma and uranium curves show unsystematic variations (e.g. flow units 6, 29, 39, 44, 51, 52, Enclosure 1). The unsystematic variations do not allow an overall profile to be established as can be done with the logs representing the physical properties. It is therefore concluded that the geochemical behaviour of the flow units is largely related to the original composition of the magma and that differ- ent flow fields may be present. The thicknesses of the flow units have been compared to those interpreted by Hald & Waagstein (1984) and are found typically to be within the range 5–25 m (of which interval flow units of a thick- ness of 10–15 m represent 23%; Fig. 3). A closer inspection of the curves through the crust and core parts of individual flow units reveals large fluctua- tions indicating that the physical properties are not con- stant in either the crust or the core. The nature of the variations in one crust is illustrated in a few examples; in flow unit 18 the porosity increases with depth and the density is rather constant, in flow unit 46 the porosity increases with depth and density decreases, in flow unit 20 the porosity decreases with depth and density increas- es. Although the porosity and density figures for the crust differ, the relationship between the crust and core seems to be constant apart from those flow units where the crust and core are of similar thickness. The idealised core as a massive part of the flow unit characterised by constant low values of transit time, low values of porosity and high density and resistivity values is only true for a few flow units, e.g. 20 and 30, whereas different situations are found in e.g. flow units 13, 19, 26, 30, 36. Most of the flow units penetrated by the Lo- pra-1/1A well are characterised by the presence of porosi- ty that may be vesicular zones or fractures within the core, though their regularity seems to indicate that they can be ascribed to vesicular zones. These zones are clearly reflect- ed by the density and porosity curves and deflect the S- sonic transit time curves to a larger degree than they do the P-sonic transit time curves. The vesicular zones with- in the core are not found at a consistent height above the base or below the top, but most of them seem to display some kind of symmetry. The variation from core to core is exemplified by e.g. flow units 26 and 40 which are most massive at the top and porous in the middle; the core of flow unit 41 is most massive in the upper part, and the vesicular zone is displaced towards the base as compared to the core of flow unit 40, the core of flow unit 58 is almost symmetrical with a porous zone in the middle, the cores of flow units 31 and 52 become more massive to- wards the base, and the core of flow unit 59 shows nu- merous vesicular zones. A minority of cores are massive in the middle and then are less massive both upwards and downwards (e.g. flow units 45 and 48). Statistical analysis shows that when vesicular zones are present, the core is at least 70% thicker than the crust, whereas the core and crust are of equal thickness when no vesicular zones are indicated by the logs. Five flow units (21, 29, 52, 59 and 75), representative of flow units of different thickness between 10 and 50 m, have been chosen to illustrate the relationship between flow thickness and mainly velocity. The thicknesses were correlated with the P-wave transit time, the calliper, the NPHI-porosity, density, gamma ray (CGR) and the thor- ium content logs (Fig. 1). It can be seen that different values of physical properties are recorded in the flow unit depending on whether crust or core is encountered. The flow units vary considerably in their geochemical expres- sion, and the bimodal behaviour of the flow units is not seen. The crust is represented by high values of the calli- per, transit time, NPHI-porosity and low values of densi- ty. Transit time in the core seems to be fairly constant and therefore independent of the thickness of the core and the flow unit. The variations in values recorded by the calli- per seem independent of the thickness of the flow unit (Fig. 1 and Enclosure 1). The relationship between the flow unit thickness and NPHI-porosity and density val- ues of the core seems to suggest that thin flows exhibit high porosity and low density whereas thick flows show 0 5 10 15 20 Fr eq ue nc y (% ) This study Hald & Waagstein (1984) 10 20 30 40 50 600 Flow unit thickness in metres Fig. 3. Frequency diagram illustrating the distribution of flow unit thicknesses interpreted in the present study and in Hald & Waag- stein (1984). The flow units have been grouped into 5 m inter- vals. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1916 17 Crossplo flow-unit Density/Soni P/dept Crossplo compoun flow Density/Soni P/dept Crossplo crust Density/Soni P/deptCrossplo core Density/Son c /dep Crossplo Dolerite Density/Soni P/dept Crossplo sediment Density/Soni Cross-plot flow unit: density/P-sonic/depth Cross-plot compound flow: density/P-sonic/depth Cross-plot dolerite: density/P-sonic/depth Cross-plot sediment: density/P-sonic/depth Cross-plot core: density/P-sonic/depth Cross-plot crust: density/P-sonic/depth 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 D en si ty 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 D en si ty D en si ty 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 D en si ty 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 D en si ty 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 D en si ty 50 70 90 110 130 50 70 90 110 130 50 70 90 110 130 P-sonic 50 70 90 110 130 P-sonic 50 70 90 110 130 P-sonic 50 70 90 110 130 P-sonic P-sonic P-sonic 2600 2200 1800 1400 1000 600 200 Depth MD Fig. 4. Cross-plot of inferred lithology: the cross-plots present the statistical behaviour of the flow units, cores, crusts, compound flows, dolerite dykes and sediments plotted for P-sonic transit time versus RHOB density and depth. It is found that the clustering of data is enhanced by plotting the inferred lithology instead of the entire basalt column. 100 80 60 40 20 0 Sonic (µs/ft) Density Flow unit Massive Porous Compound Dolerite Sediment 40 140 40 140 40 140 40 140 40 140 40 140 1.50 3.30 1.50 3.301.55 3.151.50 3.301.50 3.301.55 3.15 Avg: 55.127 Min: 40.817 Med: 52.79 Max: 103.723 Avg: 52.283 Min: 40.817 Med: 51.676 Max: 92.807 Avg: 65.0 Min: 45.989 Med: 64.4 Max: 103.723 Avg: 61.919 Min: 46.422 Med: 60.676 Max: 125.3853 Avg: 50.098 Min: 40.227 Med: 48.421 Max: 90.512 Avg: 71.606 Min: 48.286 Med: 70.878 Max: 106.435 Avg: 2.845 Min: 1.554 Med: 2.906 Max: 3.151 Avg: 2.900 Min: 1.566 Med: 2.931 Max: 3.151 Avg: 2.658 Min: 1.554 Med: 2.664 Max: 3.151 Avg: 2.739 Min: 2.033 Med: 2.757 Max: 3.011 Avg: 2.900 Min: 1.510 Med: 2.946 Max: 3.035 Avg: 2.582 Min: 1.529 Med: 2.582 Max: 2.949 100 80 60 40 20 0 Fr eq ue nc y (% ) Fr eq ue nc y (% ) Fig. 5. Histogram and statistical values of the six cross-plots presenting the inferred lithology for P-sonic transit time versus density and depth shown in Fig. 4. The average, median, minimum and maximum values are shown. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1917 18 low porosity and high density (Fig. 1) but Enclosure 1 shows that this relation is too simple. The term ‘compound flow’ is used for intervals where the cyclicity of the logs representing physical properties suggests strongly that flow units are present, but the indi- vidual parts of the flows are poorly defined and log read- ings do not reach the values normally ascribed to crust or core (Enclosure 1). Eighteen compound flow intervals have been recognised. Dolerite dykes have previously been encountered at two levels in the Lopra-1 well (Hald & Waagstein 1984). The dykes are represented by constant values of gamma ray, thorium and potassium lower than found in flow units, rather high density and the lowest values of porosity and transit time combined with the highest values of the resis- tivity measured in the well (Enclosure 1). The calliper shows that pronounced caving occurred in the upper dyke, but no such caving occurred in the lower one. The dyke intervals were expected to show rather constant values of the log curves representing the physical properties but, as is seen on Enclosure 1, a cyclic appearance of those logs is seen. This cyclic appearance is most obviously seen in the lower dolerite dyke where five symmetric intervals are re- vealed by the RHOB, NPHI, lls and lld logs. This sug- gests that the lower dolerite dyke consists of five thinner dykes. Statistics The lithological division of the volcanic sequence into flow units (consisting of porous crust and massive core), com- pound flows, dolerite dykes and sediments makes possible a statistical analysis of the properties of the various litho- logies found in the Lopra-1 well. The data obtained from especially the physical properties logs as exemplified by the values of P-, S- and Stonely-wave, porosity, density and resistivity logs (Table 3) show a much reduced scatter after being separated into the lithological divisions. A strong correlation, most pronounced for the core and to a lesser extent for the crust and sediment beds, is also found between P- and S-wave transit times, P- and S-wave tran- sit times versus density, P- and S-wave transit times versus porosity, and density versus porosity. Examples are pre- sented in Figs 4 and 5 where six cross-plots with associat- ed histograms and statistical parameters show the benefit 200 800 1400 2000 40 60 90 120 140 200 800 1400 2000 40 60 90 120 140 200 800 1400 2000 40 60 90 120 140 200 800 1400 2000 2600 2600 D ep th ( m ) 40 60 90 120 140 Flow units Composite flows Dolerites Sediments Cores Crusts 200 800 1400 2000 40 60 90 120 140 800 1400 2000 40 60 90 120 140 2600 200 2600 2600 2600 D ep th ( m ) D ep th ( m ) D ep th ( m ) D ep th ( m ) D ep th ( m ) Sonic (µs/ft) Sonic (µs/ft) Sonic (µs/ft) Sonic (µs/ft) Sonic (µs/ft) Sonic (µs/ft) Fig. 6. Cross-plot of P-sonic transit time versus depth for the inferred lithology. The figure shows that the core and sediment lithologies especially show a pronounced trend of decreasing transit time values with depth whereas crust and compound flows are ambiguous. The amount of data representing dolerite dykes is insufficient to describe their relationship with depth. The flow units show a strong increase in transit time with depth but this is due to the core part of the flow unit. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1918 19 of dividing the geological column into the suggested litho- logical units. The relationship between the physical properties of the lithological units and depth has been examined and a cor- relation is observed. The degree of correlation is strongest for the flow units (due to the core values), core, sediment/ tuff layers, dolerite dykes and to a lesser extent for crust and compound flows. P-sonic transit time curves are il- lustrated versus depth for the different lithologies as an example (Fig. 6). There it is seen that especially the sonic transit time increases with depth for the flow units and a first order regression line has been estimated as P-sonic = 58.1370–0.0022 × depth measured in metres or P-sonic = 58.0934–0.0007 × depth when depth is measured in feet. Likewise a first order regression line estimated for the cores is P-sonic = 54.0165–0.0012 × depth in metres and P-sonic = 54.7769–0.0011 × depth for depths measured in feet. This observation is in accordance with general findings made by Kern & Richter (1979). The sediment/tuff layers also show a significant corre- lation with depth and in general terms the crust does too, whereas the compound flows show no systematic varia- tion. A depth relation for dolerite dykes is not justified because of their limited numbers. The potassium and CGR gamma ray logs show a strong response to the sediment beds, but there is otherwise no correlation between the geochemical logs and the indi- vidual flow units. This probably indicates that the geo- chemical values are related to the original composition of the flow units, which varies strongly, suggesting that one magma type may represent more than one flow unit. This can be illustrated by analysing the thorium content in rela- tion to flow units and compound flows as defined by the petrophysical logs. The intervals covering the dolerite dykes are not included in this part of the study as they are repor- ted to be later intrusions (Hald & Waagstein 1984). The thorium log data have been block averaged into three groups according to their thorium content (group 1 con- sisting of 72 flow and compound flow units with a thori- um content of 0.0–1.0 ppm; group 2 containing 21 flow and compound flow units with 1.0–2.0 ppm Th; group 3 with 11 flow and compound flow units representing a content of 2.0–3.0 ppm Th). This results in 36 groups, tentatively interpreted as flow fields (Enclosure 1). Some flow fields are constructed of a number of flows and com- pound flow units, but 57% of the flow fields constitute one flow or compound flow unit (Fig. 7). The individual flow fields are separated by sediment/tuff layers in 70% of cases. Since this study is an analysis of the basalt column based primarily on a lithological division from petrophys- ical logs, further statistical treatment of the geochemical logs is outside its scope and probably should be done by dividing the lava flows into those with high, medium and low radioactivity. Discussion The flood basalt pile in the interval 200–2486 m depth in the Lopra-1 well has been divided into lithological units based on the wire-line logging (Enclosure 1). Eighty-six flow units consisting of porous crust and massive core (in- cluding basal zones), 18 compound flows and two dolerite dykes have been recognised using petrophysical logs (transit time of P-, S- and Stonely-waves, NPHI-porosity, density, resistivity and calliper). Fifty-two sedimentary/tuff layers have been identified mainly by potassium, CGS and SGR gamma logs supplemented by the petrophysical logs. Thir- ty-seven flow fields have been suggested based on the tho- rium log in combination with the petrophysical logs. The flow units have been subdivided into three parts on criteria proposed by various workers based on field evi- dence or wire-line logs from basalt covered regions (e.g. Planke 1994; Planke & Flovenz 1996; Waagstein 1999). The Lopra-1 results reveal that overall the crust, core and basal zones are similar to other basalt covered areas re- ported using wire-line log studies, but there are some dif- ferences. To describe the crusts as showing a downward increase in density and resistivity, a decrease in porosity and P- and S-wave transit times is too simple. The down- ward increase in density values is commonly replaced by zones of decreasing values, probably reflecting the pres- ence of vesicular or fractured zones. Likewise the previ- ously reported behaviour of the core as being rather uni- form although most massive in its middle part (e.g. Planke 1994) is far from true for the flows of the Lopra-1 well. Number of flow units and compound flows in flow field 11 5 1 2 3 5 64 11 14 20 Fr eq ue nc y (% ) Flow field with 0–1 ppm Th Flow field with 1–2 ppm Th Flow field with 2–3 ppm Th Fig. 7. Frequency diagram illustrating the distribution of the number of flow units or compound flows versus thorium content included in the inferred flow fields. The flow units and compound flows have been grouped into three depending on their thorium content. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1919 20 Two types of core are commonly present, one type with a ‘massive’ central part where the logs indicate that only the central part is truly massive. This type, which represents 12% of the flow units, is generally thin, and the crust is of the same thickness as the core. The other type is charac- terised by varying numbers of zones of high porosity/low density whose relative depth within the core also varies. In these flow units the thickness of the core is at least 70% of the total thickness of the flow unit, and the flow units themselves are of different thicknesses. This zonation is interpreted as vesicular or perhaps occasionally fractured zones. The presence of vesicular zones has been reported by McMillan et al. (1987) and Walker (1989). Walker (1989) described the vesicular zonation from spongy pahoe- hoe lava in Hawaii and mentioned that aa-type lava can also contain vesicular zones. This raises the question wheth- er cores of pahoehoe lavas can be distinguished from those of aa-type lava. Most papers distinguish aa from pahoe- hoe using surface morphology and perhaps the vertical thickness of the basal zone (Waagstein 1999). Unfortu- nately surface morphology is not revealed by the present logging suite, and the vertical thicknesses of the basal zone are too small to be treated on their own in the present study. However, Hald & Waagstein (1984) reported that the majority of basalt from the shallower part of the Lo- pra-1 well was aa-type lava and that pahoehoe lava in the field seems to have a similar crust and core relationship. It is therefore tentatively suggested that flow units with ves- icles may represent aa-type flows, and that the flow units where the crust and core are of the same thickness and the core is massive are pahoehoe type. However only a few studies have dealt with the question of how a flow unit is formed in relation to its distribution of vesicles and cooling history based on cross-sections. The present study has identified 86 flow units and 18 compound flows compared with the 87 reported by Waag- stein et al. (1982) and Hald & Waagstein (1984). Flow thicknesses of 10–15 m (23%) are most common followed by thicknesses of 15–20 m (18%) (Fig. 3). Thus the flow unit thicknesses reported here are slightly higher than those of Hald & Waagstein (1984). Five flow units of different thickness – tentatively in- terpreted as aa-type – were chosen by using the transit time of the P-wave, calliper, NPHI-porosity, density, thorium and gamma ray (CGR) logs to analyse the thickness of the flow units (Fig. 1). The petrophysical logs show clear- ly that different values are ascribed to the crust and core. The expected overall trend of the flow unit in relation to the petrophysical logs is valid. The P-wave transit time in the core is rather constant and not dependent on the thick- ness of the core and thus of the flow unit. In contrast to the transit time, the NPHI-porosity and density varies, and in the flow units chosen the density seems to increase and the porosity decrease with thickness of the core and thus of the flow unit. That the transit time is independent of the thickness of the core and thus of the flow unit is somewhat different from the results presented by Planke (1994) from a seaward dipping reflector sequence. The behaviour of the radioactive gamma ray and potassium logs in Lopra-1 also differs as they do not show the cyclic- ity found in the seaward dipping reflectors studied by Planke (1994). Whether this indicates different deposi- tional settings, lava types, magma type, weathering con- ditions or the fact that fresh water is encountered in the Lopra-1 well is uncertain. Eighteen compound flows, characterised by similar values of the logs reflecting the physical properties as the crust, have been interpreted. However a cyclic behaviour similar to that in the flow units can be recognised within the compound flows although the values of the P- and S- wave transit time, density, porosity and resistivity are low- er than in the flow units. It is suggested that the com- pound flows represent a number of thin flows or the toes of flow units. Two dolerite dykes were encountered and are charac- terised by the lowest values of SGR and CGR gamma ray, P- and S-wave transit time and porosity combined with the highest level of resistivity, high values of density and a constant level of the Stonely-wave transit time curve. The sonic P- and S-transit time logs and the density, porosity and resistivity curves show a cyclic behaviour resembling the pattern of a flow unit, and it is speculated that the dyke intervals represent a group of thinner dykes. Fifty-two sedimentary layers have been detected by the combination of potassium, CGR and SGR gamma log, transit time logs of P-, S- and Stonely-waves, NPHI-po- rosity, density, resistivity and calliper logs. The sediments have been identified on the basis of a locally high spike of the potassium log value in combination with a gamma ray spike larger than 25 GAPI, high values of transit time of P-, S- and Stonely-waves, high values of the NPHI- porosity and low values of the density and resistivity. The number of sedimentary layers is greater than reported by Hald & Waagstein (1984) and Nielsen et al. (1984) who reported 27 sediment layers based on the combination of cuttings and geophysical logs acquired in 1981. Of the 27 sediment layers identified by Hald & Waagstein (1984), the recently acquired geophysical logs used in the present study have confirmed 22 of them. The sedimentary layers are thought to have formed during periods of quiescent volcanism. In that case, at least 51 periods of quiescence occurred in the interval from GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1920 21 184 to 2484 m depth in the Lopra area. Mostly a single flow unit was extruded between each pause, but up to eight flows are found within one interval. The enlarged number of sediment layers indicates that the time inter- vals between eruptions taking place were shorter than pre- viously thought and thus the extrusion rate of basalt was faster than previously assumed. The amount of basalt pro- duced between the periods of quiescence, measured by the thickness of the flow unit, differs markedly at differ- ent depths in the Lopra-1 well. This may indicate that the thickness of basalt produced by a single eruption differed, but to confirm this the location and extension of each eruption needs to be known. The geochemical logs, espe- cially the thorium log, show numerous fluctuations indi- cating the presence of some highly radioactive basalt flows, similar to those reported in the Deccan Trap area by Buck- ley & Oliver (1990). To try to systematise these fluctua- tions three intervals have been defined (0–1 ppm, 1–2 ppm and 2–3 ppm), and the column was subdivided at the same time into smaller units that appear to correlate with the flow units and sedimentary beds, and on this basis it is suggested that a minimum of 36 units exists. This suggests that 36 different compositions of magma have been present during the time represented by the Lo- pra-1 well, and tentatively it is suggested that a minimum of 36 flow fields exists. Table 2 lists results of other studies of subaerially ex- truded basalts whereas Table 3 contains the statistical val- ues of the Lopra-1 well wire-line log data from this study. The Lopra-1 results, which are more detailed due to the subdivision into the suggested lithology units, the large thickness of basalt and the extensive logging suite, are in overall agreement with the range of results obtained from the literature survey. Some scattering in measured values exists; however, the amount of data is not sufficient to conclude whether it is the primary volcanic setting or post- depositional processes that are responsible for these dif- ferences. At a depth of approximately 2417 m, the cyclic nature of the wire-line log curves representing the physical prop- erties of the basalt terminates abruptly, and an increase in transit times and porosity combined with a decrease in density signal a transition into the non-subaerially extruded basalt sequence below this level. The transition generates a negative seismic reflection from this level. It is conclu- ded that the negative reflection interpreted from the VSP by Kiørboe & Pedersen (1995) corresponds to the base of the subaerially extruded basalt. The total known thick- ness of the flood basalt on the Faroe Islands can therefore be increased by 250 m. Conclusions The study shows that it is possible to use wire-line log- ging combined with results from work carried out in flood basalt regions in other parts of the world to divide the subaerial plateau basalts drilled in the Lopra-1 well into a number of lithological units. Fifty-two sedimentary lay- ers have been identified using geochemical logs supple- mented by the logs representing physical properties. Eighty-six basaltic flow units (consisting of a vesicular crust, a massive or partly vesicular core and a thin basal zone), 18 compound flows and two dolerite dykes have been identified using the logs representing physical prop- erties supplemented by the geochemical logs. Tentative aa and pahoehoe lava types are suggested, among which the aa lava is the most common. The geochemical logs appear to respond to the composition of the magma rather than reflecting individual flow units, and it is suggested that the Lopra-1 well penetrated a minimum of 36 flow fields. The present subdivision of the lithology is strongly sup- ported by statistical analysis of the logged physical pa- rameters and shows that a pronounced reduction in data scatter is obtained when the basalts are studied as individ- ual units. Thus the study presents consistent and detailed information about the lower series of the Faroe plateau basalts that has led both to new results and results in greater detail than was previous known, which is of benefit in geophysical and geological studies of basalt and pre-basalt geology. Acknowledgements Thanks are due to the Lopra consortium and the Geolo- gical Survey of Denmark and Greenland for access to the wire-line log data. Thanks are due to Dr. H. Delius for assistance in correcting the resistivity logs. Thanks are due the reviewers Professor Dr. D. Goldberg and Professor Dr. P.K. Harvey for valuable and constructive suggestions that have improved the content. J.A. Chalmers is thanked for valuable suggestions that improved the language. Land- mark is acknowledged for the University Grant issued to the Geological Institute, University of Copenhagen. References Berlitz, R. et al. 1988: Log interpretation in igneous and metamor- phic rocks with five case studies. Technical Review 36, 30–47. Berthelsen, O., Noe-Nygaard, A. & Rasmussen, J. (eds) 1984: The deep drilling project 1980–1981 in the Faeroe Islands. Annales GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1921 22 Societatis Scientiarum Faeroensis, Supplementum IX, 158 pp. Tórshavn: Føroya Fróðskaparfelag. Buckley, D.K. & Oliver, D. 1990: Geophysical logging of water exploration boreholes in the Deccan Traps, central India. In: Hurst, A., Lovell, M.A & Morton, A.C. (eds): Geological appli- cation of wire-line logs. Geological Society Special Publication (London) 48, 153–161. Cashman, K.V. & Kauahikaua, J.P. 1997: Reevaluation of vesicle distribution in basaltic lava flows. Geology 25, 419–422. Christensen, N. 1996: Poisson’s ratio and crustal seismology. Jour- nal of Geophysical Research 101(B2), 3139–3156. Goldberg, D.S., Reynolds, D.J., Williams, C.F., Witte, W.K., Olsen, P.E. & Kent, D.V. 1994: Well logging results from the Newark Rift Basin Coring Project. Scientific Drilling 4, 267–279. Hald, N. & Waagstein, R. 1984: Lithology and chemistry of a 2 km-sequence of Lower Tertiary tholeiitic lavas drilled on Su- ðuroy, Faeroe Islands (Lopra-1). In: Berthelsen, O., Noe-Ny- gaard, A. & Rasmussen, J. (eds): The deep drilling project 1980– 1981 in the Faeroe Islands. Annales Societatis Scientiarum Fae- roensis, Supplementum IX, 15–38. Tórshavn: Føroya Fróðska- parfelag. Kern, H. & Richter, A. 1979: Compressional and shear wave ve- locities at high temperature and confining pressure in basalts from the Faeroe Islands. Tectonophysics 54, 231–252. Kiørboe, L.V. & Pedersen, S.A. 1995: VSP seismic experiment in the Faeroes. In: Scrutton, R.A. et al. (eds): The tectonics, sedi- mentation and palaeoceanography of the North Atlantic region. Geological Society Special Publication (London) 90, 111–123. Lund, J. 1983: Biostratigraphy of interbasaltic coals from the Fae- roe Islands. In: Bott, M.H.P. et al. (eds): Structure and develop- ment of the Greenland–Scotland Ridge. NATO Conference Series IV: Marine Science, 417–423. New York/London: Ple- num Press. McMillan, K., Randal, C.W. & Long P.E. 1987: Two-stage vesicu- lation in the Cohassett flow of the Grande Ronde Basalt, south- central Washington. Geology 15, 809–812. Nielsen, P.H., Stefánsson, V. & Tulinius, H. 1984: Geophysical logs from Lopra-1 and Vestmanna-1. In: Berthelsen, O., Noe- Nygaard, A. & Rasmussen, J. (eds): The deep drilling project 1980–1981 in the Faeroe Islands. Annales Societatis Scientiarum Faeroensis, Supplementum IX, 115–135. Tórshavn: Føroya Fróðskaparfelag. Planke, S. 1994: Geophysical response of flood basalts from analy- sis of wire-line logs. Ocean Drilling Program Site 642, Vøring volcanic margin. Journal of Geophysical Research 99(B5), 9279– 9296. Planke, S. & Flovenz, O.G. 1996: Seismic properties of flood ba- salt (extended abstract), 6 pp. Geophysics for lithology predic- tions. Kristiansand: Norwegian Petroleum Society. Rasmussen, J. & Noe-Nygaard, A. 1969: Beskrivelse til geologisk kort over Færøerne i målestok 1:50 000. Danmarks Geologiske Undersøgelse 1. Række 24, 370 pp. + map vol. (with summaries in Faeroese and English). Rasmussen, J. & Noe-Nygaard, A. 1970: Geology of the Faeroe Islands (pre-Quaternary). Danmarks Geologiske Undersøgelse 1. Række 25, 142 pp. Rasmussen, J. & Noe-Nygaard, A. 1990: The origin of the Faeroe islands in text, pictures and on maps, 64 pp., 6 maps at 1:50 000. Copenhagen: Geological Survey of Denmark (also text in Faeroese and Danish). Rider, M. 1996: The geological interpretation of well logs, 2nd edition, 256 pp. Houston, Texas: Gulf Publishing. Sanyal, S.K., Juprasert, S. & Jubasche, M. 1980: An evaluation of a rhyolite-basalt-volcanic ash sequence from well logs. Log Ana- lyst 21, 3–9. Self, S., Thordarson, T. & Keszthelyi, L. 1997: Emplacement of continental flood basalt lava flows. In: Mahonney, J.J. & Cof- fin, M.F (eds): Large igneous provinces: continental, oceanic and planetary flood volcanism. American Geophysical Union Geophysical Monograph 100, 381–410. Singh, S.C. 1996: Delineation of basaltic lava flows through geo- physical logging. Gondwana Geological Magazine Supplement- ed 2, 539 only. Thordarson, T. & Self, S. 1998: The Roza Member, Columbia Riv- er Basalt Group. A gigantic pahoehoe lava flow field formed by endogenous processes? Journal of Geophysical Research 103, 411– 445. Waagstein, R. 1999: A geological field guide to the Palaeogene flood basalts of Suðuroy, Faroe Islands. Danmarks og Grønlands Ge- ologiske Undersøgelse Rapport 1998/130, 46 pp. Waagstein, R., Nielsen, P.H., Fine, S. & Hald, N. 1982: Lopra well no. 1, Suðuroy, Faroe Islands. Geological well completion report, 40 pp. Unpublished report, Geological Survey of Den- mark, Copenhagen. Walker, G.P.L. 1989: Spongey pahoehoe in Hawaii: a study of ves- icle-distribution patterns in basalt and their significance. Bulle- tin of Volcanology 51, 199–209. Walker, G.P.L. 1991: Structure, origin by injection of lava under surface crust, of tumuli, ‘lava rises’, ‘lava-rise pits’, and ‘lava- inflation clefts’ in Hawaii. Bulletin of Volcanology 53, 546–558. Manuscript received 20 June 2000; revision accepted 11 September 2001. GEUS Bulletin no 9 - 7 juli.pmd 07-07-2006, 14:1922