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

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.     

Composition and properties of glacigenic sediments in the southwestern Barents Sea

Abstract

The composition and properties of glacigenic sediments in the southwestern Barents Sea are described based on data from 33 shallow boreholes (< 143 m below seabed) and 11 seabed cores (<4.2m below seabed). The cores are tied into a regional seismostratigraphic framework, illustrating the relationships between different boreholes.

A massive, muddy diamicton (silty, sandy clay with scattered gravel) is found in nearly all cores. Average clay content (<2 pm) of this lithology is about 38%, but varies between about 25% and 50%. Short intervals of finely laminated, waterlain sediments or gravelly sand are cored in a few occasions. A high content of sand and gravel in the cores from near the Norwegian coast shows an influence of sediment input from the mainland, while material eroded from sedimentary rocks dominates farther offshore.

The data presented on physical properties include undisturbed and remolded undrained shear strength, natural water content, bulk density, compressional sound velocity (P waves), Atterberg consistency limits, effective preconsolidation pressure, and consolidation coefficient.

Prediction of overconsolidation from seismic mapping of erosional surfaces is confirmed by the borehole cores. High compaction is found both in Weichselian and older deposits, with a general increase in compaction toward the east as well as toward shallower water. Cores that are “underconsolidated” at their present burial depth are also reported.

The average compressional sound velocity is about 1780 m/s for the borehole cores, 1550 m/s for the seabed cores, and increases with increasing shear strength and consolidation. Both horizontal and vertical sound velocities are measured in several cores, and although the data have a considerable scatter, a slightly aniso‐tropic sound velocity is indicated.     

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.     

An experimental study of velocities under pressure for ancient oceanic rocks

Abstract

Seismic velocities of Indus suture rock types from Dras‐Sanko‐Kargil, Kashmir Himalaya, as a function of pressure up to 10 kbar were studied. The high‐pressure measurements on the rocks reflect the depthwise increase in velocity, and in general they help in better understanding and better interpretation of the regional rocks in terms of their geological observations. An attempt has been made at correlation with ocean‐dredged samples, and it was found that the ultramafics, gabbros, and dykes are compatible with oceanic rocks and other ophiolite sequences, whereas metavolcanics are incompatible, suggesting the dismembered nature of Indus ophioli‐tes.     

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.     

台湾海峡盆地九龙江凹陷火成岩发育特征及其油气地质意义

基于二维地震剖面的精细地质解释,并结合磁力和钻井资料分析,开展了九龙江凹陷火成岩发育特征、岩浆活动和区域构造与盆地演化关系及其油气地质意义的研究。分析认为,九龙江凹陷经历了始新世中晚期、中新世中晚期二期火山活动,火成岩在地震剖面上主要表现为板状、锥状及块状三种地震相类型,钻井主要揭示了喷发、侵入及火山碎屑三大类型火成岩,它们的发育与分布是区域性构造事件和盆地演化相作用的结果。始新世中晚期火山活动产生的高温有利于烃源岩在短期内经历高温的催化改造而生成油气,形成的火成岩体岩墙面可以增加和改善凹陷中央区油气向西部斜坡的运移,该火成岩附近区的一些圈闭在成岩期受热液蚀变的影响,储层孔、渗性得到改善的同时也会致使储层发生不同程度的变质,在今后的勘探开发中是一个值得关注的问题。     

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.     

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

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

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.     

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.     

Provenance identification and sedimentary analysis of the beach and bar systems in the Palaeogene of the Enping Sag,Pearl River Mouth Basin,South China Sea

During deposition in the Palaeogene, the Pearl River Mouth Basin (Enping Sag) was filled by sediments derived from northern and southern source regions. The northern sources are the fold belt Caledonian and Hercynian metamorphic rocks that are widely exposed along the South China block and stacked with Mesozoic igneous rocks, and the southern sources are Mesozoic igneous rocks from the southern Panyu lower uplift. We analysed the characteristics of the provenance and sediment-dispersal patterns of the beach and bar systems in the Palaeogene third-order sequence framework on the basis of detritus, thin section, zircon and the geophysical response test data. The source and deposition areas were distinguished based on their palaeogeomorphology. The sedimentary material composition of the deposition area was determined based on the characteristics of the detrital assemblage and cathode luminescence, and the characteristics of the age distribution were determined in detail using detrital zircon U–Pb dating. This study suggests that the provenance system in the Palaeogene mainly consisted of Palaeozoic metamorphic rocks and Mesozoic igneous rocks, with a supplement of Proterozoic metamorphic rocks. The transport capacity of three large fault troughs that delivered the northern provenances gradually increased from the Wenchang Formation (E₂w) to the Enping Formation (E₃e), whereas the southern provenances were transported by four or five incised valleys in E₂w and then submerged in E₃e. In the 3D seismic area, the features of the beach and bar sedimentary systems were described and interpreted based an analysis of the well-based sedimentary facies, seismic reflection characteristics and multiple attribute clustering. The results show that the beach and bar systems by the short-axis braided deltas mainly developed in the upper member of E₂w and that the beach and bar systems on the subaqueous high in E₃e developed during the Palaeogene. In addition, the bedrock-beach bar sedimentary system in the lower member (E₃e^L) was transferred to the braided delta–beach bar sedimentary system in the upper member (E₃e^U).     

Comparison of high-resolution seismic profiles and the geoacoustic properties of Eckernförde Bay sediments

High-resolution seismic profiles of Eckernförde Bay and the adjacent Baltic Sea were collected, and the geoacoustic properties of sediments there were measured. Bulk densities averaged ~ 1.35 g cm^–3 and ranged from ~ 1.2 to ~ 1.7 g cm^–3. Compressional wave velocities in gas-free sediments averaged ~ 1460 m s^–1 and ranged from ~ 1425 to ~ 1555 m s^–1. In nongassy sediments, bulk density variations typically controlled changes of acoustic impedance. Impedance changes were usually too small and closely spaced to be resolved seismically, although, at certain sites, significant impedance changes are far apart enough that they correlate one-to-one with seismic reflectors. Where free gas is present, velocity decreases and wave energy is scattered, causing a prominent seismic reflector.     

利用重磁震资料解释推断东海海底火成岩

利用1:100万上海幅海洋区域地质调查项目最新获得的重力、磁力和多道地震数据,根据研究区内重、磁场特征,以及已知的火成岩发育状况,在数据处理的基础上,通过剖面和平面相结合,对东海海底的火成岩进行了推断解释,从而得到火成岩在空间上的展布位置。分析认为东海海底火成岩具有沿北北东向成带发育的特点,沿深大断裂侵入于中生代地层;北北西向火成岩带发育时间较晚,规模小,与北西向断裂形成较晚相对应,火成岩分布呈现东西分带的特点,陆架盆地中低部隆起为喜马拉雅早-中期岩浆岩。     

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.     

Physical and elastic properties of marine sediments off Bombay,India

Abstract

Keeping in view the paucity of information as to the nature of the marine sediments from the continental shelf adjoining the Indian subcontinent, a number of shallow seismic surveys were carried out, nearshore and offshore Bombay between 18°45'N and 21°00'N. Representative core samples preserving their natural state were also retrieved from the region in the water depths ranging from 5 to 70 m for the determination of physical properties in the laboratory. Data on the physical, acoustic, and elastic properties of the sediment cores are reported for the first time. Useful individual least‐squares relations are presented for acoustic impedance, reflection coefficient, and bulk modulus against density; for the dependence of rigidity and bulk moduli on the constrained modulus; for the association between impedance and field sediment velocity against P‐velocity; and for rigidity against Poisson's ratio. Results indicate that the bulk modulus and Young's modulus are higher for silty clay and clayey silt samples than for the clay samples. Similarly, the acoustic impedance, reflection coefficient, and constrained modulus of silty clay and clayey silt are also higher than clay. The results are found to be comparable to the North Atlantic and Bay of Bengal sediments.     

火成岩发育区地震资料构造成图校正方法及应用

渤海海域BZ油田区新生界广泛发育火成岩。由于不同期次的高速火成岩纵横向分布极不均衡,造成其下伏构造受影响程度存在很大差异,常规的时深转换方法无法真实反映地下的构造形态。笔者提出了一种新的火成岩下构造恢复的方法：首先利用数值模拟分析不同厚度火成岩对下伏构造的影响程度,建立火成岩厚度与时间校正量的数学关系式,然后运用层速度约束反演方法定量识别火成岩的累计厚度,进而得到其引起的时间校正量,最终采用去除火成岩影响的时深关系进行时深转换得到更精确的深度域构造图。实践表明,该方法有效提高了火成岩发育区的构造成图精度,为渤海海域内首次开展新生界火成岩下油田评价提供了有力的技术支撑,具有较好的推广应用价值。     

双参数约束地震叠前同时反演在北黄海盆地某区块火成岩预测中的应用

北黄海盆地是中国东部勘探及研究程度较低的含油气盆地之一,多期构造活动伴生形成的大量火成岩进一步增加了北黄海盆地油气勘探工作深入的难度。为了掌握北黄海盆地某区块火成岩的空间分布特征,利用盆地近年最新采集处理的某区块三维地震资料,通过区内钻测井资料综合分析建立特有火成岩岩石物理模型,运用双参数约束三维地震叠前同时反演技术对研究区的火成岩分布进行预测。结果显示利用纵波阻抗和纵横波速度比双参数联合约束的地震叠前同时反演技术,可以对区内的火成岩进行较为有效的刻画,分析认为北黄海盆地某区块内火成岩由多期形成,中深部的火成岩主要以岩株、岩墙或岩脉状为主,具有规模小、分布广的特点。     

Acoustic estimations of porosity and thickness for thin soft surface layer of sea floor

Abstract

In some areas, the sea floor contains a surface layer with high porosity and a thickness of 10–30 cm. The characters of an explosive pulse reflected at such a sea floor are discussed. The thickness and porosity estimations for the surface layer by acoustic method are studied, and the estimated and measured results are listed.     

Paleogene igneous intrusion and its effect on thermal maturity of organic-rich mudstones in the Beibuwan Basin,South China Sea

The seismic, drilling and logging data reveal that a large-scale igneous intrusion with a width of 14 Km and a maximum thickness of 170 m intruded within the Paleogene Liushagang Formation in the Fushan Depression, Beibuwan Basin, South China Sea (SCS). In this study, we report the geochemistry and Sr-Nd-Pb-Hf isotopic compositions of the igneous rocks, and evaluate the thermal effect induced by this large-scale igneous intrusion on the host rocks. The analyzed igneous samples exhibit strong enrichment in light rare earth elements (LREE) and large-ion lithophile elements (LILE), having characteristics similar to intra-plate oceanic island basalts (OIB). The Sr-Nd-Pb-Hf isotopic data display narrow ranges (e.g. ⁸⁷Sr/⁸⁶Sr_i = 0.7042–0.7044, ¹⁴³Nd/¹⁴⁴Nd_i = 0.5128–0.5129, ²⁰⁶Pb/²⁰⁴Pb_i = 18.90–18.94, εHf(t) = +7.56∼+9.60). Geochemical and isotopic compositions suggest a mixed mantle source between depleted mid-ocean-ridge-basalt (MORB) mantle (DMM)-like mantle and enriched mantle type II (EMII, possibly the Hainan mantle plume). Vitrinite reflectance values, major mineralogical changes together with a maturity-related biomarker [Ts/(Ts + Tm)] all reveal significant thermal effect induced by the igneous intrusion. Vitrinite reflectance values of the host rock are up to 2.5% in the intrusion region, whereas lower reflectance values (0.62–0.73%) occur in the unaffected area of the same strata. Moreover, our results suggest that the host rocks above the igneous intrusion are characterised by higher maturity than below, which should be attributed to the different behavior of hydrothermal fluids. These observations suggest that the thermal effect of large-scale thick igneous intrusions is much more intense than that of thin igneous intrusions, and the behavior of hydrothermal fluids induced by magmatic intrusive events should be a critical impact factor during heat transfer process.     

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.