Rock physics modelling based on depositional and burial history of Barents Sea sandstones

Understanding how physical properties and seismic signatures of present day rocks are related to ancient geological processes is important for enhanced reservoir characterization. In this paper, we have studied this relationship for the Kobbe Formation sandstone in the Barents Sea. These rocks show anomalous low shear velocities and high V_P/V_S ratios, which does not agree well with conventional rock physics models for moderately to well consolidated sandstones. These sandstones have been buried relatively deeply and subsequently uplifted 1–2 km. We compared well log data of the Kobbe sandstone with velocity–depth trends modelled by integrating basin modelling principles and rock physics. We found that more accurate velocity predictions were obtained when first honouring mechanical and chemical compaction during burial, followed by generation of micro-cracks during uplift. We suspect that these micro-cracks are formed as overburden is eroded, leading to changes in the subsurface stress-field. Moreover, the Kobbe Formation is typically heterogeneous and characterized by structural clays and mica that can reduce the rigidity of grain contacts. By accounting for depositional and burial history, our velocity predictions become more consistent with geophysical observables. Our approach yields more robust velocity predictions, which are important in prospect risking and net erosion estimates.     

北京地区风沙活动及其整治的初步研究

根据两年来北京地区的野外风沙考察,半定位风沙(尘)观测,室内风洞实验以及有关区域资料的分析研究得出:风沙问题是北京区域环境存在的严重问题之一,区域风沙活动由来已久,在频繁程度上不亚于半干旱区的河北坝上和乌盟后山地区。区内各地风沙活动均以春季为主,冬春两季风沙活动约占全年的80%以上。风沙天气的年际变化较大,五十年代最为频繁,六十年代最少,七十年代初期又有回升,中期以后至八十年代中期又逐渐趋于减少。风沙活动的形成及其盛衰变化,除与自然因素有关外,还与地表植被和土壤结构等因人类活动所造成的变化息息相关。风沙活动不仅使大气遭受污染,而且使土地资源遭受破坏。全区风沙危害严重的区域(风沙区)面积达342万亩,约占平原区总面积的35.6%,其中沙质地表因长期遭受风力作用所形成的风沙化土地面积约为136万余亩,占平原区面积的14.2%,以永定河冲积平原分布面积最大。研究证明,北京地区风沙主要来源于本地,系就地起沙。因此,风沙防治要立足于本地,特别是风沙土地的治理,对根治本区风沙危害、改善首都环境具有根本性意义。     

Stochastic Modeling of Subseismic Faults Conditioned on Displacement and Orientation Maps

Mathematical Geosciences - Subseismic faults are small faults or fractures that may be difficult to determine but can have large consequences for fluid flow and pressure communication in the...     

Caprock interaction with CO2: A laboratory study of reactivity of shale with supercritical CO2 and brine

Crushed rock from two caprock samples, a carbonate-rich shale and a clay-rich shale, were reacted with a mixture of brine and supercritical CO₂ (CO₂–brine) in a laboratory batch reactor, at different temperature and pressure conditions. The samples were cored from a proposed underground CO₂ storage site near the town of Longyearbyen in Svalbard. The reacting fluid was a mixture of 1 M NaCl solution and CO₂ (110 bar) and the water/rock ratio was 20:1. Carbon dioxide was injected into the reactors after the solution had been bubbled with N₂, in order to mimic O₂-depleted natural storage conditions. A control reaction was also run on the clay-rich shale sample, where the crushed rock was reacted with brine (CO₂-free brine) at the same experimental conditions. A total of 8 batch reaction experiments were run at temperatures ranging from 80 to 250 °C and total pressures of 110 bar (∼40 bar for the control experiment). The experiments lasted 1–5 weeks.Fluid analysis showed that the aqueous concentration of major elements (i.e. Ca, Mg, Fe, K, Al) and SiO₂ increased in all experiments. Release rates of Fe and SiO₂ were more pronounced in solutions reacted with CO₂–brine as compared to those reacted with CO₂-free brine. For samples reacted with the CO₂–brine, lower temperature reactions (80 °C) released much more Fe and SiO₂ than higher temperature reactions (150–250 °C). Analysis by SEM and XRD of reacted solids also revealed changes in mineralogical compositions. The carbonate-rich shale was more reactive at 250 °C, as revealed by the dissolution of plagioclase and clay minerals (illite and chlorite), dissolution and re-precipitation of carbonates, and the formation of smectite. Carbon dioxide was also permanently sequestered as calcite in the same sample. The clay-rich shale reacted with CO₂–brine did not show major mineralogical alteration. However, a significant amount of analcime was formed in the clay-rich shale reacted with CO₂-free brine; while no trace of analcime was observed in either of the samples reacted with CO₂–brine.     

On the potential for CO2 mineral storage in continental flood basalts – PHREEQC batch- and 1D diffusion–reaction simulations

Continental flood basalts (CFB) are considered as potential CO₂ storage sites because of their high reactivity and abundant divalent metal ions that can potentially trap carbon for geological timescales. Moreover, laterally extensive CFB are found in many place in the world within reasonable distances from major CO₂ point emission sources. Based on the mineral and glass composition of the Columbia River Basalt (CRB) we estimated the potential of CFB to store CO₂ in secondary carbonates. We simulated the system using kinetic dependent dissolution of primary basalt-minerals (pyroxene, feldspar and glass) and the local equilibrium assumption for secondary phases (weathering products). The simulations were divided into closed-system batch simulations at a constant CO₂ pressure of 100?bar with sensitivity studies of temperature and reactive surface area, an evaluation of the reactivity of H₂O in scCO₂, and finally 1D reactive diffusion simulations giving reactivity at CO₂ pressures varying from 0 to 100?bar. Although the uncertainty in reactive surface area and corresponding reaction rates are large, we have estimated the potential for CO₂ mineral storage and identified factors that control the maximum extent of carbonation. The simulations showed that formation of carbonates from basalt at 40?C may be limited to the formation of siderite and possibly FeMg carbonates. Calcium was largely consumed by zeolite and oxide instead of forming carbonates. At higher temperatures (60 ?C 100?C), magnesite is suggested to form together with siderite and ankerite. The maximum potential of CO₂ stored as solid carbonates, if CO₂ is supplied to the reactions unlimited, is shown to depend on the availability of pore space as the hydration and carbonation reactions increase the solid volume and clog the pore space. For systems such as in the scCO₂ phase with limited amount of water, the total carbonation potential is limited by the amount of water present for hydration of basalt.     

Runoff modelling at two field slopes: use of in situ measurements of air permeability to characterize spatial variability of saturated hydraulic conductivity

The time required at a field site to obtain a few measurements of saturated hydraulic conductivity (K_s) will allow for many measurements of soil air permeability (k_a). This study investigates if k_a measured in situ (k_{a, in situ}) can be a substitute for measurement of K_s in relation to infiltration and surface runoff modelling. Measurements of k_{a, in situ} were carried out in two small agricultural catchments. A spatial correlation of the log‐transformed values existed having a range of approximately 100 m. A predictive relationship between K_s and k_a measured on 100‐cm³ soil samples in the laboratory was derived for one of the field slopes and showed good agreement with an earlier suggested predictive K_s–k_a relationship. In situ measurements of K_s and k_a suggested that the predictive relationships also could be used at larger scale. The K_s–k_a relationships together with the k_{a, in situ} data were applied in a distributed surface runoff (DSR) model, simulating a high‐intensity rainfall event. The DSR simulation results were highly dependent on whether the geometric average of k_{a, in situ} or kriged values of k_{a, in situ} was used as model input. When increasing the resolution of K_s in the DSR model, a limit of 30–40 m was found for both field slopes. Below this limit, the simulated runoff and hydrograph peaks were independent of resolution scale. If only a few randomly chosen values of K_s were used to represent the spatial variation within the field slope, very large deviations in repeated DSR simulation results were obtained, both with respect to peak height and hydrograph shape. In contrast, when using many predicted K_s values based on a K_s–k_a relationship and measured k_{a, in situ} data, the DSR model generally captured the correct hydrograph shape although simulations were sensitive to the chosen K_s–k_a relationship. As massive measurement efforts normally will be required to obtain a satisfactory representation of the spatial variability in K_s, the use of k_{a, in situ} to assess spatial variability in K_s appears a promising alternative. Copyright © 2004 John Wiley & Sons, Ltd.     

Neodymium isotopes in seawater from the Barents Sea and Fram Strait Arctic-Atlantic gateways

The neodymium concentration, C_Nd, and isotopic composition, ε_Nd, in seawater have been determined in the water column at five sites in the Barents Sea-Fram Strait area where most of the water exchange between the Arctic Ocean and the North Atlantic takes place. In the main Arctic Ocean inflow branch across the Barents Sea the concentration and isotopic composition (C_Nd = 15.5 pmol/kg and ε_Nd = −10.8) are similar to those reported for the northeastern Nordic Seas, which is consistent with this region being a source area for the Arctic inflow. Due to the addition of Nd from Svalbard shelf sediments, the C_Nd in the surface waters above 150 m, in the Fram Strait inflow branch is higher by a factor of 2 and the ε_Nd is shifted to lower values (−11.8).In the stratified Nansen Basin, where cold low salinity water overlies warmer Atlantic water the C_Nd and ε_Nd do not vary with the vertical temperature-salinity structure but are essentially constant and similar to those of the Atlantic inflow throughout the entire water column, down to 3700 m depth, which indicates that the Nd is to a large extent of Atlantic origin.Compared to the Atlantic inflow water, the Nd in the major Arctic Ocean outflow, the Fram Strait, show higher C_Nd in the surface waters above 150 m, and a higher ε_Nd (−9.8) throughout the entire water column down to 1300 m depth. Sources for the more radiogenic Nd isotopic composition in deep water of the Fram Strait outflow most likely involve boundary exchange with sediments on the shelf and slope as the water passes along the Canadian archipelago. River water is a possible source in the surface water but it also seems likely that Pacific water Nd, modified by interactions on the shelf, is an important component in the Fram Strait surface outflow. Changes in the relative proportions of inflow of river water and flow of Pacific water through the Arctic Ocean could thus influence the isotopic composition of Nd in the North Atlantic.     

Palaeogeographic,climatic and tectonic change in southeastern Australia: the Late Neogene evolution of the Murray Basin

The Murray Basin is a low-lying but extensive intracratonic depocentre in southeastern Australia, preserving an extraordinary record of Late Neogene sedimentation. New stratigraphic and sedimentologic data allow the long-term evolution of the basin to be re-evaluated and suggest a significant role for: (1) tectonism in controlling basin evolution, and (2) progressive and step-wise climatic change beginning in the early Pleistocene. Tectonic change is associated with regional uplift, occurring at approximately the same rate from the early Pliocene until the present day, and possibly associated with changing mantle circulation patterns or plate boundary processes. This uplift led to the defeat and re-routing of the Murray River, Australia’s major continental drainage system. Key to our interpretation is recognition of timing relationships between four prominent palaeogeographic features – the Loxton-Parilla Sands strandplain, the Gambier coastal plain, palaeo megalake Bungunnia and the Kanawinka Escarpment. Geomorphic and stratigraphic evidence suggest that during the Early Pliocene the ancestral Murray River was located in western Victoria, flowing south along the Douglas Depression. Relatively small amounts of regional uplift (<200 m) defeated this drainage system, dramatically changing the palaeogeography of southeastern Australia and forming Plio-Pleistocene megalake Bungunnia. At its maximum extent Lake Bungunnia covered more than 50,000 km², making it one of the largest known palaeo- or modern-lakes in an intracontinental setting. Magnetostratigraphic constraints suggest lake formation c. 2.4 Ma. The formation of Lake Bungunnia influenced the Pliocene coastal dynamics, depriving the coastline of a sediment source and changing the coastal system from a prograding strandline system to an erosional one. Erosion during this period formed the Kanawinka Escarpment, a palaeo sea-cliff and one of the most prominent and laterally extensive geomorphic features in southeastern Australia. Marine sediments c. 800 ka to c. 1.16 Ma represent the time of re-establishment of depositional coastal dynamics and of a permanent outlet for the Murray River. This age range is consistent with our best estimate of the age of the youngest Lake Bungunnia sediments and points towards an early Pleistocene age for the demise of the lake system. The youngest Lake Bungunnia sediment, present on a number of distinct terraces, suggests that progressive, step-wise climatic change played a role in the demise of the lake. However, in order for the ancestral Murray River system to have been able to breach the pre-existing tectonic dam, it is likely that tectonic change and/or temporarily enhanced discharge was also significant. This scenario indicates that the modern Murray River has only been in existence for at most 700 ka.