Penetration calculations and measurements for a layered soil target     

Rupert K. Byers  Albert J. Chabai 《国际地质力学数值与分析法杂志》1977,1(2):107-138

A tow-dimensional finite difference wave propagation code was used to analyse the impact of a rigid, ogival-nosed penetrator with a target consisting of thin alternating layers of silt, sand, and clay. The response of the idealized target was described with an elastic-plastic constitutive model depending upon two stress invariants and the history of plastic deformation. Interfacial friction between the penetrator and target was assumed to be negligible. Comparisons are made of calculated results and those of nominally similar experiments conducted at the Watching Hill test site in Alberta, Canada. Analysis of these comparisons reveals that numerical methods like the one employed in this study can yield insight into penomena, as well as suggesting possible improvements in the calculation technique.  相似文献   

Waterflooding optimization with the INSIM-FT data-driven model     

Zhenyu Guo  Albert C. Reynolds  Hui Zhao 《Computational Geosciences》2018,22(3):745-761

In a recent paper, we developed a physics-based data-driven model referred to as INSIM-FT and showed that it can be used for history matching and future reservoir performance predictions even when no prior geological model is available. The model requires no prior knowledge of petrophysical properties. In this work, we explore the possibility of using INSIM-FT in place of a reservoir simulation model when estimating the well controls that optimize water flooding performance where we use the net present value (NPV) of life-cycle production as our cost (objective) function. The well controls are either the flowing bottom-hole pressure (BHP) or total liquid rates at injectors and producers on the time intervals which represent the prescribed control steps. The optimal well controls that maximize NPV are estimated with an ensemble-based optimization algorithm using the history-matched INSIM-FT model as the forward model. We compare the optimal NPV obtained using INSIM-FT as the forward model with the estimate of the optimal NPV obtained using the correct full-scale reservoir simulation model when performing waterflooding optimization.  相似文献   

Simulating observations of dark matter dominated galaxies: towards the optimal halo profile     

W. J. G. de Blok  Albert Bosma  Stacy McGaugh 《Monthly notices of the Royal Astronomical Society》2003,340(2):657-678

  相似文献   

Paleoseismology of the 2016 MW 6.1 Petermann earthquake source: Implications for intraplate earthquake behaviour and the geomorphic longevity of bedrock fault scarps in a low strain-rate cratonic region     

Tamarah Rosellen King  Mark Quigley  Dan Clark  Albert Zondervan  Jan-Hendrik May  Abaz Alimanovic 《地球表面变化过程与地形》2021,46(7):1238-1256

The 20 May 2016 M_W 6.1 Petermann earthquake in central Australia generated a 21 km surface rupture with 0.1 to 1 m vertical displacements across a low-relief landscape. No paleo-scarps or potentially analogous topographic features are evident in pre-earthquake Worldview-1 and Worldview-2 satellite data. Two excavations across the surface rupture expose near-surface fault geometry and mixed aeolian-sheetwash sediment faulted only in the 2016 earthquake. A 10.6 ± 0.4 ka optically stimulated luminescence (OSL) age of sheetwash sediment provides a minimum estimate for the period of quiescence prior to 2016 rupture. Seven cosmogenic beryllium-10 (¹⁰Be) bedrock erosion rates are derived for samples < 5 km distance from the surface rupture on the hanging-wall and foot-wall, and three from samples 19 to 50 km from the surface rupture. No distinction is found between fault proximal rates (1.3 ± 0.1 to 2.6 ± 0.2 m Myr⁻¹) and distal samples (1.4 ± 0.1 to 2.3 ± 0.2 m Myr⁻¹). The thickness of rock fragments (2–5 cm) coseismically displaced in the Petermann earthquake perturbs the steady-state bedrock erosion rate by only 1 to 3%, less than the erosion rate uncertainty estimated for each sample (7–12%). Using ¹⁰Be erosion rates and scarp height measurements we estimate approximately 0.5 to 1 Myr of differential erosion is required to return to pre-earthquake topography. By inference any pre-2016 fault-related topography likely required a similar time for removal. We conclude that the Petermann earthquake was the first on this fault in the last ca. 0.5–1 Myr. Extrapolating single nuclide erosion rates across this timescale introduces large uncertainties, and we cannot resolve whether 2016 represents the first ever surface rupture on this fault, or a > 1 Myr interseismic period. Either option reinforces the importance of including distributed earthquake sources in fault displacement and seismic hazard analyses.  相似文献   

A comparison of hillslope drainage area estimation methods using high-resolution DEMs with implications for topographic studies of gullies     

Simon J. Walker  Albert I. J. M. van Dijk  Scott N. Wilkinson  Peter B. Hairsine 《地球表面变化过程与地形》2021,46(11):2229-2247

Topographic models provide a useful tool for understanding gully occurrence in the landscape but require reliable estimates of gully head drainage areas. Modern high-resolution topography data (collected using structure from motion photogrammetry or light detection and ranging) is increasingly used for topographic studies of gullies, but little work has been done to assess the variability of gully head drainage area estimates using different methods. This study evaluated alternative approaches to using high-resolution digital elevation models (DEMs) so that gully topographic models can be more readily applied to any area with suitably high-resolution data. Specifically, we investigated the impact of single- or multiple-direction flow routing algorithms, DEM hydrologic-enforcement procedures and spatial resolution on gully head drainage area estimation. We tested these methods on a 40 km² site centred on Weany Creek, a low-relief semi-arid landscape draining towards the Great Barrier Reef, Australia. Using a subroutine to separate gully heads into those with divergent or convergent flow patterns upslope, we found that divergent flow conditions occurred at half of 484 studied gullies. Drainage areas estimated by different flow routing algorithms were more variable in these divergent cases than for convergent cases. This variation caused a significant difference between topographic threshold parameters (slope b and intercept k) derived from single- or multiple-direction flow routing algorithms, respectively. Different methods of hydrologic enforcement (filling or breaching) also affected threshold analysis, resulting in estimates of the exponent b being ~188% higher if the DEM was filled than if breached. The testing of the methods to date indicates that a finer resolution (≤2 m) DEM and a multiple-direction flow routing algorithm achieve the most realistic drainage area estimates in low-relief landscapes. For Weany Creek we estimated threshold parameters k = 0.033 and b = 0.189, indicating that it is highly susceptible to gully erosion.  相似文献   

Investigation of the sampling performance of ensemble-based methods with a simple reservoir model   总被引：1，自引：0，他引：1  

Alexandre A. Emerick  Albert C. Reynolds 《Computational Geosciences》2013,17(2):325-350

The application of the ensemble Kalman filter (EnKF) for history matching petroleum reservoir models has been the subject of intense investigation during the past 10 years. Unfortunately, EnKF often fails to provide reasonable data matches for highly nonlinear problems. This fact motivated the development of several iterative ensemble-based methods in the last few years. However, there exists no study comparing the performance of these methods in the literature, especially in terms of their ability to quantify uncertainty correctly. In this paper, we compare the performance of nine ensemble-based methods in terms of the quality of the data matches, quantification of uncertainty, and computational cost. For this purpose, we use a small but highly nonlinear reservoir model so that we can generate the reference posterior distribution of reservoir properties using a very long chain generated by a Markov chain Monte Carlo sampling algorithm. We also consider one adjoint-based implementation of the randomized maximum likelihood method in the comparisons.  相似文献   

Effects of wave-wave and wave-mean flow interactions on the evolution of a baroclinic wave     

Albert Barcilon  Terrence R. Nathan 《地球物理与天体物理流体动力学》2013,107(1-4):59-79

Abstract

The weakly nonlinear evolution of a free baroclinic wave in the presence of slightly supercritical, vertically sheared zonal flow and a forced stationary wave field that consists of a single zonal scale and an arbitrary number of meridional harmonics is examined within the context of the conventional two-layer model. The presence of the (planetary-scale) stationary wave introduces zonal variations in the supercriticality and is shown to alter the growth rate and asymptotic equilibrium of the (synoptic-scale) baroclinic wave via two distinct mechanisms: The first is due to the direct interaction of the stationary wave with the shorter synoptic wave (wave-wave mechanism), and the second is due to the interaction of the synoptic wave with that portion of the mean field that is corrected by the zonally rectified stationary wave fluxes (wave-mean mechanism). These mechanisms can oppose or augment each other depending on the amplitude and spatial structure of the stationary wave field. If the stationary wave field is confined primarily to the upper (lower) layer and consists of only the gravest cross-stream mode, conditions are favorable (unfavorable) for nonzero equilibrium of the free wave.

In addition to the time dependent heat flux generated by baroclinic growth of the free wave, its interaction with a stationary wave field consisting of two or more meridional harmonics generates time dependent heat fluxes that vary with period of the free wave. However, if the stationary wave field contains several meridional harmonics of sufficiently large amplitude, the free baroclinic wave is destroyed.  相似文献   

Three-dimensional poromechanical back analysis of the pulse test accounting for transverse isotropy   总被引：3，自引：2，他引：1  

Richard Giot  Albert Giraud  Théophile Guillon  Christophe Auvray 《Acta Geotechnica》2012,7(3):151-165

The pulse test is usually considered as an efficient test for measuring the permeability of weakly permeable porous media. It is generally interpreted by comparison of experimental and theoretical curves obtained by a simplified 1D analytical solution. In a previous work, we proposed an improvement for test interpretation, through fully coupled numerical modeling in 2D axisymmetric configurations. The numerical modeling took into account the poromechanical coupling and exhibited 2D effects that cannot be assessed with classical interpretation. The modeling was coupled with an inverse method permitting the assessment of intrinsic permeability, mechanical parameters such as Young??s modulus, and poromechanical coupling parameters such as Biot coefficient. The interpretation of the tests on argillites showed anisotropy of the material and revealed the fundamental importance of this anisotropy for further works. The present paper aims at taking into account the transverse isotropy of the rock samples in the poromechanical back analysis of the test. A transverse isotropic poroelastic constitutive law has been developed and implemented in Code_Aster (Edf), and 3D modeling is required to account for all possible orientations of the samples. The poromechanical modeling is still coupled with an inverse method, also implemented in the finite element code. The method was applied to laboratory tests on Meuse/Haute-Marne argillites, on cylindrical samples with axis either parallel or perpendicular to the isotropy planes, and the results were compared with isotropic interpretation. The 3D transverse isotropic modeling provides more consistent values of the parameters than the 2D modeling, particularly for the mechanical parameters.  相似文献   

Reconstruction of hydrocarbons accumulation in sediments affected by the oil refinery industry: the case of Tehuantepec Gulf (Mexico)     

Ana Carolina Ruiz-Fernández  Mario Sprovieri  Mauro Frignani  Joan Albert Sanchez-Cabeza  Maria Luisa Feo  Luca Giorgio Bellucci  Libia Hascibe Pérez-Bernal  Michel Preda  María Luisa Machain-Castillo 《Environmental Earth Sciences》2012,67(3):727-742

The Isthmus of Tehuantepec corresponds to the shortest distance (~200?km) between the Gulf of Mexico and the Pacific Ocean in Southern Mexico, and the main economical activity of this region is oil extraction and refining. Polycyclic aromatic hydrocarbons (PAHs) and total petroleum hydrocarbons (TPHs) were determined in a ²¹⁰Pb dated sediment core collected from the continental shelf of Tehuantepec Gulf, in the vicinity of the oil refinery of Salina Cruz, Oaxaca, the main oil refining facility of the country. The sediments were mostly of coarse nature and hence PAHs and TPHs concentrations throughout the core (61?C404???g?g^?1 and 29?C154?mg?kg^?1, respectively) were below international quality benchmarks. Depth profiles of both PAHs and TPHs concentrations showed increasing trends since the early 1900s but the higher values were found from the 1950s to present. PAH congener ratios showed that these contaminants had both petrogenic and pyrolitic sources, although the former has been predominant since the 1970s. The Salina Cruz refinery started operations in 1978 but the oil industry activities in the Tehuantepec Isthmus go back to the beginning of the twentieth century with the operation of Minatitlan refinery in the Gulf of Mexico, and the Gulf of Tehuantepec being the main conduit for oil distribution in the Pacific coast. The observed changes in contaminant distributions described well the oil industry development in the area.  相似文献   

How do early successional patterns in man-made wetlands differ between cold temperate and Mediterranean regions?     

Albert Ruhí  Jan Herrmann  Stéphanie Gascón  Jordi Sala  Dani Boix 《Limnologica》2012,42(4):328-339

The early stages of succession in newly created wetlands should be impacted by the region in which the wetland is located, since climate may have filtered the dominant biological strategies of the taxa leading this process and may condition their dynamics. We studied the early successional patterns of macroinvertebrates within man made ponds, located in cold temperate (Scandinavian Peninsula) and Mediterranean regions (Iberian Peninsula), during the first three years following their creation. We predicted (1) non random subsets of the regional species pool (deterministic assembly) guiding the successional process in cold temperate wetlands, and random (stochastic) assembly in the Mediterranean region; (2) higher successional rates in Mediterranean ponds than in the cold temperate ponds, with contrary episodes of highest and lowest change throughout the year; and (3) a significant difference in the composition of biological traits between regions, due to the dominance of traits adapted to hydrological variability in the Mediterranean region. Hypotheses on community structure (1) and the composition of dominant biological traits (3) were mostly supported – deterministic assembly mechanisms in the cold temperate ponds and stochastic assembly in the Mediterranean ponds; and a dominance of different biological traits between regions, explained by the need to overcome hydrological disturbances in the Mediterranean ponds. The dynamics of succession (2) were explained by climatic factors in the cold temperate region but not in the Mediterranean ponds. We suggest that the intrinsic hydrological variability of Mediterranean lentic systems may be a major factor driving community changes in man-made wetlands in this region. In order to generalize the observed patterns, we performed a meta-analysis of the temporal trends of taxonomic distinctness parameters of other successional studies across a latitudinal gradient in Europe, which supported the differences we had observed between latitudinal extremes.  相似文献