Method for generating parameterized ecotoxicity data of dispersed oil for use in environmental modelling

The aim of the work was to establish methodology for realistic laboratory-based test exposures of organisms to oil dispersions, specifically designed to generate parameterized toxicity data. Such data are needed to improve the value of numerical models used to predict fate and effects of oil spills and different oil spill responses. A method for continuous and predictable in-line production of oil dispersions with defined size distribution of different oil qualities was successfully established. The system enables simultaneous comparison between the effects of different concentrations of dispersion and their corresponding equilibrium water soluble fractions. Thus, net effects of the oil droplet fraction may be estimated. The method provides data for comparing the toxicity of oil dispersions generated both mechanically and with the use of chemical dispersions, incorporating the toxicity of both dissolved oil and droplets of oil.     

Assessment of terrestrial laser scanning technology for obtaining high‐resolution DEMs of soils

Terrestrial Laser Scanners (TLS) provide a non‐contact method to measure soil microtopography of relatively large surface areas. The appropriateness of the technology in relatation to the derived Digital Elevation Models (DEM) however has not been reported. The suitability of TLS for soil microtopography measurements was tested on‐field for three large soil surface areas in agricultural fields. The acquired point clouds were filtered with a custom cloud import algorithm, and converted into digital elevation models (DEM) of different resolutions. To assess DEM quality, point clouds measured from different viewpoints were statistically compared. The statistical fit between point clouds from different viewpoints depends on spatial resolution of the DEM. The best results were obtained at the higher resolutions (0.02 to 0.04 cm), where less than 5 % of the grid cells showed significant differences between one viewpoint and the next (p < 0.01). Copyright © 2012 John Wiley & Sons, Ltd.     

A model for underpressure development in a glacial valley,an example from Adventdalen,Svalbard

The underpressure observed in the glacial valley Adventdalen at Svalbard is studied numerically with a basin model and analytically with a compartment model. The pressure equation used in the basin model, which accounts for underpressure generation, is derived from mass conservation of pore fluid and solid, in addition to constitutive equations. The compartment model is derived as a similar pressure equation, which is based on a simplified representation of the basin geometry. It is used to derive analytical expressions for the underpressure (overpressure) from a series of unloading (loading) intervals. The compartment model gives a characteristic time for underpressure generation of each interval, which tells when the pressure state is transient or stationary. The transient pressure is linear in time for short‐time spans compared to the characteristic time, and then it is proportional to the weight removed from the surface. We compare different contributions to the underpressure generation and find that porosity rebound from unloading is more important than the decompression of the pore fluid during unloading and the thermal contraction of the pore fluid during cooling of the subsurface. Our modelling shows that the unloading from the last deglaciation can explain the present day underpressure. The basin model simulates the subsurface pressure resulting from erosion and unloading in addition to the fluid flow driven by the topography. Basin modelling indicates that the mountains surrounding the valley are more important for the topographic‐driven flow in the aquifer than the recharging in the neighbour valley. The compartment model turns out to be useful to estimate the orders of magnitude for system properties like seal and aquifer permeabilities and decompaction coefficients, despite its geometric simplicity. We estimate that the DeGeerdalen aquifer cannot have a permeability that is higher than 1 · 10^?18 m², as otherwise, the fluid flow in the aquifer becomes dominated by topographic‐driven flow. The upper value for the seal permeability is estimated to be 1 · 10^?20 m², as higher values preclude the generation and preservation of underpressure. The porosity rebound is estimated to be <0.1% during the last deglaciation using a decompaction coefficient α_r = 1 · 10^?9 Pa^?1.     

Modeling semi-steady state near-well flow performance for horizontal wells in anisotropic reservoirs

This paper presents a novel methodology to model semi-steady state horizontal well flow performance in an anisotropic reservoir taking into account flow in the near-well region for an arbitrary well trajectory. It is based on an analytical productivity model describing coupled axial reservoir flow and radial well inflow. In order to apply this model in an anisotropic reservoir, the permeability field relative to the radial direction perpendicular to the well trajectory and the axial direction along the well trajectory must first be determined. A classical space transformation is used in concert with rotational transforms to obtain a virtual isotropic model. The transformation preserves the volumes and pressures. It is not a novel concept, but different from previous approaches in the sense that it is only applied in the near-well domain to formulate an equally isotropic media. As a result, the use of this virtual isotropic model requires the Dietz shape factor for an ellipse, transformed from the original cylindrical near-well domain. The Dietz shape factors are determined numerically in this research. The semi-steady state well/near-well model is implemented in a numerical simulator incorporating formation anisotropy and wellbore hydraulics. The specific productivity index along the well trajectory is generated using the virtual configuration. Numerical results for different anisotropy ratios and also incorporating frictional losses in the well are presented. Furthermore, the well/near-well model is applied in coupling with streamline reservoir model for a water flooding case. This appears to be the first coupling of a well hydraulics model and a streamline simulator. It presents the application of the well/near-well model in integrated reservoir simulation in an efficient and accurate manner. The results demonstrate that the coupling approach with a streamline reservoir model and the well/near-well is of great potential for advanced well simulation efficiently.     

A method for estimating sediment budgets of washover deposits using digital terrain models

Washover fans are located on small barriers in fetch-limited micro-tidal coastal environments in Denmark. These washover fans are formed during high-energy storm events and we present a method to quantify their volumes and to estimate sediment exchanges between washover fans and their adjacent morphologies. We use high resolution digital terrain models (DTMs) based on light detection and ranging (LiDAR) data. We have delineated landforms using known methods of scale analysis and geomorphometric classification. We quantified volumes of the delineated landforms and estimated the related sediment budgets. These computed volumes were compared using different pre-depositional surfaces. Finally, we assessed the sediment exchange and associated sources of sediments of the washover fans. We applied a scale analysis to determine suitable DTM resolution and focal statistics window size as input to a geomorphometric classification analysis. Landform areas and landforms were delineated using morphometric threshold values, and volumes and sediment budgets of the delineated landforms were computed using different assumptions to define the pre-depositional surface. Resulting washover fan volumes were validated against digital elevation model (DEM) of difference (DoD) derived volumes. Sediment budgets were derived from representative volumes of the washover fans and adjacent berms. We show that quantification of washover features derived from DTMs, using geomorphometric analysis is feasible and that the presented approach provides estimates of washover deposit volumes with an accuracy between 1% and 28% compared to control volumes. © 2021 John Wiley & Sons, Ltd.     

Structural geology of the Huldra Field, northern North Sea—a major tilted fault block at the eastern edge of the Horda Platform

The Huldra fault block is a rotated major fault block on the east margin of the Viking Graben in the northern North Sea. Unlike the rest of the Horda Platform area, the Jurassic section in the Huldra fault block was rotated more than 20° during slip on the listric Huldra fault, which forms a low-angle detachment beneath the Huldra fault block. The fault block is interpreted as resulting from marginal collapse of the Horda Platform after relief along the eastern margin of the Viking Graben built up in early parts of the middle to late Jurassic rifting history. The collapse resulted in NW directed transport of the Huldra fault block, consistent with a previously postulated change in extension direction from W–E to NW–SE toward the end of the Jurassic period. Minor faults within the Hulrda fault block are consistent with E–W extension and thus may have formed early during the late Jurassic rifting phase. Nevertheless, the crest (Huldra Field) seems surprisingly intact, considering its proximity to a major fault zone. Deformation bands studied from core material are non-cataclastic and concentrated in zones. Evidence for smearing along a cored fault surface indicates that minor subseismic faults may be sealing. Production data from the field indicate good communication between most wells, suggesting that the subseismic faults and deformation band zones that are present in the reservoir have relatively small influence on the flow of gas in the reservoir.     

Effects of pore fluid pressure on the seismic response of a fractured carbonate reservoir

An effective medium model for the stress-dependent seismic properties of fractured reservoirs is developed here on the basis of a combination of a general theory of viscoelastic waves in rock-like composites with recently published formulae for deformation of communicating and interacting cavities (interconnected pores/cracks and fractures at finite concentration) under drained loading. The inclusion-based model operates with spheroidal cavities at two different length scales; namely, the microscopic scale of the pores and (grain-boundary) cracks, and the mesoscopic scale of the fractures (controlling the flow of fluid). The different cavity types can in principle have any orientation and aspect ratio, but the microscopic pores/cracks and mesoscopic fractures were here assumed to be randomly and vertically oriented, respectively. By using three different aspect ratios for the relatively round pores (representing the stiff part of the pore space) and a distribution of aspect ratios for the relatively flat cracks (representing the compliant part of the pore space), we obtained a good fit between theoretical predictions and ultrasonic laboratory measurements on an unfractured rock sample under dry conditions. By using a single aspect ratio for the mesoscopic fractures, we arrived at a higher-order microstructural model of fractured porous media which represents a generalization of the first-order model developed by Chapman et al. (2002,2003). The effect of cavity size was here modelled under the assumption that the characteristic time for wave-induced (squirt) flow at the scale of a particular cavity (pore/crack vs. fracture) is proportional with the relevant scale-size. In the modelling, we investigate the effect of a decreasing pore pressure with constant confining pressure (fixed depth), and hence, increasing effective pressure. The analysis shows that the attenuation-peak due to the mesoscopic fractures in the reservoir will move downward in frequency as the effective pressure increases. In the range of seismic frequencies, our modelling indicates that the P-wave velocities may change by more than 20% perpendicular to the fractures and close to 10% parallel to the fractures. In comparison, the vertical S-wave velocities change by only about 5% for both polarization directions (perpendicular and parallel to the fractures) when the effective pressure increases from 0 to 15 MPa. This change is mainly due to the overall change in porosity with pressure. The weak pressure dependence is a consequence of the fact that the S waves will only sense if the fractures are open or not, and since all the fractures have the same aspect ratio, they will close at the same effective pressure (which is outside the analysed interval). Approximate reflection coefficients were computed for a model consisting of the fractured reservoir embedded as a layer in an isotropic shale and analysed with respect to variations in Amplitude Versus Offset and aZimuth (AVOZ) properties at seismic frequencies for increasing effective pressure. For the P-P reflections at the top of the reservoir, it is found that there is a significant dependence on effective pressure, but that the variations with azimuth and offset are small. The lack of azimuthal dependence may be explained from the approximate reflection coefficient formula as a result of cancellation of terms related to the S-wave velocity and the Thomson’s anisotropy parameter δ. For the P-S reflection, the azimuthal dependence is larger, but the pressure dependence is weaker (due to a single aspect ratio for the fractures). Finally, using the effective stiffness tensor for the fractured reservoir model with a visco-elastic finite-difference code, synthetic seismograms and hodograms were computed. From the seismograms, attenuation changes in the P wave reflected at the bottom of the reservoir can be observed as the effective pressure increases. S waves are not much affected by the fractures with respect to attenuation, but azimuthal dependence is stronger than for P waves, and S-wave splitting in the bottom reservoir P-S reflection is clearly seen both in the seismograms and hodograms. From the hodograms, some variation in the P-S reflection with effective pressure can also be observed.     

Numerical simulation of rifting in the northern Viking Graben: the mutual effect of modelling parameters

Numerous basin modelling studies have been performed on the Viking Graben in the northern North Sea during the past decades in order to understand the driving mechanisms for basin evolution and palaeo temperature estimations. In such modelling, it is important to include lithospheric flexure. The values derived for the lithospheric strength by these studies vary considerably (i.e. up to a factor of 30). In this study, which is based on new interpretation of a regional transect, we show that both the estimated value of the effective elastic thickness and the derived β-profile are dependent on the assumed value of the depth of necking. The use of models that implicitly set the level of necking at a depth of 0 km generally leads to an underestimation of the lithospheric strength, and an overestimation of the thinning factors. In the northern Viking Graben, a necking depth at intermediate crustal levels gives results comparable to the observations. Extension by faulting is modelled to be a significant factor. In conclusion, rifting in the northern Viking Graben can be explained with various models of effective elastic thicknesses (EET) varying from 1 km for a zero necking depth to the depth of the 450 °C isotherm for an intermediate level of necking.It is also shown that the development of the basin during the post-rift phase cannot be explained by pure shear/simple shear extension. Two mechanisms are proposed here to explain the post-rift subsidence pattern, namely intra-plate stress and phase boundary migration. The two extreme models for EET mentioned above (1 km for a zero necking depth to the depth of the 450 °C isotherm for an intermediate level of necking) give very different responses to compressional stress, the latter gives basically no response for realistic intra-plate stress.     

Threshold response of benthic macrofauna integrity to metal contamination in West Greenland

Sediment metal chemistry and benthic infauna surveys have been conducted over 33 years following a BACI protocol in relation to submarine tailings deposition (STD) from a lead-zinc mine in a western Greenland fjord system. We found clear predictable changes of benthic fauna composition in response to STD both temporally and spatially. Faunal re-colonization 15 years after mine closure, was slow and the impacted areas were still dominated by opportunistic species, although the most opportunistic ones (e.g. Capitella species) had decreased in importance. Concentration-response relations between sediment lead and faunal indices of benthic community integrity (e.g. the AMBI and DKI indices) indicated a threshold of ca. 200mg/kg, above which deterioration of faunal communities occurred. Above this threshold, diversity decreased dramatically and dominance of sensitive and indifferent species was substituted by tolerant or opportunistic species. Disposal of metal contaminated tailings may have long lasting effects on the biological system.