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

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

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     

国际地面沉降研究综述

文章在介绍国际地面沉降会议历史及２０００年第六届国际地面沉隆会议简况的基础上,并根据２０００年第六届国际地面沉降会议论文,对国际地面沉降研究进展情况分成如下六个方面进行了综述。地面沉降地质因素介绍了古代地面沉降、泥炭层沉降、地震砂土液化地面沉降和海平面上升研究状况。地下流体运移地面沉降方面介绍了以地下水开采为主的地面沉降问题及地面塌陷、天然气开采引起的地面沉降、均匀沉降对建筑结构的破坏和欠固结石英     

接触摩擦问题的数值模拟

无网格伽辽金法(EFGM)可脱离单元的概念,特别适合岩体裂纹面的接触摩擦分析。基于EFGM,在裂纹面引入罚参数,通过迭代计算,得到裂纹面真实的应力状态,从而模拟闭合裂纹的粘接、滑移和张开行为,数值结果表明该方法是合理可行的。     

环境模拟和GIS集成的初步研究

GIS与环境模拟在技术、研究内容、方法上的进一步集成,具有广泛的应用前景.以开展的流域土壤和水资源研究模型的集成和系统化及其应用项目的实际工作,说明了环境模拟和GIS技术集成的必然性;从GIS技术的发展,环境模拟模型的改进与完善,数据的初始化与管理等方面系统地论述了环境模拟模型和GIS集成的必要性.通过现有集成模式的分析,对二者集成的概念框架以及不同的集成方式进行了概念性的描述,结合研究项目给出了集成示例.     

天津市深层基岩地下热水系统数值模拟中若干问题的处理

天津市基岩地下热水具有密度随温度变化显著、存在水岩之间热交换的特点, 采用适合于描述深层地下热水水流和热量运移的压力场和温度场控制方程, 通过对渗流场、边界条件的合理概化和对初始条件、热储层缺失区域及倾斜断层等的合理处理, 使所建立的三维数值模型变得更为有效和实用.