Acoustoelasticity in cracked solids

A new model that accounts for the stress dependence of the phase velocity of elastodynamic waves propagating in a cracked solid under compression is presented. The phase velocities of longitudinal and shear waves are derived from the effective elastic properties of a cracked solid, which are evaluated within the framework of Kachanov's approach. Following Kachanov, the extra-compliance tensor of the cracked solid is related to the crack compliances, which display a marked non-linear behaviour when subjected to a compressive load. Such non-linear behaviour is shown to be derived from the elastic interaction between the contacting crack faces under compression. This work does not address the effect of mutual interaction among cracks and the generation of higher harmonics due to the medium non-linearity. Numerical examples are presented that illustrate the phase velocity changes occurring in a solid with a random distribution of parallel cracks as a function of an external compressive load. A distinctive feature of the acoustoelastic effect in solids with large parallel fractures and in solids with distributions of aligned microcracks is also illustrated.     

Study on Phase Behavior of CO2/Hydrocarbons in Shale Reservoirs Considering Sieving Effect and Capillary Pressure

The phase behavior of fluid is essential for predicting ultimate oil recovery and determining optimal production parameters. The pore size in shale porous media is nanopore, which causes different phase behaviors of fluid in unconventional reservoirs. Nanopores in shale media can be regard as semipermeable membrane to filter heavy components (sieving effect) in shale oil, which leads to the different distributions of fluid components and different phase behaviors. In addition, the phase behavior of fluid in nanopores can be significantly altered by large capillary pressure. In this paper, the phase behavior of fluid in shale reservoirs is investigated by a new two-phase flash algorithm considering sieving effect and capillary pressure. Firstly, membrane efficiency and capillary pressure are introduced to establish a thermodynamic equilibrium model that is solved by Rachford–Rice flash calculation and Newton–Raphson method. The capillary pressures in different pore sizes are calculated by the Young–Laplace equation. Then, the influences of sieving effect and capillary pressure on phase behavior are analyzed. The results indicate that capillary pressure can suppress the bubble point pressure of fluid in nanopores. The distributions of fluid components are different in various parts of shale media. In the unfiltered part, density and viscosity of fluid are higher. Finally, it is found that the membrane efficiency can be improved by CO₂ injection. The minimum miscibility pressure for shale oil–CO₂ system is also studied. The developed model provides a better understanding of the phase behavior of fluid in shale oil reservoirs.

     

Propagation of harmonic plane waves in a general anisotropic porous solid

Wave propagation is studied in a general anisotropic poroelastic solid. The presence of dissipation due to fluid-viscosity as well as hydraulic anisotropy of pore permeability are also considered. Biot's theory is used to derive a system of modified Christoffel equations for the propagation of plane harmonic waves in porous media. A non-trivial solution of this system is ensured by a determinantal equation. This equation is separated into two different polynomial equations. One is the quartic equation whose roots represent the complex velocities of four attenuating waves in the medium. The other is a eighth-degree polynomial whose roots represent the vertical slowness values for the four waves propagating upward and downward in a finite porous medium. Procedure is explained to associate the numerically obtained roots with the waves propagating in the medium. The slowness surfaces of waves reflected at the boundary of the medium are computed for a realistic numerical model. The behaviours of phase velocity surfaces are analysed with the help of numerical examples.     

Quadrangle-grid velocity–stress finite difference method for poroelastic wave equations

A quadrangle-grid velocity–stress finite difference method, based on a first-order hyperbolic system that is equivalent to Biot's equations, is developed for the simulation of wave propagation in 2-D heterogeneous porous media. In this method the velocity components of the solid material and of the pore fluid relative to that of the solid, and the stress components of three solid stresses and one fluid pressure are defined at different nodes for a staggered non-rectangular grid. The scheme uses non-orthogonal grids, allowing surface topography and curved interfaces to be easily modelled in the numerical simulation of seismic responses of poroelastic reservoirs. The free-surface conditions of complex geometry are achieved by using integral equilibrium equations on the surface, and the source implementations are simple. The algorithm is an extension of the quadrangle-grid finite difference method used for elastic wave equations.     

Existence of transverse waves in anisotropic poroelastic media

Out of the four waves in an anisotropic poroelastic medium, two are termed as quasi-transverse waves. The prefix 'quasi' refers to their polarizations being nearly, but not exactly, perpendicular to direction of propagation. In this composite medium, unlike perfectly elastic medium, the propagation of a longitudinal wave along a phase direction may not be accompanied by transverse waves. The existence of a transverse wave in anisotropic poroelastic media is ensured by the two equations restricting the choice of elastic coefficients of porous aggregate as well as fluid–solid coupling. Necessary and sufficient conditions for the existence of transverse waves along the coordinate axes and in the coordinate planes for general anisotropy are discussed. The discussion is extended to the case of orthotropic materials and existence for few specific phase directions is also explored. The conditions for the transverse waves decided on the basis of their apparent polarizations, that is, particle motion being perpendicular to ray direction, are also discussed. For a particular numerical model, the existence of these apparent transverse waves is solved numerically for phase directions in coordinate planes. For general directions of phase propagation, the existence of these transverse waves is checked graphically for the chosen numerical model.     

Estimability analysis of geodetic, astrometric and geodynamical quantities in very long baseline interferometry

Summary. Within the framework of Newtonian kinematics VLBI observations are analysed with respect to estimability of geodetic and astrometric quantities. An Earth model of either rigid or deformable type is designed; instrumental clock offsets and clock drifts are included. Observational patterns are studied in all detail reviewed in seven tables. Appendix A is an introduction to the set-up of the observational model for a deformable Earth both in terms of coordinate-free and coordinate-related geometry. Appendix B is a study of the invariance characteristics of VLBI observations. Interrelations of three fundamental quantities, length unit, time unit and velocity of light are discussed. An overall result of an Earth model of deformable type is the need of simultaneous observations to more than one source; VLBI time delay observations cannot distinguish between secular changes of network size (expansion or shrinking) and a common secular drift (deceleration or acceleration) of clocks used.     

基于多代理模型的城市土地利用博弈模拟

以美国达拉斯北部地区为例,通过建立基于多代理的人类系统和自然系统的耦合互动模型,研究不同土地管理策略的有效性.自然系统包含土地覆盖变化、流域水文动力学和野生动物栖息地生态系统.人类系统包括土地利益相关者(地主、屋主、开发商、政府)的影响土地利用变化价值观.系统强调自然对土地利用决策的反馈机制,不同类型代理相对空间位置的重要作用,以及同一类别中不同类型代理在特定情形下的相互转化.研究进行了两类模拟,一是模拟初始化的若干情景,包括不同类型地主和户主的分布、不同类型政府以及经济模型的假定;二是模拟主动式土地政策的实施效果.模拟结果表明,当政府购置保留地时考虑土地拥有者的价值观,较之单纯基于机会或生态因素更能导致有效的增长管理策略.     

Parallelization of the flow-path network model using a particle-set strategy

High-performance simulation of flow dynamics remains a major challenge in the use of physical-based, fully distributed hydrologic models. Parallel computing has been widely used to overcome efficiency limitation by partitioning a basin into sub-basins and executing calculations among multiple processors. However, existing partition-based parallelization strategies are still hampered by the dependency between inter-connected sub-basins. This study proposed a particle-set strategy to parallelize the flow-path network (FPN) model for achieving higher performance in the simulation of flow dynamics. The FPN model replaced the hydrological calculations on sub-basins with the movements of water packages along the upstream and downstream flow paths. Unlike previous partition-based task decomposition approaches, the proposed particle-set strategy decomposes the computational workload by randomly allocating runoff particles to concurrent computing processors. Simulation experiments of the flow routing process were undertaken to validate the developed particle-set FPN model. The outcomes of hourly outlet discharges were compared with field gauged records, and up to 128 computing processors were tested to explore its speedup capability in parallel computing. The experimental results showed that the proposed framework can achieve similar prediction accuracy and parallel efficiency to that of the Triangulated Irregular Network (TIN)-based Real-Time Integrated Basin Simulator (tRIBS).     

我国温带荒漠生物土壤结皮孔隙结构分布特征

对原状生物土壤结皮孔隙结构、孔隙率及其垂直连通性进行室内微观扫描与计算机断层摄影成像分析。结果表明,生物土壤结皮团聚体特征明显,其孔隙率的增大主要是由大孔隙数量的增加而形成,同时,土壤孔隙连通性也随大孔隙数量的增加而增大,生物土壤结皮的发育导致表层土壤具有明显的极其复杂的土壤结构演化特征。     

Seismic wave properties in time-dependent porosity homogeneous media

It is quantified the properties of seismic waves in fully saturated homogeneous porous media within the framework of Sahay's modified and reformulated poroelastic theory. The computational results comprise amplitude attenuation, velocity dispersion and seismic waveforms. They show that the behaviour of all four waves modelled as a function of offset, frequency, porosity, fluid viscosity and source bandwidth depicts realistic dissipation within the sonic–ultrasonic band. Therefore, it appears that there is no need to include material heterogeneity to model attenuation. By inference it is concluded that the fluid viscosity effects may be enhanced by dynamic porosity.     

A general-purpose parallel raster processing programming library test application using a geographic cellular automata model

A general-purpose parallel raster processing programming library (pRPL) was developed and applied to speed up a commonly used cellular automaton model with known tractability limitations. The library is suitable for use by geographic information scientists with basic programming skills, but who lack knowledge and experience of parallel computing and programming. pRPL is a general-purpose programming library that provides generic support for raster processing, including local-scope, neighborhood-scope, regional-scope, and global-scope algorithms as long as they are parallelizable. The library also supports multilayer algorithms. Besides the standard data domain decomposition methods, pRPL provides a spatially adaptive quad-tree-based decomposition to produce more evenly distributed workloads among processors. Data parallelism and task parallelism are supported, with both static and dynamic load-balancing. By grouping processors, pRPL also supports data–task hybrid parallelism, i.e., data parallelism within a processor group and task parallelism among processor groups. pSLEUTH, a parallel version of a well-known cellular automata model for simulating urban land-use change (SLEUTH), was developed to demonstrate full utilization of the advanced features of pRPL. Experiments with real-world data sets were conducted and the performance of pSLEUTH measured. We conclude not only that pRPL greatly reduces the development complexity of implementing a parallel raster-processing algorithm, it also greatly reduces the computing time of computationally intensive raster-processing algorithms, as demonstrated with pSLEUTH.     

Numerical modelling of wavefields in three-dimensional inhomogeneous media

Summary. Numerical modelling is one of the most efficient methods for an investigation of the relationship between structural features and peculiarities of observed wavefields. It is practically the only method for 2-D and 3-D inhomogeneous media.
An algorithm based on ray theory has been developed for calculations of travel times and amplitudes of seismic waves in 3-D inhomogeneous media with curved interfaces. It was applied for numerical modelling of kinematic and dynamic characteristics of seismic waves propagating in laterally inhomogeneous media.
Travel-time and amplitude patterns were studied in the 2-D and 3-D models of a geosyncline, in which velocity distribution was given by an analytical function of the coordinates. For a more complicated model representing a subducting high-velocity lithospheric plate in a transition zone between oceanic and continental upper mantle, the velocity distribution was given by discrete values on a 2-D non-rectangular grid. It was shown that when a source was placed above the lithospheric plate, a shadow zone appeared along a strike of the structure, i.e. in the direction which is perpendicular to a strong lateral velocity gradient. Travel-time residuals were calculated along the seismological profile for a 3-D velocity distribution in the upper mantle beneath Central Asia, obtained as a result of inversion of travel times by the Backus-Gilbert method. They were found to be in a good agreement with the observed data.     

Predicted present-day evolution patterns of ice thickness and bedrock elevation over Greenland and Antarctica

This paper discusses predicted evolution patterns of present-day changes of ice thickness, surface elevation, and bedrock elevation over the Greenland and Antarctic continents. These were obtained from calculations with dynamic 3-D ice sheet models which were coupled to a visco-elastic solid Earth model. The experiments were initialized over the last two glacial cycles and subsequently averaged over the last 200 years to obtain the current evolution. The calculations indicate that the Antarctic Ice Sheet is still adjusting to the last glacial-interglacial transition yielding a decreasing ice volume and a rising bedrock elevation of the order of several centimetres per year. The Greenland Ice Sheet was found to be close to a stationary state with a mean thickness change of only a few millimetres per year, but the calculations revealed large spatial differences. Predicted patterns over Greenland are characterized by a small thickening over the ice sheet interior and a general thinning of the ablation area. In Antarctica, almost all of the predicted changes are concentrated in the West Antarctic Ice Sheet, which is still retreating at both the Weddell and Ross Sea margins. Over most of both ice sheets, the model indicates that the surface elevation trend is dominated by ice thickness changes rather than by bedrock elevation changes.     

Multi-level Optimization of Reservoir Scheduling Using Multi-resolution Wavelet-Based Up-scaled Models

This paper presents a multi-level procedure for production and injection scheduling through a numerical model-based optimization of well control variables. To calculate the net present value (NPV), the objective function of optimization, this procedure uses a number of discretized systems for a reservoir model with different degrees of up-scaling prepared according to a multi-resolution wavelet technique. These up-scaled models were incorporated into optimization based on a probability function. In early optimization iterations, due to the necessity to explore the search space quickly, the coarsest grid model has a higher chance for selection than the others; however, by a selection (with a low probability) of the finest up-scaled grid model in these iterations, solutions and objective function were tuned. In the later iterations of optimization, the finest up-scaled grid model probability was the highest in order to ensure the reliability of the final solution. The optimization algorithm is an adaptive simulated annealing algorithm coupled with a polytope. This procedure was evaluated in two case studies. The first case study was a horizontal 2D oil model with water flooding. The second case study was a vertical 2D oil model with gas injection. The results show that the proposed optimization procedure provides approximately the same accuracy compared to the situation in which the fine grid model is used for all the optimization iterations. Also, the run-time for the proposed optimization procedure is comparable to the run-time of the optimization in which only the coarsest grid model is used to calculate objective function. Moreover, the superiority of the wavelet-based up-scaling over an analogous multiple grid system optimization using uniformly up-scaled models is presented.

     

The role of biogeochemical cycling for the formation and preservation of varved sediments in Soppensee (Switzerland)

We analyzed seasonally aggregated observations of temperature, conductivity, dissolved oxygen and dissolved inorganic carbon from Soppensee (District of Lucerne, Switzerland) for the yrs 1980 to 1993. Holomictic Soppensee is characterized by a strong summer stratification with a thin epilimnion separated from an anoxic hypolimnion by a strong pycnocline formed by thermal and chemical gradients. A vertical one-dimensional model was developed to simulate the observed seasonal cycles of carbon and oxygen. The processes of net community production, mineralization of organic matter, precipitation and dissolution of calcite, gas exchange, in- and outflow, sedimentation and vertical eddy diffusion are included. According to the model, the annual net community production is estimated to about 110 g C m^-2 yr^-1 and the annual net primary production to about 330 to 440 g C m^-2 yr^-1, which is a typical value for eutrophic lakes. A mass balance of the carbon cycle indicates that most of the inflow comes from groundwater which is super-saturated with respect to atmospheric CO₂. Therefore the surface waters exhibit a large capacity for calcite precipitation. The results of the model are used to constrain the conditions that favor the formation of varved sediments in Soppensee during thousands of yrs. Model calculations show that the deep waters would still turn anoxic even if the sedimentation rate of organic matter were decreased to 25%. Several physical factors such as biogenic stabilization of the deep waters due to calcite dissolution and low input of wind energy are responsible for the long term anoxia in Soppensee.     

Accessing Flexibility: A Historical Institutional Analysis of Water Use for Oil and Gas Development in Colorado

US Mountain West Water access and allocation institutions have a history of adapting policy and practice to increase flexibility for diverse water uses. We examine how flexible access has developed over time and space by operationalizing the historical institutional (HI) theoretical and methodological framework. We trace historical water access for oil and gas (OG) development in Colorado, working from contemporary water right data to examine historical critical events, policies, and political contexts. OG water use has iteratively shaped water governance institutions in the top OG producing regions of Colorado, Weld, and Garfield Counties. The analysis suggests that to more accurately capture institutional change and continuity in resource allocation systems, an analysis of informal institutions is an essential theoretical contribution to the HI framework. While increased flexibility makes multiple uses easier, policies favor the most economically lucrative beneficial uses and generate issues of transparency, an important consideration for the public’s resource. Future practices of flexibility are contingent on market structures and institutional access mechanisms shaped during previous government policy processes, illuminating the value of the HI framework to inform future water policy.     

Using adaptively coupled models and high-performance computing for enabling the computability of dust storm forecasting

Forecasting dust storms for large geographical areas with high resolution poses great challenges for scientific and computational research. Limitations of computing power and the scalability of parallel systems preclude an immediate solution to such challenges. This article reports our research on using adaptively coupled models to resolve the computational challenges and enable the computability of dust storm forecasting by dividing the large geographical domain into multiple subdomains based on spatiotemporal distributions of the dust storm. A dust storm model (Eta-8bin) performs a quick forecasting with low resolution (22 km) to identify potential hotspots with high dust concentration. A finer model, non-hydrostatic mesoscale model (NMM-dust) performs high-resolution (3 km) forecasting over the much smaller hotspots in parallel to reduce computational requirements and computing time. We also adopted spatiotemporal principles among computing resources and subdomains to optimize parallel systems and improve the performance of high-resolution NMM-dust model. This research enabled the computability of high-resolution, large-area dust storm forecasting using the adaptively coupled execution of the two models Eta-8bin and NMM-dust.     

Mineral and Lithology Mapping of Drill Core Pulps Using Visible and Infrared Spectrometry

A novel approach for using field spectrometry for determining both the mineralogy and the lithology of drill core pulps (powders) is developed and evaluated. The methodology is developed using material from a single drillhole through a mineralized sequence of rocks from central New South Wales. Mineral library spectra are used in linear unmixing routines to determine the mineral abundances in drill core pulps that represent between 1 m and 3 m of core. Comparison with X-Ray Diffraction (XRD) analyses shows that for most major constituents, spectrometry provides an estimate of quantitative mineralogy that is as reliable as that provided by XRD. Confusion between the absorption features of calcite and those of chlorite causes the calcite contents determined by spectrometry to be unreliable. Convex geometry is used to recognize the spectra of those samples that are extreme and are representative of unique lithologies. Linear unmixing is used to determine the abundance of these lithologies in each drillhole sample and these abundances are used to interpret the geology of the drillhole. The interpreted geology agrees well with conventional drillhole logs of the visible geology and photographs of the split core. The methods developed provide a quick and cost-effective way of determining the lithology and alteration mineralogy of drill core pulps.     

沙尘起动初期近地面浓度分布数值模拟研究

对沙尘事件发生初始阶段近地面空间的风沙运动进行了数值模拟研究。考虑固相的剪切黏性、总体黏性、气固相间动量交换、升力等,建立了一个二维气固两相流动数学模型。采用Fluent流体力学软件,进行数值计算。对一组风洞实验数据进行了模拟计算,验证了本模型和计算方法的正确性。研究区域高500 m、长1 000 m。近地层风速廓线采用指数式。沙尘粒度分布用R-R分布;粒径5~150 μm。模拟计算得出近地层内不同风速、沙粒直径的风沙速度及沙尘体积浓度空间分布。结果表明,沙尘起动初期分为发展壮大、回落、稳定3个阶段。沙尘卷起的发展段长度、沙尘卷扬高度及浓度分布曲线形状,均随沙粒和入口风速而明显变化。到达稳定阶段,沙尘浓度沿高度的分布表现为,近地面存在一个均匀的浓度层,即饱和层,其浓度、高度取决于风速和沙粒粒径;饱和层以上沙尘浓度呈指数衰减。     

Insights from a Simple Hotelling Model of the World Oil Market

This paper uses annual data on world oil price and consumption from 1965 to 2006 to calibrate a Hotelling model of optimal nonrenewable resource extraction. Numerical solutions are generated for various specifications of the elasticity of demand for both isoelastic demand and linear demand under each of two possible market structures: perfect competition and monopoly. Prior to the 1973 oil crisis, the model that best fits actual data is one of perfect competition with linear demand and a demand elasticity of −0.4. For the periods 1973–1981 and 1981–1990, the model that best fits actual data is one of monopoly with linear demand and demand elasticities of −0.8 and −0.7, respectively, suggesting that the market was strongly influenced by OPEC during this time. Under the model that best fits the most recent period (perfect competition with linear demand and demand elasticity −0.5), the real oil price (in 1982–1984 U.S.$) should fall in the range $60.87–$66.31/barrel over the years 2010–2030.