Reduced order models for many-query subsurface flow applications

Inverse modeling involves repeated evaluations of forward models, which can be computationally prohibitive for large numerical models. To reduce the overall computational burden of these simulations, we study the use of reduced order models (ROMs) as numerical surrogates. These ROMs usually involve using solutions to high-fidelity models at different sample points within the parameter space to construct an approximate solution at any point within the parameter space. This paper examines an input–output relational approach based on Gaussian process regression (GPR). We show that these ROMs are more accurate than the linear lookup tables with the same number of high-fidelity simulations. We describe an adaptive sampling procedure that automatically selects optimal sample points and demonstrate the use of GPR to a smooth response surface and a response surface with abrupt changes. We also describe how GPR can be used to construct ROMs for models with heterogeneous material properties. Finally, we demonstrate how the use of a GPR-based ROM in two many-query applications—uncertainty quantification and global sensitivity analysis—significantly reduces the total computational effort.     

Inversion of 2D spectral induced polarization imaging data

Laboratory measurements of various materials suggest that more information can be obtained by measuring the in‐phase and out‐of‐phase potentials at a number of frequencies. One common model used to describe the variation of the electrical properties with frequency is the Cole‐Cole model. Apart from the DC resistivity (ρ) and chargeability (m) parameters used in conventional induced‐polarization (IP) surveys, the Cole‐Cole model has two additional parameters, i.e. the time (τ) and relaxation (c) constants. Much research has been conducted on the use of the additional Cole‐Cole parameters to distinguish between different IP sources. Here, we propose a modified inversion method to recover the Cole‐Cole parameters from a 2D spectral IP (SIP) survey. In this method, an approximate inversion method is initially used to construct a non‐homogeneous starting model for the resistivity and chargeability values. The 2D model consists of a number of rectangular cells with constant resistivity (ρ), chargeability (m), time (τ) and relaxation (c) constant values in each cell. A regularized least‐squares optimization method is then used to recover the time and relaxation constant parameters as well as to refine the chargeability values in the 2D model. We present results from tests carried out with the proposed method for a synthetic data set as well as from a laboratory tank experiment.     

10 Myr evolution of sedimentation rates in a deep marine to non‐marine foreland basin system: Tectonic and sedimentary controls (Eocene,Tremp–Jaca Basin,Southern Pyrenees,NE Spain)

The propagation of the deformation front in foreland systems is typically accompanied by the incorporation of parts of the basin into wedge‐top piggy‐back basins, this process is likely producing considerable changes to sedimentation rates (SR). Here we investigate the spatial‐temporal evolution of SR for the Tremp–Jaca Basin in the Southern Pyrenees during its evolution from a wedge‐top, foreredeep, forebulge configuration to a wedge‐top stage. SR were controlled by a series of tectonic structures that influenced subsidence distribution and modified the sediment dispersal patterns. We compare the decompacted SR calculated from 12 magnetostratigraphic sections located throughout the Tremp–Jaca Basin represent the full range of depositional environment and times. While the derived long‐term SR range between 9.0 and 84.5 cm/kyr, compiled data at the scale of magnetozones (0.1–2.5 Myr) yield SR that range from 3.0 to 170 cm/kyr. From this analysis, three main types of depocenter are recognized: a regional depocenter in the foredeep depozone; depocenters related to both regional subsidence and salt tectonics in the wedge‐top depozone; and a depocenter related to clastic shelf building showing transgressive and regressive trends with graded and non‐graded episodes. From the evolution of SR we distinguish two stages. The Lutetian Stage (from 49.1–41.2 Ma) portrays a compartmentalized basin characterized by variable SR in dominantly underfilled accommodation areas. The markedly different advance of the deformation front between the Central and Western Pyrenees resulted in a complex distribution of the foreland depozones during this stage. The Bartonian–Priabonian Stage (41.2–36.9 Ma) represents the integration of the whole basin into the wedge‐top, showing a generalized reduction of SR in a mostly overfilled relatively uniform basin. The stacking of basement units in the hinterland during the whole period produced unusually high SR in the wedge‐top depozone.     

DA型脉动白矮星G255—2的双站联合观测和初步的星震学分析

介绍中法双方分别于１９９５年和１９９７年对ＤＡ型脉动白矮星Ｇ２５５—２进行的两次高精度的光电测光联合观测．在仔细的数据处理之后进行了多重频率脉动分析,得到了该星的本征脉动模式．同时,结合ＤＡ型脉动白矮星的理论模型进行了初步的星震学分析．     

Sediment Facies of a Low-Energy, Meso-Tidal, Fringing Reef, Singapore

Pulau (P.) Semakau comprises a narrow, low-lying island, surrounded by a wide fringing coral reef of late-Holocene age. The modern reef flat comprises a gently sloping surface related to modern mean low water neap tide level. Six sediment facies are recognised below the line of high water of which three, adjacent to the island, are composed of terrigenous minerals and rock fragments and three are autochthonous carbonates sediments developed since local sea-level still stand (c. 6,500 years BP). The fringing reefs of P. Semakau and Singapore differ from many of the reef forms recorded from Peninsular Malaysia waters and the Gulf of Thailand. The Singapore reefs have wide, well-developed intertidal reef flats, lack lagoons and reef crests, and have a very steep reef slope. In contrast to most other fringing reefs of the region, the Singapore reefs have developed in a low wave energy, meso-tidal environment lacking strong environmental gradients.     

Rapid inversion of data from 2D resistivity surveys with electrode displacements

Resistivity monitoring surveys are used to detect temporal changes in the subsurface using repeated measurements over the same site. The positions of the electrodes are typically measured at the start of the survey program and possibly at occasional later times. In areas with unstable ground, such as landslide‐prone slopes, the positions of the electrodes can be displaced by ground movements. If this occurs at times when the positions of the electrodes are not directly measured, they have to be estimated. This can be done by interpolation or, as in recent developments, from the resistivity data using new inverse methods. The smoothness‐constrained least squares optimisation method can be modified to include the electrode positions as additional unknown parameters. The Jacobian matrices with the sensitivity of the apparent resistivity measurements to changes in the electrode positions are then required by the optimisation method. In this paper, a fast adjoint‐equation method is used to calculate the Jacobian matrices required by the least squares method to reduce the calculation time. In areas with large near‐surface resistivity contrasts, the inversion routine sometimes cannot accurately distinguish between electrode displacements and subsurface resistivity variations. To overcome this problem, the model for the initial time‐lapse dataset (with accurately known electrode positions) is used as the starting model for the inversion of the later‐time dataset. This greatly improves the accuracy of the estimated electrode positions compared to the use of a homogeneous half‐space starting model. In areas where the movement of the electrodes is expected to occur in a fixed direction, the method of transformations can be used to include this information as an additional constraint in the optimisation routine.     

The inversion of data from very large three-dimensional electrical resistivity tomography mobile surveys

New developments in mobile resistivity meter instrumentation have made it possible to survey large areas with dense data coverage. The mobile system usually has a limited number of electrodes attached to a cable that is pulled along behind an operator so that a large area can be covered within a short time. Such surveys can produce three-dimensional datasets with hundreds of thousands of electrodes positions and data points. Similarly, the inverse model used to interpret the data can have several hundred thousand cells. It is impractical to model such large datasets within a reasonable time on microcomputers used by many small companies employing standard inversion techniques. We describe a model segmentation technique that subdivides the finite-element mesh used to calculate the apparent resistivity and Jacobian matrix values into a number of smaller meshes. A fast technique that optimizes the calculation of the Jacobian matrix values for multi-channel systems was also developed. A one-dimensional wavelet transform method was then used to compress the storage of the Jacobian matrix, in turn reducing the computer time and memory required to solve the least-squares optimization equation to determine the inverse model resistivity values. The new techniques reduce the calculation time and memory required by more than 80% while producing models that differ by less than 1% from that obtained using the standard inversion technique with a single mesh. We present results using a synthetic model and a field dataset that illustrates the effectiveness of the proposed techniques.     

Extraction of agar from Eucheuma cottonii and Gelidium amansii seaweeds with sonication pretreatment using autoclaving method

The effect of sonication pretreatment condition on Eucheuma cottonii and Gelidium amansii seaweed towards agar extraction wae studied. Four parameters were changed during sonication to investigate the effects on agar yield and quality. These parameters include the time interval, concentration ratio, frequency, and intensity. The highest amount of agar extracted from Eucheuma cottonii species could be obtained from the time interval of 30 min, seaweed weight to solvent volume ratio of 1:20, the frequency of 35 Hz, and the sonication power intensity of 30%. For Gelidium amansii species, the best agar yield also could be obtained from the time interval of 30 min, 1:20 of seaweed weight to water volume ratio, the frequency of 35 Hz, and power intensity of 30%. From the experiment, sonication pretreatment significantly influenced the yield and properties of extracted agar. The sonication with autoclaved seaweed produced agar containing less sulfate content, which is an excellent chemical property for gel electrophoresis applications.The gel strength of sonication with autoclaving for both seaweeds, Eucheuma and Gelidium species was the highest among those by sonication with direct heating, which proved that sonication pretreatment with autoclaving could enhance the physical properties of the agar.     

High-resolution simulation and characterization of density-driven flow in CO2 storage in saline aquifers

Simulations are routinely used to study the process of carbon dioxide (CO₂) sequestration in saline aquifers. In this paper, we describe the modeling and simulation of the dissolution–diffusion–convection process based on a total velocity splitting formulation for a variable-density incompressible single-phase model. A second-order accurate sequential algorithm, implemented within a block-structured adaptive mesh refinement (AMR) framework, is used to perform high-resolution studies of the process. We study both the short-term and long-term behaviors of the process. It is found that the onset time of convection follows closely the prediction of linear stability analysis. In addition, the CO₂ flux at the top boundary, which gives the rate at which CO₂ gas dissolves into a negatively buoyant aqueous phase, will reach a stabilized state at the space and time scales we are interested in. This flux is found to be proportional to permeability, and independent of porosity and effective diffusivity, indicative of a convection-dominated flow. A 3D simulation further shows that the added degrees of freedom shorten the onset time and increase the magnitude of the stabilized CO₂ flux by about 25%. Finally, our results are found to be comparable to results obtained from TOUGH2-MP.     

Provenance-related characteristics of beach sediments around the island of Menorca, Balearic Islands (western Mediterranean)

The island of Menorca, one of the Balearic Islands (Spain) located in the western Mediterranean, is characterised by a contrasting geology and landscape with two major geographic domains: (1) a southern region called Migjorn, comprised of Late Miocene calcarenites and limestones, and (2) a northern region known as Tramuntana, which is composed of folded and faulted Palaeozoic, Mesozoic and Tertiary (Oligocene) siliceous and calcareous rocks. Both domains are lined by numerous pocket beaches exhibiting a high variety of surficial sediment assemblages. Grain-size and compositional analyses revealed that cliff erosion and nearshore Posidonia oceanica meadows are the main sources of sediments consisting mostly of medium- to coarse-grained carbonate sands of marine biogenic origin, with variable amounts of terrigenous rock fragments and quartz. Based on distinctly different contributions of bioclastic material, biogenic carbonates and quartz, 320 sediment samples from 64 beaches were grouped into different facies associations dominated by either (1) biogenic sands, (2) biogenic sands with terrigenous contributions or (3) terrigenous sands with quartz. Nevertheless, there is a marked regional variability in sediment texture and composition. Thus, variable mixtures of carbonate and siliciclastic sediments characterise the beaches of the northern region, whereas the beaches of the southern region are composed mostly of carbonate sands of marine biogenic origin. An exception is the central sector of the south coast, which is enriched in quartz sand (~10 %); this can be related to outcrops of quartz-rich basement rock and also to rocks exposed in some northern drainage basins captured by southern streams since the Plio-Quaternary.