Generalized dual porosity: A numerical method for representing spatially variable sub-grid scale processes

A generalized dual porosity method (GDPM) has been developed to incorporate sub-grid scale heterogeneity into large-scale flow and transport simulations. The method is spatially variable in the sense that the method can be applied with different levels of resolution for different spatial nodes in the simulation. The method utilizes the nodal connectivity structure and linear equation solvers of unstructured grids like those used in the finite element method, and can be applied to any problem without externally modifying the numerical grid. The algorithm scales linearly in CPU time and storage with the number of GDPM nodes. We demonstrate the utility and computational efficiency of the technique with two verification problems and an example problem of a field site.     

Resolving the contributing factors to Mississippi Delta subsidence: Past and Present

To date, quantification of individual components that contribute to shallow and deep‐seated subsidence in passive margin deltas worldwide has proven problematic. A new, regional gridded chronostratigraphic dataset for the Lower Mississippi Delta region, derived from 80,928 well reports across the northern Gulf of Mexico (GOM), has bridged the disparity between geodetic mean rates measuring total land surface subsidence across annual‐to‐decadal timescales and the deep‐seated stratigraphic subsidence rates that record isostatic response over timescales of >10⁴ years. Through a quantitative assessment of gridded chronostratigraphic surfaces, sections, and subsidence rates extending from the Middle Pleistocene (0.58 Ma) to the Late Pliocene (3.85 Ma), we identify both temporal and spatial variability in deep‐seated subsidence across the northern GOM. Targeted deep‐seated subsidence data extracted across prior GOM Holocene sea‐level sample locations have revealed more than an order of magnitude greater rates of isostatic compensation in the Mississippi depocentre versus similar GOM sea‐level control sites in Florida and Alabama, casting doubt on efforts towards a representative Holocene sea‐level curve. Spatial variability in subsidence was also assessed locally in both the strike and dip directions to assess the contributions of growth faults. Fault throw displacement magnitude was discovered to decrease with depth, accounting for less than half of the total deep‐seated subsidence record of the Middle Pleistocene. Temporal subsidence complexities were also revealed including a direct, inverse logarithmic relationship between subsidence rate and time indicating variable subsidence component controls across different timescales. Despite the spatial and temporal complexities, this dataset serves as the first regional baseline for deep‐seated subsidence rates across the northern GOM.     

Effects of hull–ice friction coefficient on predictions of pack ice forces for moored offshore vessels

Designers and offshore operators frequently predict pack ice loading on offshore vessels by conducting scale model tests. One factor that can affect pack ice loading is the hull–ice friction coefficient. This research investigates the effect of hull–ice friction coefficients for a moored offshore vessel model and includes ice floe size and ice concentration as additional variables. A method of non-dimensional analysis is modified in order to deal with the multivariate nature of the new data. The resulting non-dimensional equation provides insight on relationships between the predicted pack ice force and the variables under investigation. The relationship between pack ice force and hull–ice friction coefficient is shown to be approximately a fourth root function, while the relationship for ice floe size, non-dimensionalised by the vessel beam, is approximately linear. The relationship between predicted pack ice force and ice concentration exists in a band bounded by cubic and sixth power curves. Applying the modified equation to the previous data sets shows the current analysis slightly improves the normalization of pack ice forces.     

Compositional changes in natural gas bubble plumes: observations from the Coal Oil Point marine hydrocarbon seep field

Detailed measurements of bubble composition, dissolved gas concentrations, and plume dynamics were conducted during a 9-month period at a very intense, shallow (22-m water depth) marine hydrocarbon seep in the Santa Barbara Channel, California. Methane, carbon dioxide, and heavier hydrocarbons were lost from rising seep bubbles, while nitrogen and oxygen were gained. Within the rising seawater bubble plume, dissolved methane concentrations were more than 4 orders of magnitude greater than atmospheric equilibrium concentrations. Strong upwelling flows were observed and bubble-rise times were ~40 s, demonstrating the rapid exchange of gases within the bubble plume.     

Dynamic response of a tanker moored to an articulated loading platform

An analytical model was developed for the dynamic analysis of an articulated loading platform in an operation condition, while remaining in a head seas position. The environmental excitation considered, resulting from groups of regular waves, included first- and second-order force contributions. The nylon hawser connecting the tanker to the ALP was modeled as a nonlinear spring. The hydrodynamic load on the tower was evaluated using Morison's equation, which was modified to account for the relative motion of the tower and the fluid particles. The hydrodynamic load on the tanker was calculated using linear diffraction theory based on the 2-D Helmholtz equation. The “near field” approach of Pinkster was used to evaluate the drift force.     

How processing digital elevation models can affect simulated water budgets

For regional models, the shallow water table surface is often used as a source/sink boundary condition, as model grid scale precludes simulation of the water table aquifer. This approach is appropriate when the water table surface is relatively stationary. Since water table surface maps are not readily available, the elevation of the water table used in model cells is estimated via a two-step process. First, a regression equation is developed using existing land and water table elevations from wells in the area. This equation is then used to predict the water table surface for each model cell using land surface elevation available from digital elevation models (DEM). Two methods of processing DEM for estimating the land surface for each cell are commonly used (value nearest the cell centroid or mean value in the cell). This article demonstrates how these two methods of DEM processing can affect the simulated water budget. For the example presented, approximately 20% more total flow through the aquifer system is simulated if the centroid value rather than the mean value is used. This is due to the one-third greater average ground water gradients associated with the centroid value than the mean value. The results will vary depending on the particular model area topography and cell size. The use of the mean DEM value in each model cell will result in a more conservative water budget and is more appropriate because the model cell water table value should be representative of the entire cell area, not the centroid of the model cell.     

A Practical Approach to the Design, Monitoring, and Optimization of In Situ MTBE Aerobic Biobarriers

A paradigm for the design, monitoring, and optimization of in situ methyl tert -butyl ether (MTBE) aerobic biobarriers is presented. In this technology, an oxygen-rich biologically reactive treatment zone (the "biobarrier") is established in situ and downgradient of the source of dissolved MTBE contamination in groundwater, typically gasoline-impacted soils resulting from leaks and spills at service station sites or other fuel storage and distribution facilities. The system is designed so that groundwater containing dissolved MTBE flows to, and through, the biobarrier treatment zone, ideally under natural gradient conditions so that no pumping is necessary. As the groundwater passes through the biobarrier, the MTBE is converted by microorganisms to innocuous by-products. The system also reduces concentrations of other aerobically degradable chemicals dissolved in the groundwater, such as benzene, toluene, xylenes, and tert -butyl alcohol. This design paradigm is based on experience gained while designing, monitoring, and optimizing pilot-scale and full-scale MTBE biobarrier systems. It is largely empirically based, although the design approach does rely on simple engineering calculations. The paradigm emphasizes gas injection–based oxygen delivery schemes, although many of the steps would be common to other methods of delivering oxygen to aquifers.     

龙门山褶皱逆冲带中前山逆冲带之新期构造及断层运动学(英文)

The disastrous Wenchuan, Sichuan, earthquake (MS=8.0) on 12 May 2009 ruptured several major thrust faults of the Longmenshan fold-and-thrust belt, along the western margin of the Sichuan Basin.