Seismic and Petrophysical Relationships from UK Continental Shelf Rocks at Elevated Pressures

Abstract

There are few published seismic (P- and S-waves) properties for seafloor bedrocks. At low pressures (1 to 10 MPa), velocities and attenuations are determined mainly by open microcracks. At higher pressures, the microcracks close, and the velocities and attenuations depend primarily on the matrix porosity. We have investigated both the relationships between the acoustic, petrophysical, and geological properties of the rocks at 40 MPa pressure and the effect of microcracks on the acoustic properties at 10 MPa pressure. In this paper we discuss the former; the latter will be discussed separately.

P- and S-wave velocity and attenuation measurements were carried out on a suite of seabed sedimentary and igneous rocks at effective pressures from 10 to 40 MPa at ultrasonic frequencies. The porosities and permeabilities of the rocks ranged from 0% to 32% and 0 to 110 mDarcy, respectively. Characterization of the rocks revealed that most of the sandstones have a substantial clay content (kaolinite, illite, and chlorite) and fractures. Most of the igneous rocks are chloritized.

The seismic properties of the rocks are markedly lower than those of similar continental rocks because of the microporosity formed by the alteration of feldspars, micas, and mafic minerals to clays (e.g., chloritization of pyroxenes) and the corresponding reduction of the elastic moduli. The results of this study suggest that the values of velocities and quality factors used for ocean acoustic propagation models are lower than normally assumed.     

Assessing the Validity of Seismo-Acoustic Predictor Equations for Obtaining Seabed and Sub-Surface Sediment Physical Properties

The interdependence between the seismo-acoustic properties of a marine sediment and its geotechnical/physical parameters has been known for many years, and it has been postulated that this should allow the extraction of geotechnical information from seismic data. Though in the literature many correlations have been published for the surficial layer, there is a lack of information for greater sediment depths. In this article, a desktop study on a synthetic seafloor model illustrates how the application of published near-surface prediction equations to subsurface sediments (up to several tens of meters burial depth) can lead to spurious predictions. To test this further, acoustic and geotechnical properties were measured on a number of sediment core samples, some of which were subjected to loading in acoustically-equipped consolidation cells (oedometers) to simulate greater burial depth conditions. For low effective pressures (representing small burial depths extending to around 10 meters subsurface), the general applicability of established relationships was confirmed: the prediction of porosity, bulk density, and mean grain size from acoustic velocity and impedance appears generally possible for the investigated sedimentary environments. As effective pressure increases through, the observed relationships deviate more and more from the established ones for the near-surface area. For the samples tested in this study, in some instances increasing pressure even resulted in decreasing velocities. There are several possible explanations for this abnormal behavior, including the presence of gas, overconsolidation, or bimodal grain size distribution. The results indicate that an appropriate depth correction must be introduced into the published prediction equations in order to obtain reliable estimates of physical sediment properties for greater subsurface depths.     

Microfabric and rock properties of experimentally compressed silt-clay mixtures

Oedometric mechanical compaction tests were performed on brine-saturated synthetic samples consisting of silt-clay mixtures to study changes in microfabric and rock properties as a function of effective stress. The silt consisted of crushed quartz (∼100%) with grain size range between 4 and 40 μm, whereas the clay consisted of 81% kaolinite, 14% mica/illite and 5% microcline of grain size between 0.4 and 30 μm. Five sample pairs ranging in composition from pure silt to pure clay were compacted to 5 and 50 MPa effective stress respectively. SEM studies were carried out to investigate microfabric changes in the mechanically compacted silt-clay mixtures. The degree of alignment of the different minerals present (quartz, mica/illite and kaolinite) were computed by using an image analysis software. Experimental compaction have measured the changes in the rock properties such as porosity and velocity as a function of effective stress for different mixtures of clay and silt. Clay-rich samples showed a higher degree of mineral orientation and lower porosity compared to silt-dominated samples as a function of effective stress. Pure clay sample had 11% porosity at 50 MPa effective stress whereas the pure silt sample retained about 29% porosity at the same effective stress. The experiments showed that low porosity down to 11% is possible by mechanical compaction only. A systematic increase in strain was observed in the silt-clay mixtures with increasing clay content but the porosity values found for the 50:50 silt-clay mixture were lower than that of 25:75 silt-clay mixture. No preferential mineral orientation is expected before compaction owing to the high initial porosity suggesting that the final fabric is a direct result of the effective stress. Both P- and S-wave velocities increased in all silt-clay mixtures with increasing effective stress. The maximum P- and S-wave velocities were observed in the 25:75 silt-clay mixture whereas the minimum Vp and Vs were recorded in the pure silt mixture. At 50 MPa effective stress P- wave velocities as high as 3 km/s resulted from experimental mechanical compaction alone. The results show that fine-grained sediment porosity and velocity are dependent on microfabric, which in turn is a function of grain size distribution, particle shape, sediment composition and stress. At 5 MPa effective stress, quartz orientation increased as a function of the amount of clay indicating that clay facilitate rotation of angular quartz grains. Adding clay from 25% to 75% in the silt-clay mixtures at 50 MPa effective stress decreased the quartz alignment. The clay mineral orientation increased by increasing both the amount of clay and the effective stress, the mica/illite fabric alignment being systematically higher than that of kaolinite. Even small amount of silt (25%) added to pure clay reduced the degree of clay alignment significantly. This study demonstrates that experimental compaction of well characterized synthetic mudstones can be a useful tool to understand microfabric and rock properties of shallow natural mudstones where mechanical compaction is the dominant process.     

Assessing the Validity of Seismo-Acoustic Predictor Equations for Obtaining Seabed and Sub-Surface Sediment Physical Properties

The interdependence between the seismo-acoustic properties of a marine sediment and its geotechnical/physical parameters has been known for many years, and it has been postulated that this should allow the extraction of geotechnical information from seismic data. Though in the literature many correlations have been published for the surficial layer, there is a lack of information for greater sediment depths. In this article, a desktop study on a synthetic seafloor model illustrates how the application of published near-surface prediction equations to subsurface sediments (up to several tens of meters burial depth) can lead to spurious predictions. To test this further, acoustic and geotechnical properties were measured on a number of sediment core samples, some of which were subjected to loading in acoustically-equipped consolidation cells (oedometers) to simulate greater burial depth conditions. For low effective pressures (representing small burial depths extending to around 10 meters subsurface), the general applicability of established relationships was confirmed: the prediction of porosity, bulk density, and mean grain size from acoustic velocity and impedance appears generally possible for the investigated sedimentary environments. As effective pressure increases through, the observed relationships deviate more and more from the established ones for the near-surface area. For the samples tested in this study, in some instances increasing pressure even resulted in decreasing velocities. There are several possible explanations for this abnormal behavior, including the presence of gas, overconsolidation, or bimodal grain size distribution. The results indicate that an appropriate depth correction must be introduced into the published prediction equations in order to obtain reliable estimates of physical sediment properties for greater subsurface depths.     

Q structure of the oceanic crust

Compressional wave attenuations and velocities have been measured as a function of confining pressure in ophiolite samples representing a cross-section of the oceanic crust and uppermost mantle. Data are presented for basalts, diabase dikes, gabbros and a suite of serpentinites and peridotites showing a range of serpentization. An ultrasonic pulse-echo spectral ratio technique was used to determine the attenuations to confining pressures of 500 MPa. From this data a Q profile for the oceanic crust and upper mantle is presented. Q is found to moderately increase with depth through the pillow basalts of the upper oceanic crust. The sheeted dike rocks of Layer 2C show an increase in Q with depth due to progressive metamorphism (from greenschist to amphibolite facies). Q drops abruptly from Layer 2C to Layer 3, though it is not clear why the gabbros have such low Q's. The crust-mantle boundary is a Q discontinuity; however, the Q contrast between Layer 3 and the upper mantle could be altered by upper mantle serpentinization, interlayered gabbros and peridotites at the boundary, or serpentinized peridotite diapirs intruding the gabbroic section. Q varies significantly with the percentage of serpentinization in the ultramafic samples, with the largest changes in Q being at the extremes of zero and full serpentinization. Q is sensitive to the overburden pressure for all of the samples.     

北黄海盆地东部坳陷火成岩地震相识别与分布预测

北黄海盆地东部坳陷火成岩分布复杂并具有多期性,直接影响了油气成藏作用,然而目前对其空间分布尚缺乏深入认识,制约着该区勘探开发的部署。利用全三维地震资料,根据地震反射特征识别出了不同类型的火成岩地震相模式,并对其不同特征进行了分析,通过地震相识别进行平面追踪,结合时间切片三维地震技术及属性分析对火成岩分布进行了预测。结果表明,北黄海盆地东部坳陷主要发育柱状、尖峰状、蘑菇状和丘状4种火成岩反射地震相;研究区火成岩分区性较明显,主要发育3处火成岩体,分别位于东部隆起剥蚀区、坳陷中部和坳陷的西北部。     

Unique problems associated with seismic analysis of partially gas-saturated unconsolidated sediments

Gas hydrate stability conditions restrict the occurrence of gas hydrate to unconsolidated and high water-content sediments at shallow depths. Because of these host sediments properties, seismic and well log data acquired for the detection of free gas and associated gas hydrate-bearing sediments often require nonconventional analysis. For example, a conventional method of identifying free gas using the compressional/shear-wave velocity (V_p/V_s) ratio at the logging frequency will not work, unless the free-gas saturations are more than about 40%. The P-wave velocity dispersion of partially gas-saturated sediments causes a problem in interpreting well log velocities and seismic data. Using the White, J.E. [1975. Computed seismic speeds and attenuation in rocks with partial gas saturation. Geophysics 40, 224–232] model for partially gas-saturated sediments, the difference between well log and seismic velocities can be reconciled. The inclusion of P-wave velocity dispersion in interpreting well log data is, therefore, essential to identify free gas and to tie surface seismic data to synthetic seismograms.     

From mechanical modeling to seismic imaging of faults: A synthetic workflow to study the impact of faults on seismic

Although typically interpreted as 2D surfaces, faults are 3D narrow zones of highly and heterogeneously strained rocks, with petrophysical properties differing from the host rock. Here we present a synthetic workflow to evaluate the potential of seismic data for imaging fault structure and properties. The workflow consists of discrete element modeling (DEM) of faulting, empirical relations to modify initial acoustic properties based on volumetric strain, and a ray-based algorithm simulating prestack depth migration (PSDM). We illustrate the application of the workflow in 2D to a 100 m displacement normal fault in a kilometer size sandstone-shale sequence at 1.5 km depth. To explore the effect of particle size on fault evolution, we ran two DEM simulations with particle assemblages of similar bulk mechanical behavior but different particle size, one with coarse (1–3 m particle radii) and the other with fine (0.5–1.5 m particle radii) particles. Both simulations produce realistic but different fault geometries and strain fields, with the finer particle size model displaying narrower fault zones and fault linkage at later stages. Seismic images of these models are highly influenced by illumination direction and wave frequency. Specular illumination highlights flat reflectors outside the fault zone, but fault related diffractions are still observable. Footwall directed illumination produces low amplitude images. Hanging wall directed illumination images the shale layers within the main fault segment and the lateral extent of fault related deformation. Resolution and the accuracy of the reflectors are proportional to wave frequency. Wave frequencies of 20 Hz or more are necessary to image the different fault structure of the coarse and fine models. At 30–40 Hz, there is a direct correlation between seismic amplitude variations and the input acoustic properties after faulting. At these high frequencies, seismic amplitude variations predict both the extent of faulting and the changes in rock properties in the fault zone.     

The distribution and characteristics of the igneous complexes in the northern East China Sea Shelf Basin and their implications for hydrocarbon potential

This paper presents results of two-dimensional seismic mapping of the northern East China Sea Shelf Basin. Various igneous features such as sills, volcanic edifices and stocks were identified by the geophysical exploration. The sills are most common, and are observed at more than 90 locations. Most mapped sills in the study area are characterized by high-amplitude continuous reflections with distinct terminations. Saucer- and cup-shaped sills are observed locally. The stocks are discordant (nearly vertical) igneous bodies and they are characterized by seismic transparency, with upturned host rocks and uplifted overburden. The volcanic edifices and/or necks consist of irregular mounds and peaks and are characterized by strong positive top reflections with chaotic internal facies. The oldest igneous activity in the northern East China Sea Shelf Basin is Early Cretaceous (123.3 ± 3.7). This igneous activity coincides with those observed in eastern China which has been related mainly to the subduction of the Pacific Plate beneath Eurasia Plate. The Miocene igneous activity is well constrained based on seismic stratigraphic relationships within the folded stratigraphy, age dating, and the occurrence of igneous sills. The timing of this intrusion is coincident with the intensive igneous activity as previously suggested for the eastern China. Igneous rocks can produce hydrocarbon traps, reservoirs and they can act as a seal, and therefore are of great importance in petroleum study.     

High resolution velocity analysis of ocean bottom seismometer data by theτ-p method

A method of high resolution seismic velocity analysis for ocean bottom seismometer (OBS) records is applied to the study of the shallow oceanic crust, especially sedimentary and basement layers. This method is based on the direct-p mapping and the-sum inversion. We use data obtained from a 1989 airgun-OBS experiment in the northern Yamato Basin, Japan Sea and derive P- and S-wave velocity functions that can be compared with the seismic reflection profiles. Using split-spread profile records, we obtain interface dips and true interval velocities from the OBS data. These results show good agreement with the reflection profile records, the acoustic velocities of core samples, and sonic log profiles. We also present a method for estimating errors in the derived velocity functions by calculating covariance of the derived layers' thicknesses. The estimated depth errors are about 150 m at shallow depths, which is close to the seismic wavelength used. The high resolution of this method relies on accurate determination of shot positions by GPS, spatially dense seismic observations, and the use of unsaturated reflected waves arriving after the direct water wave that are observed on low-gain component records.     

冲绳海槽南段基底组成探讨

本文根据反射地震、折射地震、磁力等资料,结合周边地质,探讨冲绳海槽南段基底组成。在海槽周边的东海陆架盆地、台湾褶皱带和琉球岛弧褶皱带,均出露不同程度变质的晚古生代、中生代和早第三纪地层。多道反射地震表明,海槽南段沉积盖层由上第三系和第四系组成,声学基底由下第三系及更老地层构成。邻近海槽的折射地震揭示,除第四系-中新统速度层之外,还存在纵波速度分别为4．7～5．3km／s和6．3km／s的下第三系和中生界速度层。磁异常分析和正反演拟合计算结果表明,海槽磁性基底主要由变质岩系构成,次为燕山期中酸性岩浆岩和喜山期中基性岩浆岩,磁性基底大部分相当声学基底。综合分析表明,海槽南段基底主要由不同程度变质的下第三系、中生界和上古生界构成;在海槽某些构造部位,已有喜山期基性岩浆岩形成。     

Prolonged Magmatic and Tectonic Development of the Caribbean Igneous Province Revealed by a Diving Submersible Survey

Using the diving submersible survey NAUTICA we investigated the central part of the Caribbean large igneous province (CLIP) to observe and sample internal portions of this proposed oceanic plateau. Most of the samples are gabbroic and doleritic rocks; basalts are scarce. Radiometric dating by ⁴⁰Ar/³⁹Ar incremental heating experiments indicate that the intrusive rocks are Campanian in age (81–75 Ma). In some places these intrusive rocks underlie older Santonian (85–83 Ma) extrusive basaltic rocks, suggesting that the Campanian rocks represent a sill injection and an underplating episode. Results of the diving program supplemented by information from ODP and DSDP drilling sites document a 20 m.y. period (94–75 Ma) of igneous activity in the submerged portion of the Caribbean large igneous province (CLIP). In the northern part of the Beata Ridge late Campanian and/or post Campanian uplift is documented by prominent Maastrichtian (71–65 Ma) erosion and the establishment of a Paleocene-middle Eocene (65–49 Ma) carbonate platform. During and after the uplift an extensional period is indicated by seismic images and the subsidence (3 km depth) of the carbonate platform. Paleocene ages (55–56 Ma) determined on some volcanic samples are attributed to localised decompression mantle melting that accompanied the extension. We document a prolonged period of magmatic and tectonic events that do not fit with the current models of short-lived plateau formation during mantle plume initiation but shares many similarities with the constructional histories of other oceanic large igneous provinces.     

Cenozoic magmatism in the northern continental margin of the South China Sea: evidence from seismic profiles

Igneous rocks in the northern margin of the South China Sea (SCS) have been identified via high resolution multi-channel seismic data in addition to other geophysical and drilling well data. This study identified intrusive and extrusive structures including seamounts and buried volcanoes, and their seismic characteristics. Intrusive features consist of piercement and implicit-piercement type structures, indicating different energy input associated with diapir formation. Extrusive structures are divided into flat-topped and conical-topped seamounts. Three main criteria (the overlying strata, the contact relationship and sills) were used to distinguish between intrusive rocks and buried volcanos. Three criteria are also used to estimate the timing of igneous rock formation: the contact relationship, the overlying sedimentary thickness and seismic reflection characteristics. These criteria are applied to recognize and distinguish between three periods of Cenozoic magmatism in the northern margin of the SCS: before seafloor spreading (Paleocene and Eocene), during seafloor spreading (Early Oligocene–Mid Miocene) and after cessation of seafloor spreading (Mid Miocene–Recent). Among them, greater attention is given to the extensive magmatism since 5.5 Ma, which is present throughout nearly all of the study area, making it a significant event in the SCS. Almost all of the Cenozoic igneous rocks were located below the 1500 m bathymetric contour. In contrast with the wide distribution of igneous rocks in the volcanic rifted margin, igneous rocks in the syn-rift stage of the northern margin of the SCS are extremely sporadic, and they could only be found in the southern Pearl River Mouth basin and NW sub-sea basin. The ocean–continent transition of the northern SCS exhibits high-angle listric faults, concentrated on the seaward side of the magmatic zone, and a sharply decreased crust, with little influence from a mantle plume. These observations provide further evidence to suggest that the northern margin of the SCS is a magma-poor rifted margin.     

Seismic Characterisation of Gas-rich Near Surface Sediments in the Arkona Basin, Baltic Sea

Gas in sediments has become an important subject of research for various reasons. It affects large areas of the sea floor where it is mainly produced. Gas and gas migration have a strong impact on the environmental situation as well as on sea floor stability. Furthermore, large research programs on gas hydrates have been initiated during the last 10 years in order to investigate their potential for future energy production and their climatic impact. These activities require the improvement of geophysical methods for reservoir investigations especially with respect to their physical properties and internal structures. Basic relationships between the physical properties and seismic parameters can be investigated in shallow marine areas as they are more easily accessible than hydrocarbon reservoirs. High-resolution seismic profiles from the Arkona Basin (SW Baltic Sea) show distinct ‘acoustic turbidity’ zones which indicate the presence of free gas in the near surface sediments. Total gas concentrations were determined from cores taken in the study area with mean concentrations of 46.5 ml/l wet sediment in non-acoustic turbidity zones and up to 106.1 ml/l in the basin centre with acoustic turbidity. The expression of gas bubbles on reflection seismic profiles has been investigated in two distinct frequency ranges using a boomer (600–2600 Hz) and an echosounder (38 kHz). A comparison of data from both seismic sources showed strong differences in displaying reflectors. Different compressional wave velocities were observed in acoustic turbidity zones between boomer and echosounder profiles. Furthermore, acoustic turbidity zones were differently characterised with respect to scattering and attenuation of seismic waves. This leads to the conclusion that seismic parameters become strongly frequency dependent due to the dynamic properties of gas bubbles.     

Appraisal of gas-hydrate/free-gas from VP/VS ratio in the Makran accretionary prism

Existence of gas-hydrate in the marine sediments elevates both the P- and S-wave seismic velocities, whereas even a small amount of underlying free-gas decreases the P-wave velocity considerably and the S-wave velocity remains almost unaffected. Study of both P- and S-wave seismic velocities or their ratio (V_P/V_S) for the hydrate-bearing sediment provides more information than that obtained by the P- or S-wave velocity alone for the quantitative assessment of gas-hydrate. We estimate the P- and S-wave seismic velocities across a BSR (interface between gas-hydrate and free-gas bearing sediments) using the travel time inversion followed by a constrained AVA modeling of multi channel seismic (MCS) data at two locations in the Makran accretionary prism. Using this V_P/V_S ratio, we then quantify the amount of gas-hydrate and free-gas based on two rock-physics models. The result shows an estimate of 12–14.5% gas-hydrate and 4.5–5.5% free-gas of the pore volume based on first model, and 13–20% gas-hydrate and 3–3.5% free-gas of the pore volume based on the second model, respectively.     

东海陆架盆地及其邻域岩浆岩时空分布特征

基于二维地震资料和重磁资料的处理、解释,研究了东海陆架盆地及其邻域岩浆岩时空分布特征,并收集整理浙闽地区相关地质资料,将海陆连为一体进行综合解释。浙闽陆区岩浆岩非常发育,并且大致沿28°N为界,北侧岩浆岩大体呈NE向带状分布,南侧则呈NNE向带状分布,且埋深较浅,甚至在有些地区出露地表,其分布特征受区内大断裂的控制。综合地震、重磁资料研究认为,海域岩浆岩多分布在隆起部位或凹陷的边缘位置,发育情况与断裂的展布密切相关。收集整理了海域钻遇岩浆岩的15口钻井资料,沿自SW至NE一线,对岩浆岩的同位素年龄进行了统计分析,并结合其侵入关系对岩浆活动进行了期次和展布特征的划分,初步划分为4期:燕山期、四川期、华北期和喜马拉雅期;4带:浙闽隆起区、中部隆起区、钓鱼岛隆褶带和琉球岛弧带。     

南海北部深水区新生代岩浆岩分布规律及其与海底地质灾害的相关性

前人对南海北部新生代陆缘盆地结构、构造样式和期次等做了大量的研究工作,但对于具体体现内动力过程的岩浆活动关注较少,尤其是北部深水区的研究鲜有涉及。本文基于二维地震剖面的解译分析,将南海北部新生代岩浆岩划分为3个集中分布区:琼东南-西沙区、神狐区和东沙区。各区域岩浆岩的产出状态有所差异,其中琼东南-西沙区岩浆岩的平面展布规律性极强,主要受右行右阶的走滑断裂控制;神狐区岩浆岩位于NW向与NE向断裂的交汇处,侵位空间受先存断裂制约;东沙区岩浆岩的产生与东沙运动期间区域伸展过程有直接关系。作为内动力的表征,岩浆作用对于海底地形地貌具有重要的改造作用,与活动断裂、海底滑坡、浅层气等海底地质灾害也具有成因上的关联。     

白云凹陷地球物理场及深部结构特征

珠江口盆地白云凹陷是南海最具代表性的第三系深水陆坡沉积区。以穿过白云凹陷中部的一条深反射地震剖面(14s)为研究基础，采用综合地球物理研究方法分析了该区地球物理场特征，根据重力异常平面等值线勾画了白云凹陷的形态，并提取该测线相对应的重磁剖面数据，利用重磁资料和地震剖面进行了综合反演。以深剖面地震资料建立了地质模型，利用所得的重力数据进行了研究深部结构的正演拟合，实测与计算值拟合较好，支持中生代俯冲洋壳存在的观点；同时结合地震资料对深部结构进行了分析，该区莫霍面由陆向海抬升，呈阶梯状变化，地壳厚度逐渐减薄，具有大陆边缘陆壳向洋壳过渡的特征。根据地质模型还进行了变密度综合反演拟合来分析基底岩性特征，该区基底主要为中酸性岩浆岩，部分为变质岩和基性火山岩，岩石密度由陆向洋逐渐减小，磁性体分布不均。     

Synthetic seismic illumination of small-scale growth faults,paralic deposits and low-angle clinoforms: A case study of the Triassic successions on Edgeøya,NW Barents Shelf

Fault activity and sandstone-body geometries and spread, all bear significant weight to understanding the potential hydrocarbon systems on the NW Barents Shelf. Synthetic seismic modelling of onshore sedimentary successions provides insight into the seismic resolution and expression of various sedimentological features on Edgeøya: (i) sandy growth basins in pro-delta mudstones (ii) paralic sedimentary deposits and (iii) low-angle tidally influenced progradational successions.Synthetic modelling suggests that the lithological contrast associated with sand-infilled growth basins in pro-deltaic shales, even with offsets as small as 50–75 m, will create distinct geometries in seismic data. Modelling suggests that while optimal 90° illumination, low frequency bands and more typical sedimentary velocities around 2000 m/s will generate very clear discernible growth faults, the angular lithological contrast should generate discernible features even with the high velocities and typical 20–30 Hz frequency band of seismic sections on the NW Barents Shelf. Comparing to actual seismic data in which multiples, noise and in places overburden are influential, it is possible to identify growth-fault geometries more confidently, and to link them to larger planar fault activity.Modelling other features identified in the paralic sedimentary system it is apparent that many of the massive channel or lenticular shaped sandstone bodies should be identifiable in the actual seismic, although their expression is less distinct and more easily misinterpreted than that of growth faults. It is apparent that features such as igneous intrusions, unless imaged in a “perfect” survey, can be difficult to properly identify, particularly near-vertical connections. The velocity contrast creates strong impedance along horizontal sections, but heavily fractured igneous intrusions with lower velocities could easily be assumed to be isolated sandstone bodies. While the modelling appreciates the overall wedging nature of the successions, the simplified lithological observations onshore cannot predict the probable erosion/condensation contrast associated with low-angle clinoforms which are visible in the seismic data, and hence not reproduce them easily in the models.     

An integrated geophysical study of Vestbakken Volcanic Province, western Barents Sea continental margin, and adjacent oceanic crust

This paper describes results from a geophysical study in the Vestbakken Volcanic Province, located on the central parts of the western Barents Sea continental margin, and adjacent oceanic crust in the Norwegian-Greenland Sea. The results are derived mainly from interpretation and modeling of multichannel seismic, ocean bottom seismometer and land station data along a regional seismic profile. The resulting model shows oceanic crust in the western parts of the profile. This crust is buried by a thick Cenozoic sedimentary package. Low velocities in the bottom of this package indicate overpressure. The igneous oceanic crust shows an average thickness of 7.2 km with the thinnest crust (5–6 km) in the southwest and the thickest crust (8–9 km) close to the continent-ocean boundary (COB). The thick oceanic crust is probably related to high mantle temperatures formed by brittle weakening and shear heating along a shear system prior to continental breakup. The COB is interpreted in the central parts of the profile where the velocity structure and Bouguer anomalies change significantly. East of the COB Moho depths increase while the vertical velocity gradient decreases. Below the assumed center for Early Eocene volcanic activity the model shows increased velocities in the crust. These increased crustal velocities are interpreted to represent Early Eocene mafic feeder dykes. East of the zone of volcanoes velocities in the crust decrease and sedimentary velocities are observed at depths of more than 10 km. The amount of crustal intrusions is much lower in this area than farther west. East of the Kn?legga Fault crystalline basement velocities are brought close to the seabed. This fault marks the eastern limit of thick Cenozoic and Mesozoic packages on central parts of the western Barents Sea continental margin.