Study of elastic wave propagation in two-phase anisotropic media by numerical modeling of pseudospectral method

When there exists anisotropy in underground media, elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media, four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.     

地质雷达正演中的频散压制和吸收边界改进方法

从麦克斯韦方程组出发,建立了地质雷达的时域有限差分法（ＦＤＴＤ）数学模型,导出了理想 散关系和超级吸收边界条件。理想频散关系 考虑了ＦＤＴＤ法的收敛性和稳定性,也考虑了高频电磁波在Ｙｅｅ氏网格中的传播特点;超吸收边界条件则用磁场分量来提高电场分量精度。数值试验表明,理想频散关系能真实地反映雷达波在地下介质中的传播规律,超吸收边界条件能有效减小截断边界的伪反射,提高正演精度。将之应用于实际计算,取得     

非均质土壤饱和稳定流中盐分运移的传递函数模拟

对室内人工构造的两种非均质土柱,以传递函数模型作为模拟手段,研究了稳定流场中饱和非均质土壤盐分优先运移的随机特征。首先计算了模型参数μ和σ2,获得了氯离子在非均质土壤中迁移时间的概率密度函数,然后应用传递函数模型对土柱中氯离子的出流浓度动态进行了随机模拟,并对非均质土壤中氯离子的均值和中值迁移时间及相应的运移体积与可动体积进行了分析和讨论,还依据质量守衡原理获得了土壤溶液中氯离子平均驻留浓度的变化。     

数字高程模型信息提取与数字水文模型研究进展

回顾了数字高程模型(DEM)数据的信息提取方法,阐述了由DEM提取的信息在水文水资源领域应用的现状,探讨了数字模型在水文科学中的作用和数字水文在数字地球所处的地位及应用前景。     

工程地震勘探数值模拟研究

用傅氏变换法求解波动方程,对工程地震勘探进行了数值模拟计算,并研究了该算法的频散特性及其稳定性。     

Object and Pixel-Based Reservoir Modeling of a Braided Fluvial Reservoir

To assess differences between object and pixel-based reservoir modeling techniques, ten realizations of a UK Continental Shelf braided fluvial reservoir were produced using Boolean Simulation (BS) and Sequential Indicator Simulation (SIS). Various sensitivities associated with geological input data as well as with technique-specific modeling parameters were analyzed for both techniques. The resulting realizations from the object-based and pixel-based modeling efforts were assessed by visual inspection and by evaluation of the values and ranges of the single-phase effective permeability tensors, obtained through upscaling. The BS method performed well for the modeling of two types of fluvial channels, yielding well-confined channels, but failed to represent the complex interaction of these with sheetflood and other deposits present in the reservoir. SIS gave less confined channels and had great difficulty in representing the large-scale geometries of one type of channel while maintaining its appropriate proportions. Adding an SIS background to the Boolean channels, as opposed to a Boolean background, resulted in an improved distribution of sheetflood bodies. The permeability results indicated that the SIS method yielded models with much higher horizontal permeability values (20–100%) and lower horizontal anisotropy than the BS versions. By widening the channel distribution and increasing the range of azimuths, however, the BS-produced models gave results approaching the SIS behavior. For this reservoir, we chose to combine the two methods by using object-based channels and a pixel-based heterogeneous background, resulting in moderate permeability and anisotropy levels.     

Physics of Atoll Groundwater Flow

Large-scale steady-state groundwater flow in atoll carbonate platforms results from temperature and salinity-induced density gradients. Atolls are built on top of a basaltic substrate that provides geothermal heating from beneath. Moreover, they are immersed in the tropical ocean where temperature decreases rapidly with depth. Groundwater circulation in these platforms has long been associated with the geothermal heat flux because it is capable of generating inward and upward flow of oceanic origin water by buoyancy effects. This study shows that hydraulic circulation occurs even in the absence of a geothermal flux because the combination of the cold subsurface ocean waters with the warm surface conditions is sufficient to maintain a convection cell within the carbonate platform. Using a one-dimensional analytical model, validated by more sophisticated two-dimensional simulations, we can investigate the interaction between these two driving forces. The flow rate inside the platform is, in fact, a function of the ratio of the geothermal flux to the temperature gradient in the ocean. It increases with the geothermal flux but decreases with the oceanic temperature gradient. This one-dimensional model also shows that taking salinity effects on density into account increases the flow rates transiting through the platform by a third.     

Development of groundwater modeling for the Azraq Basin, Jordan

 The three-dimensional groundwater flow model MODFLOW was applied to simulate water level change in the complex multi-aquifer systems (the Upper and Middle Aquifers) of the Azraq basin. The model was calibrated by matching observed and simulated drawdown for steady and transient states over the period 1970–1992. Drawdown data for the period 1993–1997 were used to test the model's ability to predict the response of the aquifers. The model performed well in representing the water level contours of the Upper and Middle Aquifers for steady state calibration. Agreement between the observed and simulated drawdowns was obtained for transient state calibration. To predict the aquifer system responses for the period of 1997–2025, four different pumping schemes (scenarios) have been investigated. The first scenario (present pumping rate) reveals that there will be approximately a 25 m drop in the water level at the well-field area in 2025. However, the worst scenario (pumping rate at 1.5 times the present rate) reveals an approximate 39 m drop in the water level at the well-field area in 2025. The safe yield for the Upper Aquifer System was found to be about 25 million cubic meters (MCM) yearly. Received: 24 June 1999 · Accepted: 30 November 1999     

Hydrogeology of a coal-seam gas exploration area, southeastern British Columbia, Canada: Part 1. Groundwater flow systems

Discovery of high contents of methane gas in coals of the Mist Mountain Formation in the Elk River valley, southeastern British Columbia, Canada, has led to increased exploration activity for coal-seam gas (CSG). CSG production requires groundwater abstraction to depressurize the coal beds and to facilitate methane flow to the production wells. Groundwater abstraction will have hydrodynamic effects on the flow system, and an understanding of the groundwater flow system is needed to evaluate these effects. The purpose of this paper is to describe the groundwater flow system in the area by means of a groundwater flow model and interpretation of hydrochemical and isotopic analyses of groundwater and surface water. Groundwater flow for the Weary Creek exploration area is modeled in two vertical sections. The model domains, based on classic upland–lowland conceptual flow models, are approximately 10,000 m long and 4,000 m deep. Each consists of a fixed water-table boundary and no-flow boundaries along the traces of major faults. Steady-state groundwater flow is calibrated to hydraulic-head, streamflow, and groundwater-recharge data. Simulated steady-state velocity fields define regional and local flow components consistent with the conceptual model. The results are consistent with regional trends in δ²H, δ¹⁸O, tritium, and TDS, which define two distinct groundwater groups (A and B) and a third of intermediate composition. An active, shallow, local flow component (group A) is recharged in beds cropping out along subdued ridges; this component discharges as seeps along lower and mid-slope positions in the southern part of the study area. The waters are tritiated, relatively enriched in δ²H and δ¹⁸O, and have low TDS. A deeper regional flow component (group B), which originates at a higher altitude and which discharges to the Elk River valley bottom, is characterized by non-tritiated groundwater with relatively depleted δ²H and δ¹⁸O, and higher TDS. Groundwater contributes less than 10% of the total direct flow to the Elk River, as indicated by flow measurements and by the absence of group A and group B characteristics in the river water. Thus it is hypothesized that groundwater extraction during CSG production will have little impact on the river. The groundwater flow model developed in this work is used in a companion paper to further test this hypothesis. Electronic Publication     

Fresh-water lenses and practical limitations of their three-dimensional simulation

Fresh-water lenses are the major sources of water supply in many atoll islands in the Pacific and Indian Oceans, particularly in dry seasons. Several two- and three-dimensional models are currently available for the simulation of atoll-island aquifers; however, 2D models cannot include 3D spatial variability of material properties, they must simplify the boundary conditions, and they cannot correctly simulate pumping wells. In an attempt to overcome these difficulties, a 3D model, SALTFLOW, was adopted for the simulation of Home Island in the Indian Ocean. This exercise required a discretisation on the order of a few metres and time steps of a few hours requiring significantly high CPU times. High CPU demand proved to be a difficult challenge but cannot be considered a serious practical limitation with today's advanced computers. The exhaustive data demands of the model (e.g., 3D distributions of hydraulic conductivity, porosity, dispersivities, and spatial and temporal variations of recharge and extraction rates) proved to be more problematical. Although the Home Island data set is unusually comprehensive by any standards, nonetheless the quality and quantity of the available data proved inadequate to meet the calibration needs of a highly karstic aquifer system. The Home Island modeling demonstrates the practical limitations of 3D models. It raises the concern that our ability to develop computer codes capable of simulating complex systems now exceeds our ability to supply the input data necessary for reliable calibration. Finally, the paper demonstrates the importance of the transient calibration in reliable simulation of various management options and emphasises that transient calibration should be considered as an integral part of any similar 2D or 3D modeling. Electronic Publication