储层水动力学数值模拟系统的研究及应用

油气在盆地中产生和运移过程中,水动力因素起着重要的控制作用［１］,如何确定这一作用的实际影响,是一个值得研究的问题。通过数学推导,得出储层中的稳态油层和非稳态油层两种情况下的水动力理论模型,并且对水动力圈闭、深度－压力系统和流体势等作出理论模拟,将其用于油气勘探中,结果与实际情况较相符。     

A very large scale GIS-based groundwater flow model for the Nubian sandstone aquifer in Eastern Sahara (Egypt,northern Sudan and eastern Libya)

A three-dimensional GIS-based groundwater flow model for the Nubian Sandstone Aquifer in the eastern Sahara was developed and calibrated under steady-state and transient conditions. The model was used to simulate the response of the aquifer to climatic changes that occurred during the last 25,000 years. The simulation results indicated that the groundwater in this aquifer was formed by infiltration during the wet periods 20,000 and 5,000 years b.p. The recharge of groundwater due to regional groundwater flow from more humid areas in the south was excluded. It also indicates that the Nubian Aquifer System is a fossil aquifer, which had been in an unsteady state condition for the last 3,000 years.

---

Resumen Fue desarrollado un modelo de flujo de agua subterránea en tres dimensiones, basado en un SIG, para el Acuífero Arenisca Nubian en el Sahara Oriental, el cual fue calibrado para condiciones de estado estacionario y transitorio. El modelo se usó para simular la respuesta del acuífero a los cambios climáticos que ocurrieron durante los últimos 25000 años. Los resultados de esta simulación indicaron que el agua subterránea en este acuífero, se formó por infiltración, durante los períodos húmedos que hubo hace 20000 y 5000 años, antes del presente. Fue excluida la recarga del acuífero debida a un flujo regional de agua subterránea proveniente de áreas con un clima más húmedo en el sur. El modelo también muestra, que el Sistema Acuífero Nubian es un acuífero fósil, el cual ha permanecido en una condición de estado no estacionario, durante los últimos tres mil años.

---

Résumé Pour laquifère gréseux Nubien de Sahara -Est on a mis au points un modèle tridimensionnel, basé sur GIS. Le modèle a été calibré tant pour lécoulement stationnaire que pour lécoulement transitoire. On a simulé après la réponse de laquifère aux changements climatiques des derniers 25000 ans. Les résultats des simulations indiquent que la nappe a été rechargée par des infiltrations pendant une période humide qui sétend 5000 et 20000 ans, dès temps actuel. On na pas pris en compte la recharge de laquifère par la zone plus humide située dans sa partie sud. Le modèle indique aussi que leau de laquifère Nubien est une eau fossile qui a eu un écoulement transitoire pendant les derniers 3000 ans.

     

后竹亭水源地地下水资源量数值模拟

后竹亭水源地是曹洼水源地的一部分.由于地处鲁西南经济发达区,工业取水项目、取水量逐年增加.为正确评价后竹亭水源地供水能力,合理开发水资源,将后竹亭水源地纳入曹洼水源地进行整体评价,并对后竹亭水源地拟订三个开采方案进行计算比较,最后选定方案二,为合理开发利用水资源提供科学依据.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

Groundwater flow and solute transport at the Mourquong saline-water disposal basin, Murray Basin, southeastern Australia

Saline groundwater and drainage effluent from irrigation are commonly stored in some 200 natural and artificial saline-water disposal basins throughout the Murray-Darling Basin of Australia. Their impact on underlying aquifers and the River Murray, one of Australia's major water supplies, is of serious concern. In one such scheme, saline groundwater is pumped into Lake Mourquong, a natural groundwater discharge complex. The disposal basin is hydrodynamically restricted by low-permeability lacustrine clays, but there are vulnerable areas in the southeast where the clay is apparently missing. The extent of vertical and lateral leakage of basin brines and the processes controlling their migration are examined using (1) analyses of chloride and stable isotopes of water (²H/¹H and ¹⁸O/¹⁶O) to infer mixing between regional groundwater and lake water, and (2) the variable-density groundwater flow and solute-transport code SUTRA. Hydrochemical results indicate that evaporated disposal water has moved at least 100 m in an easterly direction and that there is negligible movement of brines in a southerly direction towards the River Murray. The model is used to consider various management scenarios. Salt-load movement to the River Murray was highest in a "worst-case" scenario with irrigation employed between the basin and the River Murray. Present-day operating conditions lead to little, if any, direct movement of brine from the basin into the river. Electronic Publication     

ATLAS III in situ heating test in boom clay: Field data,observation and interpretation

The ATLAS III small scale in situ heating test aimed at assessing the thermo-hydro-mechanical (THM) effects on the Boom clay of the significant temperature gradients generated in this host rock as a consequence of the geological disposal of radioactive, heat-emitting wastes. This paper presents data on temperature, pore water pressure and total stress measured during the experiment and highlights several interesting observations regarding the thermal anisotropy and THM coupling in the Boom clay. The test has a simple geometry and well defined boundary conditions, which facilitates the comparison between measurement and numerical modeling studies. These studies included three dimensional coupled THM modeling of the test. The good agreement between measurement and numerical modeling of temperature and pore water pressure yields a set of THM parameters and confirms the thermo-mechanical anisotropy of the Boom clay.     

Western Alpine back-thrusting as subsidence mechanism in the Tertiary Piedmont Basin (Western Po Plain, NW Italy)

Basin formation dynamics of the Tertiary Piedmont Basin (TPB) are here investigated by means of cross-section numerical modelling. Previous works hypothesised that basin subsidence occurred due first to extension (Oligocene) and then to subsequent loading due to back-thrusting (Miocene). However, structural evidence shows that the TPB was mainly under contraction from Oligocene until post Pliocene time while extension played a minor role. Furthermore, thermal indicators strongly call for a cold (flexure-induced) mechanism but are strictly inconsistent with a hot (thermally induced) mechanism. Our new modelling shows that the TPB stratigraphic features can be reproduced by flexure of a visco-elastic plate loaded by back-thrusts active in the Western Alps in Oligo-Miocene times. Far-field compression contributed to the TPB subsidence and controlled the basin infill geometry by enhancing basin tilting, forebulge uplift and erosion of the southern margin of the basin. These results suggest that the TPB subsidence is the result of a combination of mechanisms including thrust loading and far-field compressional stresses.     

Numerical modeling of tide-surge interaction along Orissa coast of India

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.     

基于CUDA的地震相干体并行算法

相干体技术在地震勘探资料解释方面得到了广泛的应用,由于相干体技术处理的对象是三维地震数据体,所以算法运算时间较长。为了缩短解释周期,本文充分发挥GPU并行计算优势,对C3相干体算法进行并行化分析。从硬盘读取数据到GPU上计算相干值并写入硬盘的整个过程进行分析,剔除了冗余数据的读取,完成了C3相干体算法的并行化设计与实现。最后分别对串行算法与并行算法进行性能测试,结果表明本文设计的并行算法在保证精度的前提下达到了16倍左右的加速比,对加快地震资料解释具有重要意义。     

Pile driving formulas based on pile wave equation analyses

Pile driving formulas, which directly relate the pile set resulting from a hammer blow to the static load capacity of the pile, are often used to decide whether a pile will have the required design capacity. However, existing formulas do not consider soil or pile type, and do not, in general, reliably predict pile capacity. In this paper, an advanced model for dynamic pile driving analysis was used to develop accurate pile driving formulas. The proposed driving formulas are validated through well-documented case histories. Comparisons of predictions from the proposed formulas with results from static and dynamic load tests show that they produce reasonably accurate predictions of pile capacity based on pile set observations.